WRITE "HighDimNonLin"$

% Do we need to define x01 etc. as parameter??!!

% B_ IS THE VARIABLE VECTOR 
B_:={x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13,y14,y15,y16,y17,y18,y19,y20}$

FOR EACH EL_ IN B_ DO DEPEND EL_,T$

%B1_ IS THE UNKNOWN PARAMETER VECTOR
B1_:={vm,km,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16,p17,p18,p19,p20}$

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

%MODEL EQUATIONS
C_:={df(x1,t)=-vm*x1/(km+x1)-p1*x1,
     df(x2,t)=p1*x1-p2*x2,
     df(x3,t)=p2*x2-p3*x3,
     df(x4,t)=p3*x3-p4*x4,
     df(x5,t)=p4*x4-p5*x5,
     df(x6,t)=p5*x5-p6*x6,
     df(x7,t)=p6*x6-p7*x7,
     df(x8,t)=p7*x7-p8*x8,
     df(x9,t)=p8*x8-p9*x9,
     df(x10,t)=p9*x9-p10*x10,
     df(x11,t)=p10*x10-p11*x11,
     df(x12,t)=p11*x11-p12*x12,
     df(x13,t)=p12*x12-p13*x13,
     df(x14,t)=p13*x13-p14*x14,
     df(x15,t)=p14*x14-p15*x15,
     df(x16,t)=p15*x15-p16*x16,
     df(x17,t)=p16*x16-p17*x17,
     df(x18,t)=p17*x17-p18*x18,
     df(x19,t)=p18*x18-p19*x19,
     df(x20,t)=p19*x19-p20*x20,
     x1=y1,
     x2=y2,
     x3=y3,
     x4=y4,
     x5=y5,
     x6=y6,
     x7=y7,
     x8=y8,
     x9=y9,
     x10=y10,
     x11=y11,
     x12=y12,
     x13=y13,
     x14=y14,
     x15=y15,
     x16=y16,
     x17=y17,
     x18=y18,
     x19=y19,
     x20=y20}$

SEED_:=25$
DAISY()$

% INITIAL CONDITIONS
IC_:={x1=x01,x2=x02,x3=x03,x4=x04,x5=x05,x6=x06,x7=x07,x8=x08,x9=x09,x10=x010,
      x11=x011,x12=x012,x13=x013,x14=x014,x15=x015,x16=x016,x17=x017,x18=x018,x19=x019,x20=x020}$
CONDINIZ()$
END$


