WRITE "Bilirubin1"$

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={k03,k04,k13,k24,k31,k42,k43}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-k31*x1+k13*x3+u1,
     df(x2,t)=-k42*x2+k24*x4,
     df(x3,t)=k31*x1-(k03+k13+k43)*x3,
     df(x4,t)=k42*x2+k43*x3-(k04+k24)*x4,
     y1=x1,
     y2=x2}$

SEED_:=65$
DAISY()$

END$

