16 August 2022 02:33:40 PM

asa091_test():
  C++ version
  Test asa091().

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.3999988240893519
  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.01000000000000066   6.627e-16
           2    0.004988                      0.01       0.01000000000000101   1.006e-15
           1      0.1125                      0.02       0.02000000000000264   2.637e-15
           2     0.00995                      0.02       0.02000000000000399   3.993e-15
           1      0.4729                       0.4        0.3999999999997083   2.917e-13
           2      0.1813                       0.4        0.4000000000015109   1.511e-12
           3     0.05976                       0.4        0.4000000000853208   8.532e-11
           4     0.01752                       0.4        0.4000000003800369     3.8e-10
           1      0.6827                         1         0.999999999999892    1.08e-13
           2      0.3935                         1        0.9999999999992277   7.723e-13
           3      0.1987                         1        0.9999999999994041   5.959e-13
           4      0.0902                         1         1.000000000000471   4.712e-13
           5     0.03743                         1         1.000000000003392   3.392e-12
           3      0.4276                         2         1.999999999998896   1.104e-12
           3      0.6084                         3         2.999999999999514   4.858e-13
           3      0.7385                         4         4.000000000000509   5.089e-13
           3      0.8282                         5         5.000000000001442   1.442e-12
           3      0.8884                         6         6.000000000002119   2.119e-12
          10   0.0001721                         1          1.00002046236252   2.046e-05
          10     0.00366                         2         2.000000010187766   1.019e-08
          10     0.01858                         3         3.000000000006481   6.481e-12

asa091_test():
  Normal end of execution.

16 August 2022 02:33:40 PM