main()
{
double B=4095.1;
double A=B+1;
/* A-B = 1 à epsilon près !!! */
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
n=5 x = +1.0312591580864137
n=6 x = +129.0406374377594148
n=7 x = +524468.2550088063580915
n=8 x = +2148270324.2415719032287598
n=9 x = +8799530071030.8046875000000000
JFC
www.lactamme.polytechnique.fr