ρ_0 = 100 (slider from 10 to 100, step 10)
u0 = 1.0471975511965976 (slider from 0 to 6.283185307179586, step 0.1)
v0 = 0.7853981633974483 (slider from 0 to 3.141592653589793, step 0.1)
A: (ρ_0 * ((cos(u0) * sin(v0)), (sin(u0) * sin(v0)), cos(v0)))
f = Segment((0, 0, 0), A)
a = 5 (slider from 0 to 5, step 0.25)
b = 4 (slider from 0 to 5, step 0.25)
c = 3 (slider from 0 to 5, step 0.25)
Σ: x^(2) / a^(2) + y^(2) / b^(2) + z^(2) / c^(2) = 1
s = 3 (slider from 1 to 5, step 1)
Δ_{inf}: (((((((((((4 * a^(4)) * b^(2)) * c^(2)) * s^(2)) - ((((4 * a^(4)) * b^(2)) * c^(2)) * s) + ((a^(4) * b^(2)) * c^(2))) * (sin(u0))^(2)) + ((((((4 * a^(2)) * b^(4)) * c^(2)) * s^(2)) - ((((4 * a^(2)) * b^(4)) * c^(2)) * s) + ((a^(2) * b^(4)) * c^(2))) * (cos(u0))^(2))) * (sin(v0))^(2)) + ((a^(4) * b^(4)) * (cos(v0))^(2))) * z^(2)) + (((((((((((8 * a^(4)) * b^(2)) * c^(2)) * s) - ((((8 * a^(4)) * b^(2)) * c^(2)) * s^(2))) * sin(u0)) * cos(v0)) * sin(v0)) * y) + (((((((((8 * a^(2)) * b^(4)) * c^(2)) * s) - ((((8 * a^(2)) * b^(4)) * c^(2)) * s^(2))) * cos(u0)) * cos(v0)) * sin(v0)) * x)) * z) + ((((((a^(4) * c^(4)) * (sin(u0))^(2)) + ((((((4 * a^(2)) * b^(2)) * c^(4)) * s^(2)) - ((((4 * a^(2)) * b^(2)) * c^(4)) * s) + ((a^(2) * b^(2)) * c^(4))) * (cos(u0))^(2))) * (sin(v0))^(2)) + ((((((4 * a^(4)) * b^(2)) * c^(2)) * s^(2)) - ((((4 * a^(4)) * b^(2)) * c^(2)) * s) + ((a^(4) * b^(2)) * c^(2))) * (cos(v0))^(2))) * y^(2)) + ((((((((((8 * a^(2)) * b^(2)) * c^(4)) * s) - ((((8 * a^(2)) * b^(2)) * c^(4)) * s^(2))) * cos(u0)) * sin(u0)) * (sin(v0))^(2)) * x) * y) + ((((((((((4 * a^(2)) * b^(2)) * c^(4)) * s^(2)) - ((((4 * a^(2)) * b^(2)) * c^(4)) * s) + ((a^(2) * b^(2)) * c^(4))) * (sin(u0))^(2)) + ((b^(4) * c^(4)) * (cos(u0))^(2))) * (sin(v0))^(2)) + ((((((4 * a^(2)) * b^(4)) * c^(2)) * s^(2)) - ((((4 * a^(2)) * b^(4)) * c^(2)) * s) + ((a^(2) * b^(4)) * c^(2))) * (cos(v0))^(2))) * x^(2)) + ((((((((-4 * a^(4)) * b^(2)) * c^(4)) * s^(2)) + ((((4 * a^(4)) * b^(2)) * c^(4)) * s) - ((a^(4) * b^(2)) * c^(4))) * (sin(u0))^(2)) + ((((((-4 * a^(2)) * b^(4)) * c^(4)) * s^(2)) + ((((4 * a^(2)) * b^(4)) * c^(4)) * s) - ((a^(2) * b^(4)) * c^(4))) * (cos(u0))^(2))) * (sin(v0))^(2)) + ((((((-4 * a^(4)) * b^(4)) * c^(2)) * s^(2)) + ((((4 * a^(4)) * b^(4)) * c^(2)) * s) - ((a^(4) * b^(4)) * c^(2))) * (cos(v0))^(2))) / ((((((((4 * a^(2)) * c^(2)) * s^(2)) - (((4 * a^(2)) * c^(2)) * s) + (a^(2) * c^(2))) * (sin(u0))^(2)) + (((((4 * b^(2)) * c^(2)) * s^(2)) - (((4 * b^(2)) * c^(2)) * s) + (b^(2) * c^(2))) * (cos(u0))^(2))) * (sin(v0))^(2)) + (((((4 * a^(2)) * b^(2)) * s^(2)) - (((4 * a^(2)) * b^(2)) * s) + (a^(2) * b^(2))) * (cos(v0))^(2))) = 0
v_1: (5 * (1, sin(u0) / cos(u0), cos(v0) / ((cos(u0) * sin(v0)))))
v_2: (5 * (1, 0, ((-(((c^(2) * cos(u0)) * sin(v0))))) / ((a^(2) * cos(v0)))))
v_{3}: (4 * (0, 1, ((-(((c^(2) * sin(u0)) * sin(v0))))) / ((b^(2) * cos(v0)))))
w_{1,s=1}: (a, 0, 0)
w_{2,s=1}: (0, b, 0)
w_{3,s=1}: (0, 0, c)
w_{1,s=2}: (10.44947105485186 * (0.3993, 0.6248, 0.671))
w_{2,s=2}: (4.831445433742327 * (-0.9105, 0.3559, 0.2104))
w_{3,s=2}: (3.56534137263603 * (-0.1074, -0.695, 0.711))
w_{1,s=3}: (17.3051 * (0.368, 0.6172, 0.6955))
w_{2,s=3}: (4.843848603854122 * (-0.9235, 0.3295, 0.1963))
w_{3,s=3}: (3.578960237245429 * (-0.108, -0.7145, 0.6912))
w_{1,s=4}: (24.1899 * (0.3607, 0.6148, 0.7013))
w_{2,s=4}: (4.846716120843427 * (-0.9264, 0.3233, 0.193))
w_{3,s=4}: (3.582344707232289 * (-0.1081, -0.7193, 0.6862))
w_{1,s=5}: (31.08239501289134 * (0.3578, 0.6139, 0.7037))
w_{2,s=5}: (4.847838601456296 * (-0.9275, 0.3209, 0.1917))
w_{3,s=5}: (3.583695752332174 * (-0.10814, -0.72123, 0.6842))
u: v_2 ⊗ v_{3}
d = OrthogonalPlane((0, 0, 0), Vector[v_2 ⊗ v_{3}])
g = IntersectPath(d, Σ)
v: w_{2,s=5} ⊗ w_{3,s=5}
e = OrthogonalPlane((0, 0, 0), Vector[w_{2,s=5} ⊗ w_{3,s=5}])