from sympy import * x, y, z = symbols('x y z', real=True)
A = Point(-3,0,0) B = Point(3,0,0) H = Point(0,y,0) P = Point(0,y,4)
solve(Eq(P.distance(A),6), y)
[-sqrt(11), sqrt(11)]
Line(A,B).distance(H).subs(y, sqrt(11))