/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T E S T   D U   ' for(...) '   F L O T T A N T  :                                                                          */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/test_for_flottant.01$c' :                                                                                  */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  "INCLUDES.01.I"

main()
          {
          double    u;
          double    mu=+0.4,Mu=+0.5;
          double    pas;

          pas = +0.050;

          for       (u=mu ; u <= Mu ; u=u+pas)
                    {
                    printf("pas=%f     u=%f\n",pas,u);
                                        /* La derniere valeur affichee est "u=0.500000".                                             */
                    }

          printf("\n");

          pas = +0.025;

          for       (u=mu ; u <= Mu ; u=u+pas)
                    {
                    printf("pas=%f     u=%f\n",pas,u);
                                        /* La derniere valeur affichee est "u=0.475000". On notera que l'anomalie disparait si       */
                                        /* l'on ajoute une valeur superieure ou egale a 1.0e-16 a 'Mu'...                            */
                    }
          }



Copyright © Jean-François Colonna, 2021-2023.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2021-2023.