function A=trap(f,a,b,n)
format long
x=linspace(a,b,n+1);h=(b-a)/n;
A=h/2*sum(f(x(1:n))+f(x(2:n+1)));