% si "ripulisce" il workspace e lo schermo
clear;
clc;

NP = 5;

ek = 0:NP;

nk = 10.^ek;

for ik = 1:length(nk)
    pi_app(ik) = 4.*partsum_arctg( 1, nk(ik) );
end

Ek = abs(pi_app - (pi));

Erk = Ek/pi;


barra = '----------------------------------------------------------------';
titolo = 'Tabella degli errori nel calcolo approssimato di pi greco';

disp(titolo);
disp(barra);
fprintf('  ek  \t     n  \t  pi \t\t |E_k| \t\t E_r \n');
disp(barra);
fprintf('  %d \t %7d \t %8.7g \t %4.3e  \t %4.3g     \n',...
    [ek; nk; pi_app; Ek; Erk]);

