V'=-gl*(V-Vl) - (gampa*fa(t) + gnmda/(1+.15*exp(-.08*V)))*V-ggaba*fg(t)*(V-vgaba)
par gnmda=2,gl=.1,vl=-70,vgaba=-80,gampa=.1,ggaba=.35
par ta1=20,ta2=30,tg1=150,tg2=160
fa(t)=heav(t-ta1)*heav(ta2-t)
fg(t)=1+.25*heav(t-tg1)*heav(tg2-t)
init v=-70.55
@ ylo=-90,yhi=0,xhi=200,total=200,dt=.1
done