Created two unit tests showing the target_functions.relax_fit.jacobian() function is correct.
This compares the calculated Jacobian to the numerically integrated values from the
test_suite/shared_data/curve_fitting/numeric_gradient/jacobian.py script.
