6 May       2023   4:20:53.121 PM      
  
 quad2d():
   FORTRAN77 version
   Estimate the integral of f(x,y) over [0,1]x[0,1].
   f(x,y) = 1 / ( 1 - x * y ).
  
   A        =    0.0000000000000000     
   B        =    1.0000000000000000     
   NX       =        32768
   NY       =        32768
   N        =   1073741824
   Exact    =    1.6449340668482264     
  
   Estimate =    1.6449188090126949     
   Error    =    1.5257835531468089E-005
   Time     =    17.067444999999999     
  
 quad2d():
   Normal end of execution.
  
 6 May       2023   4:21:10.189 PM