A: (3, -1, 1)
u = Vector((3, -2, 7))
b = OrthogonalPlane((3, -1, 1), u)