科學(xué)計(jì)算應(yīng)該要安裝numpy,你可以去下Anaconda,全套自動(dòng)打包好了的。

創(chuàng)新互聯(lián)公司專(zhuān)注于企業(yè)成都全網(wǎng)營(yíng)銷(xiāo)推廣、網(wǎng)站重做改版、若羌網(wǎng)站定制設(shè)計(jì)、自適應(yīng)品牌網(wǎng)站建設(shè)、H5技術(shù)、成都商城網(wǎng)站開(kāi)發(fā)、集團(tuán)公司官網(wǎng)建設(shè)、成都外貿(mào)網(wǎng)站建設(shè)、高端網(wǎng)站制作、響應(yīng)式網(wǎng)頁(yè)設(shè)計(jì)等建站業(yè)務(wù),價(jià)格優(yōu)惠性?xún)r(jià)比高,為若羌等各大城市提供網(wǎng)站開(kāi)發(fā)制作服務(wù)。
把代碼改了一下:
#?-*-?coding:?utf-8?-*-
import?math
import?sys
import?time
import?numpy?as?np
def?area():
x1=1????????????#起始區(qū)間
x2=4.5????????????#結(jié)束區(qū)間
dx=0.0001????????#步長(zhǎng),步長(zhǎng)越小越精確
y1=0
y2=0
y3=0????????????#結(jié)果
for?x?in?np.arange(x1,?x2,?dx):
y1=y1+dx*x*x????????????#細(xì)小矩形取左邊為高
y2=y2+?dx*(x+0.5*dx)*(x+0.5*dx)?????????????#取中間
y3=y3+dx*(x+dx)*(x+dx)
print?(y1,?y2,?y3)
start?=?time.clock()????????
area()
end?=?time.clock()
print?(end-start)
你好,目前python有一個(gè)sympy的包。他提供一個(gè)積分計(jì)算的函數(shù):
integrate(函數(shù),(變量,下限, 上限))
所以你上面的那個(gè)積分計(jì)算可以如下:
from sympy import *
x = symbols('x')
print(integrate(abs(x-a), (x, b1, b2)))
用蒙特卡洛方法計(jì)算定積分
計(jì)算定積分
利用蒙特卡洛計(jì)算方法,核心步驟是求取隨機(jī)的 g(X1),………,g(Xn),n∈[a,b],由數(shù)學(xué)期望和大數(shù)定理可以近似計(jì)算定積分,公式為
原函數(shù):
導(dǎo)函數(shù):
計(jì)算導(dǎo)函數(shù)在[10,15]上的定積分;
Python
用蒙特卡洛方法計(jì)算的定積分:
直接用原函數(shù)計(jì)算的定積分:
偏差程度為:
本文名稱(chēng):python中求積分函數(shù) python 積分函數(shù)
文章分享:http://chinadenli.net/article38/hijgsp.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供商城網(wǎng)站、App開(kāi)發(fā)、網(wǎng)站設(shè)計(jì)、網(wǎng)頁(yè)設(shè)計(jì)公司、網(wǎng)站收錄、企業(yè)網(wǎng)站制作
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶(hù)投稿、用戶(hù)轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話(huà):028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)