У меня есть следующий тройной интеграл:
интеграл
данные: t 1D массив (размер 124), pr 1D массив (размер 10),Lat 1D массив (241) и LON 1D массив (размер 480) V 4D массив (124, 10, 241, 480)
Я хочу интегрировать v по t, pr и lon в каждом лат. Итак, я использовалследующий код:
def M(T, lam, P, V, phi):
return integrate.tplquad(V*R*np.cos(phi), 0, T,
lambda T: 0, lambda T: lam,
lambda T,lam: 0, lambda T,lam: P)[0]
for i in range(lat.shape[0]):
a = M(t, lon, pr, v[:, :, i, :], lat[i])
Но я получил ошибку:
ValueError: Значение истинности массива с более чем одним элементом неоднозначно.Используйте a.any () или a.all ()
То, как работает tplquad, до сих пор сбивает меня с толку.Любая помощь в выполнении интеграла с tplquad или других функций?