1、请编一个函数fun，它的功能是：根据以下公式求π的值(要求满足精度0.0005，即，某项小于0.0005时停止迭代)：
　　
　　程序运行后，如果输入精度0.0005，则程序输出多少。
　　注意：部分源程序存在PROG1.C中，请勿改动主函数和其他函数中的内容，仅在函数fun的指定的部位填入你编写的若干语句。
　　试题源程序如下：
　　#include ＜stdio.h＞
　　#include ＜math.h＞
　　double fun(double eps)
　　{
　　}
　　main()
　　{ double x;
　　　printf("Input eps:");
　　　scanf("%lf",&x); printf("\neps=%lf, PI=%lf\n",x,fun(x));
　　}

分析：
(1)本题所用基本算法应为累加。假设累加值放在变量s中，累加项放在变量t中，累加操作由语句s=s+t;来实现。
　　(2)若称 为第1累加项，则其前的1为第0累加项，其后的一项为第2累加项，按给定的公式可知，从第1累加项开始，后一项的累加项是前一项的值乘以 。所以当前的累加项应当是：t=t*n/(2.0*n+1.0)。表达式右边的t中是前一项的值，表达式左边的t为当前的每累加累加项。请注意，不要写成：t*n/(2*n+1)而进行整除。
　　(3)若第0累加项的1作为s的初值，语句：s=s+t;执行n次，就把第1到第n项累加到了s中。每进行一次累加，n值增1。
　　(4)把以上操作放在循环中。按本题规定，当某项小于eps(0.0005)时停止迭代，因此若用while循环，可用t＞=eps作为控制循环进行的条件：
　　while( t＞=eps ){ s+=t; n++; t=t*n/(2.0*n+1); }
　　(5)注意应给所用变量赋适当的初值。
　　(6)退出循环后，函数的返回值应是：2*s。