sol : rk([yp,6.0*y*y + t],[y,yp],[1.0,0.0],[t,0.0,1.1,0.01]);