function r =sol2(x1,x2,t)
  r = x1*exp(-x2*t);