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

#include  <stdio.h>

#define   ITER      100

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

     for       (n=0 ; n<ITER ; n++)
               {
               pi = pi + (1/puissance)*(- 32.0/(4*n+1) - 1.0/(4*n+3)
                                        + 256.0/(10*n+1) - 64.0/(10*n+3) - 4.0/(10*n+5) - 4.0/(10*n+7) + 1.0/(10*n+9)
                                        );
               puissance = -(1024.0 * puissance);
                                        /* Cette formule est donnee a la page 53 du journal Hypotheses (avril ou mai 1997) de        */
                                        /* Texas Instrument.                                                                         */
               }
     pi = pi/64.0;

     printf("\n pi = %f",pi);
     }