P = (2.2357107570376065, 1.5276353777825877, 2)
O = (0, 0, 0)
P'_1 = Mirror(P, O)
a_1, c, b_1 = Cone(P'_1, P, 2)
ellipse = IntersectPath(se, a_1)
e = OrthogonalLine(P, se)
P' = Intersect(e, se)
u = Vector(P, P')
Hauptachse = Line(P', O)
A, B = Intersect(Hauptachse, ellipse)
M: (A + B) / 2
Nebenachse = OrthogonalLine(M, Hauptachse, xOyPlane)
C, D = Intersect(Nebenachse, ellipse)
a = Segment(M, B)
b = Segment(M, C)