The BCF model. Speed variations

Excitation cI(t) with period T is added to the system at the left border
    ∂t u = DΔu - (Jfi + Jso + Jsi ) + cI(t) ,     I(t) ~ [cos(πt/T)]500,     -δ ∫ δI(t)dt = 1.
Δtms Δxcm D T c

delay Ymin Ymax it fps
This script uses simple explicit algorithm. It makes it time steps per frame. x length is 1024Δx.
Heart rhythms     updated 8 Dec 2011