f = Function(x^(2.0), -2.0, 2.0)
O = (0, 0)
I = (1, 0)
J = (0, 1)
K = (1, 1)
poly1, o, i, k, j = Polygon(O, I, K, J)
Ω = (0.5, 0.5)
c = Circle(O, Ω)
L2, L1 = Intersect(c, xAxis)
a = OrthogonalLine(L1, xAxis)
M = Intersect(f, a)
b = Segment(M, L1)
d = OrthogonalLine(L2, xAxis)
N = Intersect(f, d)
e = Segment(N, L2)