如何用SymPy和SciPy.quad正确积分指定函数？代码报错求解

你的SymPy代码返回错误结果，核心原因是运算优先级误解。Python中乘法与除法为左结合，你定义的z = (((x**2)-4*x+5*y)/4*x )实际等价于((x² - 4x + 5y) * x)/4；若你原本想表达的是(x² - 4x + 5y)/(4x)，写法错误会直接导致积分结果偏差。

一、SymPy 符号积分实现 #

正确代码示例 #

import sympy as sp

# 定义符号变量
x, y = sp.symbols('x y')

# 按左结合逻辑的原函数（对应你的代码写法）
z = ((x**2 - 4*x + 5*y)/4)*x
# 若你实际想表达分式形式，替换为：
# z = (x**2 - 4*x + 5*y)/(4*x)

# 计算积分并化简
integral = sp.integrate(z, x).simplify()
print(integral)

对应积分结果 #

• 左结合形式((x²-4x+5y)*x)/4的积分结果：

  x⁴/16 - x³/3 + (5x²y)/8 + C（C为积分常数）

• 分式形式(x²-4x+5y)/(4x)的积分结果：

  x/4 - log(x) + (5y log(x))/4 + C（C为积分常数）

二、SciPy.quad 数值积分实现 #

SciPy.quad是数值积分工具，需将y视为固定参数，传入具体数值计算：

左结合形式代码示例 #

from scipy.integrate import quad

# 定义被积函数
def integrand(x, y):
    return ((x**2 - 4*x + 5*y)*x)/4

# 设定y的具体值与积分区间
y_val = 2
result, error = quad(integrand, 0, 1, args=(y_val,))

print(f"数值积分结果：{result}")
print(f"积分误差估计：{error}")

分式形式代码示例 #

注意分式形式中x不能取0，需避开该点设置积分区间：

from scipy.integrate import quad

def integrand(x, y):
    return (x**2 - 4*x + 5*y)/(4*x)

y_val = 2
# 积分区间避开x=0，示例取[1,2]
result, error = quad(integrand, 1, 2, args=(y_val,))

print(f"数值积分结果：{result}")
print(f"积分误差估计：{error}")

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