13 March 2024 3:56:40.588 PM toms708_test(): Fortran77 version Test toms708(). TEST01 BRATIO computes the Beta ratio function. X Y W W1 0.01 0.99 0.6475472957845341E-07 0.9999999403953552 0.02 0.98 0.2361296083108755E-05 0.9999976158142090 0.03 0.97 0.1873401197372004E-04 0.9999812841415405 0.04 0.96 0.7957524940138683E-04 0.9999204277992249 0.05 0.95 0.2399659133516252E-03 0.9997600317001343 0.06 0.94 0.5825495463795960E-03 0.9994174242019653 0.07 0.93 0.1217441866174340E-02 0.9987825155258179 0.08 0.92 0.2279699314385653E-02 0.9977203011512756 0.09 0.91 0.3924931865185499E-02 0.9960750341415405 0.10 0.90 0.6323639769107103E-02 0.9936763644218445 0.11 0.89 0.9654808789491653E-02 0.9903451800346375 0.12 0.88 0.1409919280558825E-01 0.9859008193016052 0.13 0.87 0.1983261480927467E-01 0.9801673889160156 0.14 0.86 0.2701987139880657E-01 0.9729801416397095 0.15 0.85 0.3580895066261292E-01 0.9641910791397095 0.16 0.84 0.4632619768381119E-01 0.9536738395690918 0.17 0.83 0.5867260321974754E-01 0.9413273930549622 0.18 0.82 0.7292070239782333E-01 0.9270793199539185 0.19 0.81 0.8911225199699402E-01 0.9108877182006836 0.20 0.80 0.1072581112384796 0.8927419185638428 0.21 0.79 0.1273372471332550 0.8726627826690674 0.22 0.78 0.1492973715066910 0.8507026433944702 0.23 0.77 0.1730579435825348 0.8269420862197876 0.24 0.76 0.1985097825527191 0.8014901876449585 0.25 0.75 0.2255204021930695 0.7744796276092529 0.26 0.74 0.2539350986480713 0.7460649013519287 0.27 0.73 0.2835809290409088 0.7164191007614136 0.28 0.72 0.3142708241939545 0.6857291460037231 0.29 0.71 0.3458069562911987 0.6541930437088013 0.30 0.70 0.3779838383197784 0.6220161914825439 0.31 0.69 0.4105930626392365 0.5894069671630859 0.32 0.68 0.4434260129928589 0.5565739870071411 0.33 0.67 0.4762778282165527 0.5237221717834473 0.34 0.66 0.5089496970176697 0.4910503029823303 0.35 0.65 0.5412513613700867 0.4587486386299133 0.36 0.64 0.5730050206184387 0.4269949793815613 0.37 0.63 0.6040445566177368 0.3959554731845856 0.38 0.62 0.6342210769653320 0.3657789230346680 0.39 0.61 0.6634005308151245 0.3365994393825531 0.40 0.60 0.6914667487144470 0.3085332512855530 0.41 0.59 0.7183209657669067 0.2816790044307709 0.42 0.58 0.7438825368881226 0.2561174929141998 0.43 0.57 0.7680886387825012 0.2319113761186600 0.44 0.56 0.7908935546875000 0.2091064453125000 0.45 0.55 0.8122686743736267 0.1877313107252121 0.46 0.54 0.8322014808654785 0.1677985340356827 0.47 0.53 0.8506937026977539 0.1493062973022461 0.48 0.52 0.8677610754966736 0.1322389096021652 0.49 0.51 0.8834316730499268 0.1165682971477509 0.50 0.50 0.8977444171905518 0.1022555902600288 TEST02: BRATIO evaluates the normalized incomplete Beta function BETA_CDF(A,B,X). BETA_CDF_VALUES returns some exact values. A B X Exact F BRATIO(A,B,X) DIFF 0.50 0.50 0.01 0.6376855820417404E-01 0.6376855075359344E-01 0.75E-08 0.50 0.50 0.10 0.2048327624797821 0.2048327475786209 0.15E-07 0.50 0.50 1.00 1.000000000000000 1.000000000000000 0.0 1.00 0.50 0.00 0.000000000000000 0.000000000000000 0.0 1.00 0.50 0.01 0.5012562964111567E-02 0.5012562964111567E-02 0.0 1.00 0.50 0.10 0.5131670087575912E-01 0.5131670832633972E-01 0.75E-08 1.00 0.50 0.50 0.2928932309150696 0.2928932011127472 0.30E-07 1.00 1.00 0.50 0.5000000000000000 0.5000000000000000 0.0 2.00 2.00 0.10 0.2800000086426735E-01 0.2799999900162220E-01 0.19E-08 2.00 2.00 0.20 0.1040000021457672 0.1040000021457672 0.0 2.00 2.00 0.30 0.2160000056028366 0.2160000056028366 0.0 2.00 2.00 0.40 0.3519999980926514 0.3520000576972961 0.60E-07 2.00 2.00 0.50 0.5000000000000000 0.5000000000000000 0.0 2.00 2.00 0.60 0.6480000019073486 0.6480000615119934 0.60E-07 2.00 2.00 0.70 0.7839999794960022 0.7839999794960022 0.0 2.00 2.00 0.80 0.8960000276565552 0.8960000276565552 0.0 2.00 2.00 0.90 0.9720000028610229 0.9720000028610229 0.0 5.50 5.00 0.50 0.4361908733844757 0.4361908435821533 0.30E-07 10.00 0.50 0.90 0.1516409069299698 0.1516408324241638 0.75E-07 10.00 5.00 0.50 0.8978271484375000E-01 0.8978272229433060E-01 0.75E-08 10.00 5.00 1.00 1.000000000000000 1.000000000000000 0.0 10.00 10.00 0.50 0.5000000000000000 0.5000000000000000 0.0 20.00 5.00 0.80 0.4598773419857025 0.4598772525787354 0.89E-07 20.00 10.00 0.60 0.2146816104650497 0.2146816104650497 0.0 20.00 10.00 0.80 0.9507364630699158 0.9507364034652710 0.60E-07 20.00 20.00 0.50 0.5000000000000000 0.4999997913837433 0.21E-06 20.00 20.00 0.60 0.8979413509368896 0.8979415893554688 0.24E-06 30.00 10.00 0.70 0.2241297513246536 0.2241300046443939 0.25E-06 30.00 10.00 0.80 0.7586405277252197 0.7586401700973511 0.36E-06 40.00 20.00 0.70 0.7001783251762390 0.7001782655715942 0.60E-07 1.00 0.50 0.10 0.5131670087575912E-01 0.5131670832633972E-01 0.75E-08 1.00 0.50 0.20 0.1055728122591972 0.1055728197097778 0.75E-08 1.00 0.50 0.30 0.1633399724960327 0.1633399724960327 0.0 1.00 0.50 0.40 0.2254033237695694 0.2254033237695694 0.0 1.00 2.00 0.20 0.3600000143051147 0.3599999547004700 0.60E-07 1.00 3.00 0.20 0.4880000054836273 0.4879999756813049 0.30E-07 1.00 4.00 0.20 0.5903999805450439 0.5903999805450439 0.0 1.00 5.00 0.20 0.6723200082778931 0.6723200082778931 0.0 2.00 2.00 0.30 0.2160000056028366 0.2160000056028366 0.0 3.00 2.00 0.30 0.8370000123977661E-01 0.8370000123977661E-01 0.0 4.00 2.00 0.30 0.3078000061213970E-01 0.3077999129891396E-01 0.93E-08 5.00 2.00 0.30 0.1093500014394522E-01 0.1093499734997749E-01 0.28E-08 TEST03: BETALN evaluates the logarithm of the Gamma function. BETA_LOG_VALUES returns some exact values. X Y Exact F BETALN(X,Y) DIFF 0.2000 1.0000 1.609437942504883 1.609437942504883 0.000 0.4000 1.0000 0.9162907600402832 0.9162907004356384 0.5960E-07 0.6000 1.0000 0.5108256340026855 0.5108255743980408 0.5960E-07 0.8000 1.0000 0.2231435477733612 0.2231435477733612 0.000 1.0000 0.2000 1.609437942504883 1.609437942504883 0.000 1.0000 0.4000 0.9162907600402832 0.9162907004356384 0.5960E-07 1.0000 1.0000 0.000000000000000 -0.000000000000000 0.000 2.0000 2.0000 -1.791759490966797 -1.791759490966797 0.000 3.0000 3.0000 -3.401197433471680 -3.401197433471680 0.000 4.0000 4.0000 -4.941642284393311 -4.941642284393311 0.000 5.0000 5.0000 -6.445719718933105 -6.445719718933105 0.000 6.0000 2.0000 -3.737669706344604 -3.737669467926025 0.2384E-06 6.0000 3.0000 -5.123963832855225 -5.123964309692383 0.4768E-06 6.0000 4.0000 -6.222576141357422 -6.222576141357422 0.000 6.0000 5.0000 -7.138866901397705 -7.138866901397705 0.000 6.0000 6.0000 -7.927324295043945 -7.927324295043945 0.000 7.0000 7.0000 -9.393661499023438 -9.393661499023438 0.000 TEST04: ERF evaluates the error function. ERF_VALUES returns some exact values. X Exact F ERF(X) DIFF 0.0000 0.000000000000000 0.000000000000000 0.000 0.1000 0.1124629154801369 0.1124629154801369 0.000 0.2000 0.2227025926113129 0.2227025926113129 0.000 0.3000 0.3286267518997192 0.3286267817020416 0.2980E-07 0.4000 0.4283923506736755 0.4283923506736755 0.000 0.5000 0.5204998850822449 0.5204998850822449 0.000 0.6000 0.6038560867309570 0.6038560867309570 0.000 0.7000 0.6778011918067932 0.6778011918067932 0.000 0.8000 0.7421009540557861 0.7421009540557861 0.000 0.9000 0.7969081997871399 0.7969081997871399 0.000 1.0000 0.8427007794380188 0.8427007794380188 0.000 1.1000 0.8802050948143005 0.8802050948143005 0.000 1.2000 0.9103139638900757 0.9103139638900757 0.000 1.3000 0.9340079426765442 0.9340079426765442 0.000 1.4000 0.9522851109504700 0.9522851109504700 0.000 1.5000 0.9661051630973816 0.9661051630973816 0.000 1.6000 0.9763484001159668 0.9763484001159668 0.000 1.7000 0.9837904572486877 0.9837904572486877 0.000 1.8000 0.9890905022621155 0.9890905022621155 0.000 1.9000 0.9927904009819031 0.9927904009819031 0.000 2.0000 0.9953222870826721 0.9953222870826721 0.000 TEST05: GRAT1 evaluates the incomplete Gamma function for A <= 1. GAMMA_INC_VALUES returns some exact values. A X Exact F GAMMA_INC(A,X) DIFF 0.10 0.03 0.7382350564002991 0.7382350564002991 0.000 0.10 0.30 0.9083579778671265 0.9083579778671265 0.000 0.10 1.50 0.9886559844017029 0.9886559844017029 0.000 0.50 0.08 0.3014646470546722 0.3014646470546722 0.000 0.50 0.75 0.7793286442756653 0.7793286442756653 0.000 0.50 3.50 0.9918490052223206 0.9918490648269653 0.5960E-07 1.00 0.10 0.9516258537769318E-01 0.9516257792711258E-01 0.7451E-08 1.00 1.00 0.6321205496788025 0.6321205496788025 0.000 1.00 5.00 0.9932620525360107 0.9932620525360107 0.000 1.10 0.10 0.7205974310636520E-01 Unavailable 1.10 1.00 0.5891809463500977 Unavailable 1.10 5.00 0.9915367960929871 Unavailable 2.00 0.15 0.1018582750111818E-01 Unavailable 2.00 1.50 0.4421746134757996 Unavailable 2.00 7.00 0.9927049279212952 Unavailable 6.00 2.50 0.4202103987336159E-01 Unavailable 6.00 12.00 0.9796589612960815 Unavailable 11.00 16.00 0.9226039648056030 Unavailable 26.00 25.00 0.4470785856246948 Unavailable 41.00 45.00 0.7444549202919006 Unavailable TEST06: GAMLN evaluates the logarithm of the Gamma function. GAMMA_LOG_VALUES returns some exact values. X Exact F GAMLN(X) DIFF 0.20 1.524064183235168 1.524063825607300 0.3576E-06 0.40 0.7966780066490173 0.7966777682304382 0.2384E-06 0.60 0.3982337117195129 0.3982338309288025 0.1192E-06 0.80 0.1520599126815796 0.1520596742630005 0.2384E-06 1.00 0.000000000000000 -0.000000000000000 0.000 1.10 -0.4987246543169022E-01 -0.4987245053052902E-01 0.1490E-07 1.20 -0.8537410944700241E-01 -0.8537410944700241E-01 0.000 1.30 -0.1081747934222221 -0.1081748157739639 0.2235E-07 1.40 -0.1196128949522972 -0.1196129247546196 0.2980E-07 1.50 -0.1207822039723396 -0.1207822412252426 0.3725E-07 1.60 -0.1125917658209801 -0.1125917583703995 0.7451E-08 1.70 -0.9580771625041962E-01 -0.9580769389867783E-01 0.2235E-07 1.80 -0.7108385115861893E-01 -0.7108388096094131E-01 0.2980E-07 1.90 -0.3898428380489349E-01 -0.3898428380489349E-01 0.000 2.00 0.000000000000000 0.000000000000000 0.000 10.00 12.80182743072510 12.80182743072510 0.000 20.00 39.33988571166992 39.33988189697266 0.3815E-05 30.00 71.25704193115234 71.25704193115234 0.000 TEST07: PSI evaluates the PSI function. PSI_VALUES returns some exact values. X Exact F PSI(X) DIFF 1.00 -0.5772156715393066 -0.5772156715393066 0.000 1.10 -0.4237549304962158 -0.4237548708915710 0.5960E-07 1.20 -0.2890399098396301 -0.2890397906303406 0.1192E-06 1.30 -0.1691908836364746 -0.1691909581422806 0.7451E-07 1.40 -0.6138454377651215E-01 -0.6138456240296364E-01 0.1863E-07 1.50 -0.3648997470736504E-01 0.3648997843265533E-01 0.7298E-01 1.60 0.1260474473237991 0.1260474771261215 0.2980E-07 1.70 0.2085478752851486 0.2085478901863098 0.1490E-07 1.80 0.2849914431571960 0.2849914133548737 0.2980E-07 1.90 0.3561841547489166 0.3561841547489166 0.000 2.00 0.4227843284606934 0.4227842986583710 0.2980E-07 toms708_test(): Normal end of execution. 13 March 2024 3:56:40.589 PM