O1 = Intersect(xAxis, yAxis)
A = Point(xAxis)
c1 = Circle(O1, A)
P1 = Point(c1)
C3 = Point(c1)
H1', H1 = Intersect(c1, xAxis)
H2', H2 = Intersect(c1, Circle[H1, Distance[O1, A]])
h1 = Ray(H2, O1)
D2, D1 = Intersect(c1, yAxis)
p1 = Line(P1, D2)
G2', G2 = Intersect(c1, LineBisector[H2, D1])
g2 = Ray(G2, O1)
G1, G1' = Intersect(c1, LineBisector[D1, H1])
g1 = Ray(G1, D2)
G0 = Intersect(g2, g1)
h2 = Ray(H1, D2)
H0 = Intersect(h2, h1)
C2, C4 = Intersect(c1, Circle[C3, (Distance[D1, C3] * 2)])
n2 = Ray(C2, O1)
C1, C1' = Intersect(c1, Circle[C2, (Distance[D1, C3] * 2)])
n1 = Ray(C1, D2)
C0 = Intersect(n1, n2)
O2 = Intersect(LineBisector[C0, G0], LineBisector[G0, H0])
c2 = Circle(O2, Distance[O2, H0])
Q4 = Mirror(P1, yAxis)
P0, P0' = Intersect(c2, p1)
p2 = Ray(P0, O1)
P2', Q2_{1} = Intersect(c1, p2)
Q3', Q3 = Intersect(c1, LineBisector[Q2_{1}, Q4])
Q2', Q2 = Intersect(c1, LineBisector[P1, Q3])
α = Angle(Vector[O1, Q4], Vector[O1, P1])
f = Segment(O1, P1)
g = Segment(O1, Q2)
h = Segment(O1, Q3)
i = Segment(O1, Q4)
arc1 = If(α ≤ 180°, Arc[c1, Q2, P1])
arc2 = If(α ≤ 180°, Arc[c1, Q3, Q2])
arc3 = If(α ≤ 180°, Arc[c1, Q4, Q3])
Diffarc: arc2 - arc3
α1 = If(arc1 > 0, Angle[Vector[O1, Q2], Vector[O1, P1]])
α2 = If(arc1 > 0, Angle[Vector[O1, Q3], Vector[O1, Q2]])
α3 = If(arc1 > 0, Angle[Vector[O1, Q4], Vector[O1, Q3]])
Diffangle: α2 - α3
β = Angle(Vector[O1, C4], Vector[O1, C1])