# Gray-Scott

u0=u1

v0=v1

u[1..500]'=-u[j]*(v[j])^2+f*(1-u[j])+D1*(u[j+1]-2*u[j]+u[j-1])/(h*h)

v[1..500]'=u[j]*(v[j])^2-(k+f)*v[j]+D2*(v[j+1]-2*v[j]+v[j-1])/(h*h)

u501=u500

v501=v500

init u[1..100]=0

init u[101..150]=0.2

init u[151..200]=0

init u[201..300]=1

init u[301..500]=0

init v[1..100]=0

init v[101..150]=0.3

init v[151..500]=0

par f=0.02,k=0.035,D1=0.0002,D2=0.0001,h=0.01

@ total=500, maxstor=10000

done

# Gray-Scott ODE

u'=-u*v^2+f*(1-u)

v'=u*v^2-(k+f)*v

par f=.02,k=.035

@ XP=u, YP=v, XLO=0, XHI=1, YLO=0, YHI=1

done