dsolve(y'=t^3*y^2+2*t*y);