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