用公式π/4=1-1/3+1/5-1/7...求π的近似值,直到发现某一项的绝对值小于10^6为止(该项不累加)

解:程序:

#include

#include

int main()

{

int sign = 1;

double pi = 0.0, n = 1.0, term = 1.0;//term表示当前项

while (fabs(term) >= 1e-6)

{

pi += term;

n += 2;

sign = -sign;

term = sign / n;

}

pi *= 4;

printf("pi=%10.8f\n", pi);

return 0;

}

结果:

pi=3.14159065

请按任意键继续. . .

本程序输出的结果是pi=3.14159065,虽然输出了8位小数,但是只有前5位小数3,14159是准确的,因为第7位已经小于10^-6,后面的项没有累加。

再看如下两个精度不同的程序:

程序1:

#include

#include

int main()

{

int sign=1;

int count = 0;

double pi = 0.0, n = 1.0, term = 1.0;//term表示当前项

while(fabs(term)>=1e-6)

{

count++;

pi += term;

n += 2;

sign = -sign;

term = sign / n;

}

pi *= 4;

printf("pi=%10.8f\n",pi);

printf("count=%d\n",count);

return 0;

}

结果:

pi=3.14159065

count=500000

请按任意键继续. . .

程序2:

#include

#include

int main()

{

int sign=1;

int count = 0;

double pi = 0.0, n = 1.0, term = 1.0;//term表示当前项

while(fabs(term)>=1e-8)//变化部分

{

count++;

pi += term;

n += 2;

sign = -sign;

term = sign / n;

}

pi *= 4;

printf("pi=%10.8f\n",pi);

printf("count=%d\n",count);

return 0;

}

结果:

pi=3.14159263

count=50000000

请按任意键继续. . .

精度不同,运行时间不同,程序2精度更高,但是运行次数是程序1的100倍。