WRITE "Transfection2"$

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={d1,d2,d3,b,kTL,m0,e0}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-d1*x1-d2*x1*x3,
     df(x2,t)=kTL*x1-b*x2,
     df(x3,t)=d3*x4-d2*x1*x3,
     df(x4,t)=-d3*x4+d2*x1*x3,
     y1=x2}$

SEED_:=65$
DAISY()$

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

