I have the following concern with a nonlinear equation system in python :

```
from scipy.optimize import fsolve as fs
def nonlinear(z):
'''Solve a sistems 2x2 not linear
'''
x, y = z[0], z[1]
f1 = (x**2)+(y**2)-2*(4.41*x+2.68*y)+25.59
f2 = (x**2)+(y**2)-2*(3.23*x+2.1*y)+14.49
return [f1, f2]
```

How can I find the multiple solutions of the previous system with the code I made?

Because if I execute the previous function with initial values, I would throw only one solution:

```
x, y = fs(nonlinear, [1,1])
print(x,y)
3.38998152293 2.67210655679
```

All the solutions that the system has are the following: