06 February 2008 08:19:12 AM

ASA091_PRB:
  C++ version

  Test routines in the ASA091 library.

TEST01:
  Perform a simple sample calculation using
  PPCHI2 to invert the Chi-Squared CDF.

  P =                                  0.017523
  V =                                         4
  G Log(Gamma(V/2)) =                         0
  VALUE =                    0.3999988241105069
  VALUE (correct) =                         0.4

  Error flag IFAULT = 0

TEST02:
  PPCHI2 computes percentage points of the Chi-Square CDF.
  Compare to tabulated values.

         N        CDF       X                         X2                    DIFF
                           (tabulated)               (PPCHI2)

           1     0.07966                      0.01      0.009999999999567905   4.321e-13
           2    0.004988                      0.01       0.01000000000000101   1.006e-15
           1      0.1125                      0.02       0.01999999999911529   8.847e-13
           2     0.00995                      0.02       0.02000000000000399   3.993e-15
           1      0.4729                       0.4        0.3999999999845615   1.544e-11
           2      0.1813                       0.4        0.4000000000015109   1.511e-12
           3     0.05976                       0.4        0.4000000000743913   7.439e-11
           4     0.01752                       0.4        0.4000000004011919   4.012e-10
           1      0.6827                         1        0.9999999999419344   5.807e-11
           2      0.3935                         1        0.9999999999992277   7.723e-13
           3      0.1987                         1        0.9999999999665969    3.34e-11
           4      0.0902                         1         1.000000000066364   6.636e-11
           5     0.03743                         1         0.999999999912167   8.783e-11
           3      0.4276                         2         1.999999999916729   8.327e-11
           3      0.6084                         3         3.000000000053088   5.309e-11
           3      0.7385                         4         4.000000000051578   5.158e-11
           3      0.8282                         5         5.000000000050925   5.093e-11
           3      0.8884                         6         6.000000000050505   5.051e-11
          10   0.0001721                         1         1.000020461888895   2.046e-05
          10     0.00366                         2         2.000000008541654   8.542e-09
          10     0.01858                         3         2.999999996729284   3.271e-09

ASA091_PRB:
  Normal end of execution.

06 February 2008 08:19:12 AM