(* *) (* Numerical Mathematics and Computing, Fifth Edition *) (* Ward Cheney & David Kincaid *) (* Brooks/Cole Publ. Co. *) (* (c) 2004 *) (* No warranties implied or expressed. *) (* File: ode2 *) (* numerical solution ODE & plot curve *) NDSolve[{x'[t] == 1 + x[t]^2 + t^3, x[1] == -4}, x, {t, 1, 2}] x[2] /. % Plot[Evaluate[x[t] /. %%], {t, 1, 2}]