//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