u1'=-u1+f(w11*u1+w12*u2+w13*u3-te1)
u2'=(-u2+f(w21*u1+w22*u2+w23*u3-te2))/tau
u3'=(-u3+f(w31*u1+w32*u2+w33*u3-te3))/tau3
par te1=-2,te2=4.5,tau=1.35,w11=8,w12=-6,w21=16,w22=-2
par te3=8,w13=-10,w23=0,tau3=60,w31=8,w32=0,w33=0
f(u)=1/(1+exp(-u))
init u1=.8,u2=.99,u3=.098
@ dt=.05,njmp=5,total=1000,xhi=1000,yhi=1,ylo=0
done