Voir le sujet précédent :: Voir le sujet suivant |
Auteur |
Message |
michel
Inscrit le: 23 Juin 2006 Messages: 72
|
Posté le: 02 Mar 2007 17:27 Sujet du message: série harmonique |
|
|
Bonjour,
Peut-on majorer la série harmonique autrement que par tatonnements ?
C'est-à-dire pour un entier S donné, peut-on trouver n tel que :
1+1/2+1/3....+1/n < S < 1+1/2+1/3...1/(n+1) ?
Si oui peut-on écrire une procédure qui en fonction de S par exemple S=8 donnerait n ? Merci. |
|
Revenir en haut de page |
|
 |
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 02 Mar 2007 19:00 Sujet du message: |
|
|
Bonjour, appelons S(n)=1+1/2+1/3....+1/n. La suite (S(n)) est strictement croissante et de limite +infini.
Donc si S>0 est donné il existe n tel que S(n)<S<S(n+1)=S(n)+1/(n+1)
Code: |
harmonique:=proc(S)
> local n,s;
> n:=0; s:=0;
> while s<=S do
> n:=n+1;
> s:=evalf(s+1/n)
> end do;
> n,s
> end proc:
> harmonique(8);
1674, 8.000485592
|
Pour la première partie de la question, il existe la méthode de comparaison avec des intégrales : cf ici http://fr.wikipedia.org/wiki/Comparaison_s%C3%A9rie-int%C3%A9grale
Au revoir. |
|
Revenir en haut de page |
|
 |
michel
Inscrit le: 23 Juin 2006 Messages: 72
|
Posté le: 04 Mar 2007 11:08 Sujet du message: Nieme terme de série harmonique |
|
|
Bonjour,
Pour la rapidité d'exécution de la procédure "harmonique", je propose de remplacer en tête n:=0:s:=0 par n:=floor(evalf(exp(S-gamma))): et
s:=evalf(sum(1/k,k=1..n)):
Qu'en pensez-vous ? |
|
Revenir en haut de page |
|
 |
ALS
Inscrit le: 11 Sep 2006 Messages: 647
|
Posté le: 04 Mar 2007 12:07 Sujet du message: |
|
|
Bonjour,
oui d'accord, on peut aussi utiliser l'option remember dans une procédure récursive, ce qui accélérera les calculs.
Bon dimanche. |
|
Revenir en haut de page |
|
 |
|