pl(xz)

u = Vector((1, 1, 2))
v = Vector((-2, 5, 1))
n: u ⊗ v
P_{0} = Point({1, -2, 3})
P = OrthogonalPlane(P_{0}, n)
pl = OrthogonalPlane(Point[{0, 0, 0}], Vector[(0, 1, 0)])
r = IntersectPath(pl, P)