WRITE "BIOMD03"$

% B_ IS THE VARIABLE VECTOR
B_:={x1,x2,x3,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}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=p4-x1*p5-x1*x3*p6*(x1+p7)^-1,
     df(x2,t)=-x2*p9*(p10+x2)^-1-x1*p1*(x1+p3)^-1*(x2-1)*(p8-x2+1)^-1,
     df(x3,t)=-p13*x3*(p12+x3)^-1-x2*p2*(x3-1)*(p11-x3+1)^-1,
     y1=x1,
     y2=x2,
     y3=x3}$

SEED_:=65$
DAISY()$

% INITIAL CONDITIONS
IC_:={x1=1/100,x2=1/100,x3=1/100}$
CONDINIZ()$
END$


