write a functioni = csimpson(f,h)that takes as inputsfan array of functionvalues andhthe distance between grid values and computesithe compositesimpson’s approximation of the integral in matlab

main.m
function I = csimpson(f, h)
% CSIMPSON  Composite Simpson's rule for numerical integration
% Inputs:
%   f: array of function values
%   h: distance between grid values
% Output:
%   I: composite Simpson's approximation of the integral

% Check that length(f) is odd
if mod(length(f), 2) == 0
    error('Length of f must be odd')
end

% Compute the weights
n = length(f) - 1;
w = ones(1, n+1);
for i = 2:n
    if mod(i, 2) == 0
        w(i+1) = 2;
    else
        w(i+1) = 4;
    end
end

% Compute the approximation
I = h/3 * sum(w .* f);
end
545 chars
28 lines

main.m
>> f = [1, 5, 3, 8, 2];
>> h = 1;
>> I = csimpson(f, h)
I =
   14
66 chars
6 lines