Newton-Raphson implementation, SymPy linsolve producing empty set?

Question

I wrote am implementation of the Newton-Raphson method (a multivariate generalization of Newton\'s method) to solve the system:

sqrt( (x - x_i) ** 2 + (y - y_