Bonus : (programmé en Python)
amplitude = 0.001
init : (x,t,s,s-t,n)= (1.5, -1, -1, 0, 3)
étape 1 : (x,t,s,s-t,n)= (0.5, 0.375, 0.5, 0.125, 3)
étape 2 : (x,t,s,s-t,n)= (0.5, 0.4010416666666667, 0.4166666666666667, 0.015625, 5)
étape 3 : (x,t,s,s-t,n)= (0.5, 0.4046875, 0.40729166666666666, 0.002604166666666685, 7)
étape 4 : (x,t,s,s-t,n)= (0.5, 0.40531529017857143, 0.40580357142857143, 0.00048828125, 9)
---> remarque : t et s ont été intervertis dans l'étape 4
(t, s, n) = 0.405315290179 0.405803571429 9
à comparer à ln(1.5) 0.405465108108
amplitude = 1e-05
(t, s, n) = 0.40546372234 0.405468081994 15
à comparer à ln(1.5) 0.405465108108
n = 5 :
P(x) = (x-1) - (x-1)^2 / 2 + (x-1)^3 / 3 - (x-1)^4 / 4
(x-1) > 0
majorant s : g(x) = (x-1) - (x-1)^2 / 2 + (x-1)^3 / 3
minorant t : f(x) = (x-1) - (x-1)^2 / 2 + (x-1)^3 / 3 - (x-1)^4 / 4