PTA 的题近似求PI

本题要求编写程序，根据下式求π的近似值，直到最后一项小于给定精度eps。

输入格式：
输入在一行中给出精度eps，可以使用以下语句来读输入：

scanf("%le", &eps);
输出格式：
在一行内，按照以下格式输出π的近似值（保留小数点后5位）：

PI = 近似值
输入样例：

1E-5

输出样例：

PI = 3.14158

代码

#include <stdio.h>
double fenzi(int m);//这是一个求(分子)阶层 
double fenmu(int m);//这是一个求分母 
int main()
{
     double eps,i;
     double s=0;double d;
     scanf("%le",&eps);

     for(i=0;1;i++) 
     {
      d=fenzi(i)/fenmu(2*i+1);
      s=s+d;
      if(d<eps)
      break;
     }

 printf("PI = %.5lf\n",2.0*s);
 return 0;
 } 
 double fenzi(int m)//阶层 
 {
      double s=1.0;
      int i;
      for(i=1;i<=m;i++) 
       	s=s*i;

  return s;
 }

 double fenmu(int m)//单个数分母（1x3x5......）
 {  
      double s=1;
      int i;
      for(i=1;i<=m;i=i+2)
	      s=s*i;
      
  return s;
 }