S = (1.7145055188062944E-17, -1.3800000000000001, 0.3999999999999999)
O = (0, 0, 0)
M = Mirror(S, O)
a_1, c, b_1 = Cone(M, S, 2)
kegelschnitt = IntersectPath(se, a_1)
e = OrthogonalLine(S, se)
S' = Intersect(e, se)
u = Vector(S, S')
Hauptachse = Line(S', O)
A, B = Intersect(Hauptachse, kegelschnitt)
M_1: (A + B) / 2
Nebenachse = OrthogonalLine(M_1, Hauptachse, xOyPlane)
C, D = Intersect(Nebenachse, kegelschnitt)
a = Segment(M_1, B)
b = Segment(M_1, C)
d = Mirror(a_1, S)
f = IntersectPath(se, d)