n = Vector((0, -1, 0))
P = Point({3, -1, 2})
A = OrthogonalPlane(P, n)