编写复合梯形公式Tn的程序

下面是一个编写复合梯形公式Tn的Python程序示例：

def f(x):
    # 定义被积函数
    return x**2

def Trapezoidal(f, a, b, n):
    # 复合梯形公式Tn
    h = (b - a) / n
    x = [a + i*h for i in range(n+1)]
    s = sum([f(x[i]) + f(x[i+1]) for i in range(n)])
    return h/2 * s

# 测试
print(Trapezoidal(f, 0, 1, 100))

该程序定义了被积函数 `f(x)`，并使用复合梯形公式Tn计算该函数在区间 `[a, b]` 上的积分值。其中，`n` 为等分数，`h` 为区间长度除以等分数，`x` 为等分点的列表，`s` 为被积函数在等分点上的取值之和。最后，程序输出积分近似值。

相关问题

用matlab编写复合梯形公式Tn=h/2*[f(a)+f(b)+2∑f(xk)]的程序

可以按照以下方式编写复合梯形公式Tn的Matlab程序：

function T = compositeTrapezoidal(f, a, b, n)
% f: 被积函数
% a, b: 积分区间
% n: 子区间个数

h = (b - a) / n; % 子区间宽度
x = linspace(a, b, n+1); % 子区间节点
y = f(x); % 子区间节点处的函数值
T = h/2 * (y(1) + 2*sum(y(2:end-1)) + y(end)); % 复合梯形公式
end

其中，`f`是被积函数，`a`和`b`是积分区间的端点，`n`是子区间的个数。函数中使用`linspace`函数生成了节点向量`x`，并且通过`f`函数计算出节点处的函数值`y`。最后，根据复合梯形公式计算出积分值`T`。