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