main() { double B=4095.1; double A=B+1; double x=1; int n; printf("initialisation x=%+.16f\n",x); for (n=1 ; n<=9 ; n++) { x = (A*x) - B; printf("iteration %01d x=%+.16f\n",n,x); } }Lorsque je me suis intéressé il y a de nombreuses années à ces problèmes dits d'erreurs d'arrondi, je n'imaginai pas qu'un jour je réussirai à écrire quelque chose qui soit à la fois aussi élementaire et aussi "monstrueux" !