function y = int_rect1(a,b,n)
% numerical integration
% composite rectangle rule
h = (b-a)/n ;
x = a:h:b ;
s = a^2 ;
for i = 1:n
    s = s+x(i)^2;
end
y = s*h;


end

