22 December 2011 09:10:00 AM
TOMS179_PRB:
C++ version
Test 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
0.2 1.524063822430784 1.524063822415266 1.552e-11
0.4 0.7966778177017837 0.7966778176896321 1.215e-11
0.6 0.3982338580692348 0.3982338580596541 9.581e-12
0.8 0.1520596783998375 0.1520596783922423 7.595e-12
1 0 -1.995014743683132e-11 1.995e-11
1.1 -0.04987244125983972 -0.04987244127741915 1.758e-11
1.2 -0.08537409000331583 -0.08537409001883486 1.552e-11
1.3 -0.1081748095078604 -0.1081748095215821 1.372e-11
1.4 -0.1196129141723712 -0.1196129141845221 1.215e-11
1.5 -0.1207822376352452 -0.1207822376460281 1.078e-11
1.6 -0.1125917656967557 -0.1125917657063365 9.581e-12
1.7 -0.09580769740706586 -0.09580769741558978 8.524e-12
1.8 -0.07108387291437215 -0.07108387292196741 7.595e-12
1.9 -0.03898427592308333 -0.03898427592986082 6.777e-12
2 0 -1.995014743683132e-11 1.995e-11
3 0.6931471805599453 0.6931471805399952 1.995e-11
4 1.791759469228055 1.791759469208105 1.995e-11
10 12.80182748008147 12.80182748008065 8.225e-13
20 39.33988418719949 39.33988418719949 0
30 71.25703896716801 71.257038967168 1.421e-14
TEST02:
MDBETA estimates the value of th modified Beta function.
Compare with tabulated values.
X P Q Beta Beta DIFF
(Tabulated) (MDBETA)
0.01 0.5 0.5 0.06376856085851985 0.0637685608586229 1.03e-13
0.1 0.5 0.5 0.2048327646991335 0.2048327646994644 3.309e-13
1 0.5 0.5 1 1 0
0 1 0.5 0 0 0
0.01 1 0.5 0.005012562893380045 0.005012562893480049 1e-13
0.1 1 0.5 0.0513167019494862 0.05131670195051003 1.024e-12
0.5 1 0.5 0.2928932188134525 0.2928932188192955 5.843e-12
0.5 1 1 0.5 0.5000000000099751 9.975e-12
0.1 2 2 0.028 0.02800000000055863 5.586e-13
0.2 2 2 0.104 0.1040000000020749 2.075e-12
0.3 2 2 0.216 0.2160000000043092 4.309e-12
0.4 2 2 0.352 0.3520000000070226 7.023e-12
0.5 2 2 0.5 0.5000000000099752 9.975e-12
0.6 2 2 0.648 0.6479999999929774 7.023e-12
0.7 2 2 0.784 0.7839999999956906 4.309e-12
0.8 2 2 0.896 0.8959999999979251 2.075e-12
0.9 2 2 0.972 0.9719999999994414 5.586e-13
0.5 5.5 5 0.4361908850559777 0.4361908850689236 1.295e-11
0.9 10 0.5 0.1516409096347099 0.15164090962251 1.22e-11
0.5 10 5 0.08978271484375 0.08978271484561308 1.863e-12
1 10 5 1 1 0
0.5 10 10 0.5 0.5000000000008402 8.402e-13
0.8 20 5 0.4598773297575791 0.4598773297468096 1.077e-11
0.6 20 10 0.2146816102371739 0.2146816102365308 6.431e-13
0.8 20 10 0.9507364826957875 0.9507364826957472 4.03e-14
0.5 20 20 0.5 0.4999999999999913 8.66e-15
0.6 20 20 0.8979413687105918 0.8979413687105925 6.661e-16
0.7 30 10 0.2241297491808366 0.2241297491802104 6.262e-13
0.8 30 10 0.7586405487192086 0.7586405487190115 1.971e-13
0.7 40 20 0.7001783247477069 0.7001783247477009 5.995e-15
0.1 1 0.5 0.0513167019494862 0.05131670195051003 1.024e-12
0.2 1 0.5 0.1055728090000841 0.1055728090021903 2.106e-12
0.3 1 0.5 0.1633399734659245 0.1633399734691831 3.259e-12
0.4 1 0.5 0.2254033307585166 0.2254033307630134 4.497e-12
0.2 1 2 0.36 0.360000000007182 7.182e-12
0.2 1 3 0.488 0.4880000000097358 9.736e-12
0.2 1 4 0.5904 0.5904000000117785 1.178e-11
0.2 1 5 0.67232 0.6723200000134134 1.341e-11
0.3 2 2 0.216 0.2160000000043092 4.309e-12
0.3 3 2 0.0837 0.08370000000166977 1.67e-12
0.3 4 2 0.03078 0.03078000000061405 6.14e-13
0.3 5 2 0.010935 0.01093500000021816 2.182e-13
TOMS179_PRB:
Normal end of execution.
22 December 2011 09:10:00 AM