Python中使用scipy.integrate.quad积分函数报错的解决求助

我来帮你梳理清楚问题的根源，以及对应的解决办法：

首先拆解你遇到的两个错误：

错误原因分析 #

1. 不加self时的报错：TypeError: integral() takes 0 positional arguments but 1 was given
   scipy的quad函数在调用积分函数时，会自动把**积分变量（也就是这里的E）**作为第一个参数传递给你的自定义函数。但你最初的函数没有定义任何参数，导致quad传参时出现不匹配的错误。

2. 加self时的报错：TypeError: only length-1 arrays can be converted to Python scalars
self是Python类方法的专属参数，你的函数并不是类的一部分，完全不需要加它。另外，你在函数里用arange(1, 5000, 1)生成了E数组，但quad是逐个传入单个标量E值的，数组和标量混合运算就触发了类型转换错误。

修正后的完整代码 #

核心思路是：让积分函数接收单个标量参数E，移除不必要的数组生成，同时删掉多余的self参数：

from numpy import radians, cos
from scipy.integrate import quad

def integrand(E):
    theta = 0
    costh = cos(radians(theta))
    a = 18 / (E * costh + 145)
    b = (E + 2.7 / costh)**-2.7
    c = (E + 5) / (E + 5 / costh)
    return a * b * c * (1**4)  # 1的4次方就是1，可直接简化为1

# 调用quad计算500到1000区间的积分
result, error_estimate = quad(integrand, 500, 1000)
print("积分结果：", result)
print("误差估计值：", error_estimate)

额外简化建议 #

因为theta=0，radians(0)等于0，cos(0)的结果是1，我们可以直接把costh替换成1，让代码更简洁高效：

from scipy.integrate import quad

def integrand(E):
    a = 18 / (E + 145)
    b = (E + 2.7)**-2.7
    c = (E + 5) / (E + 5)  # 这里分子分母相同，结果就是1，也可以直接省略c的计算
    return a * b * c

result, error_estimate = quad(integrand, 500, 1000)
print("积分结果：", result)
print("误差估计值：", error_estimate)

修改后，quad会自动遍历积分区间内的E值，逐个传入integrand函数计算，最终返回积分结果和对应的误差估计值。

内容的提问来源于stack exchange，提问作者Pedro Paxton