probleme de precision

bonjour à tous j'ai un ptit soucis pour un exo

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

#define N 69   #include <stdio.h>   int main() { int n; float sigma=0; printf("N vaut %d\n",N); for (n=1; n <= N; n++) { sigma = sigma + 1.0; } printf("somme calculee = %f\n",sigma); return 0; }

---

quand N est plus petit que 16777216 , ca marche ... quand N est plus grand ca marche plus , quelqu'un saurait il à quoi c'est du ? à la précision , sans aucun doute

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

#include <stdio.h>   /*   calculer :   U(n+1) <- A*U(n) -B   avec A=64.1 B=63.1 et U0=1   Normalement pour tout n, Un=1 */   int main() { double A=64.1; double B=63.1; double x=1; int n;   for (n=1; n <= 14; n++) { printf("n=%02d x=%+.16f\n",n,x); x = (A*x) - B; } return 0; }

---

j'ai cet autre soucis ... quelqu'un peut m'aider ? :oops:

Un est en effet supposé égal à 1.0.
Mais avec les erreurs d'approximation, je suppose qu'on obtient rapidement une différence, très vite accentuée par la multiplication par A...

Citation:

---

Envoyé par Andrey

quand N est plus petit que 16777216 , ca marche ... quand N est plus grand ca marche plus , quelqu'un saurait il à quoi c'est du ? à la précision , sans aucun doute

---

Les nombres en virgules flottantes ont pour valeur:

-1^s m 2^e

avec
s qui vaut 0 ou 1,
m entier positif
e entier positif

Tu dépasses ce qui est stockable dans m pour un float.

Citation:

---

Envoyé par Andrey

Code:

---

1 2	double A=64.1; double B=63.1;

---

---

Voir mon message précédant. 64.1 et 63.1 ne sont pas représentables en virgule flottante, A et B contiennent donc des approximations. Ce qui fait que par après tu n'as pas le résultat auquel tu t'attends.