scipy.integrate 计算的值不匹配

Question

在 scipy 模块中使用积分时，我注意到积分的值与预期值（我手动计算的）不匹配。

我有一个简单的函数如下：

import math
def integrand(x, L):
    return (math.cos(sc.pi * x / L))**2
     

L=10
I = quad(integrand, 0, 1, args=(L))

I = 0.9677446418943196 # from scipy

但理论上，结果应该是 0.5519233944327677。我很困惑为什么存在差异。

我试图找到的积分如下：

Answer 1

scipy 是正确的。你的计算不正确：

请注意，给定的 Decimal 与 scipy 获得的一致。这个问题最适合stackexchange。