
function somma=seriebernoulli1(n)

%--------------------------------------------------------------------------
% Input:
% n: numero naturale positivo
%--------------------------------------------------------------------------
% Output:
% somma: numero reale pari a sum_{k=1}^n k^2/2^k che e' circa "6".
%--------------------------------------------------------------------------

if n <= 0
    fprintf('\n \t La somma non puo'' essere calcolata');
end

v=1:n;
numeratore=v.^2;
denominatore=2.^v;

argomento_somma=numeratore./denominatore;
somma=sum(argomento_somma);

fid=fopen('risultati.txt','w');
fprintf(fid,'\n \t %1.25f', somma);
fprintf(fid,'\n \t %1.1e', 6-somma);
s=fclose(fid);
