PopravekSmeri_C4

a: 10
e = 0.6 (slider from 0 to 1, step 0.05)
fi1 = 0.209439510239319 (slider from ((-pi)) / 2 to pi / 2, step pi / 30)
la1 = 0.10471975511966 (slider from (-pi) to pi, step pi / 30)
h1 = 1.7 (slider from 0 to 8, step 0.1)
fi2 = 1.0471975511966 (slider from ((-pi)) / 2 to pi / 2, step pi / 30)
la2 = 1.25663706143592 (slider from (-pi) to pi, step pi / 30)
h2 = 3.5 (slider from 0 to 8)
rN1: a / sqrt(1 - (e^(2) * (sin(fi1))^(2)))
rN2: a / sqrt(1 - (e^(2) * (sin(fi2))^(2)))
b: (a * sqrt(1 - e^(2)))
EL: (x / a)^(2) + (y / a)^(2) + (z / b)^(2) = 1
EQ = CurveCartesian((a * cos(t)), (a * sin(t)), 0, t, 0, (2 * pi))
MX = CurveCartesian((a * cos(t)), 0, (b * sin(t)), t, 0, (2 * pi))
LM1 = CurveCartesian(((a * cos(t)) * cos(la1)), ((a * cos(t)) * sin(la1)), (b * sin(t)), t, 0, (2 * pi))
LM2 = CurveCartesian(((a * cos(t)) * cos(la2)), ((a * cos(t)) * sin(la2)), (b * sin(t)), t, 0, (2 * pi))
T1: ((((rN1 + h1) * cos(fi1)) * cos(la1)), (((rN1 + h1) * cos(fi1)) * sin(la1)), (((rN1 * (1 - e^(2))) + h1) * sin(fi1)))
T2: ((((rN2 + h2) * cos(fi2)) * cos(la2)), (((rN2 + h2) * cos(fi2)) * sin(la2)), (((rN2 * (1 - e^(2))) + h2) * sin(fi2)))
eh1: ((cos(fi1) * cos(la1)), (cos(fi1) * sin(la1)), sin(fi1))
eh2: ((cos(fi2) * cos(la2)), (cos(fi2) * sin(la2)), sin(fi2))
pNor1 = Segment(T1 - ((rN1 + h1 + 2) * eh1), T1 + (2 * eh1))
pNor2 = Segment(T2 - ((rN2 + h2 + 2) * eh2), T2 + (2 * eh2))
T2E: T2 - (h2 * eh2)
T1E: T1 - (h1 * eh1)
NP1 = Plane(T1, T1E, T2)
NP1E = Plane(T1, T1E, T2E)