14 June 2018 09:52:17 AM

ASA091_TEST:
  C version
  Test the ASA091 library.

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

  P =                        0.0175230000000000
  V =                        4.0000000000000000
  G Log(Gamma(V/2)) =        0.0000000000000000
  VALUE =                    0.3999988240893519
  VALUE (correct) =          0.4000000000000000

  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.0000      0.0797        0.0100000000000000        0.0100000000000007  6.6266e-16
      2.0000      0.0050        0.0100000000000000        0.0100000000000010  1.0061e-15
      1.0000      0.1125        0.0200000000000000        0.0200000000000026  2.6368e-15
      2.0000      0.0100        0.0200000000000000        0.0200000000000040  3.9933e-15
      1.0000      0.4729        0.4000000000000000        0.3999999999997083  2.9171e-13
      2.0000      0.1813        0.4000000000000000        0.4000000000015109  1.5109e-12
      3.0000      0.0598        0.4000000000000000        0.4000000000853208  8.5321e-11
      4.0000      0.0175        0.4000000000000000        0.4000000003800369  3.8004e-10
      1.0000      0.6827        1.0000000000000000        0.9999999999998920  1.0802e-13
      2.0000      0.3935        1.0000000000000000        0.9999999999992277  7.7227e-13
      3.0000      0.1987        1.0000000000000000        0.9999999999994041  5.9586e-13
      4.0000      0.0902        1.0000000000000000        1.0000000000004712  4.7118e-13
      5.0000      0.0374        1.0000000000000000        1.0000000000033920  3.3920e-12
      3.0000      0.4276        2.0000000000000000        1.9999999999988956  1.1044e-12
      3.0000      0.6084        3.0000000000000000        2.9999999999995142  4.8583e-13
      3.0000      0.7385        4.0000000000000000        4.0000000000005089  5.0893e-13
      3.0000      0.8282        5.0000000000000000        5.0000000000014424  1.4424e-12
      3.0000      0.8884        6.0000000000000000        6.0000000000021192  2.1192e-12
     10.0000      0.0002        1.0000000000000000        1.0000204623625195  2.0462e-05
     10.0000      0.0037        2.0000000000000000        2.0000000101877657  1.0188e-08
     10.0000      0.0186        3.0000000000000000        3.0000000000064806  6.4806e-12

ASA091_TEST:
  Normal end of execution.

14 June 2018 09:52:17 AM