Problème avec les solutions renvoyé par fsolve

Bonjour à tous,

Je dois réaliser un algorithme de triangulation et résoudre pour cela un système de 3 équations à 3 inconnues de la forme: (x-a)^2+(y-b)^2+(z-c)^2 = r^2. (intersection de 3 sphères donc)
Le système n'étant pas linéaire j'utilise, pour le resoudre, la fonction fsolve de scipy qui doit logiquement me renvoyer 2 couples de valeurs pour x,y,z.
Cependant, l'algorithme ne me renvoie qu'un seul couple de valeur. En comparant avec Matlab (qui me renvoie bien 2 solutions), je me suis rendu compte qu'en fait fsolve réalise la moyenne des solutions sur x,y,z.
Cela m'ammène donc à vous poser la question suivante: Comment faire pour que fsolve me sorte bien les 2 couples de valeur dont j'ai besoin ?

Merci par avance :D