WRITE "GoodwinPoly"$

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

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

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

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

%MODEL EQUATIONS
C_:={df(x1,t)=p1*xi1 - p4*x1,
     df(x2,t)=p5*x1 - p6*x2,
     df(x3,t)=p7*x2 - p8*x3,
     df(xi1,t)=-p3*x3^(p3 - 1)*xi1^2*(p7*x2 - p8*x3),
     y1=x1,
     y2=x2,
     y3=x3}$

SEED_:=25$
DAISY()$

% INITIAL CONDITIONS
IC_:={x1=3/10,
      x2=9/10,
      x3=13/10,
      xi1=1/(p2 + (13/10)^p3)}$
CONDINIZ()$
END$

