Kaidah Trapesium dengan Matlab

function hasil=trapesium(a,b,n,f)
awal=0;
h=(b-a)/n;
for i=1:n-1,
awal=awal+ f(a+(i*h));
end
hasil=((h/2)*(f(a)+f(b)))+ (h*awal);