The simp38 function uses Simpson's 3/8 rule to calculate the
integral of the function f over the interval from a
to b. The parameter m sets the number of intervals
to use when evaluating. Additional options are passed to the
function f when evaluating.