function xdot=mlsynfun(t,x,flag,alpha,beta,ton,bswitch)
%%  Carol Lucas
% Morris-Lecar dynamics with modified parameters.
% Notably: vc=-5; vd=10; phi=0.5  (originally: 2; 30; 0.04; respectively)
%
% Use weak alpa beta= 3,1 for exitatory, alpha beta = 3,.1 for inhibatory
%  For Fig 6.7, beta swithes from 1 to .1 at tbeta=800
%%  init=[-44;0;0;.3;0;.3] Fig 6.6 or init=[-55;0;-20;0;0;0] Fig 6.7

V1=x(1,1);
w1=x(2,1);
V2=x(3,1);
w2=x(4,1);
s1=x(5,1);
s2=x(6,1);

if and(bswitch>0,t>800)
	beta=bswitch;
	end
vk=-84;vl=-60;vca=120;
i=15;gk=8;gl=1;c=20;
va=-1.2;vb=18;
vc=-5;vd=10;phi=0.1;gca=8;
thetasyn=20;ksyn=2;gsyn=1;vsyn=0;

minf1=.5*(1+tanh((V1-va)/vb));
minf2=.5*(1+tanh((V2-va)/vb));
winf1=.5*(1+tanh((V1-vc)/vd));
tauw1=1/cosh((V1-vc)/(2*vd));
winf2=.5*(1+tanh((V2-vc)/vd));
tauw2=1/cosh((V2-vc)/(2*vd));
sinf1=1/(1+exp(-(V1-thetasyn)/ksyn));
sinf2=1/(1+exp(-(V2-thetasyn)/ksyn));

xdot(1,1)=(i-gca*minf1*(V1-vca)-gk*w1*(V1-vk)-gl*(V1-vl)+(t>=ton)*gsyn*s2*(vsyn-V1))/c;
xdot(2,1) = phi*(winf1-w1)/tauw1;
xdot(3,1) = (i-gca*minf2*(V2-vca)-gk*w2*(V2-vk)-gl*(V2-vl)+(t>=ton)*gsyn*s1*(vsyn-V2))/c;
xdot(4,1) = phi*(winf2-w2)/tauw2;
xdot(5,1) = alpha*sinf1*(1-s1)-beta*s1;
xdot(6,1) = alpha*sinf2*(1-s2)-beta*s2;