Composite trapezoidal rule algorithm in Matlab

Composite trapezoidal rule algorithm in Matlab :

function trapez(a, b, n) %trapez(a, b, n) approximates the integral of a function f(x) in the %interval [a;b], by the composite trapezoidal rule %n is the number of subintervals %the user needs to specify the function f(x) at the bottom %Author: Alain G. Kapitho %Date : Jan 2006 h = (b-a)/n; sum = 0; for i = 1:n-1 x(i) = a + i*h; sum = sum + f(x(i)); end integral = h*(f(a) + 2*sum + f(b))/2 %this needs to be changed accordingly with the specific problem you have at %hand, before proceeding to the command line function y = f(x) y = 2 + sin(2*sqrt(x));