pipprox.c

   1 #include <stdio.h>
   2 #include <math.h>
   3
   4 int main(void)
   5 {
   6   double delta = 1e30;
   7
   8   for (unsigned y=1; y < ~0U/M_PI; y++)
   9     {
  10       unsigned x = y*M_PI + 0.5;
  11       double a = (double)x / y;
  12       double d = fabs(a-M_PI);
  13       if (d < delta)
  14         {
  15           delta = d;
  16           printf("%u/%u = %.20f (%.20f)\n", x, y, a, d);
  17         }
  18     }
  19
  20   return 0;
  21 }