The LR model. 1D travelling waves
Again we include diffusion term D Δu =
D ∂2u / ∂x2 to get running waves
(D = 0.001 cm2/ms).
For Δx = 0.01 cm you can see alternans.
For Δx = 0.008 cm waves can not propagate.
This script uses simple explicit algorithm. It makes it time steps
per frame. x length is 1024Δx.
Set it = 100 to test wave speed.
updated 9 Feb 2012