function I=trap(f,a,b,n)
I=(f(a)+f(b))/2;
x=a;
h=(b-a)/n;
for i=1:n-1
    x=x+h;
    I=I+f(x);
end
I=I/n;