Sun Oct 23 15:59:51 2022

quad2d_test():
  Python version: 3.6.9
  Test quad2d().
  Estimate the integral of f(x,y) over [a,b]x[c,d].
  f(x,y) = 1 / ( 1 - x * y ).
  using a  16384  x  16384  grid of  268435456  points.
   0.0  <= x <=  1.0
   0.0  <= y <=  1.0

  Exact =     1.6449340668482264
  Estimate =  1.6449035528760592
  Error    =  3.0513972167245385e-05
  Time     =  85.367237

quad2d_test():
  Normal end of execution.
Sun Oct 23 16:01:16 2022