P = (-3, 1)
b = OrthogonalLine(P, r)
c = Line(P, r)