%set interval
h=0.1;
%initialize nodes
x=[1:h:1.6];
%evaluate fn values at nodes
fx=exp(x.^2);
n=length(fx);

%evaluate interior nodes' sum
fxint=0;
for ii=2:2:(n-1) 
fxint=fxint+fx(ii);
end

%evaluate interior end-nodes' sum
fxintend=0;
for iii=3:2:(n-2) 
fxintend=fxintend+fx(iii);
end

%apply composite simpson's 1/3 rule
integral=h*(fx(1)+4*fxint+2*fxintend+fx(n))/3