v'= -gl*(v-vl)-gca*m*(v-vca)-gk*n*(v-vk)+i+gc1*(vp-v)+ips*pulse(t)
vp'=-gl*(vp-vl)+gc2*(v-vp)+ipd*pulse(t)
m'=am(v)*(1-m)-bm(v)*m
#m=am(v)/(am(v)+bm(v))
n'=an(v)*(1-n)-bn(v)*n
pulse(t)=1/((1+exp(-(t-ton)*4))*(1+exp(4*(t-toff))))
init v=-70,vp=-70
aux minf=am(v)/(am(v)+bm(v))
am(v)=exp((v-v1)/v2)
bm(v)=exp(-(v-v3)/v4)
an(v)=exp((v-v5)/v6)
bn(v)=exp(-(v-v7)/v8)
par vl=-70,vk=-90,vca=150,i=0,gca=1,gk=2.5,gl=.1
par gc1=1,gc2=1
par ips=0,ipd=0,ton=10,toff=20
par v1=-8,v2=8,v3=-11,v4=8
par v5=16,v6=8,v7=-26,v8=8
@ meth=cvode,toler=1e-5,atol=1e-4
@ ylo=-80,yhi=30,total=50,xhi=50
done