Exercice n°51 page 159
Le théorème de Van Aubel

A = (6.94, 0.96)
B = (9.66, 0.76)
a = Segment(A, B)
C = (10.36, -1.28)
b = Segment(C, B)
D = (6.08, -0.96)
c = Segment(C, D)
d = Segment(D, A)
e = OrthogonalLine(A, a)
f = OrthogonalLine(B, a)
g = OrthogonalLine(A, a)
h = Circle(A, B)
k = Circle(B, A)
E = Intersect(k, f, 2)
F = Intersect(h, e, 2)
poly1, f_1, e_1, b_1, a_1 = Polygon(F, E, B, A)
i = OrthogonalLine(B, b)
j = OrthogonalLine(C, b)
p = Circle(C, B)
q = Circle(B, C)
G = Intersect(p, j, 1)
H = Intersect(q, i, 1)
poly2, b_2, c_1, g_1, h_1 = Polygon(B, C, G, H)
l = OrthogonalLine(A, d)
m = OrthogonalLine(D, d)
r = Circle(D, A)
s = Circle(A, D)
I = Intersect(r, m, 2)
J = Intersect(s, l, 2)
poly3, j_1, a_2, d_1, i_1 = Polygon(J, A, D, I)
n = OrthogonalLine(D, c)
t = OrthogonalLine(C, c)
k_1 = Circle(D, C)
p_1 = Circle(C, D)
K = Intersect(k_1, n, 2)
L = Intersect(p_1, t, 2)
poly4, k_2, d_2, c_2, l_1 = Polygon(K, D, C, L)
poly5, d_3, a_3, b_3, c_3 = Polygon(D, A, B, C)