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

k1p=0.7; k1m=0.2; k2p=0.05; k2m=0.005; cad=0.3;
K1=k1m/k1p;
tau=(1+K1)/(k2p*cad+k2m*(K1+1));
tauact=1/(k1p+k1m);
yinf=k2m*tau;

xc=x(1,1);
xo=x(2,1);
xoalt=x(3,1);
y=x(4,1);

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

%% First order models
xdot(3,1)=-(xoalt-(1/(1+K1)))/tauact;
xdot(4,1)=(yinf-y)/tau;