The BCF model. 1D travelling waves

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 rhythms     updated 8 Dec 2011