WRITE "JAK-STAT"$

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,x5,x6,x7,x8,x9,u1,y1,y2}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={p1,p2,p3,p4,p5,p6,p7}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-u1+(p7*p4*x9)/p6,
     df(x2,t)=u1-2*p2*x2^2,
     df(x3,t)=p2*x2^2-p3*x3,
     df(x4,t)=-(p7*p4*x4-p6*p3*x3)/p7,
     df(x5,t)=-p4*(x5-2*x4),
     df(x6,t)=p4*(x5-x6),
     df(x7,t)=p4*(x6-x7),
     df(x8,t)=p4*(x7-x8),
     df(x9,t)=p4*(x8-x9),
     u1=p1*x1,
     y1=(x2+2*x3)/p5,
     y2=(x1+x2+2*x3)/p5}$

SEED_:=65$
DAISY()$

% INITIAL CONDITIONS
%IC_:={x1=p5,x2=0,x3=0,x4=0,x5=0,x6=0,x7=0,x8=0,x9=0}$
%CONDINIZ()$
END$
