dx1/dt=p1*x6/(p3+x6)-p5*x1/(p12+x1)+(p26*x7)*u1;
dx2/dt=p19*x1-p22*x2+p23*x3-p6*x2/(p13+x2);
dx3/dt=p22*x2-p23*x2-p7*x3/(p14+x3);
dx4/dt=p2*p4^2/(p4^2+x3^2)-p8*x4/(p15+x4);
dx5/dt=p20*x4-p24*x5+p25*x6-p9*x6/(p16+x5);
dx6/dt=p24*x5-p25*x6-p10*x6/(p17+x6);
dx7/dt=p21-p11*x7/(p18+x7)+-(p21+p27*x7)*u2;
y1=x1;
y2=x2;
x1(0)=0;
x2(0)=0;
x3(0)=0;
x4(0)=0;
x5(0)=0;
x6(0)=0;
x7(0)=0
