C_f(x) = 1.0 / x
A = (-0.5, -2)
B = (3, 0.3333333333333333)
M = Point(C_f)
poly1, b, m, a = Polygon(M, A, B)
c = OrthogonalLine(M, m)
d = OrthogonalLine(A, a)
C = Intersect(d, c)