function [I]=trapez(f,a,b,m)

 

% integral funkcije f v mejah od a do b z m podintervalov

 

h=(b-a)/m;

 

x=linspace(a,b,m+1);

y=eval(vectorize(f));

 

I=h*(sum(y)-y(1)/2-y(m+1)/2);

 

return;