P = (2.0291075850330924, 1.4832508469072743, 2)
O = (0, 0, 0)
A'_1 = Mirror(P, O)
a_1, c, d, b_1 = Cylinder(P, A'_1, 2)
ellipse = IntersectPath(se, a_1)
f = OrthogonalLine(P, se)
P' = Intersect(se, f)
u = Vector(P, P')
Hauptachse = Line(O, P')
Nebenachse = OrthogonalLine(O, Hauptachse, xOyPlane)
B, A = Intersect(Hauptachse, ellipse)
D, C = Intersect(Nebenachse, ellipse)
b = Segment(O, C)
a = Segment(O, B)