main()
{
double B=4095.1;
double A=B+1;
/* A-B = 1 */
double x=1;
int n;
printf("initialisation x=%+.16f\n",x);
for (n=1 ; n<=9 ; n++)
{
x = (A*x) - B;
/* La première fois : */
/* */
/* (A*x)-B --> (A*1)-B --> A-B --> 1 */
/* */
/* donc 'x' est un invariant (x = 1)... */
printf("iteration %01d x = %+.16f\n",n,x);
}
}
initialisation x = +1.0000000000000000
iteration 1 x = +1.0000000000004547
iteration 2 x = +1.0000000018630999
JFC
www.lactamme.polytechnique.fr