分类：龙贝格法

2020
09-23

Python龙贝格法求积分实例

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