main() { double B=4095.1; double A=B+1; /* A-B = 1 */ double x=1; /* 'x' est invariant (x = 1)... */ int n; for (n=1 ; n<=9 ; n++) { printf("n=%02d x = %+.16f\n",n,x); x = (A*x) - B; /* 'x' est invariant (x = 1)... */ } } n=1 x = +1.0000000000000000 n=2 x = +1.0000000000004547 JFC www.lactamme.polytechnique.fr