Я точно не знаю вашего ожидаемого результата, потому что я не знаю, каковы ваши переменные. Как правило, используйте sympy для определения вашей алгебры:
import sympy as smp
from sympy.solvers.polysys import solve_poly_system
# Define variables
x, x2, x3, y, y2, y3, d, d3 = smp.symbols("x x_2 x_3 y y_2 y_3 d_1 d_3")
# Define equations
eq1 = smp.Eq((x - x2)**2 + (y - y2)**2, d**2)
eq2 = smp.Eq((x - x3)**2 + (y - y3)**2, d3**2)
Затем можно решить систему уравнений с одним из sympys solvers . Я думаю, вы хотите solve_poly_system