13 August    2014   9:53:47.154 PM      
 
ELLIPSOID_MONTE_CARLO_PRB
  FORTRAN77 version
  Test the ELLIPSOID_MONTE_CARLO library.
 
TEST01
  Use ELLIPSOID_SAMPLE to estimate integrals
  in a 2D ellipse x' * A * x <= r^2.

  Ellipsoid radius R =    2.00000    
 
  Ellipsoid center V:
 
         1:    0.0000000    
         2:    0.0000000    
 
  Ellipsoid matrix A:
 
  Col         1             2       
  Row
 
    1:   9.00000       1.00000    
    2:   1.00000       4.00000    

  Ellipsoid volume =    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
 
TEST02
  Use ELLIPSOID_SAMPLE to estimate integrals
  in a 2D ellipse (x-v)' * A * (x-v) <= r^2.

  Ellipsoid radius R =    2.00000    
 
  Ellipsoid center V:
 
         1:    2.0000000    
         2:    3.0000000    
 
  Ellipsoid matrix A:
 
  Col         1             2       
  Row
 
    1:   9.00000       1.00000    
    2:   1.00000       4.00000    

  Ellipsoid volume =    2.12410    
 
         N        1              X               Y               X^2               XY             Y^2             X^3
 
         1     2.12410         4.70381         5.99625         10.4165         13.2786         16.9271         23.0673    
         2     2.12410         4.21757         6.52009         8.37837         12.9452         20.0141         16.6520    
         4     2.12410         4.52408         6.06031         9.64167         12.9251         17.3858         20.5609    
         8     2.12410         4.25154         6.12128         8.53494         12.2345         17.6933         17.1837    
        16     2.12410         4.25356         6.44421         8.54085         12.8946         19.6659         17.1944    
        32     2.12410         4.19055         6.18815         8.29749         12.1838         18.1612         16.4898    
        64     2.12410         4.22319         6.39639         8.42422         12.6943         19.4073         16.8599    
       128     2.12410         4.25242         6.30937         8.54177         12.6086         18.8869         17.2145    
       256     2.12410         4.24985         6.43045         8.53131         12.8462         19.6011         17.1824    
       512     2.12410         4.25578         6.38932         8.55485         12.7789         19.3604         17.2529    
      1024     2.12410         4.24053         6.35772         8.49555         12.6693         19.1792         17.0796    
      2048     2.12410         4.24327         6.39310         8.50502         12.7501         19.3817         17.1036    
      4096     2.12410         4.24053         6.37095         8.49515         12.6962         19.2514         17.0772    
      8192     2.12410         4.24754         6.38084         8.52337         12.7375         19.3093         17.1625    
     16384     2.12410         4.24736         6.37775         8.52256         12.7302         19.2915         17.1599    
     32768     2.12410         4.24699         6.37182         8.52116         12.7172         19.2559         17.1559    
     65536     2.12410         4.24950         6.37085         8.53125         12.7230         19.2485         17.1864    
 
TEST03
  Use ELLIPSOID_SAMPLE to estimate integrals
  in a 3D ellipse (x-v)' * A * (x-v) <= r^2.

  Ellipsoid radius R =   0.500000    
 
  Ellipsoid center V:
 
         1:    1.0000000    
         2:    2.0000000    
         3:    3.0000000    
 
  Ellipsoid matrix A:
 
  Col         1             2             3       
  Row
 
    1:   9.00000       6.00000       3.00000    
    2:   6.00000       5.00000       4.00000    
    3:   3.00000       4.00000       9.00000    

  Ellipsoid volume =   0.872665E-01
 
         N        1              X               Y                Z                X^2            YZ              Z^3
 
         1    0.872665E-01    0.118556        0.127999        0.270935        0.161063         1.80968         2.61156    
         2    0.872665E-01    0.505427E-01    0.232324        0.249758        0.301574E-01     5.06528         2.04716    
         4    0.872665E-01    0.917626E-01    0.165736        0.265086        0.984727E-01     2.91490         2.45045    
         8    0.872665E-01    0.885784E-01    0.173441        0.261981        0.972310E-01     3.17841         2.37401    
        16    0.872665E-01    0.954258E-01    0.160572        0.265292        0.114340         2.81944         2.46932    
        32    0.872665E-01    0.815371E-01    0.184104        0.258630        0.844829E-01     3.48828         2.28260    
        64    0.872665E-01    0.873417E-01    0.173498        0.262541        0.947942E-01     3.18914         2.38842    
       128    0.872665E-01    0.908308E-01    0.169208        0.262950        0.103113         3.05812         2.40171    
       256    0.872665E-01    0.844516E-01    0.179415        0.260577        0.905178E-01     3.36271         2.33783    
       512    0.872665E-01    0.863446E-01    0.176134        0.261319        0.954505E-01     3.27468         2.35974    
      1024    0.872665E-01    0.868068E-01    0.175389        0.261589        0.954095E-01     3.24681         2.36496    
      2048    0.872665E-01    0.880433E-01    0.173281        0.262051        0.977763E-01     3.18198         2.37747    
      4096    0.872665E-01    0.870301E-01    0.174981        0.261639        0.957777E-01     3.23227         2.36634    
      8192    0.872665E-01    0.867427E-01    0.175376        0.261592        0.951869E-01     3.24487         2.36523    
     16384    0.872665E-01    0.875100E-01    0.174143        0.261904        0.967846E-01     3.20988         2.37371    
     32768    0.872665E-01    0.870833E-01    0.174799        0.261744        0.959327E-01     3.22908         2.36940    
     65536    0.872665E-01    0.873715E-01    0.174369        0.261825        0.965099E-01     3.21578         2.37158    
 
ELLIPSOID_MONTE_CARLO_PRB
  Normal end of execution.
 
13 August    2014   9:53:47.408 PM