python积分求助

Eternal03

1. from mpmath import *
   
2. import math
   
3. 
   
4. mp.dps = 53
   
5. 
   
6. def f(x):
   
7.     return 3 * x ** 2 + 181 * x - 312
   
8. 
   
9. def g(x):
   
10.     return exp(2 * x) + log10(x ** 2 + 5 * x) + 32 * x - 201
    
11. 
    
12. def product_fg(x):
    
13.     return f(x) * g(x)
    
14. a = mp.mpf(input("Please enter the lower limit of integration (0.0001 < a < 10000): "))
    
15. b = mp.mpf(input("Please enter the upper limit of integration (0.0001 < b < 10000): "))
    
16. 
    
17. result = quad(product_fg, [a, b], method='tanh-sinh', epsabs=1e-15, epsrel=1e-15)
    
18. integral, error = result
    
19. print(f"The result of integration is：{integral}")  
    
20. print(f"The error is estimated as：{error}")

复制代码

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

FishC

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

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

复制代码

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

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

1. pip install mpmath

复制代码

希望对你有帮助！

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

小甲鱼的二师兄

1. from mpmath import *
   
2. import math
   
3. 
   
4. # 设置高精度
   
5. mp.dps = 53
   
6. 
   
7. # 定义函数f
   
8. def f(x):
   
9.     return 3 * x ** 2 + 181 * x - 312
   
10. 
    
11. # 定义函数g
    
12. def g(x):
    
13.     return exp(2 * x) + log10(x ** 2 + 5 * x) + 32 * x - 201
    
14. 
    
15. # 定义函数product_fg
    
16. def product_fg(x):
    
17.     return f(x) * g(x)
    
18. 
    
19. # 获取用户输入并检查范围
    
20. a = mp.mpf(input("Please enter the lower limit of integration (0.0001 < a < 10000): "))
    
21. b = mp.mpf(input("Please enter the upper limit of integration (0.0001 < b < 10000): "))
    
22. 
    
23. if not (0.0001 < a < 10000 and 0.0001 < b < 10000):
    
24.     print("Error: The limits must be in the range (0.0001, 10000).")
    
25. else:
    
26.     # 使用mpmath的quad函数进行积分
    
27.     result = quad(product_fg, [a, b], method='tanh-sinh', epsabs=1e-15, epsrel=1e-15)
    
28.    
    
29.     # 打印结果
    
30.     print(f"The result of integration is: {nstr(result, 10)}")

复制代码

主要修正点：

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