我就废话不多说了,直接上代码吧!
# 龙贝格法求积分
import math
a=0 # 积分下限
b=1 # 积分上限
eps=10**-5 # 精度
T=[] # 复化梯形序列
S=[] # Simpson序列
C=[] # Cotes序列
R=[] # Romberg序列
def func(x): # 被积函数
y=math.exp(-x)
return y
def Romberg(a,b,eps,func):
h = b - a
T.append(h * (func(a) + func(b)) / 2)
ep=eps+1
m=0
while
1