u_0

u_1

u_2

r = 2 (slider from 0 to 4, step 0.01)
f = Function(r * x * abs(1 - x), 0, 1)
u_0 = 1.1 (slider from -1 to 5, step 0.01)
n = 43 (slider from 1 to 60, step 1)
A: (u_0, 0)
B: (u_0, f(u_0))
L1 = IterationList(f, u_0, n)
list1 = Sequence(Segment[(Element[L1, i], Element[L1, i + 1]), (Element[L1, i + 1], Element[L1, i + 1])], i, 1, n - 1)
list2 = Sequence(Segment[(Element[L1, i], Element[L1, i]), (Element[L1, i], Element[L1, i + 1])], i, 2, n - 1)
b = Segment(B, A)
c = OrthogonalLine(B, yAxis)
C = Intersect(c, a)
d = OrthogonalLine(C, xAxis)
D = Intersect(d, xAxis)
e = Segment(C, D)
E = Intersect(f, d, (0.22000000000000022, 0.3432000000000003))
g = OrthogonalLine(E, yAxis)
F = Intersect(g, a)
h = OrthogonalLine(F, xAxis)
G = Intersect(h, xAxis)
i = Segment(F, G)