function y=f(t,x)

y = 0.03134*x - 1.5887e-10*x^2;
