create an integral function of trapezoid in matlab

main.m
function result = integral_trapezoid(f, a, b, n)
    % f: function to integrate
    % a: lower limit of integration
    % b: upper limit of integration
    % n: number of subintervals
    
    x = linspace(a, b, n+1);
    y = f(x);
    
    h = (b - a) / n;
    
    result = h * (sum(y) - 0.5 * (y(1) + y(end)));
end
318 chars
14 lines