function Oct7()


   function Q = AdaptSimpson1(f,a,b,tol,fa,fm,fb)
 % Q is an estimate the integral of f from a to b.
 % tol is a heuristic error bound on the error.
 % Estimate the error two ways:
 if nargin==4
     fa = f(a); fm = f((a+b)/2); fb = f(b);
 end
 