使用牛顿迭代法求圆周率的近似值。 已经知道无穷级数4/1 - 4/3 + 4/5 - 4/7 + ……的和是圆周率π,通过这一无穷级数前n项的和即可估算圆周率的值。当计算的圆周率与math库中的常数pi误差小于自然数e的(-m)次方的时候停止计算。
输入格式:
输入一个数m,表示误差参数,即当计算的圆周率与math库中的常数pi误差小于自然数e的(-m)次方的时候停止计算。
输出格式:
输出求取的无穷级数项数
输出估算的圆周率值,保留6位小数
输入样例:
在这里给出一组输入。例如:
5
输出样例:
在这里给出相应的输出。例如:
149
3.148304
解决方案:此题用的Python。①需要导入math库,先输入m,表示误差参数。定义sum,为估算的圆周率。定义n,用来求每一项的正负判断。②根据题目要求,当计算的圆周率与math库中的常数pi误差小于自然数e的(-m)次方的时候停止计算,则大于等于时,计算一直进行,圆周率sum更新,且n++。③最后输出项数n-1和sum。
核心代码:
import math
m = int(input()) #输入一个数m,表示误差参数
sum = 0 #估算的圆周率值
n = 1 #用来求每一项的正负判断的
#根据题目要求 小于时停止计算 所以进行计算的条件是大于等于
while math.fabs(sum - math.pi) >= math.pow(math.e, -1*m):
sum += 4/(2*n-1) * pow(-1,n+1)#计算圆周率,注意符号
n += 1#计算一次 n累加
print("%d\n%.6f"%(n-1,sum),end="")#注意输出的项数是n-1,不是n