create a script with simpsons rule to find area under the curve in matlab

main.m
% Define the function
f = @(x) x.^2 + 2*x + 1;

% Define the interval and the number of subdivisions
a = 0;
b = 4;
n = 20;

% Calculate the step size
h = (b-a)/n;

% Use Simpson's Rule
x = a:h:b;
y = f(x);
area = h/3*(y(1) + 4*sum(y(2:2:end-1)) + 2*sum(y(3:2:end-2)) + y(end));

% Display the area
disp(['The area under the curve is approximately ' num2str(area)]);
366 chars
19 lines