11 August 2022 8:35:29.868 PM asa109_test(): FORTRAN90 version Test asa109(). test01(): xinbta() inverts the incomplete Beta function. Compare with tabulated values. A B FX X X DIFF (tabulated) (XINBTA) 0.50 0.50 0.063769 0.1000000000000000E-01 0.1000000000000001E-01 0.6939E-17 0.50 0.50 0.204833 0.1000000000000000 0.9999999999999999E-01 0.1388E-16 0.50 0.50 1.000000 1.000000000000000 1.000000000000000 0.000 1.00 0.50 0.000000 0.000000000000000 0.000000000000000 0.000 1.00 0.50 0.005013 0.1000000000000000E-01 0.9999999999999998E-02 0.1735E-17 1.00 0.50 0.051317 0.1000000000000000 0.1000000000000000 0.2776E-16 1.00 0.50 0.292893 0.5000000000000000 0.5000000000000003 0.3331E-15 1.00 1.00 0.500000 0.5000000000000000 0.5000000000000000 0.000 2.00 2.00 0.028000 0.1000000000000000 0.9999999999999999E-01 0.1388E-16 2.00 2.00 0.104000 0.2000000000000000 0.2000000000000000 0.2776E-16 2.00 2.00 0.216000 0.3000000000000000 0.3000000000000000 0.000 2.00 2.00 0.352000 0.4000000000000000 0.4000000000000000 0.000 2.00 2.00 0.500000 0.5000000000000000 0.4999999999999999 0.5551E-16 2.00 2.00 0.648000 0.6000000000000000 0.6000000000000000 0.000 2.00 2.00 0.784000 0.7000000000000000 0.7000000000000000 0.000 2.00 2.00 0.896000 0.8000000000000000 0.8000000000000000 0.000 2.00 2.00 0.972000 0.9000000000000000 0.8999999999999999 0.1110E-15 5.50 5.00 0.436191 0.5000000000000000 0.5000000000000006 0.5551E-15 10.00 0.50 0.151641 0.9000000000000000 0.9000000000000010 0.9992E-15 10.00 5.00 0.089783 0.5000000000000000 0.5000000000000009 0.8882E-15 10.00 5.00 1.000000 1.000000000000000 1.000000000000000 0.000 10.00 10.00 0.500000 0.5000000000000000 0.4999999999999998 0.1665E-15 20.00 5.00 0.459877 0.8000000000000000 0.7999999999999988 0.1221E-14 20.00 10.00 0.214682 0.6000000000000000 0.5999999999999998 0.2220E-15 20.00 10.00 0.950736 0.8000000000000000 0.8000000000000000 0.000 20.00 20.00 0.500000 0.5000000000000000 0.5000000000000002 0.2220E-15 20.00 20.00 0.897941 0.6000000000000000 0.5999999999999999 0.1110E-15 30.00 10.00 0.224130 0.7000000000000000 0.7000000000000051 0.5107E-14 30.00 10.00 0.758641 0.8000000000000000 0.7999999999999998 0.2220E-15 40.00 20.00 0.700178 0.7000000000000000 0.6999999999999980 0.1998E-14 1.00 0.50 0.051317 0.1000000000000000 0.1000000000000000 0.2776E-16 1.00 0.50 0.105573 0.2000000000000000 0.2000000000000000 0.5551E-16 1.00 0.50 0.163340 0.3000000000000000 0.3000000000000000 0.5551E-16 1.00 0.50 0.225403 0.4000000000000000 0.4000000000000000 0.000 1.00 2.00 0.360000 0.2000000000000000 0.2000000000000000 0.000 1.00 3.00 0.488000 0.2000000000000000 0.2000000000000000 0.2776E-16 1.00 4.00 0.590400 0.2000000000000000 0.2000000000000000 0.5551E-16 1.00 5.00 0.672320 0.2000000000000000 0.2000000000000000 0.5551E-16 2.00 2.00 0.216000 0.3000000000000000 0.3000000000000000 0.000 3.00 2.00 0.083700 0.3000000000000000 0.3000000000000000 0.5551E-16 4.00 2.00 0.030780 0.3000000000000000 0.3000000000000001 0.1110E-15 5.00 2.00 0.010935 0.3000000000000000 0.3000000000000000 0.000 1.31 11.76 0.918885 0.2256090000000000 0.2256089999999995 0.4996E-15 1.31 11.76 0.210530 0.3355680000000000E-01 0.3355680000000000E-01 0.000 1.31 11.76 0.182413 0.2952220000000000E-01 0.2952220000000001E-01 0.1388E-16 test02(): beta_inc_values() returns values of the incomplete Beta function. A B X BETA_INC(A,B)(X) 0.50000000 0.50000000 0.010000 0.6376856085851985E-01 0.6376856085851985E-01 0.0000E+00 0.50000000 0.50000000 0.100000 0.2048327646991335 0.2048327646991334 0.8327E-16 0.50000000 0.50000000 1.000000 1.000000000000000 1.000000000000000 0.0000E+00 1.00000000 0.50000000 0.000000 0.000000000000000 0.000000000000000 0.0000E+00 1.00000000 0.50000000 0.010000 0.5012562893380045E-02 0.5012562893380046E-02 0.8674E-18 1.00000000 0.50000000 0.100000 0.5131670194948620E-01 0.5131670194948618E-01 0.1388E-16 1.00000000 0.50000000 0.500000 0.2928932188134525 0.2928932188134523 0.2220E-15 1.00000000 1.00000000 0.500000 0.5000000000000000 0.5000000000000000 0.0000E+00 2.00000000 2.00000000 0.100000 0.2800000000000000E-01 0.2800000000000002E-01 0.2082E-16 2.00000000 2.00000000 0.200000 0.1040000000000000 0.1040000000000000 0.1388E-16 2.00000000 2.00000000 0.300000 0.2160000000000000 0.2159999999999999 0.1110E-15 2.00000000 2.00000000 0.400000 0.3520000000000000 0.3520000000000001 0.1110E-15 2.00000000 2.00000000 0.500000 0.5000000000000000 0.5000000000000000 0.0000E+00 2.00000000 2.00000000 0.600000 0.6480000000000000 0.6479999999999999 0.1110E-15 2.00000000 2.00000000 0.700000 0.7840000000000000 0.7840000000000000 0.0000E+00 2.00000000 2.00000000 0.800000 0.8960000000000000 0.8960000000000000 0.0000E+00 2.00000000 2.00000000 0.900000 0.9720000000000000 0.9720000000000000 0.0000E+00 5.50000000 5.00000000 0.500000 0.4361908850559777 0.4361908850559766 0.1110E-14 10.00000000 0.50000000 0.900000 0.1516409096347099 0.1516409096347084 0.1527E-14 10.00000000 5.00000000 0.500000 0.8978271484375000E-01 0.8978271484374982E-01 0.1804E-15 10.00000000 5.00000000 1.000000 1.000000000000000 1.000000000000000 0.0000E+00 10.00000000 10.00000000 0.500000 0.5000000000000000 0.4999999999999992 0.7772E-15 20.00000000 5.00000000 0.800000 0.4598773297575791 0.4598773297575741 0.4996E-14 20.00000000 10.00000000 0.600000 0.2146816102371739 0.2146816102371752 0.1332E-14 20.00000000 10.00000000 0.800000 0.9507364826957875 0.9507364826957873 0.2220E-15 20.00000000 20.00000000 0.500000 0.5000000000000000 0.4999999999999992 0.7772E-15 20.00000000 20.00000000 0.600000 0.8979413687105918 0.8979413687105914 0.4441E-15 30.00000000 10.00000000 0.700000 0.2241297491808366 0.2241297491808347 0.1915E-14 30.00000000 10.00000000 0.800000 0.7586405487192086 0.7586405487192094 0.7772E-15 40.00000000 20.00000000 0.700000 0.7001783247477069 0.7001783247477168 0.9881E-14 1.00000000 0.50000000 0.100000 0.5131670194948620E-01 0.5131670194948618E-01 0.1388E-16 1.00000000 0.50000000 0.200000 0.1055728090000841 0.1055728090000841 0.4163E-16 1.00000000 0.50000000 0.300000 0.1633399734659245 0.1633399734659245 0.2776E-16 1.00000000 0.50000000 0.400000 0.2254033307585166 0.2254033307585167 0.5551E-16 1.00000000 2.00000000 0.200000 0.3600000000000000 0.3600000000000000 0.0000E+00 1.00000000 3.00000000 0.200000 0.4880000000000000 0.4879999999999999 0.5551E-16 1.00000000 4.00000000 0.200000 0.5904000000000000 0.5904000000000001 0.1110E-15 1.00000000 5.00000000 0.200000 0.6723200000000000 0.6723200000000000 0.0000E+00 2.00000000 2.00000000 0.300000 0.2160000000000000 0.2159999999999999 0.1110E-15 3.00000000 2.00000000 0.300000 0.8370000000000000E-01 0.8369999999999995E-01 0.4163E-16 4.00000000 2.00000000 0.300000 0.3078000000000000E-01 0.3077999999999997E-01 0.2429E-16 5.00000000 2.00000000 0.300000 0.1093500000000000E-01 0.1093500000000000E-01 0.5204E-17 1.30625000 11.75620000 0.225609 0.9188846846205180 0.9188846846205182 0.2220E-15 1.30625000 11.75620000 0.033557 0.2105297748941900 0.2105297748941901 0.5551E-16 1.30625000 11.75620000 0.029522 0.1824130512500673 0.1824130512500674 0.8327E-16 asa109_test(): Normal end of execution. 11 August 2022 8:35:29.868 PM