We include diffusion term D Δu =
D ∂2u / ∂x2 to get running waves (D = 0.001 cm2/ms)
∂t u = DΔu -
(Jfi + Jso + Jsi ).
Δtms
Δxcm
D
speed
delay
Ymin
Ymax
it
fps
This script uses simple explicit algorithm. It makes it time steps
per frame. x length is 1024Δx.
Set it = 100, delay = 100 to test wave speed.
Heart rhythmsupdated 8 Dec 2011