Scipy 集成模块缺少 quad_vec

Question

Win 10、Python 3.7、Scipy 1.3.2、Anaconda

我在这里查看this question。

我试用了答案中的代码，

import numpy as np
from scipy import integrate
from scipy.special import erf
from scipy.special import j0
import quadpy

q = np.linspace(0.03, 1.0, 500)


def f(t):
    return t * 0.5 * (erf((t - 40) / 3) - 1) * j0(np.multiply.outer(q, t))


y, _ = integrate.quad_vec(f, 0, 50)
y1, _ = quadpy.quad(f, 0, 50)

print(y - y1)

但出现以下错误...

  File "D:/Python codes/temp.py", line 14, in <module>
    y, _ = integrate.quad_vec(f, 0, 50)

AttributeError: module 'scipy.integrate' has no attribute 'quad_vec'

鉴于提问者或回答者都没有遇到此问题，是什么导致我的系统出现此问题？

Answer 1

根据release notes，scipy.integrate.quad_vec 与 scipy 版本 1.4.0 一起发布。因此你需要升级你的 scipy。