function y=iterazione(phi,x0,k)
    %esegue k passi di iterazione di punto fisso su phi

    y=zeros(k,1);
    y(1)=phi(x0);
    for i=2:k
        y(i)=phi(y(i-1));
    end
end