f(x) = If[(x ≥ -2) ∧ (x ≤ 2),x^(2) - (4 * x) + 1]
A1 = -2
C1 = 2
B1: (A1 + C1) / 2
A2 = If((f(A1) * f(B1)) > 0, B1, A1)
C2 = If((f(B1) * f(C1)) > 0, B1, C1)
B2: (A2 + C2) / 2
A3 = If((f(A2) * f(B2)) > 0, B2, A2)
C3 = If((f(B2) * f(C2)) > 0, B2, C2)
B3: (A3 + C3) / 2
A4 = If((f(A3) * f(B3)) > 0, B3, A3)
C4 = If((f(B3) * f(C3)) > 0, B3, C3)
B4: (A4 + C4) / 2
A5 = If((f(A4) * f(B4)) > 0, B4, A4)
C5 = If((f(B4) * f(C4)) > 0, B4, C4)
B5: (A5 + C5) / 2
A6 = If((f(A5) * f(B5)) > 0, B5, A5)
C6 = If((f(B5) * f(C5)) > 0, B5, C5)
B6: (A6 + C6) / 2
A7 = If((f(A6) * f(B6)) > 0, B6, A6)
C7 = If((f(B6) * f(C6)) > 0, B6, C6)
B7: (A7 + C7) / 2
A8 = If((f(A7) * f(B7)) > 0, B7, A7)
C8 = If((f(B7) * f(C7)) > 0, B7, C7)
B8: (A8 + C8) / 2
A9 = If((f(A8) * f(B8)) > 0, B8, A8)
C9 = If((f(B8) * f(C8)) > 0, B8, C8)
B9: (A9 + C9) / 2
A10 = If((f(A9) * f(B9)) > 0, B9, A9)
C10 = If((f(B9) * f(C9)) > 0, B9, C9)
B10: (A10 + C10) / 2
A11 = If((f(A10) * f(B10)) > 0, B10, A10)
C11 = If((f(B10) * f(C10)) > 0, B10, C10)
B11: (A11 + C11) / 2
L_1: {A1, A2, A3, A4, A5, A6, A7, A8, A9, A10, A11}
L_2: {B1, B2, B3, B4, B5, B6, B7, B8, B9, B10, B11}
L_3: {C1, C2, C3, C4, C5, C6, C7, C8, C9, C10, C11}
milieu = Sequence((Element[L_2, k], 0), k, 1, Length[L_2])
G = (2, -3)
B = (2, 0)
k = Segment(B, G)
F = (-2, 13)
A = Intersect(f, a, (-10.108148148148151, 3.478518518518515))
D = Intersect(a, xAxis)
E = Intersect(a, xAxis)
I = Intersect(f, a, (0.26794919243118864, -2.284838984678572E-13))
C = (-2, 0)
b = Segment(F, C)
Etape = 4 (slider from 0 to 11, step 1)
x_0 = Element(milieu, Etape + 1)
pointsA = Sequence((Element[L_1, k], 0), k, 1, Length[L_1])
pointsB = Sequence((Element[L_3, k], 0), k, 1, Length[L_3])
c = OrthogonalLine(x_0, a)
H = Intersect(f, c, (0.375, -0.359375))
d = Segment(x_0, H)
N = Element(pointsB, Etape + 1)
S = Element(pointsA, Etape + 1)
e = OrthogonalLine(S, a)
g = OrthogonalLine(N, a)
J = Intersect(e, b)
K = Intersect(e, b)
T = Intersect(f, e, (0.25, 0.0625))
h = Segment(S, T)
L = Intersect(f, g, (0.5, -0.75))
i = Segment(N, L)