B = (-2, -2.12)
C = (2, -2.12)
D = (4, 1)
A = (0, 1)
v1, a, b, c, d = Polygon(B, C, D, A)
d1, c_1, a_1, b_1 = Polygon(A, B, C)
f = OrthogonalLine(C, d)
g = OrthogonalLine(A, a)
h = OrthogonalLine(B, b_1)
H = Intersect(f, g)
i = Line(H, d)
E = Intersect(a, a_1)
P = Intersect(i, a)
Q = Intersect(c, i)
j = Line(H, a)
R = Intersect(j, d)
S = Intersect(b, j)
e = Circle(P, R, S)