C语言：圆周率的计算Description根据下面关系式,求圆周率pi的值,直到最后一项的值小于epsilon.epsi

C语言：圆周率的计算
Description
根据下面关系式,求圆周率pi的值,直到最后一项的值小于epsilon.epsilon的值由键盘输入.
pi/2=1+1/3+2!/(3*5)+3!/(3*5*7)+...+n!/(3*537*...*(2n+1))
Input
输入一个浮点数表示epsilon的值
Output
输出圆周率Pi的值,保留6位小数.
Sample Input
0.000001
Sample Output
3.141590
HINT
(参考：epsilon为0.01时,圆周率3.121501；epsioon为0.0001时,圆周率为3.141358)

#include<stdio.h>
int main()
{
int n;
double e,f,fz,fm,sum,pi;
n=0;
fz=1;
fm=1;
sum=0;
scanf("%lf",&e);
for(f=1;f>=e;n++)
{
if(n==0)
fz=1;
else
fz=fz*n;
fm=fm*(2*n+1);
f=fz/fm;
sum=sum+f;
}
pi=2*sum;
printf("%.6f\n",pi);
return 0;
}
可是交上去就是wrong answer,而且我的答案小数点最后一位和题目的答案差1.希望帮我看看哪里错了.

在for循环后加一句: sum -= f;
你多加了最后那个小于epsioon的值.