WRITE "Thyroid2"$

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={k02,k03,k12,k13,k21,k31,V1}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=k13*x3+k12*x2-(k21+k31)*x1+u1,
     df(x2,t)=k21*x1-(k12+k02)*x2,
     df(x3,t)=k31*x1-(k13+k03)*x3,
     y1=x1/V1}$

SEED_:=65$
DAISY()$

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

