python梯形积分

以下是Python中梯形积分的实现方法：

from sympy import *

# 定义被积函数
def f(x):
    return 2000*log(140000/(140000-2100*x))-9.8*x

# 定义积分区间
a = 8
b = 30

# 真值
truth = integrate(f(x), (x, a, b)).evalf()
print("真值：", truth)

# 梯形积分算法
n = 10 # 步长，将(a,b)区间分为n个块
h = (b-a)/n # 每个块的宽度
tra_result = 0 # 梯形积分结果
for i in range(n):
    tra_result += 1/2*h*(f(a+i*h)+f(a+(i+1)*h))

print("梯形积分值：", tra_result)