A = (1, 2)
B = (4, 2)
poly1, f, g, h, i, C, D = Polygon(A, B, 4)
P = Point(f)
j = OrthogonalLine(P, f)
c = Circle(D, Distance[A, P])
d = Circle(D, Distance[A, P])
R = Intersect(d, i, 1)
k = OrthogonalLine(R, i)
l = OrthogonalLine(P, f)
Q = Intersect(k, l)
q1, a, p, f_1, e = Polygon(A, P, Q, R)
n = Segment(R, P)
r = Line(R, P)
m = Line(C, Q)
E = Intersect(r, m)
α = Angle(Q, E, R)