A = (-2, -3, 4)
B = (3, 3, 1)
s = Segment(A, B)
M: (A + B) / 2
s_{AB} = OrthogonalPlane(M, s)
P = PointIn(s_{AB})
a = Segment(P, A)
b = Segment(P, B)