WRITE "Bilirubin2"$

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,u1,y1}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={k01,k12,k13,k14,k21,k31,k41}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-(k21+k31+k41+k01)*x1+k12*x2+k13*x3+k14*x4+u1,
     df(x2,t)=k21*x1-k12*x2,
     df(x3,t)=k31*x1-k13*x3,
     df(x4,t)=k41*x1-k14*x4,
     y1=x1}$

SEED_:=65$
DAISY()$

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


