
function demoI

% definizione della funzione
f=@(x) exp(-x.^2);
a=0;
b=10000;

% ciclo-for
for k=1:16
    N(k)=2^k;
    S(k) = simpson (f,a,b,N(k));
    fprintf('\n \t N(k): %5.0f S(k): %1.15e',N(k),S(k));
end

% valore corretto dell'integrale
I=(sqrt(pi)/2)*erf(b);

fprintf('\n \n \t I: %1.15e \n \n',I)

% errore assoluto
err=abs(S-I);

% grafico in scala semilogaritmica
semilogy(N,err,'k-');

% salvataggio su file
fid=fopen('errori.txt','w');
for k=1:16
    fprintf(fid,'\n \t | %5.0f | %1.5e |',N(k),err(k));
end
fclose(fid);
