使用matlab编写好梯形公式的函数M文件后如何使用它去计算积分?

使用matlab编写好梯形公式的函数M文件后如何使用它去计算积分?
已知
复化梯形公式M文件：
function y=trap(‘f’,a,b,n)
h=(b-a)/n;
k=n-1;
sum=0;
for i=1:k
sum=sum+feval(f,a+i*h);
end
y=h/2*( feval (f,a)+2*sum+ feval (f,b));