A_1 = (-4, 5)
B_1 = (-0.28, 4.98)
s = Segment(A_1, B_1)
C = (-4, 4)
D = (-2, 4)
p = Segment(C, D)
O = Point(yAxis)
c = Circle(O, Segment[A_1, B_1])
G, S = Intersect(c, xAxis)
a = OrthogonalLine(S, xAxis)
d = Circle(S, Segment[C, D])
E, B = Intersect(d, a)
A = Point(yAxis)
b = Segment(A, B)
M = Midpoint(A, B)
e = Circle(M, B)
X_1, X_2 = Intersect(e, xAxis)
f = Circle(O, Segment[C, D])
P, I = Intersect(f, yAxis)