Functions for creating Gradient and Hessian matrices by numerical differentiation (Richardson's method) of the partial derivatives
Useful links