function xdot=mlgapfun(t,x,flag,gc)
%%  Carol Lucas
% Morris-Lecar dynamics with modified parameters.
% Notably: vc=-5; vd=10; phi=0.5  (originally: 2; 30; 0.04; respectively)
%
% Use weak gc (=1) for antiphase and strong gc (=2) for inphase.
%%  init=[-20;.2;-30;0]

V1=x(1,1);
w1=x(2,1);
V2=x(3,1);
w2=x(4,1);

vk=-84;vl=-60;vca=120;
i=-10;gk=8;gl=1;c=20;
va=-1.2;vb=18;ton=100;
vc=-5;vd=10;phi=0.5;gca=8;

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));

xdot(1,1)=(i-gca*minf1*(V1-vca)-gk*w1*(V1-vk)-gl*(V1-vl)+gc*(t>=ton)*(V2-V1))/c;
xdot(2,1) = phi*(winf1-w1)/tauw1;
xdot(3,1) = (i-gca*minf2*(V2-vca)-gk*w2*(V2-vk)-gl*(V2-vl)+gc*(t>=ton)*(V1-V2))/c;
xdot(4,1) = phi*(winf2-w2)/tauw2;
%