根据如图,求π的近似值:
题目描述
用如下公式
求圆周率PI的近似值,直到发现某一项的绝对值小于10-6为止(该项不累加)。
要求输出的结果总宽度占10位,其中小数部分为8位。
程序中使用浮点型数据时,请定义为双精度double类型。
如果需要计算绝对值,可以使用C语言数学库提供的函数fabs,如求x的绝对值,则为fabs(x).
输入
无
输出
PI=圆周率的近似值
输出的结果总宽度占10位,其中小数部分为8位。
末尾输出换行。
样例输入
无
1
样例输出
PI=3.14159065
这一题做的比较坎坷,首先网速不好图片没加载出来,没懂为什么前几位跟acos(-1)截取的不同= =
后来网上搜看到图片之后,先写了一个错误版本,调试很久总是不停打印输出,到底不知道哪里的问题。。。
这是错误的代码
#include <stdio.h>
#include <math.h>
int main(){
double a,b;
double pi1 = 0,pi2 = 0;
int n = 0;
for(n = 0;;n++){
if(n % 2 != 0){
a = (-1.0)* (1.0/(2.0 * n + 1.0));
while(fabs(a) < pow(10.0,-6.0)){
int N = n;
for(n = 0;n < N;n++){
if(n % 2 != 0){
a = (-1.0)* (1.0/(2.0 * n + 1.0));
pi1 = pi1 + a;
}
else{
b = 1.0/(2.0 * n + 1.0);
pi2 = pi2 + b;
}
printf("PI=%10.8f\n",4 * (pi1 + pi2));
}
break;
}
//break;
}
//break;
else{
a = 1.0/(2.0 * n + 1.0);
//printf(%)
while(fabs(a) < pow(10.0,-6.0)){
int N = n;
for(n = 0;n < N;n++){
if(n % 2 != 0){
a = (-1.0)* (1.0/(2.0 * n + 1.0));
pi1 = pi1 + a;
}
else{
b = 1.0/(2.0 * n + 1.0);
pi2 = pi2 + b;
}
printf("PI=%10.8f\n",4 * (pi1 + pi2));
}
break;
}
//break;
}
}
}
处理不好break的位置,等以后有时间再调试;
下面是正确的版本,而且正确的更简单容易理解:
#include <stdio.h>
#include <math.h>
int main(){
double a,b;
double pi1 = 0,pi2 = 0;
int n = 0;
a = 1.0 * (1.0/(2.0 * n + 1.0));
while(fabs(a) > pow(10.0,-6.0)){
n++;
a = 1.0 * (1.0/(2.0 * n + 1.0));
}
int N = n;
for(n = 0;n < N;n++){
if(n % 2 != 0){
a = (-1.0)* (1.0/(2.0 * n + 1.0));
pi1 = pi1 + a;
}
else{
b = 1.0/(2.0 * n + 1.0);
pi2 = pi2 + b;
}
}
printf("PI=%10.8f\n",4 * (pi1 + pi2));
}
思路就是先找到绝对值<10的-6的那个n,将其赋值给N,再用个for循环打印输出就OK了!
那个错误的版本想一步到位,一边用for先无限尝试,找到绝对值小的数再计算,break,反而更复杂了。