% comparing linear to quadratic convergence 
% of iterative methods:
epsilon=0.5;
e_lin=0.5;
e_quad=0.5
e_secant=0.5
fprintf(' iter       linear                secant        quadratic \n')
for i=1:10
    fprintf(' %2d        %14.12f          %14.12f        %14.12f\n' ...
        ,i,e_lin,e_secant,e_quad);
    pause
    e_lin=e_lin*epsilon;
    e_quad=e_quad^2;
    e_secant=e_secant^1.6;
end
