7-4 近似求PI (10分)
本题要求编写程序,根据下式求π的近似值,直到最后一项小于给定精度eps。
输入格式:
输入在一行中给出精度eps,可以使用以下语句来读输入:
scanf("%le", &eps);
输出格式:
在一行内,按照以下格式输出π的近似值(保留小数点后5位):
PI = 近似值
输入样例:
1E-5
输出样例:
PI = 3.14158
#include <stdio.h>
#include<string.h>
double Fact( int n)//用来求分子
{
int i;
double result = 1;
for (i=1; i<=n; i++){
result = result * i * 1.0;
}
return result;
}
double PR(int m)//用来求分母
{
double js=1;
int j;
for(j=1; j<=m; j=j+2)
{
js=js*j;
}
return js;
}
int main()
{
int i=0;
double PI=0,eps,x;
scanf("%le",&eps);
while(i>=0)//建立循环
{
x=Fact(i)/PR(i*2+1);//0的阶乘为1
PI=PI+x;
i++;
if(x<eps)//结束循环的条件
break;
}
printf("PI = %.5lf",PI*2.0);
return 0;
}