Problème de Sangaku
Author: Yvan Monka - m@ths et tiques 2005

A = (-3.2666666666666666, 2.7333333333333334)
B = (5.933333333333334, 5.266666666666667)
C = (10.3, -3.8666666666666667)
c = Circle(A, B, C)
D = Point(c)
P, a, b, c_1, d = Polygon(A, B, C, D)
e = Segment(A, C)
f = Segment(D, B)
g = AngularBisector(B, A, C)
h = AngularBisector(B, C, A)
H = Intersect(g, h)
i = OrthogonalLine(H, b)
F = Intersect(b, i)
k = Circle(H, F)
j = AngularBisector(C, A, D)
l = AngularBisector(D, C, A)
J = Intersect(j, l)
m = OrthogonalLine(J, e)
H_1 = Intersect(m, e)
p = Circle(J, H_1)
n = AngularBisector(B, D, C)
q = AngularBisector(D, B, C)
I = Intersect(n, q)
r = OrthogonalLine(I, c_1)
J_1 = Intersect(r, c_1)
s = Circle(I, J_1)
t = AngularBisector(A, B, D)
a_1 = AngularBisector(A, D, B)
K = Intersect(a_1, t)
b_1 = OrthogonalLine(K, d)
L = Intersect(b_1, d)
d_1 = Circle(K, L)
Q, k_1, g_1, i_1, e_1 = Polygon(K, J, I, H)