D:\OneDrive - Univerza v Ljubljani\GeoGebra\RefElipsoid_GLGinLASistem

ae = 10 (slider from 8 to 20, step 0.1)
ee = 0.7 (slider from 0 to 1, step 0.05)
b: (ae * sqrt(1 - ee^(2)))
pfi = 0.502654824574366 (slider from ((-pi)) / 2 to pi / 2, step pi / 50)
pla = 0.75398223686155 (slider from 0 to (2 * pi), step pi / 25)
ph = 1 (slider from 0 to 5, step 0.1)
pski = 0.1 (slider from -0.2 to 0.2, step 0.01)
peta = 0.1 (slider from -0.2 to 0.2, step 0.01)
rN: ae / sqrt(1 - (ee^(2) * (sin(pfi))^(2)))
T: ((((rN + ph) * cos(pfi)) * cos(pla)), (((rN + ph) * cos(pfi)) * sin(pla)), (((rN * (1 - ee^(2))) + ph) * sin(pfi)))
eX = Vector((1, 0, 0))
eY = Vector((0, 1, 0))
eZ = Vector((0, 0, 1))
eFI = Vector((((-sin(pfi)) * cos(pla)), ((-sin(pfi)) * sin(pla)), cos(pfi)))
eLA = Vector(((-sin(pla)), cos(pla), 0))
eH = Vector(((cos(pfi) * cos(pla)), (cos(pfi) * sin(pla)), sin(pfi)))
EL = Surface(((ae / sqrt(1 - (ee^(2) * (sin(ff))^(2))) * cos(ff)) * cos(ll)), ((ae / sqrt(1 - (ee^(2) * (sin(ff))^(2))) * cos(ff)) * sin(ll)), ((ae / sqrt(1 - (ee^(2) * (sin(ff))^(2))) * (1 - ee^(2))) * sin(ff)), ff, (-(pi / 2)), pi / 2, ll, 0, (2 * pi))
MX = CurveCartesian((ae * cos(t)), 0, (b * sin(t)), t, 0, (2 * pi))
EQ = CurveCartesian((ae * cos(t)), (ae * sin(t)), 0, t, 0, (2 * pi))
LM = CurveCartesian(((ae * cos(t)) * cos(pla)), ((ae * cos(t)) * sin(pla)), (b * sin(t)), t, 0, (2 * pi))
xLG = Vector(T, T + (3 * eFI))
zLG = Vector(T, T + (3 * eH))
yLG = Vector(T, T + (3 * eLA))
pNorm = Segment(T - ((rN + ph + 1) * eH), T + (2 * eH))
TLG = Vector(T, T + (0.75 * xLG) + (0.8 * yLG) + (0.5 * zLG))
zLA = Vector(T, T + (0.1 * xLG) + (0.05 * yLG) + zLG)
xLA = Vector(T, T + xLG + (0.05 * yLG) - (0.1 * zLG))
yLA = Vector(T, T - (0.1 * xLG) + yLG - (0.1 * zLG))