% Find a root of polynomial x^5 - 3x + 1 = 0 on the interval [0,1]
x = 1/2;
for n=1:10
    x = 1/3*(x^5+1)
end