# Morris Lecar Model as in book and my notes # # p Iapp=-10 p gc=0 p C=20 p Vk=-84 p gk=8 p Vca=120 p gca=8 p Vlk=-60 p glk=1 p phi=.5 p v1=-1.2 p v2=18 p v3=-5 p v4=10 i v=-20 i w=.2 i vv=-30 i ww=0 minf(v)=.5*(1+tanh((v-v1)/v2)) winf(v)=.5*(1+tanh((v-v3)/v4)) tau(v)=1/cosh((v-v3)/(2*v4)) v'=(-gca*minf(v)*(v-Vca)-gk*w*(v-Vk)-glk*(v-Vlk)+Iapp+gc*(vv-v))/C w'=phi*(winf(v)-w)/tau(v) vv'=(-gca*minf(vv)*(vv-Vca)-gk*ww*(vv-Vk)-glk*(vv-Vlk)+Iapp+gc*(v-vv))/C ww'=phi*(winf(vv)-ww)/tau(vv) @ total=300,dt=.2,xhi=300,ylo=-30,yhi=30 done