Discussion :

[Andrey]

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

Code :

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

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

[Andrey]

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
#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 ?

[Médinoc]

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...

[Jean-Marc.Bourguet]

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.

[Jean-Marc.Bourguet]

Code :

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

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.