u1'=-u1+f(w11*u1+w12*u2-te1+bee*u1p-bie*u2p)
u2'=(-u2+f(w21*u1+w22*u2-te2+bei*u1p-bii*u1p))/tau
u1p'=-u1p+f(w11*u1p+w12*u2p-te1+bee*u1-bie*u2)
u2p'=(-u2p+f(w21*u1p+w22*u2p-te2+bei*u1-bii*u1))/tau
init u1=.5,u2=.5,u1p=.4,u2p=.6
par bee=0,bie=0,bii=0,bei=0
par te1=1,te2=6,tau=1,w11=12,w12=-8,w21=16,w22=-2
f(u)=1/(1+exp(-u))
@ total=100,xlo=80,xhi=100,ylo=0,yhi=1
done