t = 0.5235987755982988 (slider from -1.5707963267948966 to 1.5707963267948966, step π / 6)
r: 3
β = 0.8 (slider from 0 to 1, step 0.1)
a = CurveCartesian(((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t)), s, 0, (2 * π))
s: ((2 * π) * β)
A: (((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t)))
gam_s = Vector(A, A + ((((-r) * cos(t)) * sin(s)), ((r * cos(t)) * cos(s)), 0))
gam_t = Vector(A, A + ((((-r) * sin(t)) * cos(s)), (((-r) * sin(t)) * sin(s)), (r * cos(t))))
b = Surface(((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t)), s, 0, (2 * π), t, ((-π)) / 2, π / 2)
B: (((r * cos((-t))) * cos(s)), ((r * cos((-t))) * sin(s)), (r * sin((-t))))
gm_s = Vector(B, B + ((((-r) * cos((-t))) * sin(s)), ((r * cos((-t))) * cos(s)), 0))
gm_t = Vector(B, B + ((((-r) * sin((-t))) * cos(s)), (((-r) * sin((-t))) * sin(s)), (r * cos((-t)))))
c = CurveCartesian(((r * cos((-t))) * cos(s)), ((r * cos((-t))) * sin(s)), (r * sin((-t))), s, 0, (2 * π))
ae = CurveCartesian((r * cos(s)), (r * sin(s)), (r * 0), s, 0, (2 * π))
la = CurveCartesian(((r * cos(t)) * cos(((2 * π) * β))), ((r * cos(t)) * sin(((2 * π) * β))), (r * sin(t)), t, ((-π)) / 2, π / 2)
δ: 0.2
C: (((r * cos(t)) * cos(((2 * π) * δ))), ((r * cos(t)) * sin(((2 * π) * δ))), (r * sin(t)))
rho: 2
n = Vector(C, C + ((cos(t) * cos(((2 * π) * δ))), (cos(t) * sin(((2 * π) * δ))), sin(t)))
v = Vector(C, C + ((((rho * rho) * cos(t)) * cos(((2 * π) * δ))), (cos(t) * sin(((2 * π) * δ))), 0))
f = Segment((0, 0, 4), (0, 0, -4))
D = (-0.026164874358128025, 0.7915289054042907)
E = (0.13631708737288428, 0.7915289054042907)