% Script to test various composite quadrature rules
clear all;
fcn='oneOverX'; a=1; b=2; Itrapezoid=trapezoid(fcn,a,b,4);