%%LCtsfun
%% ODEs for normalized LCmodel (F4_3)
%% Equations 4.13 and 4.14

function xdot=LCtsfun(t,x)
k1p=0.7; k1m=0.2; k2p=0.05; k2m=0.005; cad=0.3;

xcts=x(1,1);
xots=x(2,1);

eps=k2m/k1m;
xdot(1,1)=(-(k1p/k1m)*xcts+xots)/eps;
xdot(2,1)=((k1p/k1m)*xcts-xots-eps*(k2p/k2m)*cad*xots+eps*(1-xots-xcts))/eps;