creating a function for repeated simpsons rule using a 2 column data set in matlab

main.m
function I = repeated_simpsons_rule(data)
% This function performs repeated Simpsons rule using a 2-column dataset
% - data: a n-by-2 matrix, where the first column is the x values and the
%         second column is the y values
% - I: the calculated integral value

x = data(:, 1);
y = data(:, 2);

n = length(x) - 1;
h = (x(n+1) - x(1)) / n;

f = y(1) + y(n+1) + 4*sum(y(2:2:n)) + 2*sum(y(3:2:n-1));
I = h * f / 3;
end
421 chars
16 lines

main.m
data = [x_values, y_values];
integral_value = repeated_simpsons_rule(data);
76 chars
3 lines