Implementation of the fitting function Jacobian matrix (work in progress)