参照例4-24编程,用梯形公式求(1+x^3)^(1/2) 从[1,3]的定积分

首先,自定义函数traint（）
traint.m
function y = traint(a,b,n,f)
%a,b为区间值
%n为a,b间的等分数
%f为计算积分的原函数
h = (b - a) / n;
x = linspace(a,b,n+1);
y1 = h * feval(f,x);
y1(1) = y1(1) / 2;
y1(n+1) = y1(n+1) / 2;
y = sum(y1);
然后,在命令窗口下输入
f=inline('(1+x.^3).^0.5');
traint(1,3,10,f)
ans =
6.2349