python實現數值積分的Simpson方法實例分析

 　　本文實例講述了python實現數值積分的Simpson方法。分享給大家供大家參考。具體如下：

　　?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 #coding = utf-8 #simpson 法計算積分，數值積分，效果非常理想 from math import * def func(x): """ 定義被積分函數 """ return x*sin(x) def Get_N(a,b,width): # width為步長 N=int((b-a)/width + 1) if N%2 == 0: N=N+1 return N def GenerateData(a,b,n,width): datas = [] r=a for i in range(0,n): datas.append(func(r)) r = r+width return datas def simpson_integral(datas,width,n): sum = datas[0]+datas[n-1] for i in range(2,n): if i%2== 0: sum = sum +4*datas[i-1] else: sum = sum +2*datas[i-1] return sum*width/3.0 if __name__ == "__main__": a=1.0 #積分上限 b=3.0 #積分下限 width=0.0625 #步長 N=Get_N(a,b,width) datas = GenerateData(a,b,N,width) print simpson_integral(datas,width,N)

　　希望本文所述對大家的Python程序設計有所幫助。