使用scipy.integrate.quadrature时max函数报错的原因咨询

执行以下Python脚本时出现错误：

from scipy import integrate
integrate.quadrature(lambda t: max(1,t), -2, 2)[0]

报错信息：

ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

发现integrate.quadrature和max函数不兼容，但不清楚具体原因。

scipy.integrate.quadrature为了提升计算效率，会向被积函数传入数组类型的参数（而非单个标量值）。但Python内置的max()函数只能处理单个数值的比较，无法对数组执行元素级的最大值计算。当数组传入max(1, t)时，函数会尝试判断整个数组的布尔值，这就触发了"数组真值不明确"的报错。

有几种可行的修复方式：

• 改用NumPy的np.maximum()，它支持元素级的最大值运算：

from scipy import integrate
import numpy as np
integrate.quadrature(lambda t: np.maximum(1, t), -2, 2)[0]

• 用条件表达式手动实现元素级判断：

from scipy import integrate
integrate.quadrature(lambda t: 1 if t < 1 else t, -2, 2)[0]

• 替换为scipy.integrate.quad，它默认处理标量输入，原逻辑可以直接运行：

from scipy import integrate
integrate.quad(lambda t: max(1, t), -2, 2)[0]

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