
function somma=seriebernoulli2(n)

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

if n <= 0
    somma=0; return;
end

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

argomento_somma=numeratore./denominatore;
somma=0;
for k=1:length(argomento_somma)
    t=argomento_somma(k);
    somma=somma+t;
end
