# Gray-Scott
u0=u1
v0=v1
u[1..100]'=-u[j]*(v[j])^2+F*(1-u[j])+D1*(u[j+1]-2*u[j]+u[j-1])/(h*h)
v[1..100]'=u[j]*(v[j])^2-(k+F)*v[j]+D2*(u[j+1]-2*u[j]+u[j-1])/(h*h)
u101=u100
v101=v100
init u[1..100]=0.25
init v[1..100]=0.6 
par F=.12,k=.03,D1=0.00002,D2=0.00001,h=0.01
@ total=20, maxstor=1000000
done