Write "ArabidopsisPoly"

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,x5,x6,x7,xi1,xi2,xi3,xi4,xi5,xi6,xi7,xi8,xi9,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 p8 p9 p10 p11 p12 p13 p14 p15 p16 p17 p18 p19 p20 p21 p22 p23 p24 p25 p26 p27}$

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


%MODEL EQUATIONS
C_:={df(x1,t)=p1*x6*x9 - p5*x1*x10,
     df(x2,t)=p19*x1 - p22*x2 + p23*x3 - p6*x2*x11,
     df(x3,t)=p22*x2 - p23*x2 - p7*x3*x12,
     df(x4,t)=p2*p4^2*x8 - p8*x4*x13,
     df(x5,t)=p20*x4 - p24*x5 + p25*x6 - p9*x6*x14,
     df(x6,t)=p24*x5 - p25*x6 - p10*x6*x15,
     df(x7,t)=p21 - p11*x7*x16,
     df(x8,t)=2*x3*x8^2*(p23*x2 - p22*x2 + p7*x3*x12),
     df(x9,t)=x9^2*(p25*x6 - p24*x5 + p10*x6*x15),
     df(x10,t)=-x10^2*(p1*x6*x9 - p5*x1*x10),
     df(x11,t)=-x11^2*(p19*x1 - p22*x2 + p23*x3 - p6*x2*x11),
     df(x12,t)=x12^2*(p23*x2 - p22*x2 + p7*x3*x12),
     df(x13,t)=-x13^2*(p2*p4^2*x8 - p8*x4*x13),
     df(x14,t)=-x14^2*(p20*x4 - p24*x5 + p25*x6 - p9*x6*x14),
     df(x15,t)=x15^2*(p25*x6 - p24*x5 + p10*x6*x15),
     df(x16,t)=-x16^2*(p21 - p11*x7*x16),
     x1=y1,
     x4=y2}$

SEED_:=25$
DAISY()$

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




