python實現(xiàn)數(shù)值積分的Simpson方法實例分析
更新時間:2015年06月05日 10:12:16 作者:小蘿莉
這篇文章主要介紹了python實現(xiàn)數(shù)值積分的Simpson方法,實例分析了Python實現(xiàn)積分運算的相關技巧,需要的朋友可以參考下
本文實例講述了python實現(xiàn)數(shù)值積分的Simpson方法。分享給大家供大家參考。具體如下:
#coding = utf-8 #simpson 法計算積分,數(shù)值積分,效果非常理想 from math import * def func(x): """ 定義被積分函數(shù) """ 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程序設計有所幫助。
相關文章
python?使用第三方庫requests-toolbelt?上傳文件流的示例
這篇文章主要介紹了python?使用第三方庫requests-toolbelt?上傳文件流,本文通過實例代碼給大家介紹的非常詳細,對大家的學習或工作具有一定的參考借鑒價值,需要的朋友可以參考下2022-09-09
Python Opencv實戰(zhàn)之印章提取的實現(xiàn)
很多時候我們需要電子版的章,但有些同學并不會通過Photoshop摳圖獲取。因此本文將利用Python OpenCV來實現(xiàn)印章的提取,感興趣的可以了解一下2022-03-03

