We include diffusion term D ΔE =
D ∂2E / ∂x2 to get running waves
∂t E = DΔE + F(E, n) =
DΔE + [A - (n/nb)M][1 - tanh(E - 3)]E2/2 - E,
∂t n = H(E, n) = ε [Θ(E - 1) - n].
M
nb
ε
Δt
Δx
D
speed
Ymin
Ymax
delay
it
fps
This script uses simple explicit algorithm.
It makes it time steps per frame.