用数字解决使用Python的集成限制

Question

我是Python新手，我正在尝试解决集成的上限：

def integrand(t):
    return (Dou)*np.exp((-Eau)/(8.314*(996.74 + (1037.1*np.exp(-0.2696*t)))))

def func(x):
    y, err=quad(integrand, 0, x)
    return y

sol = fsolve(func, 1.0)

print(sol[b])

其中Dou，Eau和b是随机定义的变量。

我想解决积分的上限，使积分等于b。我无法弄清楚这段代码有什么问题 - 任何明显的错误？

Answer 1

方法fsolve找到方程func = 0的根。如果等式右边有b而不是0，则从{{1}的输出中减去b }。

要创建一个自包含的示例，我将为b添加一个数值，以及代码中的其他未定义常量。

func