P = (0.9759468824994504, 4.315540049700682)
B = (-0.9149115553924083, 1.6956759490071427)
h = Segment(P, B)
D = (3.459122421417675, 1.2628288367186449)
f = Segment(P, D)
g = OrthogonalLine(B, h)
i = OrthogonalLine(D, f)
j = AngularBisector(B, P, D)
C_{1} = Point(g)
k = Line(P, C_{1})
k' = Mirror(k, j)
c = Circle(P, B, C_{1})
D_{1} = Intersect(i, k')
d = Circle(P, D_{1}, D)
A = Intersect(c, d, 2)
loc1 = Locus(A, C_{1})
A' = Mirror(A, j)
n = Line(P, A')
m = Segment(P, A)
α: Angle[P, D, A] + Angle[A, B, P] - Angle[B, P, D]
p = Line(B, D)
p' = Rotate(p, pi / 2 - α, D)
q = OrthogonalLine(B, p)
E = Intersect(p', q)
e = Circle(B, E, D)
C = Intersect(e, n, 2)
l = Segment(B, A)
r = Segment(A, D)
s = Segment(P, C)
t = Segment(C, B)
a = Segment(C, D)