python實(shí)現(xiàn)數(shù)值積分的Simpson方法實(shí)例分析
本文實(shí)例講述了python實(shí)現(xiàn)數(shù)值積分的Simpson方法。分享給大家供大家參考。具體如下:
#coding = utf-8 #simpson 法計(jì)算積分,數(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)
希望本文所述對(duì)大家的Python程序設(shè)計(jì)有所幫助。
相關(guān)文章
python操作手機(jī)app的實(shí)現(xiàn)步驟
本文主要介紹了python操作手機(jī)app的實(shí)現(xiàn)步驟,本文將結(jié)合實(shí)例代碼,具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2021-07-07
Python比較set的規(guī)則及簡(jiǎn)單例子
在Python中,集合可以通過比較運(yùn)算符進(jìn)行比較,檢查子集、超集、相等性等關(guān)系,文中通過代碼介紹的非常詳細(xì),對(duì)大家學(xué)習(xí)或者使用python具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2024-11-11
python?使用第三方庫requests-toolbelt?上傳文件流的示例
這篇文章主要介紹了python?使用第三方庫requests-toolbelt?上傳文件流,本文通過實(shí)例代碼給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2022-09-09
python實(shí)現(xiàn)簡(jiǎn)易五子棋游戲(控制臺(tái)版)
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)簡(jiǎn)易五子棋游戲,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2022-05-05
Python Opencv實(shí)戰(zhàn)之印章提取的實(shí)現(xiàn)
很多時(shí)候我們需要電子版的章,但有些同學(xué)并不會(huì)通過Photoshop摳圖獲取。因此本文將利用Python OpenCV來實(shí)現(xiàn)印章的提取,感興趣的可以了解一下2022-03-03
python實(shí)現(xiàn)簡(jiǎn)單倒計(jì)時(shí)功能
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)簡(jiǎn)單倒計(jì)時(shí)功能,文中示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2021-04-04
跟老齊學(xué)Python之關(guān)于循環(huán)的小伎倆
不管是while還是for,所發(fā)起的循環(huán),在python編程中是經(jīng)常被用到的。特別是for,一般認(rèn)為,它要比while快,而且也容易寫(是否容易,可能因人而異,但是,執(zhí)行時(shí)間快,是的確的),因此在實(shí)踐中,for用的比較多點(diǎn)。2014-10-10

