Question

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

Answer 1

首先,自定义函数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