p gamma=1,eps=.1,a=.2,d=1,d2=1
i v[1..10]=0
i v[11..20]=1
i v[21..160]=0
i w[1..10]=.2
i w[10..160]=0

f(v)=v*(1-v)*(v-a)

v0'=f(v0)-w0+d*(v1-v0)
v[1..100]'=f(v[j])-w[j]+d*(v[j-1]-2*v[j]+v[j+1])
v[101..159]'=f(v[j])-w[j]+d2*(v[j-1]-2*v[j]+v[j+1])
v160'=f(v160)-w160+d2*(v159-v160)
w[0..160]'=eps*(v[j]-gamma*w[j])
@ meth=qualrk,dt=.25,total=300,xhi=150
done