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.希望帮我看看哪里错了。

谢谢！

相关了解……