利用公式求π的值。

现有一公式：π/2 = 1+1/3+(1+2)/(3+5)+...+(1+2+3+...+n)/(3+5+7+...+(2n+1)),现在要求利用此公式，用户输入准确率求得π的值，比如准确率0.0005，输出3.14，即公式中任意一项比0.0005小的话，就停止迭代。

double fun(double eps)

{ 

  double s= 0.0,s1=1.0;

  int n =1;

while(s1>=eps)

{

  s= s+s1;

  s1=s1*n/(2*n+1);

  n++;

}

return 2*s;

}