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                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
                                                                                     
 JFC                                                                                 
 www.lactamme.polytechnique.fr