Discussion :

[emprex]

Bonjour,
Je débute en C et je me trouve face à une petite diffculté. Voila je dois faire un code qui calcule une expression avec des valeurs assez grosses. Le résultat de cette expression est 1. Lorque je code avec la fonction pow(x,y),
avec double sous visual c++ 6 le résultat vaut 2. Mais lorsque je le code sous dev c++ avec int le reultat est 1. Peut on m'expliquer.

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
/**************************************************/
/* programme qui calcule 9x^4-y^4+2y^2                                 */
/* avec x=10864 et y=18817                                                   */
/**************************************************/
 
#include <stdio.h>
#include <math.h>
 
 
 
int resultat,x ,y;
 
int main(void)
{
 
	x=10864;
	y=18817;
    resultat=9*pow(x,4)-pow(y,4)+2*pow(y,2);
	printf("Si x vaut %d et si y vaut %d alors le resultat de l'equation 9x^4-y^4+2y^2 vaut %d\n",x,y,resultat);
    return 0;
}

[souviron34]

man pow...

pow est une fonction qui prend des arguments DOUBLE et retourne un DOUBLE.

Et le format pour imprimer des doubles est "%lf" ou "%g"...

^^

[nicolas.sitbon]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
#include <stdio.h>
#include <math.h>
 
int main (void)
{
   long double resultat;
   long double x = 10864;
   long double y = 18817;
   resultat = 9 * powl (x, 4.0L) - powl (y, 4.0L) + 2 * powl (y, 2.0L);
   printf ("Soit x = %Lg et y = %Lg alors\n"
           "9x^4-y^4+2y^2 = %Lg\n", x, y, resultat);
   return 0;
}

[emprex]

Merci pour vos solutions. Mais lorsque j'éxecute le code de Nicolas le résultat est 2 or je dois trouver 1. Mon problème est la.

[nicolas.sitbon]

Voila j'ai corrigé mon erreur.
Cordialement.

[ToTo13]

Bonjour,

petite question toute bète : est ce que pow(x, 4.0) n'est pas plus grand que la valeur maximale qu'un double puisse contenir ?

[aoyou]

Tu ne serais pas en train de tester les polynômes qui font vraiment mal, ceux qui mettent particulièrement en évidence la perte de précision dans les calculs. Apparemment, ici ça passe avec un long double mais si tu en as encore du même genre, tu peux regarder la librairie GMP http://gmplib.org/.

[aoyou]

Bonjour,

petite question toute bète : est ce que pow(x, 4.0) n'est pas plus grand que la valeur maximale qu'un double puisse contenir ?

Oh non encore heureux, voici la plage des valeurs d'un double : -1.7*10^-308 à 1.7*10^308 (à confirmer). Ce programme illustre juste la perte de précision avec des doubles.

Par contre, je suis étonné qu'avec un int sous Dev C++, il trouve 1, car là effectivement il y a dépassement. Sauf si notre ami est en 64 bits...