WRITE "BIOMD10"$

% B_ IS THE VARIABLE VECTOR
B_:={x1,x2,x3,x4,x5,x6,x7,x8,y1,y2,y3}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16,p17,p18,p19,p20,p21,p22}$

%NUMBER OF STATES 
NX_:=8$
%NUMBER OF INPUTS 
NU_:=0$
%NUMBER OF OUTPUTS 
NY_:=3$

%MODEL EQUATIONS
C_:={df(x1,t)=(x2*p5)/(p6+x2)-(x1*p1)/(((x8/p2)^p3+1)*(p4+x1)),
     df(x2,t)=(x1*p1)/(((x8/p2)^p3+1)*(p4+x1))-(x2*p5)/(p6+x2),
     df(x3,t)=(x4*p13)/(p14+x4)-(x3*x2*p7)/(p8+x3),
     df(x4,t)=(x5*p11)/(p12+x5)-(x4*p13)/(p14+x4)+(x3*x2*p7)/(p8+x3)-(x4*x2*p9)/(p10+x4),
     df(x5,t)=(x4*x2*p9)/(p10+x4)-(x5*p11)/(p12+x5),
     df(x6,t)=(x7*p21)/(p22+x7)-(x6*x5*p15)/(p16+x6),
     df(x7,t)=(x8*p19)/(p20+x8)-(x7*p21)/(p22+x7)+(x6*x5*p15)/(p16+x6)-(x7*x5*p17)/(p18+x7),
     df(x8,t)=(x7*x5*p17)/(p18+x7)-(x8*p19)/(p20+x8),
     y1=x6,
     y2=x7,
     y3=x8}$

SEED_:=65$
DAISY()$

% INITIAL CONDITIONS
IC_:={x1=90,x2=10,x3=280,x4=10,x5=10,x6=280,x7=10,x8=10}$
CONDINIZ()$
END$

