14 August 2022 05:04:53 PM

asa147_test():
  C version
  Test asa147().

test01():
  gammds() computes the incomplete Gamma function.
  Compare to tabulated values.

             A             X          FX                        FX2
                                      (Tabulated)               (GAMMDS)            DIFF

        0.1000        0.0300        0.7382350532339351        0.7382350532333595  5.7565e-13
        0.1000        0.3000        0.9083579897300343        0.9083579897007394  2.9295e-11
        0.1000        1.5000        0.9886559833621947        0.9886559832970216  6.5173e-11
        0.5000        0.0750        0.3014646416966613        0.3014646416964176  2.4375e-13
        0.5000        0.7500        0.7793286380801532        0.7793286380721339  8.0193e-12
        0.5000        3.5000        0.9918490284064972        0.9918490283318239  7.4673e-11
        1.0000        0.1000       0.09516258196404043       0.09516258196381351  2.2692e-13
        1.0000        1.0000        0.6321205588285577        0.6321205588240379  4.5197e-12
        1.0000        5.0000        0.9932620530009145        0.9932620528409565  1.5996e-10
        1.1000        0.1000       0.07205974576054322       0.07205974576039777  1.4545e-13
        1.1000        1.0000        0.5891809618706485        0.5891809618671922  3.4563e-12
        1.1000        5.0000        0.9915368159845525        0.9915368158487908  1.3576e-10
        2.0000        0.1500       0.01018582711118352       0.01018582711109095  9.2572e-14
        2.0000        1.5000        0.4421745996289254        0.4421745996212465  7.6789e-12
        2.0000        7.0000        0.9927049442755639        0.9927049441757336  9.9830e-11
        6.0000        2.5000       0.04202103819530612       0.04202103819182567  3.4805e-12
        6.0000       12.0000        0.9796589705830716        0.9796589704272012  1.5587e-10
       11.0000       16.0000        0.9226039842296428        0.9226039839748174  2.5483e-10
       26.0000       25.0000        0.4470785799755852        0.4470785798412643  1.3432e-10
       41.0000       45.0000        0.7444549220718699        0.7444549215211238  5.5075e-10

asa147_test():
  Normal end of execution.

14 August 2022 05:04:53 PM