IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

 C Discussion :

comment calculer la somme d'une série


Sujet :

C

  1. #1
    Membre à l'essai
    Inscrit en
    Novembre 2009
    Messages
    16
    Détails du profil
    Informations forums :
    Inscription : Novembre 2009
    Messages : 16
    Points : 12
    Points
    12
    Par défaut comment calculer la somme d'une série
    comment je peut calculer la somme de cette série:
    1-1/2²+1/3²-1/4²+...+-1/n²+-...
    avec précision :
    epsilon=10e-8(10e-10,10e-12,10e-14), utilisant l'approximation suivante;
    lim k-->inf Xk=Xn, où le critère du fin est |Xn-Xn+1|<= epsilon
    la valeur n n'est pas connue en avance, sa valeur est donné du epsilon.
    la somme partielle suivante doit étre calculé du somme précedente, (not recalculated from the degining)

    the sum calculated analytacally is equal to pi2/12.

    j'ai essayé de résoudre cette exercice mais le programme que j'ai écrit ça marche pas, est ce vous avez une idée SVP
    merci

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    #include <stdio.h>
    #include<math.h>
     
    int main()
    {
    int n,i;
    double xprev,xnew,s,eps;
     
    eps=1e-8;
    s=0;
    for(i=1;;i=i+1){
    	xnew=(-1)pow(n+1)*1/n*n;
    	if (abs(xprev-xnew)<=eps)
    		break;
    	s=s+xnew;
     
    }
    n=i-1;
    	printf("s=");
    	printf("n=");
     
     
    return 0;
    }

  2. #2
    Invité(e)
    Invité(e)
    Par défaut
    Bonjour,

    je pense que l'expression de xnew est à revoir :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    xnew=(-1)pow(n+1)*1/n*n;
    • C'est quoi pow ?
    • Il faut faire attention quand on utilise des entier dans un calcul (ex double x = 3/5 donnera x = 0.)

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    if(n%2) {
       /* n impair */
       xnew = 1./(double)(n*n);
    } else {
       /* pair */
       xnew = -1./(double)(n*n);
    }
    Et n'oublie pas d'assigner xprev.

    EDIT :

    Attention abs ne s'applique qu'aux entiers, il faut utiliser fabs pour les flottant.

    Et plutôt qu'une boucle for, on peut utiliser une boucle do while qui évite d'utiliser i:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    n = 1;
    xnew = 0;
    do {  
       /* mise à jour xprev*/
       xprev = xnew
       /* calcul de xnew*/
       ...
     
      /* calcul du prochain n */
      ++n;
    }while(fabs(xprev-xnew)>eps)
    Dernière modification par Invité(e) ; 08/03/2010 à 15h21.

  3. #3
    Membre à l'essai
    Inscrit en
    Novembre 2009
    Messages
    16
    Détails du profil
    Informations forums :
    Inscription : Novembre 2009
    Messages : 16
    Points : 12
    Points
    12
    Par défaut re
    Citation Envoyé par mabu Voir le message
    Bonjour,

    je pense que l'expression de xnew est à revoir :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    xnew=(-1)pow(n+1)*1/n*n;
    • C'est quoi pow ?
    • Il faut faire attention quand on utilise des entier dans un calcul (ex double x = 3/5 donnera x = 0.)

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    if(n%2) {
       /* n impair */
       xnew = 1./(double)(n*n);
    } else {
       /* pair */
       xnew = -1./(double)(n*n);
    }
    Et n'oublie pas d'assigner xprev.

    EDIT :

    Attention abs ne s'applique qu'aux entiers, il faut utiliser fabs pour les flottant.

    Et plutôt qu'une boucle for, on peut utiliser une boucle do while qui évite d'utiliser i:

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    n = 1;
    xnew = 0;
    do {  
       /* mise à jour xprev*/
       xprev = xnew
       /* calcul de xnew*/
       ...
     
      /* calcul du prochain n */
      ++n;
    }while(fabs(xprev-xnew)>eps)
    pour le reste c'est juste, est ce que vous pouvez me donner un programme correcte pour ce task

  4. #4
    Invité(e)
    Invité(e)
    Par défaut
    Citation Envoyé par osabri Voir le message
    pour le reste c'est juste, est ce que vous pouvez me donner un programme correcte pour ce task
    Voici une boucle de calcul correcte, mais pour l'intégration, je te laisse faire.

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
        s = 0;
        n = 1;
        xnew = 0;
        do {
            /* mise à jour xprev */
            xprev = xnew;
            /* calcul de xnew */
            if (n % 2) {
                /* n impair */
                xnew = 1. / (double) (n * n);
            } else {
                /* pair */
                xnew = -1. / (double) (n * n);
            }
            /* ajout a la somme */
            s += xnew;
            /* calcul du prochain n */
            ++n;
        } while (fabs(xprev - xnew) > eps);

  5. #5
    Membre à l'essai
    Inscrit en
    Novembre 2009
    Messages
    16
    Détails du profil
    Informations forums :
    Inscription : Novembre 2009
    Messages : 16
    Points : 12
    Points
    12
    Par défaut re
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    #include <stdio.h>
    #include<math.h>
     
      int main()
    {
        int n;
     
        double xprev,xnew,s,eps;
     
        s = 0;
        n = 1;
        xnew = 0;
        do {
     
            xprev = xnew;
            /* calculating xnew */
            if (n % 2) {
                /* n odd */
                xnew = 1. / (double) (n * n);
            } else {
                /* even */
                xnew = -1. / (double) (n * n);
            }
            /* adding the sum */
            s += xnew;
            /* calculating next n */
            ++n;
        } while (fabs(xprev - xnew) > eps);
     
    	printf("s= %lf\n",s);
    	printf("n= %d\n",n);
    return 0;
    }

    mon programme affiche rien ( j'ai besoin de ça cette semaine

  6. #6
    Invité(e)
    Invité(e)
    Par défaut
    Il ne faut pas ignorer les warning à la compilation.
    SI ton compilateur ne t'en donne pas, configure le correctement : http://emmanuel-delahaye.developpez....-codage-c/#LVI
    Code gcc : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    10-03-09-1.c: In function 'main':
    10-03-09-1.c:28: warning: 'eps' is used uninitialized in this function

    N'oublie pas de donner une valeur à eps.

    EDIT :
    Et les balises code (bouton #) font toujours joli.
    Dernière modification par Deepin ; 25/07/2011 à 10h29.

  7. #7
    Membre à l'essai
    Inscrit en
    Novembre 2009
    Messages
    16
    Détails du profil
    Informations forums :
    Inscription : Novembre 2009
    Messages : 16
    Points : 12
    Points
    12
    Par défaut je peut pas changer car c'est la console du laboratoire
    s'il vous plait une correction de mon programme

  8. #8
    Expert éminent sénior
    Avatar de diogene
    Homme Profil pro
    Enseignant Chercheur
    Inscrit en
    Juin 2005
    Messages
    5 761
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Essonne (Île de France)

    Informations professionnelles :
    Activité : Enseignant Chercheur
    Secteur : Enseignement

    Informations forums :
    Inscription : Juin 2005
    Messages : 5 761
    Points : 13 926
    Points
    13 926
    Par défaut
    Ton programme marche à condition de donner une valeur à eps comme déjà signalé par mabu.

  9. #9
    Membre à l'essai
    Inscrit en
    Novembre 2009
    Messages
    16
    Détails du profil
    Informations forums :
    Inscription : Novembre 2009
    Messages : 16
    Points : 12
    Points
    12
    Par défaut re
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    #include <stdio.h>
    #include<math.h>
    #define eps 10e-8
     
    int main()
    {
        int n;
     
        double xprev,xnew,s,eps;
     
        s = 0;
        n = 1;
        xnew = 0;
        do {
     
            xprev = xnew;
            /* calculating xnew */
            if (n % 2) {
                /* n odd */
                xnew = 1. / (double) (n * n);
            } else {
                /* even */
                xnew = -1. / (double) (n * n);
            }
            /* adding the sum */
            s += xnew;
            /* calculating next n */
            ++n;
        } while (fabs(xprev - xnew) > eps);
     
    	printf("s= %lf\n",s);
    	printf("n= %d\n",n);
    return 0;
    }
    s'il vous plait corriger mes fautes pour que mon programme marche : ((

  10. #10
    Expert éminent sénior
    Avatar de diogene
    Homme Profil pro
    Enseignant Chercheur
    Inscrit en
    Juin 2005
    Messages
    5 761
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Essonne (Île de France)

    Informations professionnelles :
    Activité : Enseignant Chercheur
    Secteur : Enseignement

    Informations forums :
    Inscription : Juin 2005
    Messages : 5 761
    Points : 13 926
    Points
    13 926
    Par défaut
    On a déjà répondu deux fois à ta question

    Est-ce que tu lis les réponses ?

  11. #11
    Invité(e)
    Invité(e)
    Par défaut
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    #define eps 10e-8
    ...
    double eps;
    eps ne peut pas être à la fois une macro et une variable...

Discussions similaires

  1. Comment calculer la taille d'une base de données ?
    Par say dans le forum Décisions SGBD
    Réponses: 8
    Dernier message: 01/04/2011, 16h48
  2. comment calculer la somme d'une colonne d'un dbgrid
    Par kirty dans le forum Bases de données
    Réponses: 4
    Dernier message: 24/07/2006, 20h55
  3. [Matrices] Comment calculer le Déterminant d'une matrice 4x4
    Par cyber_N dans le forum Algorithmes et structures de données
    Réponses: 70
    Dernier message: 19/08/2005, 15h47
  4. [CR 8.5] Calculer la somme d'une somme
    Par Frederic Vincent dans le forum Formules
    Réponses: 4
    Dernier message: 12/02/2004, 17h53

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo