eq1: (13 * x) + (4 * y) - (7 * z) - 45 = 0
P: (-5, 3, 8)
Vector[(3, -1, 5)] ⊗ Vector[(2, -3, 2)]
p = OrthogonalPlane((2, 3, -1), Vector[(3, -1, 5)] ⊗ Vector[(2, -3, 2)])
f = OrthogonalLine(P, eq1)
A = Intersect(f, eq1)