Approximates the Jacobian of m functions in n unknowns using divided differences.