Follow one nerve impulse from a dendritic input to the next cell and measure its speed and delay.
Educational model, not medical advice and not a diagnostic tool. Channel sliders are abstract parameters.
Threshold: press the button (it bisects the dendritic conductance that makes the axon hillock fire).
The science
Membrane (Hodgkin-Huxley 1952). Cm dV/dt = I - gNa m³h (V-ENa) - gK n⁴ (V-EK) - gL (V-EL), gating dx/dt = φ [αx(1-x) - βx x], classic squid values (Cm = 1 uF/cm², gNa = 120, gK = 36, gL = 0.3 mS/cm², ENa = 50, EK = -77, EL = -54.4 mV) and a temperature factor φ = 3(T - 6.3)/10. Above about 32 C these squid parameters stop producing propagating spikes on a bare axon (a known limit of the classic parameter set), so the temperature slider ends there.
Cable. The neuron is a chain of compartments: C dV/dt = ∑ Gij(Vj - Vi) - Iion + Isyn, axial conductance G = π a² / (Ri L) (Ri = 35.4 ohm cm). Sub-threshold voltage decays as V(x) = V0 e-x/λ, λ = √(a Rm / 2Ri) (376 um for the 2 um dendrite). Solved with backward Euler in V (stable for the stiff axial term), exponential (Rush-Larsen) gate updates, tridiagonal sweep; the textbook point neuron is also integrated with RK4 for the checks.
Myelin. Active channels sit only at the nodes of Ranvier (10x the squid density, same g ratios so the rest potential stays -65 mV). Each node carries the capacitance and leak of its internode divided by the number of myelin lamellae (25 per um of diameter, constant g-ratio). Then speed scales with diameter (about 6 m/s per um in the literature; here roughly 5-6 m/s per um), while a bare axon scales with the square root of the diameter. Both laws are verified numerically below.
Spike initiation. The axon initial segment has a high Na channel density, so the hillock reaches threshold first; a dendritic input only decays toward it (attenuation by e-x/λ). Inputs are a dual-exponential conductance gs(t) = gmax(e-t/τd - e-t/τr)/norm, τr = 0.5 ms, τd = 3 ms, Es = 0 mV.
Synapse. The presynaptic spike amplitude sets the calcium signal, Prel = Ca⁴ / (Ca⁴ + K⁴) (K = 10 uM, Ca = 20 uM x amplitude²), 10 release sites give a binomial number of quanta, and the delay is a calcium/fusion lag (Q10 = 2) plus 0.2 ms of diffusion, about 0.9 ms here (0.5 to 1 ms is typical for warm mammalian synapses; this one is at 24 C). Cells B and C are leaky integrate-and-fire units (tau 10 ms, threshold -50 mV, 2 ms refractory); a negative D weight is an inhibitory input (E = -80 mV).
Honest simplifications. One-dimensional cable, fixed geometry, lumped synaptic input, no calcium-activated K channels, no channel noise, no ion concentration changes, deterministic except for the seeded release. The 3D scene is not to scale: the 5 mm axon is drawn as a path of about 17 scene units, and the foreground axon is decorative. All parameters are model values, check a physiology text before quoting numbers. Models were built with Blender and rendered with Three.js.
Photosensitivity: the glow follows the simulated voltage locally and smoothly; the background sparkle twinkles slower than 1 Hz and auto-fire is limited to one input per 2 s.
Educational model, not medical advice.