08 November 2016  10:43:22.560 PM
 
ELLIPSE_MONTE_CARLO_PRB
  FORTRAN90 version
  Test the ELLIPSE_MONTE_CARLO library.

ELLIPSE_AREA1_TEST
  FORTRAN90 version
  ELLIPSE_AREA1 computes the area of an ellipse.

  R =    10.0000    
 
  Matrix A in ellipse definition x*A*x=r^2
 
  Col          1             2      
  Row
 
    1:      5.            1.      
    2:      1.            2.      
  Area =    104.720    

ELLIPSE_AREA1_TEST
  Normal end of execution.

ELLIPSE_AREA2_TEST
  FORTRAN90 version
  ELLIPSE_AREA2 computes the area of an ellipse.

  Ellipse:    5.00000     * x^2 +    2.00000     * xy +    2.00000     * y^2 =    10.0000    
  Area =    104.720    

ELLIPSE_AREA2_TEST
  Normal end of execution.
 
TEST01
  Use ELLIPSE01_SAMPLE to estimate integrals
  in the ellipse x' * A * x <= r^2.
  Ellipse area =    2.12410    
 
         N        1              X               Y               X^2               XY             Y^2             X^3
 
         1     2.12410        0.455600       -0.376066        0.977220E-01   -0.806627E-01    0.665814E-01    0.209604E-01
         2     2.12410       -0.306358E-01    0.147780        0.449434E-02   -0.307640E-02    0.105018E-01   -0.181719E-03
         4     2.12410        0.275870       -0.312000        0.417762E-01   -0.231694E-01    0.140882        0.700437E-02
         8     2.12410        0.333166E-02   -0.251031        0.251927E-01   -0.180584E-01    0.825508E-01   -0.312286E-03
        16     2.12410        0.534942E-02    0.718986E-01    0.230367E-01   -0.985476E-02    0.117591       -0.846337E-03
        32     2.12410       -0.576574E-01   -0.184163        0.317003E-01   -0.195772E-01    0.149220       -0.135138E-02
        64     2.12410       -0.250151E-01    0.240724E-01    0.278674E-01   -0.233873E-01    0.145951       -0.251368E-05
       128     2.12410        0.421042E-02   -0.629404E-01    0.285123E-01   -0.227768E-01    0.147640        0.484888E-04
       256     2.12410        0.164242E-02    0.581420E-01    0.283226E-01   -0.196689E-01    0.135295       -0.829320E-04
       512     2.12410        0.756784E-02    0.170093E-01    0.281605E-01   -0.224068E-01    0.141451        0.239302E-03
      1024     2.12410       -0.767481E-02   -0.145921E-01    0.298340E-01   -0.230815E-01    0.149816       -0.185755E-03
      2048     2.12410       -0.493584E-02    0.207855E-01    0.283452E-01   -0.212838E-01    0.140061       -0.484776E-04
      4096     2.12410       -0.768025E-02   -0.135906E-02    0.294489E-01   -0.226947E-01    0.142611       -0.155983E-03
      8192     2.12410       -0.663898E-03    0.852250E-02    0.296095E-01   -0.221527E-01    0.141227       -0.336923E-04
     16384     2.12410       -0.849378E-03    0.544156E-02    0.295450E-01   -0.227768E-01    0.141865        0.318350E-04
     32768     2.12410       -0.121951E-02   -0.488945E-03    0.296223E-01   -0.228369E-01    0.141868       -0.455365E-04
     65536     2.12410        0.128864E-02   -0.146042E-02    0.296739E-01   -0.225971E-01    0.140368        0.210398E-04
 
ELLIPSE_MONTE_CARLO_PRB
  Normal end of execution.
 
08 November 2016  10:43:22.660 PM