P: (6, -20, 41)
eq1: x - (4 * y) + (8 * z) - 9 = 0
f = OrthogonalLine(P, eq1)
A = Intersect(f, eq1)