<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def jac(x,a):
 """ Jacobian matrix given Christophe's suggestion of f """

    return (x-a) / np.sqrt(((x-a)**2).sum(1))[:,np.newaxis]