快捷導(dǎo)航

python實(shí)現(xiàn)數(shù)值積分的Simpson方法實(shí)例分析

更新時(shí)間：2015年06月05日 10:12:16 作者：小蘿莉

這篇文章主要介紹了python實(shí)現(xiàn)數(shù)值積分的Simpson方法,實(shí)例分析了Python實(shí)現(xiàn)積分運(yùn)算的相關(guān)技巧,需要的朋友可以參考下

本文實(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為步長(zhǎng)
 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 #步長(zhǎng)
 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常見(jiàn)錯(cuò)誤類(lèi)型及解決方法
如果說(shuō)寫(xiě)代碼最害怕什么，那無(wú)疑是Bug。而對(duì)于新手來(lái)說(shuō)，剛剛接觸編程，在享受寫(xiě)代碼的成就感時(shí)，往往也會(huì)被各式各樣的Bug弄得暈頭轉(zhuǎn)向。Python常見(jiàn)錯(cuò)誤類(lèi)型及解決方法，拯救你的代碼！消滅Bug的過(guò)程就是技能經(jīng)驗(yàn)的累積過(guò)程。
2023-06-06
python操作手機(jī)app的實(shí)現(xiàn)步驟
本文主要介紹了python操作手機(jī)app的實(shí)現(xiàn)步驟，本文將結(jié)合實(shí)例代碼，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
2021-07-07
Python比較set的規(guī)則及簡(jiǎn)單例子
在Python中,集合可以通過(guò)比較運(yùn)算符進(jìn)行比較,檢查子集、超集、相等性等關(guān)系,文中通過(guò)代碼介紹的非常詳細(xì),對(duì)大家學(xué)習(xí)或者使用python具有一定的參考借鑒價(jià)值,需要的朋友可以參考下
2024-11-11
python?使用第三方庫(kù)requests-toolbelt?上傳文件流的示例
這篇文章主要介紹了python?使用第三方庫(kù)requests-toolbelt?上傳文件流,本文通過(guò)實(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實(shí)現(xiàn)多層感知器
這篇文章主要為大家詳細(xì)介紹了python實(shí)現(xiàn)多層感知器的相關(guān)資料，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下
2019-01-01
Python Opencv實(shí)戰(zhàn)之印章提取的實(shí)現(xiàn)
很多時(shí)候我們需要電子版的章，但有些同學(xué)并不會(huì)通過(guò)Photoshop摳圖獲取。因此本文將利用Python OpenCV來(lái)實(shí)現(xiàn)印章的提取，感興趣的可以了解一下
2022-03-03
Python中三種花式打印的示例詳解
在Python中有很多好玩的花式打印，我們今天就來(lái)挑戰(zhàn)下面三個(gè)常見(jiàn)的花式打印。文中的示例代碼講解詳細(xì)，感興趣的小伙伴快跟隨小編一起學(xué)習(xí)一下吧
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快，而且也容易寫(xiě)（是否容易，可能因人而異，但是，執(zhí)行時(shí)間快，是的確的），因此在實(shí)踐中，for用的比較多點(diǎn)。
2014-10-10