
function R=calcola_R(xn,L)

% Calcolo dell'ordine di convergenza di una sequenza "xn" convergente a "L"
% con lim_n E_{n+1}/E^R_n=C, C non uguale a 0.

metodo=2;

switch metodo
    
    case 1 % vettoriale
        
        En=abs(xn-L);
        gamma_n=En(1:end-1)./En(2:end);
        R=log(gamma_n(2:end))./log(gamma_n(1:end-1));
        
    case 2 % ciclo for
        
        En=abs(xn-L);
        N=length(xn);
        
        % calcolo gamma_n
        for i=1:N-1
            gamma_n(i)=En(i)/En(i+1);
        end
        
        % calcolo R
        for j=1:N-2
            R(j)=log(gamma_n(j+1))./log(gamma_n(j));
        end
        
end

