function xdot=LCfullE5fun(t,x,flag,Vfinal)
% Carol Lucas
% Called function for F4_3.m

V=-100;
if t>=50 & t<=300 
    V=Vfinal;
end;
k1p=0.78/(1+exp(-(3+V)/10)); 
k1m=0.78-k1p;
caout=3;
k2pcad=-7.56e-4*caout*V/(1-exp(V/13.4));
k2m=0.002;

xc=x(1,1);
xo=x(2,1);

%% Second order model
V1=k1p*xc-k1m*xo;
V2=k2pcad*xo-k2m*(1-xc-xo);
xdot(1,1)=-V1;
xdot(2,1)=V1-V2;