'Evaluate jacobian at point x.
Parameters
x : ndarray, shape (n,)
Vector of residuals f(x).
Returns
ndarray, shape (m, n)
Jacobian matrix of `self.fun` at point `x`.'

def jac(self, x):

    raise NotImplementedError
