simpsonsRule.m

function area = simpsonsRule(f,a,b,n)
% simpsonsRule takes four arguments
% f is a function handle
% a and b are the endpoints of the integral
% n is the number of intervals to use for the approximation
%
% area is the value computed by using a three-point Simpson's rule
% approximation on each subinterval

% Fill in the details of the implementation here so your 
% code will return a correct approximation of the area
area = 0;

end