Как выполнить нелинейную оптимизацию с помощью scipy / numpy или sympy?

Я пытаюсь найти оптимальное решение для следующей системы уравнений в Python:

(x-x1)^2 + (y-y1)^2 - r1^2 = 0
(x-x2)^2 + (y-y2)^2 - r2^2 = 0
(x-x3)^2 + (y-y3)^2 - r3^2 = 0

Даны значения точки (x, y) и радиуса (r):

x1, y1, r1 = (0, 0, 0.88)
x2, y2, r2 = (2, 0, 1)
x3, y3, r3 = (0, 2, 0.75)

Как лучше всего найти оптимальное решение для точки (x, y) Используя приведенный выше пример, это будет:

~ (1, 1)

Ответы на вопрос(4)

Ваш ответ на вопрос