Added code for finite difference matrices generated by the Fornberg algorithm

Added the functions fdweights and fdmatrix. fdweights generates finite difference weights using the Fornberg algorithm. fdmatrix calls fdweights to construct central difference differentiation matrices using stencils no wider than 9 points.