//Fonction par boucles pour calculer la somme pour k variant de 1 à n des 1/k
function S = s1(n)
S=1;
for i=2:n
S=S+1/i;
end
endfunction
//Fonction vectorisée pour calculer la somme pour k variant de 1 à n des 1/k
function S = s2(n)
S = sum((1:n).^(-1))';
endfunction