10 January 2011 3:00:06.972 PM TOMS708_PRB: FORTRAN90 version Test the TOMS708 library. TEST01 BRATIO computes the Beta ratio function. X Y W W1 = 1-W 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_INC(A,B,X). BETA_INC_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.7451E-08 0.50 0.50 0.10 0.2048327624797821 0.2048327475786209 0.1490E-07 0.50 0.50 1.00 1.000000000000000 1.000000000000000 0.000 1.00 0.50 0.00 0.000000000000000 0.000000000000000 0.000 1.00 0.50 0.01 0.5012562964111567E-02 0.5012562964111567E-02 0.000 1.00 0.50 0.10 0.5131670087575912E-01 0.5131670832633972E-01 0.7451E-08 1.00 0.50 0.50 0.2928932309150696 0.2928932011127472 0.2980E-07 1.00 1.00 0.50 0.5000000000000000 0.5000000000000000 0.000 2.00 2.00 0.10 0.2800000086426735E-01 0.2799999900162220E-01 0.1863E-08 2.00 2.00 0.20 0.1040000021457672 0.1040000021457672 0.000 2.00 2.00 0.30 0.2160000056028366 0.2160000056028366 0.000 2.00 2.00 0.40 0.3519999980926514 0.3520000576972961 0.5960E-07 2.00 2.00 0.50 0.5000000000000000 0.5000000000000000 0.000 2.00 2.00 0.60 0.6480000019073486 0.6480000615119934 0.5960E-07 2.00 2.00 0.70 0.7839999794960022 0.7839999794960022 0.000 2.00 2.00 0.80 0.8960000276565552 0.8960000276565552 0.000 2.00 2.00 0.90 0.9720000028610229 0.9720000028610229 0.000 5.50 5.00 0.50 0.4361908733844757 0.4361908435821533 0.2980E-07 10.00 0.50 0.90 0.1516409069299698 0.1516408324241638 0.7451E-07 10.00 5.00 0.50 0.8978271484375000E-01 0.8978272229433060E-01 0.7451E-08 10.00 5.00 1.00 1.000000000000000 1.000000000000000 0.000 10.00 10.00 0.50 0.5000000000000000 0.5000000000000000 0.000 20.00 5.00 0.80 0.4598773419857025 0.4598772525787354 0.8941E-07 20.00 10.00 0.60 0.2146816104650497 0.2146816104650497 0.000 20.00 10.00 0.80 0.9507364630699158 0.9507364034652710 0.5960E-07 20.00 20.00 0.50 0.5000000000000000 0.4999997913837433 0.2086E-06 20.00 20.00 0.60 0.8979413509368896 0.8979415893554688 0.2384E-06 30.00 10.00 0.70 0.2241297513246536 0.2241300046443939 0.2533E-06 30.00 10.00 0.80 0.7586405277252197 0.7586401700973511 0.3576E-06 40.00 20.00 0.70 0.7001783251762390 0.7001782655715942 0.5960E-07 1.00 0.50 0.10 0.5131670087575912E-01 0.5131670832633972E-01 0.7451E-08 1.00 0.50 0.20 0.1055728122591972 0.1055728197097778 0.7451E-08 1.00 0.50 0.30 0.1633399724960327 0.1633399724960327 0.000 1.00 0.50 0.40 0.2254033237695694 0.2254033237695694 0.000 1.00 2.00 0.20 0.3600000143051147 0.3599999547004700 0.5960E-07 1.00 3.00 0.20 0.4880000054836273 0.4879999756813049 0.2980E-07 1.00 4.00 0.20 0.5903999805450439 0.5903999805450439 0.000 1.00 5.00 0.20 0.6723200082778931 0.6723200082778931 0.000 2.00 2.00 0.30 0.2160000056028366 0.2160000056028366 0.000 3.00 2.00 0.30 0.8370000123977661E-01 0.8370000123977661E-01 0.000 4.00 2.00 0.30 0.3078000061213970E-01 0.3077999129891396E-01 0.9313E-08 5.00 2.00 0.30 0.1093500014394522E-01 0.1093499734997749E-01 0.2794E-08 TEST03: BETALN evaluates the logarithm of the Beta function. BETA_LOG_VALUES returns some exact values. X Y Exact F BETALN(X,Y) DIFF 0.20 1.00 1.609437942504883 1.609437942504883 0.000 0.40 1.00 0.9162907600402832 0.9162907004356384 0.5960E-07 0.60 1.00 0.5108256340026855 0.5108255743980408 0.5960E-07 0.80 1.00 0.2231435477733612 0.2231435477733612 0.000 1.00 0.20 1.609437942504883 1.609437942504883 0.000 1.00 0.40 0.9162907600402832 0.9162907004356384 0.5960E-07 1.00 1.00 0.000000000000000 -0.000000000000000 0.000 2.00 2.00 -1.791759490966797 -1.791759490966797 0.000 3.00 3.00 -3.401197433471680 -3.401197433471680 0.000 4.00 4.00 -4.941642284393311 -4.941642284393311 0.000 5.00 5.00 -6.445719718933105 -6.445719718933105 0.000 6.00 2.00 -3.737669706344604 -3.737669467926025 0.2384E-06 6.00 3.00 -5.123963832855225 -5.123964309692383 0.4768E-06 6.00 4.00 -6.222576141357422 -6.222576141357422 0.000 6.00 5.00 -7.138866901397705 -7.138866901397705 0.000 6.00 6.00 -7.927324295043945 -7.927324295043945 0.000 7.00 7.00 -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.00 0.000000000000000 0.000000000000000 0.000 0.10 0.1124629154801369 0.1124629154801369 0.000 0.20 0.2227025926113129 0.2227025926113129 0.000 0.30 0.3286267518997192 0.3286267817020416 0.2980E-07 0.40 0.4283923506736755 0.4283923506736755 0.000 0.50 0.5204998850822449 0.5204998850822449 0.000 0.60 0.6038560867309570 0.6038560867309570 0.000 0.70 0.6778011918067932 0.6778011918067932 0.000 0.80 0.7421009540557861 0.7421009540557861 0.000 0.90 0.7969081997871399 0.7969081997871399 0.000 1.00 0.8427007794380188 0.8427007794380188 0.000 1.10 0.8802050948143005 0.8802050948143005 0.000 1.20 0.9103139638900757 0.9103139638900757 0.000 1.30 0.9340079426765442 0.9340079426765442 0.000 1.40 0.9522851109504700 0.9522851109504700 0.000 1.50 0.9661051630973816 0.9661051630973816 0.000 1.60 0.9763484001159668 0.9763484001159668 0.000 1.70 0.9837904572486877 0.9837904572486877 0.000 1.80 0.9890905022621155 0.9890905022621155 0.000 1.90 0.9927904009819031 0.9927904009819031 0.000 2.00 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.1392920166254044 Unavailable 2.00 1.50 0.7768698334693909 Unavailable 2.00 7.00 0.9990881085395813 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_PRB: Normal end of execution. 10 January 2011 3:00:06.974 PM