 function xprime = exode10_1_5(t,x)
     xprime = t+x+x.^2;
 end