E = (1, -4)
D = Point(Circle[E, 4])
a = Segment(E, D)
B' = Rotate(D, 20*°, E)
α = Angle(D, E, B')
b = Line(E, B')
c = Circle(E, 5)
A, L = Intersect(c, b)
d = Segment(D, L)
e = Segment(L, E)
poly1, d_1, e_1, l = Polygon(E, L, D)
I_L = Midpoint(a)
I_E = Midpoint(d)
f = OrthogonalLine(I_L, a)
g = OrthogonalLine(I_E, d)
O = Intersect(f, g)
h = Circle(O, L)