disp('On va comparer les algorithmes s1 et s2 numériquement')
exec s1_s2.sci;    //appel des fonctions s1 et s2F

//entrée simple de la valeur de n, ou avec fenêtre de dialogue :
n = input('\nvaleur de l''indice n>100 de la somme partielle des 1/k ?\n');
//n = evstr(x_dialog('nombre de points ?','100'))F

//Première résolution
timer();
S1=s1(n);
t1=timer();
printf('\n Le résultat de la première méthode est :\n..
la somme pour k variant de 1 a %g des 1/k vaut  S1 = %g\n ',n,S1)
//voir l'aide sur printf_conversion pour "%g"
printf('Le temps mis par la première méthode est :\n..
t1 = %.4g millisecondes\n ',t1*1000)F

//Seconde résolution
timer();
S2=s2(n);
t2=timer();
printf('\n Le résultat de la seconde méthode est :\n..
la somme pour k variant de 1 a %g des 1/k vaut  S2 = %g\n ',n,S2)
printf('Le temps mis par la seconde méthode est :\n..
t2 = %.4g millisecondes\n ',t2*1000)F

//On sait que la somme pour k variant de 1 a n des 1/k est équivalente
//quand n tend vers l'infini a ln(n), qui est note log(n) par Scilab.
printf('\n Noter que  ln(%3d) = %g !\n',n,log(n))