dx1/dt=k13*x3+k12*x2-(k21+k31)*x1+u1;
dx2/dt=k21*x1-(k12+k02)*x2;
dx3/dt=k31*x1-(k13+k03)*x3;
y1=x1/V1
