python积分求助

Eternal03

from mpmath import *
import math

mp.dps = 53

def f(x):
    return 3 * x ** 2 + 181 * x - 312

def g(x):
    return exp(2 * x) + log10(x ** 2 + 5 * x) + 32 * x - 201

def product_fg(x):
    return f(x) * g(x)
a = mp.mpf(input("Please enter the lower limit of integration (0.0001 < a < 10000): "))
b = mp.mpf(input("Please enter the upper limit of integration (0.0001 < b < 10000): "))

result = quad(product_fg, [a, b], method='tanh-sinh', epsabs=1e-15, epsrel=1e-15)
integral, error = result
print(f"The result of integration is：{integral}")  
print(f"The error is estimated as：{error}")

一直报错显示TypeError: cannot unpack non-iterable mpf object，搞不懂

FishC

这个报错是因为 `quad` 函数返回的结果是一个元组，你需要使用适当的变量来接收返回的值。在这里，`result` 是一个包含积分结果和误差的元组。你可以使用下划线 `_` 来表示你不关心的变量。请将最后几行代码修改为以下方式：

result, _ = quad(product_fg, a, b, method='tanh-sinh', epsabs=1e-15, epsrel=1e-15)

现在，积分结果将被正确地赋值给 `result`，而误差将被忽略。然后你可以打印积分结果，并可选地打印误差。

另外，请确保你已经安装了 `mpmath` 这个第三方库，如果没有安装的话可以通过运行以下命令进行安装：

pip install mpmath

希望对你有帮助！

以上回复来自 FishC AI 小助理，如未能正确解答您的问题，请继续追问。

小甲鱼的二师兄

from mpmath import *
import math

# 设置高精度
mp.dps = 53

# 定义函数f
def f(x):
    return 3 * x ** 2 + 181 * x - 312

# 定义函数g
def g(x):
    return exp(2 * x) + log10(x ** 2 + 5 * x) + 32 * x - 201

# 定义函数product_fg
def product_fg(x):
    return f(x) * g(x)

# 获取用户输入并检查范围
a = mp.mpf(input("Please enter the lower limit of integration (0.0001 < a < 10000): "))
b = mp.mpf(input("Please enter the upper limit of integration (0.0001 < b < 10000): "))

if not (0.0001 < a < 10000 and 0.0001 < b < 10000):
    print("Error: The limits must be in the range (0.0001, 10000).")
else:
    # 使用mpmath的quad函数进行积分
    result = quad(product_fg, [a, b], method='tanh-sinh', epsabs=1e-15, epsrel=1e-15)
    
    # 打印结果
    print(f"The result of integration is: {nstr(result, 10)}")

主要修正点：

修正quad函数的返回值：quad函数返回的是单个值，所以不需要解包为integral, error。
增加输入范围检查：确保输入的上下限在规定范围内。
使用nstr函数打印结果：nstr函数可以格式化高精度数字以便于阅读。