13 January   2008   1:26:49.722 PM      
 
ASA241_PRB:
  FORTRAN77 version
 
  Test the routines in ASA241.
 
TEST01:
  Let FX = NormalCDF ( X ).
 
  NORMAL_01_CDF_VALUES returns some values of ( X, FX ).
 
  PPND7 takes the value of FX, and
    computes an estimate X2, of the corresponding input,
    argument, accurate to about 7 decimal places.
 
      FX                         X                         X2                     DIFF
 
    0.5000000000000000         0.000000000000000         0.000000000000000       0.000    
    0.5398278372770290        0.1000000000000000        0.9999996423721313E-01  0.3576E-07
    0.5792597094391030        0.2000000000000000        0.1999999582767487      0.4172E-07
    0.6179114221889526        0.3000000000000000        0.2999999225139618      0.7749E-07
    0.6554217416103242        0.4000000000000000        0.3999999761581421      0.2384E-07
    0.6914624612740131        0.5000000000000000        0.5000000000000000       0.000    
    0.7257468822499270        0.6000000000000000        0.5999999046325684      0.9537E-07
    0.7580363477769270        0.7000000000000000        0.7000000476837158      0.4768E-07
    0.7881446014166033        0.8000000000000000        0.7999999523162842      0.4768E-07
    0.8159398746532405        0.9000000000000000        0.9000000357627869      0.3576E-07
    0.8413447460685429        1.0000000000000000        1.0000000000000000       0.000    
    0.9331927987311419         1.500000000000000         1.499999761581421      0.2384E-06
    0.9772498680518208         2.000000000000000         2.000000238418579      0.2384E-06
    0.9937903346742240         2.500000000000000         2.499999284744263      0.7153E-06
    0.9986501019683699         3.000000000000000         2.999992847442627      0.7153E-05
    0.9997673709209645         3.500000000000000         3.499990940093994      0.9060E-05
    0.9999683287581669         4.000000000000000         4.000158309936523      0.1583E-03
 
TEST02:
  Let FX = NormalCDF ( X ).
 
  NORMAL_01_CDF_VALUES returns some values of ( X, FX ).
 
  PPND16 takes the value of FX, and
    computes an estimate X2, of the corresponding input, 
    argument, accurate to about 16 decimal places.
 
      FX                         X                         X2                     DIFF
 
    0.5000000000000000         0.000000000000000         0.000000000000000       0.000    
    0.5398278372770290        0.1000000000000000        0.9999999999999999E-01  0.1388E-16
    0.5792597094391030        0.2000000000000000        0.1999999999999999      0.1110E-15
    0.6179114221889526        0.3000000000000000        0.2999999999999998      0.2220E-15
    0.6554217416103242        0.4000000000000000        0.4000000000000000      0.5551E-16
    0.6914624612740131        0.5000000000000000        0.4999999999999998      0.1665E-15
    0.7257468822499270        0.6000000000000000        0.6000000000000016      0.1665E-14
    0.7580363477769270        0.7000000000000000        0.6999999999999998      0.1110E-15
    0.7881446014166033        0.8000000000000000        0.7999999999999998      0.2220E-15
    0.8159398746532405        0.9000000000000000        0.9000000000000000       0.000    
    0.8413447460685429        1.0000000000000000        1.0000000000000000       0.000    
    0.9331927987311419         1.500000000000000         1.500000000000000      0.2220E-15
    0.9772498680518208         2.000000000000000         2.000000000000000       0.000    
    0.9937903346742240         2.500000000000000         2.500000000000004      0.3997E-14
    0.9986501019683699         3.000000000000000         2.999999999999997      0.3109E-14
    0.9997673709209645         3.500000000000000         3.499999999999983      0.1688E-13
    0.9999683287581669         4.000000000000000         4.000000000000000      0.4441E-15
 
ASA241_PRB:
  Normal end of execution.
 
13 January   2008   1:26:49.725 PM