c - 用阿基米德方程逼近 Pi

Question

我已经根据阿基米德方程创建了一个 pi 近似值。我已经用 for 模式和递归模式这样做了。在 for 模式下，我创建了类似上面的内容:

  double Pi_approximation(double r, double L){

   int i;
   double fin;
   double y;

   for(i=1; i<4; i++){

           y =sqrt(2*((r*r) - r*(sqrt(4*((r*r) - (L)))))); 

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

   return y;
}

我的问题是在 for 的第二个循环中。在第一个 y 中正常计算但在第二个循环中当我打印 y 和 L 时，它会打印出它们的指针!有什么想法吗？

Answer 1

很难确定发生了什么，但无论如何我都会猜一猜!您传递给其中一个 sqrt() 调用的值很可能是负数。发生这种情况时，y(因此 L)将为 NaN。然后输出将取决于编译器。例如在我的编译器上输出是:

-1.#IND00