M = Point(xAxis)
A = Point(xAxis)
c = Circle(M, A)
C = Point(yAxis)
d = Circle(M, C)
P = Point(c)
a = Segment(M, P)
b = OrthogonalLine(P, xAxis)
B = Intersect(d, a, 1)
e = OrthogonalLine(B, yAxis)
E = Intersect(e, b)