C语言入门：轻松掌握圆周率计算的技巧——莱布尼茨公式

在给出代码前，先讲解一下计算圆周率的方法，即“莱布尼茨公式”。

       莱布尼茨公式是一种用来近似计算圆周率π的数学公式。它基于级数的思想，级数是一种将无穷多个数相加的运算。

       莱布尼茨公式的形式如下：

在这个公式中，每一项的分母都是奇数，且逐渐递增，并且它们的符号（正负）交替出现，而分子总是为1。通过将越来越多的项相加，我们可以逐渐接近圆周率π的值。

例如，当我们计算前几项时（下面的例子，将等式左边的”π/4“的4乘到右边去，便于表示π的值的变化）：
π ≈ 4/1 = 4
π ≈ 4/1 - 4/3 = 4 - 1.333... ≈ 2.666...
π ≈ 4/1 - 4/3 + 4/5 = 4 - 1.333... + 0.8 ≈ 3.466...

随着我们计算更多的项，所得到的结果会越来越接近圆周率π的真实值。因此，通过使用莱布尼茨公式，我们可以近似计算圆周率π的值。

下面直接“上代码”，虽然下面的代码已经有了详细的注释，但是结合上面的公式可以更好地理解下面的代码：

#include <stdio.h>
#include <math.h>

int main()
{ 
    int i=1, sign=1; // 定义整型变量i和sign，分别表示序列项和符号
    double t=1, pi=0; // 定义双精度浮点型变量t和pi，分别表示序列项和π的近似值

    while(fabs(t)>=1e-4) // 当序列项的绝对值大于等于 1e-4 时循环执行以下操作
    {
        pi=pi+t; // 将当前序列项累加到pi中，更新pi的值
        sign=-sign; // 更改符号，用于计算下一个序列项
        i=i+2; // 更新序列项，每次递增2
        t=1.0*sign/i; // 计算下一个序列项的值，并赋给t
    }
    
    printf("圆周率约为%f\n",4*pi); // 输出π的近似值，乘以4得到圆周率的近似值
}