dx1/dt=kprod-kdeg*x1-k1*x1*u1;
dx2/dt=-k3*x2-kdeg*x2-a2*x2*x10+t1*x4-a3*x2*x13+t2*x5(k1*x1-k2*x2*x8)*u2;
dx3/dt=k3*x2-kdeg*x3k2*x2*x8*u3;
dx4/dt=a2*x2*x10-t1*x4;
dx5/dt=a3*x2*x13-t2*x5;
dx6/dt=c6a*x13-a1*x6*x10+t2*x5-i1*x6;
dx7/dt=i1*kv*x6-a1*x11*x7;
dx8/dt=c4*x9-c5*x8;
dx9/dt=c2+c1*x7-c3*x9;
dx10/dt=-a2*x2*x10-a1*x10*x6+c4a*x12-c5a*x10-i1a*x10+e1a*x11;
dx11/dt=-a1*x11*x7+i1a*kv*x10-e1a*kv*x11;
dx12/dt=c2a+c1a*x7-c3a*x12;
dx13/dt=a1*x10*x6-c6a*x13-a3*x2*x13+e2a*x14;
dx14/dt=a1*x11*x7-e2a*kv*x14;
dx15/dt=c2c+c1c*x7-c3c*x15;
y1=x7;
y2=x10+x13;
y3=x9;
y4=x1+x2+x3;
y5=x2;
y6=x12;
x1(0)=x01;
x2(0)=0;
x3(0)=0;
x4(0)=0;
x5(0)=0;
x6(0)=x06;
x7(0)=x07;
x8(0)=x08;
x9(0)=x09;
x10(0)=x010;
x11(0)=x011;
x12(0)=x012;
x13(0)=NF;
x14(0)=x014;
x15(0)=x015
