08-Jan-2022 07:59:40
quad_trapezoid_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  quad_trapezoid() applies
  an n-interval trapezoid rule to estimate the integral of a
  function f(x) over the interval [a,b].

  Estimate integral ( 0.000000 < x < 1.000000 ) x^2 dx
  quad_trapezoid: n = 10, integral estimate q =                    0.335

  Estimate integral ( 0.0 < x < 2.0 ) sqrt(x) dx
  quad_trapezoid: n = 10, integral estimate q =        1.868202538231812

  Estimate integral ( 0.0 < x < 2.0 ) humps(x) dx
  quad_trapezoid: n = 1, integral estimate q =       0.3212981744421901
  quad_trapezoid: n = 2, integral estimate q =        16.16064908722109
  quad_trapezoid: n = 4, integral estimate q =        16.17515212981744
  quad_trapezoid: n = 8, integral estimate q =        30.06323338261111
  quad_trapezoid: n = 16, integral estimate q =        28.91205380707738
  quad_trapezoid: n = 32, integral estimate q =        29.30623828354212
  quad_trapezoid: n = 64, integral estimate q =        29.32098832348961
  quad_trapezoid: n = 128, integral estimate q =        29.32490602380703
  quad_trapezoid: n = 256, integral estimate q =        29.32588676924805

quad_trapezoid_test:
  Normal end of execution.
08-Jan-2022 07:59:40