%set interval
h=0.01;
%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:(n-1) 
%fxint=fxint+fx(ii);
%end

%apply composite trap. rule
%integral=0.5*h*(fx(1) + 2*fxint + fx(n))
integrall=0.5*h*(fx(1) + 2*sum(fx(2:(n-1))) +fx(n))