eq1: x + (2 * y) + z = 1
f = Line((2, 3, 7), Vector[(3, 1, 3)])
A = Intersect(f, eq1)