O = Point({0, 0, 0})
P = Point({2, 2, 0})
n = Vector((1, 2, 2))
r = Line(O, n)
p = OrthogonalPlane(P, n)
C = Intersect(r, p)
OC = Segment(O, C)