# wilson-cowan-1972

se(x) = 1/(1+exp(-ae*(x-the))) - 1/(1+exp(ae*the))
si(x) = 1/(1+exp(-ai*(x-thi))) - 1/(1+exp(ai*thi))
dE/dt = (-E+(ke-re*E)*se(c1*E-c2*I+P))/taue
dI/dt = (-I+(ki-ri*I)*si(c3*E-c4*I+Q))/taui

Param P=0,Q=0
param c1=12,c2=4,c3=13,c4=11
param ae=1.2,the=2.8,ai=1,thi=4
param re=1,ri=1
param taue=1,taui=1
param ke=1,ki=1
init E=.1,I=.05
@ total=40
@ xp=I,yp=E,xlo=-.1,xhi=.5,ylo=-.1,yhi=.5
done