matlab复合梯形法,matlab实现复合梯形法则

复合梯形法则：

function int_f = CompoundEchelon( f, a, b, m )

% input : f : function handler

% a : the lower limit of integral

% b : the upper limit of integral

% m : cut integral area into m peace

% output : int_f : the answer of the integral

h = (b - a) / m;

int_f = 0;

if m >= 2

for i = 1 : m-1

int_f = int_f + 2 * f(a + h * i);

end

end

int_f = int_f + f(a) + f(b);

int_f = int_f * h / 2;

end

例子：

clear all

format long

clc

%% (a)

fprintf(' (a) \n')

f = @(x) x./((x.^2+9).^0.5);

int1_16 = CompoundEchelon(f,