20 May 2008 01:50:23 PM TOMS179_PRB: C version Test routines in the TOMS179 library. TEST01: ALOGAM computes the logarithm of the Gamma function. We compare the result to tabulated values. X FX FX2 (Tabulated) (ALOGAM) DIFF 2.0000000000000001e-01 1.5240638224307841e+00 1.5240638224152658e+00 1.5518e-11 4.0000000000000002e-01 7.9667781770178370e-01 7.9667781768963208e-01 1.2152e-11 5.9999999999999998e-01 3.9823385806923478e-01 3.9823385805965406e-01 9.5807e-12 8.0000000000000004e-01 1.5205967839983750e-01 1.5205967839224230e-01 7.5952e-12 1.0000000000000000e+00 0.0000000000000000e+00 -1.9950147436831323e-11 1.9950e-11 1.1000000000000001e+00 -4.9872441259839723e-02 -4.9872441277419147e-02 1.7579e-11 1.2000000000000000e+00 -8.5374090003315833e-02 -8.5374090018834864e-02 1.5519e-11 1.3000000000000000e+00 -1.0817480950786040e-01 -1.0817480952158205e-01 1.3722e-11 1.3999999999999999e+00 -1.1961291417237120e-01 -1.1961291418452211e-01 1.2151e-11 1.5000000000000000e+00 -1.2078223763524520e-01 -1.2078223764602811e-01 1.0783e-11 1.6000000000000001e+00 -1.1259176569675571e-01 -1.1259176570633646e-01 9.5808e-12 1.7000000000000000e+00 -9.5807697407065864e-02 -9.5807697415589782e-02 8.5239e-12 1.8000000000000000e+00 -7.1083872914372154e-02 -7.1083872921967411e-02 7.5953e-12 1.8999999999999999e+00 -3.8984275923083331e-02 -3.8984275929860819e-02 6.7775e-12 2.0000000000000000e+00 0.0000000000000000e+00 -1.9950147436831323e-11 1.9950e-11 3.0000000000000000e+00 6.9314718055994529e-01 6.9314718053999524e-01 1.9950e-11 4.0000000000000000e+00 1.7917594692280550e+00 1.7917594692081054e+00 1.9950e-11 1.0000000000000000e+01 1.2801827480081469e+01 1.2801827480080647e+01 8.2245e-13 2.0000000000000000e+01 3.9339884187199495e+01 3.9339884187199495e+01 0.0000e+00 3.0000000000000000e+01 7.1257038967168015e+01 7.1257038967168000e+01 1.4211e-14 TEST02: MDBETA estimates the value of th modified Beta function. Compare with tabulated values. X P Q Beta Beta DIFF (Tabulated) (MDBETA) 0.0100 0.5000 0.5000 6.3768560858519854e-02 6.3768560858622897e-02 1.0304e-13 0.1000 0.5000 0.5000 2.0483276469913350e-01 2.0483276469946440e-01 3.3090e-13 1.0000 0.5000 0.5000 1.0000000000000000e+00 1.0000000000000000e+00 0.0000e+00 0.0000 1.0000 0.5000 0.0000000000000000e+00 0.0000000000000000e+00 0.0000e+00 0.0100 1.0000 0.5000 5.0125628933800452e-03 5.0125628934800494e-03 1.0000e-13 0.1000 1.0000 0.5000 5.1316701949486197e-02 5.1316701950510031e-02 1.0238e-12 0.5000 1.0000 0.5000 2.9289321881345248e-01 2.9289321881929548e-01 5.8430e-12 0.5000 1.0000 1.0000 5.0000000000000000e-01 5.0000000000997513e-01 9.9751e-12 0.1000 2.0000 2.0000 2.8000000000000001e-02 2.8000000000558627e-02 5.5863e-13 0.2000 2.0000 2.0000 1.0400000000000000e-01 1.0400000000207488e-01 2.0749e-12 0.3000 2.0000 2.0000 2.1600000000000000e-01 2.1600000000430924e-01 4.3092e-12 0.4000 2.0000 2.0000 3.5199999999999998e-01 3.5200000000702258e-01 7.0226e-12 0.5000 2.0000 2.0000 5.0000000000000000e-01 5.0000000000997524e-01 9.9752e-12 0.6000 2.0000 2.0000 6.4800000000000002e-01 6.4799999999297742e-01 7.0226e-12 0.7000 2.0000 2.0000 7.8400000000000003e-01 7.8399999999569059e-01 4.3094e-12 0.8000 2.0000 2.0000 8.9600000000000002e-01 8.9599999999792512e-01 2.0749e-12 0.9000 2.0000 2.0000 9.7199999999999998e-01 9.7199999999944142e-01 5.5855e-13 0.5000 5.5000 5.0000 4.3619088505597770e-01 4.3619088506892362e-01 1.2946e-11 0.9000 10.0000 0.5000 1.5164090963470989e-01 1.5164090962250998e-01 1.2200e-11 0.5000 10.0000 5.0000 8.9782714843750000e-02 8.9782714845613079e-02 1.8631e-12 1.0000 10.0000 5.0000 1.0000000000000000e+00 1.0000000000000000e+00 0.0000e+00 0.5000 10.0000 10.0000 5.0000000000000000e-01 5.0000000000084022e-01 8.4022e-13 0.8000 20.0000 5.0000 4.5987732975757911e-01 4.5987732974680962e-01 1.0769e-11 0.6000 20.0000 10.0000 2.1468161023717389e-01 2.1468161023653076e-01 6.4312e-13 0.8000 20.0000 10.0000 9.5073648269578748e-01 9.5073648269574718e-01 4.0301e-14 0.5000 20.0000 20.0000 5.0000000000000000e-01 4.9999999999999134e-01 8.6597e-15 0.6000 20.0000 20.0000 8.9794136871059183e-01 8.9794136871059249e-01 6.6613e-16 0.7000 30.0000 10.0000 2.2412974918083661e-01 2.2412974918021045e-01 6.2617e-13 0.8000 30.0000 10.0000 7.5864054871920861e-01 7.5864054871901154e-01 1.9706e-13 0.7000 40.0000 20.0000 7.0017832474770691e-01 7.0017832474770092e-01 5.9952e-15 0.1000 1.0000 0.5000 5.1316701949486197e-02 5.1316701950510031e-02 1.0238e-12 0.2000 1.0000 0.5000 1.0557280900008410e-01 1.0557280900219028e-01 2.1062e-12 0.3000 1.0000 0.5000 1.6333997346592449e-01 1.6333997346918311e-01 3.2586e-12 0.4000 1.0000 0.5000 2.2540333075851660e-01 2.2540333076301342e-01 4.4968e-12 0.2000 1.0000 2.0000 3.5999999999999999e-01 3.6000000000718202e-01 7.1820e-12 0.2000 1.0000 3.0000 4.8799999999999999e-01 4.8800000000973576e-01 9.7358e-12 0.2000 1.0000 4.0000 5.9040000000000004e-01 5.9040000001177850e-01 1.1778e-11 0.2000 1.0000 5.0000 6.7232000000000003e-01 6.7232000001341341e-01 1.3413e-11 0.3000 2.0000 2.0000 2.1600000000000000e-01 2.1600000000430924e-01 4.3092e-12 0.3000 3.0000 2.0000 8.3699999999999997e-02 8.3700000001669772e-02 1.6698e-12 0.3000 4.0000 2.0000 3.0779999999999998e-02 3.0780000000614045e-02 6.1405e-13 0.3000 5.0000 2.0000 1.0935000000000000e-02 1.0935000000218157e-02 2.1816e-13 TOMS179_PRB: Normal end of execution. 20 May 2008 01:50:23 PM