Excitation cI(t) with period T at the left border is added
to the system
∂t E = DΔE + F(E, n) + c I(t)δ(x)
= DΔE + [A - (n/nb)M][1 - tanh(E - 3)]E2/2 - E +
c I(t)δ(x),
∂t n = H(E, n) = ε [Θ(E - 1) -
n],
I(t) ~ [cos(πt/T)]500,
-δ ∫ δI(t)dt = 1.
M
nb
ε
Δt
Δx
D
T
c
Ymin
Ymax
delay
it
fps
The script makes it time steps per frame.
Really for large M wave speed doesn't depend on recovery state.
But e.g. for M = 1 the first impulse is much faster.