dsolve(z'=-(2*z*t+t^3)