dx1/dt=k1-k2*x1*x2;
dx2/dt=-x2^2*(k1-k2*x1*x2);
y1=x1;
x1(0)=0;
