求助何洁月C++中的一个问题

笑尽英雄

//π= 16arctan(1/5)-4*arctan(1/239)
//arctan(x) = x - x^3/3 + x^5/5 - x^7/7......
#include <iostream.h>
#include <math.h>
double arctan(double x);
void main()
{
        double a,b,c;
        a = 16.0*arctan(1/5.0);
        b = 4.0*arctan(1/239.0);
        c = a-b;
        cout<<"π = "<<a-b<<endl;
}
double arctan(double x)
{
        int i;
        int flag = -1;
        double r,e,f,sqr;
        sqr = x*x;
        r = 0;
        e = x;
        i = 1;
        f = e/i;
        while (abs(f)>1e-15)
        {
               
                r = r+f;
                e = e*sqr*flag;
                i += 2;
                f = e/i;
        }
        return r;
}

笑尽英雄

看下什么问题哦，照何洁月说的敲的程序，为什么求的值就不对呢

笑尽英雄

额原来应该用 while (fabs(f)>1e-15)
。。。。