This computation of 
does not explain how to compute a derivative
approximation. But if the point b is extended to a series of
uniformly spaced points
each a multiple of h from a and we write the TS expansion as an
expression, we get
Assume the difference operator is centered at the origin. Now formulating this as a matrix equation using terms up to the second derivative, we have
