C = (3.0539854992254845, 5.029952996094672)
i = Line(C, yAxis)
D = (0, 1)
A = (0, 0)
B = Point(i)
h = AngularBisector(C, B, A)
m = AngularBisector(D, A, B)
a = Line(A, B)
b = Line(B, C)
c = Line(A, D)
d = Line(D, C)
E = Point(a)
e = Circle(B, E)
k = Circle(E, Segment[E, B])
H = Intersect(e, b, 2)
p = Circle(H, Segment[H, B])
I = Intersect(k, p, 1)
j = AngularBisector(C, D, A)
l = AngularBisector(B, C, D)
poly1, a_1, b_1, c_1, d_1 = Polygon(A, B, C, D)
F = Intersect(j, m)
J = Intersect(l, h)
K = Intersect(j, l)
G = Intersect(m, h)
poly2, f, g, j_1, k_1 = Polygon(F, G, J, K)
L = (2, -0.5)
M = Point(yAxis)
N = Point(xAxis)
O = (3.5, 0.5)
P = (3, 7)