a = 0;
b = 1;
disp(['           n                         sum                       a+n*h '])
format long e
for n = 10000:10000:100000
   sum1 = a;
   h = (b-a)/n;
   for j = 1:n
      sum1 = sum1 + h;
   end 
   disp([n  sum1  a+n*h])
end