c语言：求π的近似值

Posted 2020-06-13

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

解：程序：

#include<stdio.h>

#include<math.h>

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<stdio.h>

#include<math.h>

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<stdio.h>

#include<math.h>

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倍。

本文出自 “岩枭” 博客，请务必保留此出处http://yaoyaolx.blog.51cto.com/10732111/1741580