g1 = inline('x.^(1/3) + 2');
g2 = inline('(x-2).^3');
g3 = inline('(6+2*(x.^(1/3)))./(3 - x^(-2/3))');

x1 = 3;
x2 = x1;
x3 = x2;

strh='k       g1          g2            g3';
strl='------------------------------------';
disp(strh);
disp(strl);
for k=2:10
    x1(k) = g1(x1(k-1));
    x2(k) = g2(x2(k-1));
    x3(k) = g3(x3(k-1));
    str=sprintf('%d  %0.10g  %0.10g  %0.10g',k,x1(k),x2(k),x3(k));
    disp(str);
end
plot(1:9,x1(2:10)./x1(1:9),1:9,x3(2:10)./x3(1:9));
legend('g1','g3')