r = 3
rho: 2
c_1 = 0
l1 = Sequence(Sequence[Vector[(((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t))), ((((r + 1) * cos(t)) * cos(s)), (((r + 1) * cos(t)) * sin(s)), ((r + 1) * sin(t)))], s, 0, (2 * π), π / 7], t, ((-π)) / 2, π / 2, π / 6)
l2 = Sequence(Sequence[Vector[(((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t))), ((((r + rho) * cos(t)) * cos(s)), (((r + rho) * cos(t)) * sin(s)), (r * sin(t)))], s, 0, (2 * π), π / 7], t, ((-π)) / 2, π / 2, π / 6)
a = Surface(((r * cos(t)) * cos(s)), ((r * cos(t)) * sin(s)), (r * sin(t)), t, 0, (2 * π), s, ((-π)) / 2, π / 2)
c_{1} = 0
A = (-0.025400253361746797, 0.8083505673246775)
B = (0.13708170836926542, 0.8083505673246775)