#S Coombes Jan 2007 #A cortical neuron model #Based around FitzHugh-Nagumo model with quadratic recovery variable #Exhibits saddle-node on a limit cycle bifurcation and illustrates type I firing par J=0.1 mu=0.01 beta=.5 v1=0 v2=0.2 a=0.1 v(0)=0 w(0)=0 f(x)=x*(x-a)*(1-x) v'=1.0/mu*(f(v) - w + J) w'=beta*(v-v1)*(v-v2)-0.5*w @ meth=cvode,dt=0.001,total=50,maxstore=10000000,xp=t,xlo=0,xhi=50,yp=v,ylo=-0.3,yhi=1.1 done done #saddle node on limit cycle at J=0.009396