u1'=-u1+f(w11*u1+w12*u2-te1)
u2'=(-u2+f(w21*u1+w22*u2-te2))/tau
par te1=2,te2=6,tau=1,w11=12,w12=-8,w21=16,w22=-2
f(u)=1/(1+exp(-u))
done