/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        C A L C U L   D E   ' P I '   P A R   L A   F O R M U L E   D E   P L O U F F E  :                                         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/pi_plouff.01$c' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, AAAAMMJJhhmmss).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  <stdio.h>

#define   ITER      100

main()
     {
     double    n;
     double    seize=1.0;
     double    pi=0.0;

     for       (n=0 ; n<ITER ; n++)
               {
               pi = pi + (1/seize)*((4/((8*n)+1)) - (2/((8*n)+4)) - (1/((8*n)+5)) - (1/((8*n)+6)));
               seize = 16.0 * seize;
                                        /* Cette formule est donnee a la page 104 du numero 231 de Pour La Science (01/1997).        */
               }
     printf("\n pi = %f",pi);
     }