WRITE "pharmacokineticsPoly"$

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={a1,a2,b1,b2,ka,kc,vm}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-a1*(x1-x2)-ka*vm*x1*xi1,
     df(x2,t)=a2*(x1-x2),
     df(x3,t)=-b1*(x3-x4)-kc*vm*x3*xi1,
     df(x4,t)=b2*(x3-x4),
     df(xi1,t)=ka*xi1^2*(a1*(x1-x2)+ka*vm*x1*xi1)+kc*xi1^2*(b1*(x3-x4)+kc*vm*x3*xi1),
     y1=x1,
     y2=x4}$

SEED_:=25$
DAISY()$

% INITIAL CONDITIONS
%IC_:={x1=c0,
%      x2=0,
%      x3=c0*g,
%      x4=0,
%      xi1=1/(c0*ka+ka*kc+c0*g*kc)}$
%CONDINIZ()$
END$
