Numerical Differentiation using Mathematica: Redux

While I was demonstrating this Mathematica script in class, I noticed there was an obvious error that had escaped into the wild. In particular, when I tried to get the center-difference approximation to the second derivative via AppDeriv[3, 2], I got the following:

Notice the "zero" truncation error for the second formula. That should have tipped me off (if I had done any testing, duh). The other errors are also incorrect. It can be fixed relatively easily by the following code.

AppDeriv[m_, n_] := Module[{points, IP, nthDeriv, e1},

  points = Table[{x0 + i h, Subscript[f, i]}, {i, 0, m - 1}];

  IP = InterpolatingPolynomial[points, x];

  nthDeriv = D[IP, {x, n}];

  e1 = D[F[X[x]]* Apply[Times,

            (x - Table[{x0 + i h}, {i, 0, m - 1}])], {x, n}];

  Formula = TableForm[

    Table[{Subsuperscript[f, i, n], nthDeriv /. x -> (x0 + i h), 

       O Superscript[h, Exponent[e1[[1]] /. 

       x -> (x0 + i h), h]]}, {i,0, m - 1}] // Simplify ]]

Here is my Mathematica screenshot with the new code. It at least seems to fix the problem above.