S = (3.1840816777831187E-17, -0.06, 0.78)
M = PointIn(e1)
h = Segment(S, M)
b = OrthogonalPlane(M, h)
A = PointIn(b)
a, c, d = Cone(M, S, Distance[M, A])
e = IntersectPath(se, a)
a' = Mirror(a, S)
hyperbel = IntersectPath(se, d)
A' = Mirror(A, M)
α = Angle(A', S, A)
β: 90° - α / 2
a_1(x) = (tan(β) * x) + y(S)
a_2(x) = ((-tan(β)) * x) + y(S)