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
n=3 x = +1.0000000018630999
n=4 x = +1.0000076314440776
JFC
www.lactamme.polytechnique.fr