dx1/dt=k1-k2*x1/(k3+x1);
y1=x1;

