/*
* pi.c
* (c) Neil Gershenfeld 2/10/08
* calculation of pi by a scalar sum
* pi = 3.14159265358979323846
*/

#include <stdio.h>
#define NPTS 2000000
#define NLOOP 10
double a,b,c,x;
int i,loop;
main() {
    a = 0.5;
    b = 0.75;
    c = 0.25;
    for (loop = 1; loop <= NLOOP; ++loop) {
       x = 0;
       for (i = 1; i <= NPTS; ++i)
          x += a/((i-b)*(i-c)); }
    printf("pi.c: NPTS, nloop = %d, %d\n",NPTS,NLOOP);
    printf("pi = %.20f\n",x);
    }