WRITE "Glycolysis"$

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,x5,u1,u2,u3,u4,y1,y2,y3,y4,y5}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR    
B1_:={k1,k2,k3,k4,kM}$

%NUMBER OF STATES 
NX_:=5$
%NUMBER OF INPUTS 
NU_:=4$
%NUMBER OF OUTPUTS 
NY_:=5$ 

%MODEL EQUATIONS
C_:={df(x1,t)=(-k1*x1/(x1+kM))*u1,
     df(x2,t)=(k1*x1/(x1+kM))*u1-(k2*x2/(x2+kM))*u2,
     df(x3,t)=(k2*x2/(x2+kM))*u2-(k3*x3/(x3+kM))*u3,
     df(x4,t)=(k2*x2/(x2+kM))*u2+(k3*x3/(x3+kM))*u3-(k4*x4/(x4+kM))*u4,
     df(x5,t)=(k4*x4/(x4+kM)*u5,
     u1=1,
     u2=1,
     u3=1,
     u4=1,
     u5=1,
     x1=y1,
     x2=y2,
     x3=y3,
     x4=y4,
     x5=y5}$

SEED_:=25$
DAISY()$

% INITIAL CONDITIONS
END$


