WRITE "Arabidopsis"$

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,x5,x6,x7,u1,u2,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}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=p1*x6/(p3+x6)-p5*x1/(p12+x1)+u1,
     df(x2,t)=p19*x1-p22*x2+p23*x3-p6*x2/(p13+x2),
     df(x3,t)=p22*x2-p23*x2-p7*x3/(p14+x3),
     df(x4,t)=p2*p4^2/(p4^2+x3^2)-p8*x4/(p15+x4),
     df(x5,t)=p20*x4-p24*x5+p25*x6-p9*x6/(p16+x5),
     df(x6,t)=p24*x5-p25*x6-p10*x6/(p17+x6),
     df(x7,t)=p21-p11*x7/(p18+x7)+u2,
     u1=p26*x7,
     u2=-p21-p27*x7,
     y1=x1,
     y2=x4}$

SEED_:=65$
DAISY()$

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

