实验6-1 近似求PI (15分)
本题要求编写程序,根据下式求π的近似值,直到最后一项小于给定精度eps。
输入格式:
输入在一行中给出精度eps,可以使用以下语句来读输入:
scanf("%le", &eps);
输出格式:
在一行内,按照以下格式输出π的近似值(保留小数点后5位):
PI = 近似值
输入样例:
1E-5
输出样例:
PI = 3.14158
代码:
#include <stdio.h>
double fact(int n);
double oddfact(int n);
int main()
{
int i;
double sum=1,eps,pi,item=1;
scanf("%le", &eps);
for(i=1;item>eps;i++)
{
item=fact(i)/oddfact(i);
//printf("PI = %.5f",item);
sum=sum+item;
//printf("PI = %.5f",sum*2);
}
printf("PI = %.5f",sum*2);
return 0;
}
double fact(int n)
{
int i;
double sum=1;
for(i=2;i<=n;i++)
sum=sum*i;
return sum;
}
double oddfact(int n)
{
int i;
double sum=1;
for(i=1;i<=n;i++)
sum=sum*(2*i+1);
return sum;
}