我深深需要計算this積分。我一直試圖這樣做幾個月,使用Python中的numpy包,特別是integrate.tplquad函數。 from __future__ import division
from math import *
import numpy as np
import scipy.special as special
import scipy.integrate as i
我試着同時計算積分,但是我的程序最終比用普通的for循環計算積分慢。我究竟做錯了什麼? package main
import (
"fmt"
"math"
"sync"
"time"
)
type Result struct {
result float64
lock sync.RWMutex
}
var wg sync.
我需要計算一個雙積分,其中內函數依賴於x和y(我將它整合到y上,其中一個限制取決於x),然後我將剩餘在x上。這是代碼的例子: [email protected](x,y) x.^2+y; %inner function to be integrated on y
[email protected](x) sqrt(x).*integral(@(y)f1(x,y),x,3); %second f
我已經定義了一個功能: def probability(x,t,i):
return np.real(np.conjugate(TD_Psi(x,t,i))*TD_Psi(x,t,i))
我到目前爲止已經使用整合的以下方法在Python for t in range(0,10):
PD = (np.trapz(probability(x,t,initial_state)