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