# fhnfor.ode
# Fitzhugh-Nagumo equations with sinusoidal forcing
v'=I+v*(1-v)*(v-a) -w+c*u
w'=eps*(v-gamma*w)
u'=u*(1-u^2-z^2)-f*z
z'=z*(1-u^2-z^2)+f*u
init u=1,v=.0505,w=.1911
par c=.05,f=.8
par I=0.2,a=.1,eps=.1,gamma=.25
@ xp=V,yp=w,xlo=-.25,xhi=1.25,ylo=-.5,yhi=1,total=7.86
@ dsmax=.05,parmin=-.5,parmax=1
@ autoxmax=1,autoxmin=0,autoymin=-1.5,autoymax=1.5
done