% Name:  Autar Kaw
% Section: Both
function [value] = my_own_integral(f,a,b,n)
sum1=0;
h=(b-a)/n;
for i=0:n
    sum1=sum1+f(a+i*h);
end
value=sum1*(b-a)/(n+1);