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。