使用格雷戈里公式求pi的近似值

问题背景：

• 输入一个精度e,使用格雷戈里公式计算pi的近似值。

• 格雷戈里公式：

  

解决方案：

• 首先定义函数funpi(e)用来求pi的近似值。
• 使精确度精确到最后一项的绝对值小于e。
• 使用do-while循环控制是否需要持续输入。

程序源代码：

#include<stdio.h>
#include<math.h>
double funpi(double e);//函数声明
int main(void)
{
	char ch;
	double e,pi;
	do{
		printf("Enter e:");
		scanf("%lf",&e);
		pi=funpi(e);//调用函数
		printf("pi=%f\n",pi);

		printf("Do you want to continue?(Y/N):");
		scanf(" %c",&ch);
	} while((ch=='y')||(ch == 'Y'));

	printf("Thank you.\n");
	while(1);//显示
	return 0;
}

double funpi(double e)//定义计算pi的函数
{
	int denominator,flag;
	double item,sum;

	flag=1;
	denominator=1;
	item=1.0;
	sum=0;

	while(fabs(item)>=e){
		item=flag*1.0/denominator;
		sum=sum+item;//累加
		flag=-flag;
		denominator=denominator+2;
	}

	return sum*4;
}

编译运行结果：

Enter e:0.001
pi=3.143589
Do you want to continue?(Y/N):Y
Enter e:0.0001
pi=3.141793
Do you want to continue?(Y/N):Y
Enter e:0.00001
pi=3.141613
Do you want to continue?(Y/N):n
Thank you.