# fitzhugh-nagumo action potential
f(v,w)=v*(1-v)*(v-a)-w+i
g(v,w)=eps*(v-gam*w)
par a=.05,i=0,eps=.01,gam=.2
par d=.5
v0'=f(v0,w0)+d*(v1-v0)
v[1..49]'=f(v[j],w[j])+d*(v[j+1]-2*v[j]+v[j-1])
v50'=f(v50,w50)+d*(v49-v50)
w[0..50]'=g(v[j],w[j])
@ total=200,dt=.25,meth=qualrk,tol=1e-6
@ xhi=200,yp=v20
done