为何无法用变量定义数值积分代码的积分上限？

嘿，作为Python新手能自己搭起输入框架已经很棒了！你提到的“无法使用变量定义积分上限”的问题，其实是因为还没把b这个变量融入到具体的积分计算逻辑里。咱们一步步来完善代码：

第一步：先明确核心逻辑 #

你已经成功获取了积分上下限a、b和分割数N，接下来的关键是让这些变量参与到积分计算中——b会通过步长h=(b-a)/N，拆解成每个小区间的端点，最终完成近似积分的计算。

第二步：定义被积函数 #

首先得明确你要计算的积分对应的函数，这里拿示例函数f(x)=x²演示，你之后可以直接替换成自己需要的函数：

def f(x):
    return x ** 2  # 替换成你的目标被积函数即可

第三步：实现矩形法则（左/右矩形两种方式） #

矩形法则是把每个小区间近似为矩形，取区间左端点或右端点的函数值作为高：

左矩形法则 #

def left_rectangle_rule(a, b, N, f):
    h = (b - a) / N
    integral = 0.0
    for i in range(N):
        x = a + i * h  # 取每个区间的左端点
        integral += f(x) * h
    return integral

右矩形法则 #

def right_rectangle_rule(a, b, N, f):
    h = (b - a) / N
    integral = 0.0
    for i in range(1, N+1):
        x = a + i * h  # 取每个区间的右端点
        integral += f(x) * h
    return integral

第四步：实现梯形法则 #

梯形法则是把每个小区间近似为梯形，取区间两端点函数值的平均值作为高：

def trapezoidal_rule(a, b, N, f):
    h = (b - a) / N
    integral = 0.5 * (f(a) + f(b))  # 首尾端点各算一半
    for i in range(1, N):
        x = a + i * h
        integral += f(x)
    integral *= h
    return integral

第五步：整合到你的输入逻辑中 #

把上面的函数和你已有的输入代码结合，就能完整运行啦（还加了分割数的输入验证，避免非整数或负数输入）：

def f(x):
    return x ** 2  # 替换为你的被积函数

# 获取积分上下限
a = float(input('Lower limit ---> '))
while True:
    b = float(input('Upper limit ---> '))
    if b > a:
        break
    elif a == b:
        print('Integral = 0.')
        exit()  # 上下限相等时直接结束程序
    else:
        print('Invalid input. Upper limit must be greater than lower limit.')

# 获取正整数分割数
while True:
    try:
        N = int(input('Number of integral divisions ---> '))
        if N > 0:
            break
        else:
            print('Number of divisions must be a positive integer.')
    except ValueError:
        print('Please enter an integer.')

h = (b - a) / N
print(f'For the integral in the range {a:.2f} to {b:.2f} with {N} divisions:')

# 计算并输出三种方法的结果
left_rect = left_rectangle_rule(a, b, N, f)
right_rect = right_rectangle_rule(a, b, N, f)
trapezoidal = trapezoidal_rule(a, b, N, f)

print(f'Left Rectangle Rule Result: {left_rect:.4f}')
print(f'Right Rectangle Rule Result: {right_rect:.4f}')
print(f'Trapezoidal Rule Result: {trapezoidal:.4f}')

关键说明 #

• 你之前的困惑点“无法使用变量定义积分上限”，本质是**b**需要通过步长h拆解到每个小区间的计算中，上面的代码已经把b完全融入到积分逻辑里了。
• 分割数**N**必须是正整数，所以增加了输入验证逻辑，避免程序出错。
• 只需要修改**f(x)**函数的内容，就能快速切换到你需要计算的任意积分。

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