Nota : please, click inside the following rectangle in order to access the next page [cliquez à l'intérieur du rectangle ci-dessous afin de passer à la page suivante].

                                                                                     
     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                                                       



[Go back to AVirtualMachineForExploringSpaceTimeAndBeyond [Retour à AVirtualMachineForExploringSpaceTimeAndBeyond]]
[Site Map, Help and Search [Plan du Site, Aide et Recherche]]
[Mail [Courrier]]


Copyright (c) Jean-François Colonna, 2019.
Copyright (c) CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / Ecole Polytechnique, 2019.