function I=QuadComposite(f,A,B,n,method)
h=(B-A)/n;
x=A:h:B;
I=0;
for k=1:n
  I=I+method(f,x(k),x(k+1));
end
