(* *) (* Numerical Mathematics and Computing, Fifth Edition *) (* Ward Cheney & David Kincaid *) (* Brooks/Cole Publ. Co. *) (* (c) 2003 *) (* No warranties implied or expressed. *) (* File: bvp1 *) NDSolve[{x''[t] == Exp[t] - 3*Sin[t] + x'[t] - x[t], x[1] == 1.09737491, x[2] == 8.63749661}, x, {t, 1.0, 2.0}, WorkingPrecision->15] Table[N[{{t, x[t]} /. %},9], {t, 1., 2., 1./11.}] TableForm[%]