P: (1, -2, 1)
Q: (-4, 2, -1)
f = Segment(P, Q)
M = Midpoint(P, Q)
a = Vector((0, 0, 0), M)
n = Vector(M, Q)
E: (((x, y, z) - a) * n) = 0