Computing the Jacobian matrix, using (at most) Numpy

Question

I\'m attempting to write a simple implementation of the Newton-Raphson method, in Python. I\'ve already done so using the SymPy library, however wh