% Inizializzo x1 e N
x1 = 1; N = 70;

% Stampo a schermo i valori
fprintf('\n \t x1 : %4.0f \t N : %4.0f \n', x1, N);

% Calcolo la sequenza fino al valore N
sequenza = sequenza_log(x1,N);

% Stampo a schermo l'ultima iterata
fprintf('\n \t Ultima iterata: %1.3e \n', sequenza(end));

% Calcolo il valore esatto
x_exact = log(x1+1);

% Vettore degli errori assoluti
err_ass = abs(sequenza-x_exact);

% Grafico
figure(1)
semilogy(1:N, err_ass,'ko-');
hold on;

% Sequenza senza cancellazione numerica
sequenza2 = sequenza_log2(x1,N);

% Vettore degli errori assoluti
err_ass2 = abs(sequenza2-x_exact);

semilogy(1:N,err_ass2,'ro-');