Python 龙贝格/Romberg算法

Python 龙贝格/Romberg算法

龙贝格求积公式也称为逐次分半加速法。是数值计算方法之一,用以求解数值积分。是在梯形公式、辛普森公式和柯特斯公式之间关系的基础上,构造出一种加速计算积分的方法。 作为一种外推算法,在不增加计算量的前提下提高了误差的精度。 在等距基点的情况下,用计算机计算积分值通常都采用把区间逐次分半的方法进行。这样,前一次分割得到的函数值在分半以后仍可被利用,且易于编程。 # 用于储存函数,返回函数值 def fun(x): return 1 / (x + 1) # 存放求积分范围 Min = float(input('请输入求积分下限:')) Max = float(input('请输入求积分上限:')) Mid = Max - Min # 存放积分结果的精度要求 precision = 10 ** int('-'+input('请输入求积精度10^-(输入值):')) # 用于存放 T S C R 的计算结果 T = [[j for j in...