Mon Nov 14 18:18:42 2022 prob_test(): Python version: 3.6.9 Test prob(). angle_cdf_test(): angle_cdf() evaluates the Angle CDF PDF parameter N = 5 PDF argument X = 0.5 CDF value = 0.0107809 angle_mean_test(): angle_mean() computes the Angle mean PDF parameter N = 5 PDF mean = 1.5708 angle_pdf_test(): angle_pdf() evaluates the Angle PDF PDF parameter N = 5 PDF argument X = 0.5 PDF value = 0.0826466 anglit_cdf_test(): anglit_cdf() evaluates the Anglit CDF anglit_cdf_inv() inverts the Anglit CDF. anglit_pdf() evaluates the Anglit PDF X PDF CDF CDf_inv -0.226441 0.326423 0.281221 -0.226441 -0.544113 -0.298221 0.0570973 -0.544113 0.177702 0.908969 0.673985 0.177702 0.206213 0.931248 0.700416 0.206213 0.251653 0.960475 0.741162 0.251653 -0.266559 0.249613 0.24589 -0.266559 -0.270769 0.24145 0.242273 -0.270769 -0.557967 -0.32455 0.0508393 -0.557967 0.30811 0.985724 0.788977 0.30811 -0.35763 0.070081 0.172093 -0.35763 anglit_sample_test(): anglit_mean() computes the Anglit mean anglit_sample() samples the Anglit distribution anglit_variance() computes the Anglit variance. PDF mean = 0 PDF variance = 0.11685 Sample size = 1000 Sample mean = 0.00270002 Sample variance = 0.111213 Sample maximum = 0.72363 Sample minimum = -0.776124 arcsin_cdf_test(): arcsin_cdf() evaluates the Arcsin CDF arcsin_cdf_inv() inverts the Arcsin CDF. arcsin_pdf() evaluates the Arcsin PDF PDF parameter A = 1 X PDF CDF CDf_inv 0.655214 0.421355 0.727422 0.655214 0.982799 1.72358 0.940875 0.982799 0.278500 0.331422 0.589837 0.2785 0.823356 0.560874 0.807901 0.823356 0.906529 0.754031 0.861277 0.906529 -0.313465 0.335204 0.39851 -0.313465 0.416817 0.35018 0.636854 0.416817 0.344184 0.339023 0.611844 0.344184 -0.508854 0.369761 0.330069 -0.508854 0.354699 0.340446 0.615417 0.354699 arcsin_sample_test(): arcsin_mean() computes the Arcsin mean arcsin_sample() samples the Arcsin distribution arcsin_variance() computes the Arcsin variance. PDF parameter A = 1 PDF mean = 0 PDF variance = 0.5 Sample size = 1000 Sample mean = 0.0016445 Sample variance = 0.493507 Sample maximum = 1 Sample minimum = -0.999996 PDF parameter A = 16 PDF mean = 0 PDF variance = 128 Sample size = 1000 Sample mean = -0.230706 Sample variance = 126.924 Sample maximum = 16 Sample minimum = -16 benford_cdf_test(): benford_cdf() evaluates the Benford CDF. N CDF(N) CDF(N) by summing 1 0.30103 0.30103 2 0.477121 0.477121 3 0.60206 0.60206 4 0.69897 0.69897 5 0.778151 0.778151 6 0.845098 0.845098 7 0.90309 0.90309 8 0.954243 0.954243 9 1 1 N CDF(N) CDF(N) by summing 10 0.0413927 0.0413927 11 0.0791812 0.0791812 12 0.113943 0.113943 13 0.146128 0.146128 14 0.176091 0.176091 15 0.20412 0.20412 16 0.230449 0.230449 17 0.255273 0.255273 18 0.278754 0.278754 19 0.30103 0.30103 20 0.322219 0.322219 21 0.342423 0.342423 22 0.361728 0.361728 23 0.380211 0.380211 24 0.39794 0.39794 25 0.414973 0.414973 26 0.431364 0.431364 27 0.447158 0.447158 28 0.462398 0.462398 29 0.477121 0.477121 30 0.491362 0.491362 31 0.50515 0.50515 32 0.518514 0.518514 33 0.531479 0.531479 34 0.544068 0.544068 35 0.556303 0.556303 36 0.568202 0.568202 37 0.579784 0.579784 38 0.591065 0.591065 39 0.60206 0.60206 40 0.612784 0.612784 41 0.623249 0.623249 42 0.633468 0.633468 43 0.643453 0.643453 44 0.653213 0.653213 45 0.662758 0.662758 46 0.672098 0.672098 47 0.681241 0.681241 48 0.690196 0.690196 49 0.69897 0.69897 50 0.70757 0.70757 51 0.716003 0.716003 52 0.724276 0.724276 53 0.732394 0.732394 54 0.740363 0.740363 55 0.748188 0.748188 56 0.755875 0.755875 57 0.763428 0.763428 58 0.770852 0.770852 59 0.778151 0.778151 60 0.78533 0.78533 61 0.792392 0.792392 62 0.799341 0.799341 63 0.80618 0.80618 64 0.812913 0.812913 65 0.819544 0.819544 66 0.826075 0.826075 67 0.832509 0.832509 68 0.838849 0.838849 69 0.845098 0.845098 70 0.851258 0.851258 71 0.857332 0.857332 72 0.863323 0.863323 73 0.869232 0.869232 74 0.875061 0.875061 75 0.880814 0.880814 76 0.886491 0.886491 77 0.892095 0.892095 78 0.897627 0.897627 79 0.90309 0.90309 80 0.908485 0.908485 81 0.913814 0.913814 82 0.919078 0.919078 83 0.924279 0.924279 84 0.929419 0.929419 85 0.934498 0.934498 86 0.939519 0.939519 87 0.944483 0.944483 88 0.94939 0.94939 89 0.954243 0.954243 90 0.959041 0.959041 91 0.963788 0.963788 92 0.968483 0.968483 93 0.973128 0.973128 94 0.977724 0.977724 95 0.982271 0.982271 96 0.986772 0.986772 97 0.991226 0.991226 98 0.995635 0.995635 99 1 1 X PDF CDF CDF_INV 8 0.0511525 0.954243 8 5 0.0791812 0.778151 5 1 0.30103 0.30103 1 1 0.30103 0.30103 1 5 0.0791812 0.778151 5 1 0.30103 0.30103 1 4 0.09691 0.69897 4 8 0.0511525 0.954243 8 1 0.30103 0.30103 1 5 0.0791812 0.778151 5 benford_pdf_test(): benford_pdf() evaluates the Benford PDF. N PDF(N) 1 0.30103 2 0.176091 3 0.124939 4 0.09691 5 0.0791812 6 0.0669468 7 0.0579919 8 0.0511525 9 0.0457575 N PDF(N) 10 0.0413927 11 0.0377886 12 0.0347621 13 0.0321847 14 0.0299632 15 0.0280287 16 0.0263289 17 0.0248236 18 0.0234811 19 0.0222764 20 0.0211893 21 0.0202034 22 0.0193052 23 0.0184834 24 0.0177288 25 0.0170333 26 0.0163904 27 0.0157943 28 0.01524 29 0.0147233 30 0.0142404 31 0.0137883 32 0.013364 33 0.012965 34 0.0125891 35 0.0122345 36 0.0118992 37 0.0115819 38 0.011281 39 0.0109954 40 0.0107239 41 0.0104654 42 0.0102192 43 0.00998422 44 0.00975984 45 0.00954532 46 0.00934003 47 0.00914338 48 0.00895484 49 0.00877392 50 0.00860017 51 0.00843317 52 0.00827253 53 0.00811789 54 0.00796893 55 0.00782534 56 0.00768683 57 0.00755314 58 0.00742402 59 0.00729924 60 0.00717858 61 0.00706185 62 0.00694886 63 0.00683942 64 0.00673338 65 0.00663058 66 0.00653087 67 0.00643411 68 0.00634018 69 0.00624895 70 0.00616031 71 0.00607415 72 0.00599036 73 0.00590886 74 0.00582954 75 0.00575233 76 0.00567713 77 0.00560388 78 0.00553249 79 0.0054629 80 0.00539503 81 0.00532883 82 0.00526424 83 0.00520119 84 0.00513964 85 0.00507953 86 0.0050208 87 0.00496342 88 0.00490733 89 0.0048525 90 0.00479888 91 0.00474644 92 0.00469512 93 0.00464491 94 0.00459575 95 0.00454763 96 0.0045005 97 0.00445434 98 0.00440912 99 0.00436481 benford_sample_test(): benford_mean() computes the mean; benford_sample() samples the distribution; benford_variance() computes the variance. PDF mean = 3.440236967123206 PDF variance = 6.056512631375667 Sample size = 10000 Sample mean = 3.4166 Sample variance = 5.950044439999999 Sample maximum = 9.0 Sample minimum = 0.0 bernoulli_cdf_test(): bernoulli_cdf() evaluates the Bernoulli CDF bernoulli_cdf_inv() inverts the Bernoulli CDF. bernoulli_pdf() evaluates the Bernoulli PDF PDF parameter A = 0.75 X PDF CDF CDf_inv 1 0.75 1 1 1 0.75 1 1 1 0.75 1 1 1 0.75 1 1 1 0.75 1 1 0 0.25 0.25 0 1 0.75 1 1 1 0.75 1 1 1 0.75 1 1 0 0.25 0.25 0 bernoulli_sample_test(): bernoulli_mean() computes the Bernoulli mean bernoulli_sample() samples the Bernoulli distribution bernoulli_variance() computes the Bernoulli variance. PDF parameter A = 0.75 PDF mean = 0.75 PDF variance = 0.1875 Sample size = 1000 Sample mean = 0.755 Sample variance = 0.184975 Sample maximum = 1 Sample minimum = 0 bessel_i0_test(): bessel_i0() evaluates the Bessel function I0(X) X Exact F I0(X) 0 1 1 0.2 1.010025027795146 1.010025027795146 0.4 1.040401782229341 1.040401782229341 0.6 1.09204536431734 1.092045364317339 0.8 1.166514922869803 1.166514922869803 1 1.266065877752008 1.266065877752008 1.2 1.393725584134064 1.393725584134064 1.4 1.553395099731217 1.553395099731216 1.6 1.749980639738909 1.749980639738909 1.8 1.989559356618051 1.989559356618051 2 2.279585302336067 2.279585302336067 2.5 3.289839144050123 3.289839144050123 3 4.880792585865024 4.880792585865024 3.5 7.37820343222548 7.37820343222548 4 11.30192195213633 11.30192195213633 4.5 17.48117185560928 17.48117185560928 5 27.23987182360445 27.23987182360445 6 67.23440697647798 67.23440697647796 8 427.5641157218048 427.5641157218047 10 2815.716628466254 2815.716628466254 bessel_i0_values_test(): bessel_i0_values() stores values of the Bessel I function. of order 0. X I(0,X) 0.000000 1 0.200000 1.010025027795146 0.400000 1.040401782229341 0.600000 1.09204536431734 0.800000 1.166514922869803 1.000000 1.266065877752008 1.200000 1.393725584134064 1.400000 1.553395099731217 1.600000 1.749980639738909 1.800000 1.989559356618051 2.000000 2.279585302336067 2.500000 3.289839144050123 3.000000 4.880792585865024 3.500000 7.37820343222548 4.000000 11.30192195213633 4.500000 17.48117185560928 5.000000 27.23987182360445 6.000000 67.23440697647798 8.000000 427.5641157218048 10.000000 2815.716628466254 bessel_i1_test(): bessel_i1() evaluates the Bessel function I1(X) X Exact F I1(X) 0.000000 0 0 0.200000 0.1005008340281251 0.1005008340281251 0.400000 0.2040267557335706 0.2040267557335706 0.600000 0.3137040256049221 0.3137040256049221 0.800000 0.4328648026206398 0.4328648026206398 1.000000 0.565159103992485 0.5651591039924849 1.200000 0.7146779415526431 0.7146779415526432 1.400000 0.8860919814143274 0.8860919814143273 1.600000 1.08481063512988 1.08481063512988 1.800000 1.317167230391899 1.317167230391899 2.000000 1.590636854637329 1.590636854637329 2.500000 2.516716245288698 2.516716245288698 3.000000 3.953370217402609 3.953370217402608 3.500000 6.205834922258365 6.205834922258364 4.000000 9.759465153704451 9.759465153704447 4.500000 15.38922275373592 15.38922275373592 5.000000 24.33564214245053 24.33564214245052 6.000000 61.34193677764024 61.34193677764024 8.000000 399.8731367825601 399.8731367825602 10.000000 2670.988303701255 2670.988303701254 bessel_i1_values_test(): bessel_i1_values() stores values of the Bessel I function. of order 1. X I(1,X) 0.000000 0 0.200000 0.1005008340281251 0.400000 0.2040267557335706 0.600000 0.3137040256049221 0.800000 0.4328648026206398 1.000000 0.565159103992485 1.200000 0.7146779415526431 1.400000 0.8860919814143274 1.600000 1.08481063512988 1.800000 1.317167230391899 2.000000 1.590636854637329 2.500000 2.516716245288698 3.000000 3.953370217402609 3.500000 6.205834922258365 4.000000 9.759465153704451 4.500000 15.38922275373592 5.000000 24.33564214245053 6.000000 61.34193677764024 8.000000 399.8731367825601 10.000000 2670.988303701255 beta_binomial_cdf_test(): beta_binomial_cdf() evaluates the Beta Binomial CDF beta_binomial_cdf_inv() inverts the Beta Binomial CDF. beta_binomial_pdf() evaluates the Beta Binomial PDF PDF parameter A = 2 PDF parameter B = 3 PDF parameter C = 4 X PDF CDF CDf_inv 2 0.257143 0.757143 2 2 0.257143 0.757143 2 0 0.214286 0.214286 0 3 0.171429 0.928571 3 2 0.257143 0.757143 2 0 0.214286 0.214286 0 1 0.285714 0.5 1 3 0.171429 0.928571 3 2 0.257143 0.757143 2 4 0.0714286 1 4 beta_binomial_sample_test(): beta_binomial_mean() computes the Beta Binomial mean beta_binomial_sample() samples the Beta Binomial distribution beta_binomial_variance() computes the Beta Binomial variance. PDF parameter A = 2 PDF parameter B = 3 PDF parameter C = 4 PDF mean = 1.6 PDF variance = 1.44 Sample size = 1000 Sample mean = 1.637 Sample variance = 1.45323 Sample maximum = 4 Sample minimum = 0 beta_cdf_test(): beta_cdf() evaluates the Beta CDF beta_cdf_inv() inverts the Beta CDF. beta_pdf() evaluates the Beta PDF PDF parameter A = 12 PDF parameter B = 12 A B X PDF CDF CDf_inv 12 12 0.687948 0.724094 0.970788 0.687948 12 12 0.434089 3.18992 0.260465 0.434089 12 12 0.610761 2.22407 0.861608 0.610761 12 12 0.506538 3.86106 0.525273 0.506538 12 12 0.412013 2.73682 0.19492 0.412013 12 12 0.788918 0.04431 0.999004 0.788918 12 12 0.588781 2.71945 0.807244 0.588781 12 12 0.336401 1.1132 0.0513808 0.336401 12 12 0.692561 0.661562 0.973983 0.692561 12 12 0.487025 3.83977 0.449932 0.487025 beta_cdf_values_test(): beta_cdf_values() stores values of the Beta function. A B X beta_cdf(A,B,X) 0.500000 0.500000 0.010000 0.06376856085851985 0.500000 0.500000 0.100000 0.2048327646991335 0.500000 0.500000 1.000000 1 1.000000 0.500000 0.000000 0 1.000000 0.500000 0.010000 0.005012562893380045 1.000000 0.500000 0.100000 0.0513167019494862 1.000000 0.500000 0.500000 0.2928932188134525 1.000000 1.000000 0.500000 0.5 2.000000 2.000000 0.100000 0.028 2.000000 2.000000 0.200000 0.104 2.000000 2.000000 0.300000 0.216 2.000000 2.000000 0.400000 0.352 2.000000 2.000000 0.500000 0.5 2.000000 2.000000 0.600000 0.648 2.000000 2.000000 0.700000 0.784 2.000000 2.000000 0.800000 0.896 2.000000 2.000000 0.900000 0.972 5.500000 5.000000 0.500000 0.4361908850559777 10.000000 0.500000 0.900000 0.1516409096347099 10.000000 5.000000 0.500000 0.08978271484375 10.000000 5.000000 1.000000 1 10.000000 10.000000 0.500000 0.5 20.000000 5.000000 0.800000 0.4598773297575791 20.000000 10.000000 0.600000 0.2146816102371739 20.000000 10.000000 0.800000 0.9507364826957875 20.000000 20.000000 0.500000 0.5 20.000000 20.000000 0.600000 0.8979413687105918 30.000000 10.000000 0.700000 0.2241297491808366 30.000000 10.000000 0.800000 0.7586405487192086 40.000000 20.000000 0.700000 0.7001783247477069 1.000000 0.500000 0.100000 0.0513167019494862 1.000000 0.500000 0.200000 0.1055728090000841 1.000000 0.500000 0.300000 0.1633399734659245 1.000000 0.500000 0.400000 0.2254033307585166 1.000000 2.000000 0.200000 0.36 1.000000 3.000000 0.200000 0.488 1.000000 4.000000 0.200000 0.5904 1.000000 5.000000 0.200000 0.67232 2.000000 2.000000 0.300000 0.216 3.000000 2.000000 0.300000 0.0837 4.000000 2.000000 0.300000 0.03078 5.000000 2.000000 0.300000 0.010935 1.306250 11.756200 0.225609 0.918884684620518 1.306250 11.756200 0.033557 0.21052977489419 1.306250 11.756200 0.029522 0.1824130512500673 beta_inc_test(): beta_inc() evaluates the normalized incomplete Beta function beta_inc(A,B,X). A B X Exact F beta_inc(A,B,X) 0.5 0.5 0.01 0.0637686 0.0637686 0.5 0.5 0.1 0.204833 0.204833 0.5 0.5 1 1 1 1 0.5 0 0 0 1 0.5 0.01 0.00501256 0.00501256 1 0.5 0.1 0.0513167 0.0513167 1 0.5 0.5 0.292893 0.292893 1 1 0.5 0.5 0.5 2 2 0.1 0.028 0.028 2 2 0.2 0.104 0.104 2 2 0.3 0.216 0.216 2 2 0.4 0.352 0.352 2 2 0.5 0.5 0.5 2 2 0.6 0.648 0.648 2 2 0.7 0.784 0.784 2 2 0.8 0.896 0.896 2 2 0.9 0.972 0.972 5.5 5 0.5 0.436191 0.436191 10 0.5 0.9 0.151641 0.151641 10 5 0.5 0.0897827 0.0897827 10 5 1 1 1 10 10 0.5 0.5 0.5 20 5 0.8 0.459877 0.459877 20 10 0.6 0.214682 0.214682 20 10 0.8 0.950736 0.950736 20 20 0.5 0.5 0.5 20 20 0.6 0.897941 0.897941 30 10 0.7 0.22413 0.22413 30 10 0.8 0.758641 0.758641 40 20 0.7 0.700178 0.700178 1 0.5 0.1 0.0513167 0.0513167 1 0.5 0.2 0.105573 0.105573 1 0.5 0.3 0.16334 0.16334 1 0.5 0.4 0.225403 0.225403 1 2 0.2 0.36 0.36 1 3 0.2 0.488 0.488 1 4 0.2 0.5904 0.5904 1 5 0.2 0.67232 0.67232 2 2 0.3 0.216 0.216 3 2 0.3 0.0837 0.0837 4 2 0.3 0.03078 0.03078 5 2 0.3 0.010935 0.010935 1.30625 11.7562 0.225609 0.918885 0.918885 1.30625 11.7562 0.0335568 0.21053 0.21053 1.30625 11.7562 0.0295222 0.182413 0.182413 beta_inc_values_test(): beta_inc_values() stores values of the BETA function. A B X beta_inc(A,B,X) 0.500000 0.500000 0.010000 0.06376856085851985 0.500000 0.500000 0.100000 0.2048327646991335 0.500000 0.500000 1.000000 1 1.000000 0.500000 0.000000 0 1.000000 0.500000 0.010000 0.005012562893380045 1.000000 0.500000 0.100000 0.0513167019494862 1.000000 0.500000 0.500000 0.2928932188134525 1.000000 1.000000 0.500000 0.5 2.000000 2.000000 0.100000 0.028 2.000000 2.000000 0.200000 0.104 2.000000 2.000000 0.300000 0.216 2.000000 2.000000 0.400000 0.352 2.000000 2.000000 0.500000 0.5 2.000000 2.000000 0.600000 0.648 2.000000 2.000000 0.700000 0.784 2.000000 2.000000 0.800000 0.896 2.000000 2.000000 0.900000 0.972 5.500000 5.000000 0.500000 0.4361908850559777 10.000000 0.500000 0.900000 0.1516409096347099 10.000000 5.000000 0.500000 0.08978271484375 10.000000 5.000000 1.000000 1 10.000000 10.000000 0.500000 0.5 20.000000 5.000000 0.800000 0.4598773297575791 20.000000 10.000000 0.600000 0.2146816102371739 20.000000 10.000000 0.800000 0.9507364826957875 20.000000 20.000000 0.500000 0.5 20.000000 20.000000 0.600000 0.8979413687105918 30.000000 10.000000 0.700000 0.2241297491808366 30.000000 10.000000 0.800000 0.7586405487192086 40.000000 20.000000 0.700000 0.7001783247477069 1.000000 0.500000 0.100000 0.0513167019494862 1.000000 0.500000 0.200000 0.1055728090000841 1.000000 0.500000 0.300000 0.1633399734659245 1.000000 0.500000 0.400000 0.2254033307585166 1.000000 2.000000 0.200000 0.36 1.000000 3.000000 0.200000 0.488 1.000000 4.000000 0.200000 0.5904 1.000000 5.000000 0.200000 0.67232 2.000000 2.000000 0.300000 0.216 3.000000 2.000000 0.300000 0.0837 4.000000 2.000000 0.300000 0.03078 5.000000 2.000000 0.300000 0.010935 1.306250 11.756200 0.225609 0.918884684620518 1.306250 11.756200 0.033557 0.21052977489419 1.306250 11.756200 0.029522 0.1824130512500673 beta_sample_test(): beta_mean() computes the Beta mean beta_sample() samples the Beta distribution beta_variance() computes the Beta variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 0.4 PDF variance = 0.04 Sample size = 1000 Sample mean = 0.391272 Sample variance = 0.0373467 Sample maximum = 0.939416 Sample minimum = 0.00650115 beta_values_test(): beta_values() stores values of the Beta function. X Y BETA(X,Y) 0.200000 1.000000 5 0.400000 1.000000 2.5 0.600000 1.000000 1.666666666666667 0.800000 1.000000 1.25 1.000000 0.200000 5 1.000000 0.400000 2.5 1.000000 1.000000 1 2.000000 2.000000 0.1666666666666667 3.000000 3.000000 0.03333333333333333 4.000000 4.000000 0.007142857142857143 5.000000 5.000000 0.001587301587301587 6.000000 2.000000 0.02380952380952381 6.000000 3.000000 0.005952380952380952 6.000000 4.000000 0.001984126984126984 6.000000 5.000000 0.0007936507936507937 6.000000 6.000000 0.0003607503607503608 7.000000 7.000000 8.325008325008325e-05 binomial_cdf_test(): binomial_cdf() evaluates the Binomial CDF binomial_cdf_inv() inverts the Binomial CDF. binomial_pdf() evaluates the Binomial PDF PDF parameter A = 5 PDF parameter B = 0.65 X PDF CDF CDf_inv 4 0.312386 0.883971 4 1 0.0487703 0.0540225 1 3 0.336416 0.571585 3 4 0.312386 0.883971 4 4 0.312386 0.883971 4 3 0.336416 0.571585 3 5 0.116029 1 5 3 0.336416 0.571585 3 4 0.312386 0.883971 4 3 0.336416 0.571585 3 binomial_sample_test(): binomial_mean() computes the Binomial mean binomial_sample() samples the Binomial distribution binomial_variance() computes the Binomial variance. PDF parameter A = 5 PDF parameter B = 0.3 PDF mean = 1.5 PDF variance = 1.05 Sample size = 1000 Sample mean = 1.489 Sample variance = 1.10188 Sample maximum = 5 Sample minimum = 0 birthday_cdf_test(): birthday_cdf() evaluates the Birthday CDF birthday_cdf_inv() inverts the Birthday CDF. birthday_pdf() evaluates the Birthday PDF N PDF CDF CDF_inv 1 0 0 1 2 0.00273973 0.00273973 2 3 0.00546444 0.00820417 3 4 0.00815175 0.0163559 4 5 0.0107797 0.0271356 5 6 0.0133269 0.0404625 6 7 0.0157732 0.0562357 7 8 0.0180996 0.0743353 8 9 0.0202885 0.0946238 9 10 0.0223243 0.116948 10 11 0.0241932 0.141141 11 12 0.0258834 0.167025 12 13 0.0273855 0.19441 13 14 0.0286922 0.223103 14 15 0.0297988 0.252901 15 16 0.0307027 0.283604 16 17 0.0314037 0.315008 17 18 0.0319038 0.346911 18 19 0.0322071 0.379119 19 20 0.0323199 0.411438 20 21 0.03225 0.443688 21 22 0.032007 0.475695 22 23 0.0316019 0.507297 23 24 0.031047 0.538344 24 25 0.0303554 0.5687 25 26 0.0295411 0.598241 26 27 0.0286185 0.626859 27 28 0.0276022 0.654461 28 29 0.0265071 0.680969 29 30 0.0253477 0.706316 30 birthday_sample_test(): birthday_sample() samples the Birthday distribution. N Mean PDF 10 0.021 0.0223243 11 0.021 0.0241932 12 0.025 0.0258834 13 0.019 0.0273855 14 0.038 0.0286922 15 0.023 0.0297988 16 0.029 0.0307027 17 0.044 0.0314037 18 0.028 0.0319038 19 0.036 0.0322071 20 0.032 0.0323199 21 0.04 0.03225 22 0.041 0.032007 23 0.031 0.0316019 24 0.027 0.031047 25 0.027 0.0303554 26 0.026 0.0295411 27 0.036 0.0286185 28 0.031 0.0276022 29 0.034 0.0265071 30 0.018 0.0253477 31 0.021 0.0241384 32 0.021 0.0228929 33 0.023 0.0216243 34 0.018 0.020345 35 0.011 0.0190664 36 0.015 0.0177989 37 0.012 0.0165519 38 0.014 0.0153338 39 0.013 0.0141518 40 0.005 0.0130121 bradford_cdf_test(): bradford_cdf() evaluates the Bradford CDF bradford_cdf_inv() inverts the Bradford CDF. bradford_pdf() evaluates the Bradford PDF PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 X PDF CDF CDf_inv 1.48318 0.883453 0.646252 1.48318 1.77978 0.648045 0.869781 1.77978 1.16325 1.45263 0.287532 1.16325 1.0224 2.02779 0.0469116 1.0224 1.54926 0.817304 0.702392 1.54926 1.51711 0.848201 0.675625 1.51711 1.58319 0.787048 0.729603 1.58319 1.75424 0.663264 0.853037 1.75424 1.22694 1.28749 0.374586 1.22694 1.21289 1.32062 0.356259 1.21289 bradford_sample_test(): bradford_mean() computes the Bradford mean bradford_sample() samples the Bradford distribution bradford_variance() computes the Bradford variance. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF mean = 1.38801 PDF variance = 0.0807807 Sample size = 1000 Sample mean = 1.38162 Sample variance = 0.0837054 Sample maximum = 1.99892 Sample minimum = 1.00035 buffon_box_pdf_test(): buffon_box_pdf() evaluates the Buffon-Laplace PDF, the probability that, on a grid of cells of width A and height B, a needle of length L, dropped at random, will cross at least one grid line. A B L PDF 1 1 0 0 1 1 0.2 0.241916 1 1 0.4 0.458366 1 1 0.6 0.649352 1 1 0.8 0.814873 1 1 1 0.95493 1 2 0 0 1 2 0.2 0.18462 1 2 0.4 0.356507 1 2 0.6 0.515662 1 2 0.8 0.662085 1 2 1 0.795775 1 3 0 0 1 3 0.2 0.165521 1 3 0.4 0.322554 1 3 0.6 0.471099 1 3 0.8 0.611155 1 3 1 0.742723 1 4 0 0 1 4 0.2 0.155972 1 4 0.4 0.305577 1 4 0.6 0.448817 1 4 0.8 0.58569 1 4 1 0.716197 1 5 0 0 1 5 0.2 0.150242 1 5 0.4 0.295392 1 5 0.6 0.435448 1 5 0.8 0.570411 1 5 1 0.700282 2 1 0 0 2 1 0.2 0.18462 2 1 0.4 0.356507 2 1 0.6 0.515662 2 1 0.8 0.662085 2 1 1 0.795775 2 2 0 0 2 2 0.4 0.241916 2 2 0.8 0.458366 2 2 1.2 0.649352 2 2 1.6 0.814873 2 2 2 0.95493 2 3 0 0 2 3 0.4 0.203718 2 3 0.8 0.39046 2 3 1.2 0.560225 2 3 1.6 0.713014 2 3 2 0.848826 2 4 0 0 2 4 0.4 0.18462 2 4 0.8 0.356507 2 4 1.2 0.515662 2 4 1.6 0.662085 2 4 2 0.795775 2 5 0 0 2 5 0.4 0.173161 2 5 0.8 0.336135 2 5 1.2 0.488924 2 5 1.6 0.631527 2 5 2 0.763944 3 1 0 0 3 1 0.2 0.165521 3 1 0.4 0.322554 3 1 0.6 0.471099 3 1 0.8 0.611155 3 1 1 0.742723 3 2 0 0 3 2 0.4 0.203718 3 2 0.8 0.39046 3 2 1.2 0.560225 3 2 1.6 0.713014 3 2 2 0.848826 3 3 0 0 3 3 0.6 0.241916 3 3 1.2 0.458366 3 3 1.8 0.649352 3 3 2.4 0.814873 3 3 3 0.95493 3 4 0 0 3 4 0.6 0.213268 3 4 1.2 0.407437 3 4 1.8 0.582507 3 4 2.4 0.738479 3 4 3 0.875352 3 5 0 0 3 5 0.6 0.196079 3 5 1.2 0.376879 3 5 1.8 0.5424 3 5 2.4 0.692642 3 5 3 0.827606 4 1 0 0 4 1 0.2 0.155972 4 1 0.4 0.305577 4 1 0.6 0.448817 4 1 0.8 0.58569 4 1 1 0.716197 4 2 0 0 4 2 0.4 0.18462 4 2 0.8 0.356507 4 2 1.2 0.515662 4 2 1.6 0.662085 4 2 2 0.795775 4 3 0 0 4 3 0.6 0.213268 4 3 1.2 0.407437 4 3 1.8 0.582507 4 3 2.4 0.738479 4 3 3 0.875352 4 4 0 0 4 4 0.8 0.241916 4 4 1.6 0.458366 4 4 2.4 0.649352 4 4 3.2 0.814873 4 4 4 0.95493 4 5 0 0 4 5 0.8 0.218997 4 5 1.6 0.417623 4 5 2.4 0.595876 4 5 3.2 0.753758 4 5 4 0.891268 5 1 0 0 5 1 0.2 0.150242 5 1 0.4 0.295392 5 1 0.6 0.435448 5 1 0.8 0.570411 5 1 1 0.700282 5 2 0 0 5 2 0.4 0.173161 5 2 0.8 0.336135 5 2 1.2 0.488924 5 2 1.6 0.631527 5 2 2 0.763944 5 3 0 0 5 3 0.6 0.196079 5 3 1.2 0.376879 5 3 1.8 0.5424 5 3 2.4 0.692642 5 3 3 0.827606 5 4 0 0 5 4 0.8 0.218997 5 4 1.6 0.417623 5 4 2.4 0.595876 5 4 3.2 0.753758 5 4 4 0.891268 5 5 0 0 5 5 1 0.241916 5 5 2 0.458366 5 5 3 0.649352 5 5 4 0.814873 5 5 5 0.95493 buffon_box_sample_test(): buffon_box_sample() simulates a Buffon-Laplace needle dropping experiment. On a grid of cells of width A and height B a needle of length L is dropped at random. We count the number of times it crosses at least one grid line, and use this to estimate the value of PI. Cell width A = 1 Cell height B = 1 Needle length L = 1 Trials Hits Est(Pi) Err 10 10 3 0.141593 100 99 3.0303 0.11129 10000 9555 3.13972 0.00187523 1000000 955198 3.14071 0.00088256 buffon_pdf_test(): buffon_pdf() evaluates the Buffon PDF, the probability that, on a grid of cells of width A, a needle of length L, dropped at random, will cross at least one grid line. A L PDF 1 0 0 1 0.2 0.127324 1 0.4 0.254648 1 0.6 0.381972 1 0.8 0.509296 1 1 0.63662 2 0 0 2 0.4 0.127324 2 0.8 0.254648 2 1.2 0.381972 2 1.6 0.509296 2 2 0.63662 3 0 0 3 0.6 0.127324 3 1.2 0.254648 3 1.8 0.381972 3 2.4 0.509296 3 3 0.63662 4 0 0 4 0.8 0.127324 4 1.6 0.254648 4 2.4 0.381972 4 3.2 0.509296 4 4 0.63662 5 0 0 5 1 0.127324 5 2 0.254648 5 3 0.381972 5 4 0.509296 5 5 0.63662 buffon_sample_test(): buffon_sample() simulates a Buffon-Laplace needle dropping experiment. On a grid of cells of width A, a needle of length L is dropped at random. We count the number of times it crosses at least one grid line, and use this to estimate the value of PI. Cell width A = 1 Needle length L = 1 Trials Hits Est(Pi) Err 10 8 2.5 0.641593 100 59 3.38983 0.248238 10000 6377 3.13627 0.00532168 1000000 636960 3.13991 0.00167806 burr_cdf_test(): burr_cdf() evaluates the Burr CDF burr_cdf_inv() inverts the Burr CDF. burr_pdf() evaluates the Burr PDF PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF parameter D = 2 X PDF CDF CDf_inv 2.544 0.574393 0.53093 2.544 3.83056 0.106553 0.932 3.83056 2.42728 0.602791 0.462073 2.42728 3.51548 0.177602 0.888117 3.51548 2.03083 0.542306 0.22636 2.03083 3.96231 0.0857721 0.94462 3.96231 2.95446 0.396517 0.732434 2.95446 3.19526 0.288545 0.814596 3.19526 2.45507 0.59758 0.478751 2.45507 2.92448 0.41082 0.720333 2.92448 burr_sample_test(): burr_mean() computes the Burr mean burr_variance() computes the Burr variance burr_sample() samples the Burr distribution PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF parameter D = 2 PDF mean = 2.61227 PDF variance = 0.62513 Sample size = 1000 Sample mean = 2.662326 Sample variance = 0.737534 Sample maximum = 12.487652 Sample minimum = 1.190114 cardioid_cdf_test(): cardioid_cdf() evaluates the Cardioid CDF cardioid_cdf_inv() inverts the Cardioid CDF. cardioid_pdf() evaluates the Cardioid PDF PDF parameter A = 0 PDF parameter B = 0.25 X PDF CDF CDf_inv 0.232098 0.236599 0.555244 0.232098 0.639636 0.223001 0.649301 0.639636 0.469226 0.230132 0.610664 0.469223 0.87832 0.209961 0.701037 0.87832 0.383282 0.232958 0.590761 0.383282 1.3041 0.180127 0.784318 1.3041 1.77554 0.142975 0.860502 1.77554 -2.11074 0.118245 0.0958086 -2.11074 -0.139714 0.237957 0.466682 -0.139714 0.515405 0.228395 0.621252 0.515405 cardioid_sample_test(): cardioid_mean() computes the Cardioid mean cardioid_sample() samples the Cardioid distribution cardioid_variance() computes the Cardioid variance. PDF parameter A = 0 PDF parameter B = 0.25 PDF mean = 0 PDF variance = 0 Sample size = 1000 Sample mean = 0.0294218 Sample variance = 2.17995 Sample maximum = 3.12763 Sample minimum = -3.12797 cauchy_cdf_test(): cauchy_cdf() evaluates the Cauchy CDF cauchy_cdf_inv() inverts the Cauchy CDF. cauchy_pdf() evaluates the Cauchy PDF PDF parameter A = 2 PDF parameter B = 3 X PDF CDF CDf_inv 3.88865 0.075987 0.678847 3.88865 8.61666 0.0180926 0.864502 8.61666 12.4323 0.00810417 0.910869 12.4323 1.98421 0.1061 0.498325 1.98421 -0.148942 0.0701229 0.302141 -0.148942 1.79784 0.105624 0.478582 1.79784 -3.32673 0.0255506 0.163267 -3.32673 -0.206314 0.0688594 0.298155 -0.206314 5.89118 0.0395559 0.790937 5.89118 2.82172 0.0986985 0.5851 2.82172 cauchy_sample_test(): cauchy_mean() computes the Cauchy mean cauchy_variance() computes the Cauchy variance cauchy_sample() samples the Cauchy distribution. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 2 PDF variance = 1.79769e+308 Sample size = 1000 Sample mean = 4.32306 Sample variance = 128152 Sample maximum = 8098.97 Sample minimum = -7682.83 chebyshev1_cdf_test(): chebyshev1_cdf() evaluates the Chebyshev1 CDF chebyshev1_cdf_inv() inverts the Chebyshev1 CDF. chebyshev1_pdf() evaluates the Chebyshev1 PDF X PDF CDF CDF_inv 0.399692 0.347254 0.630883 0.399692 -0.571467 0.387887 0.306375 -0.571467 0.766742 0.495844 0.778123 0.766742 -0.888892 0.694823 0.151476 -0.888892 0.834453 0.577613 0.814216 0.834453 0.0846427 0.319456 0.526975 0.0846427 0.115486 0.320454 0.536843 0.115486 -0.997875 4.88493 0.0207563 -0.997875 -0.861297 0.626468 0.169653 -0.861297 -0.628785 0.40936 0.283553 -0.628785 chebyshev1_sample_test(): chebyshev1_mean() computes the Chebyshev1 mean chebyshev1_sample() samples the Chebyshev1 distribution chebyshev1_variance() computes the Chebyshev1 variance. PDF mean = 0 PDF variance = 0.5 Sample size = 1000 Sample mean = -0.0013102 Sample variance = 0.495622 Sample maximum = 1 Sample minimum = -1 chi_cdf_test(): chi_cdf() evaluates the Chi CDF. chi_cdf_inv() inverts the Chi CDF. chi_pdf() evaluates the Chi PDF. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 X PDF CDF CDf_inv 2.3278 0.14106 0.0683007 2.32812 4.98041 0.218074 0.734284 4.98047 4.04242 0.290264 0.49017 4.04248 3.35515 0.27655 0.291341 3.35547 6.76608 0.0519641 0.960013 6.76562 6.34945 0.0797928 0.932857 6.34961 5.3852 0.173338 0.813552 5.38477 3.6069 0.289847 0.362838 3.60693 3.73692 0.292905 0.400751 3.73682 4.06372 0.289595 0.496346 4.06348 chi_sample_test(): chi_mean() computes the Chi mean chi_variance() computes the Chi variance chi_sample() samples the Chi distribution. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF mean = 4.19154 PDF variance = 1.81408 Sample size = 1000 Sample mean = 4.17029 Sample variance = 1.72413 Sample maximum = 8.53516 Sample minimum = 1.4563 chi_square_cdf_test(): chi_square_cdf() evaluates the Chi Square CDF chi_square_cdf_inv() inverts the Chi Square CDF. chi_square_pdf() evaluates the Chi Square PDF PDF parameter A = 4 X PDF CDF CDf_inv 3.65276 0.147019 0.455034 3.07735 1.6779 0.181282 0.794728 5.9191 9.95646 0.0171406 0.0411675 0.637278 5.07203 0.100402 0.279986 2.08562 3.82892 0.141116 0.429653 2.92563 1.70956 0.181802 0.788981 5.8451 10.6742 0.0128353 0.0304805 0.539733 6.69021 0.0589683 0.153193 1.38495 7.9438 0.0374105 0.0936586 1.02266 5.38708 0.091097 0.249835 1.92166 chi_square_sample_test(): chi_square_mean() computes the Chi Square mean chi_square_sample() samples the Chi Square distribution chi_square_variance() computes the Chi Square variance. PDF parameter A = 10 PDF mean = 10 PDF variance = 20 Sample size = 1000 Sample mean = 10.0199 Sample variance = 18.6666 Sample maximum = 29.374 Sample minimum = 1.74437 chi_square_noncentral_sample_test(): chi_square_noncentral_mean() computes the Chi Square Noncentral mean. chi_square_noncentral_sample() samples the Chi Square Noncentral PDF. chi_square_noncentral_variance() computes the Chi Square Noncentral variance. PDF parameter A = 3 PDF parameter B = 2 PDF mean = 5 PDF variance = 14 Sample size = 1000 Sample mean = 5.05636 Sample variance = 13.283 Sample maximum = 21.3861 Sample minimum = 0.0332598 circular_normal_sample_test(): circular_normal_mean() computes the Circular Normal mean circular_normal_sample() samples the Circular Normal distribution circular_normal_variance() computes the Circular Normal variance. PDF means = 1 5 PDF variances = 0.5625 0.5625 Sample size = 1000 Sample mean = 0.977814 4.98644 Sample variance = 0.573 0.535419 Sample maximum = 3.36568 7.38804 Sample minimum = -2.22462 2.73736 circular_normal_01_sample_test(): circular_normal_01_mean() computes the Circular Normal 01 mean circular_normal_01_sample() samples the Circular Normal 01 distribution circular_normal_01_variance() computes the Circular Normal 01 variance. PDF means = 0 0 PDF variances = 1 1 Sample size = 1000 Sample mean = -0.0133557 0.00592203 Sample variance = 0.942748 0.879608 Sample maximum = 2.80773 2.62021 Sample minimum = -3.07971 -3.02239 cosine_cdf_test(): cosine_cdf() evaluates the Cosine CDF. cosine_cdf_inv() inverts the Cosine CDF. cosine_pdf() evaluates the Cosine PDF. PDF parameter A = 2 PDF parameter B = 1 X PDF CDF CDf_inv 1.32083 0.123838 0.291934 1.32083 1.57585 0.145052 0.366996 1.57585 1.93711 0.15884 0.479987 1.93711 2.72557 0.119067 0.721088 2.72557 2.1626 0.157056 0.551644 2.1626 0.177631 -0.0396181 0.0558159 0.177631 3.52171 0.00780936 0.901151 3.52171 2.84599 0.105518 0.753791 2.84599 2.0115 0.159144 0.503662 2.0115 2.23623 0.154735 0.574847 2.23623 cosine_sample_test(): cosine_mean() computes the Cosine mean cosine_sample() samples the Cosine distribution cosine_variance() computes the Cosine variance. PDF parameter A = 2 PDF parameter B = 1 PDF mean = 2 PDF variance = 1.28987 Sample size = 1000 Sample mean = 2.00321 Sample variance = 1.27621 Sample maximum = 4.87161 Sample minimum = -0.650719 coupon_sample_test(): coupon_sample() samples the coupon PDF. Number of coupon types is 5 Expected wait is about 8.04719 0 24 1 5 2 8 3 7 4 21 5 20 6 7 7 9 8 25 9 14 Average wait was 14 Number of coupon types is 10 Expected wait is about 23.0259 0 23 1 37 2 44 3 32 4 16 5 30 6 39 7 40 8 38 9 31 Average wait was 33 Number of coupon types is 15 Expected wait is about 40.6208 0 52 1 28 2 36 3 78 4 44 5 35 6 35 7 44 8 55 9 63 Average wait was 47 Number of coupon types is 20 Expected wait is about 59.9146 0 54 1 101 2 113 3 71 4 63 5 119 6 130 7 50 8 92 9 67 Average wait was 86 Number of coupon types is 25 Expected wait is about 80.4719 0 90 1 89 2 109 3 126 4 81 5 219 6 45 7 59 8 76 9 87 Average wait was 98.1 coupon_complete_pdf_test(): coupon_complete_pdf() evaluates the Coupon Complete PDF. Number of coupon types is 2 BOX_NUM PDF CDF 1 0 0 2 0.5 0.5 3 0.25 0.75 4 0.125 0.875 5 0.0625 0.9375 6 0.03125 0.96875 7 0.015625 0.984375 8 0.0078125 0.992188 9 0.00390625 0.996094 10 0.00195312 0.998047 11 0.000976562 0.999023 12 0.000488281 0.999512 13 0.000244141 0.999756 14 0.00012207 0.999878 15 6.10352e-05 0.999939 16 3.05176e-05 0.999969 17 1.52588e-05 0.999985 18 7.62939e-06 0.999992 19 3.8147e-06 0.999996 20 1.90735e-06 0.999998 Number of coupon types is 3 BOX_NUM PDF CDF 1 0 0 2 0 0 3 0.222222 0.222222 4 0.222222 0.444444 5 0.17284 0.617284 6 0.123457 0.740741 7 0.085048 0.825789 8 0.0576132 0.883402 9 0.0387136 0.922116 10 0.0259107 0.948026 11 0.0173077 0.965334 12 0.0115497 0.976884 13 0.00770358 0.984587 14 0.00513698 0.989724 15 0.00342507 0.993149 16 0.00228352 0.995433 17 0.00152239 0.996955 18 0.00101494 0.99797 19 0.000676634 0.998647 20 0.000451091 0.999098 Number of coupon types is 4 BOX_NUM PDF CDF 1 0 0 2 0 0 3 0 0 4 0.09375 0.09375 5 0.140625 0.234375 6 0.146484 0.380859 7 0.131836 0.512695 8 0.110229 0.622925 9 0.0884399 0.711365 10 0.0692368 0.780602 11 0.0533867 0.833988 12 0.040771 0.874759 13 0.0309441 0.905703 14 0.0233911 0.929094 15 0.0176349 0.946729 16 0.0132719 0.960001 17 0.00997682 0.969978 18 0.00749406 0.977472 19 0.00562627 0.983098 20 0.00422256 0.987321 deranged_cdf_test(): deranged_cdf() evaluates the Deranged CDF deranged_cdf_inv() inverts the Deranged CDF. deranged_pdf() evaluates the Deranged PDF PDF parameter A = 7 X PDF CDF CDf_inv 0 0.367857 0.367857 0 1 0.368056 0.735913 1 2 0.183333 0.919246 2 3 0.0625 0.981746 3 4 0.0138889 0.995635 4 5 0.00416667 0.999802 5 6 0 0.999802 5 7 0.000198413 1 7 deranged_sample_test(): deranged_mean() computes the Deranged mean. deranged_variance() computes the Deranged variance. deranged_sample() samples the Deranged distribution. PDF parameter A = 7 PDF mean = 1 PDF variance = 1 Sample size = 1000 Sample mean = 1.016 Sample variance = 0.991744 Sample maximum = 5 Sample minimum = 0 digamma_test(): digamma() computes the Digamma or Psi function. Compare the result to tabulated values. X FX FX2 (Tabulated) (DIGAMMA) DIFF 0.1 -10.42375494041108 -10.42375494041114 5.862e-14 0.2 -5.289039896592188 -5.289039896592243 5.507e-14 0.3 -3.502524222200133 -3.502524222200181 4.796e-14 0.4 -2.561384544585116 -2.561384544585158 4.174e-14 0.5 -1.963510026021423 -1.963510026021564 1.406e-13 0.6 -1.54061921389319 -1.540619213893313 1.232e-13 0.7 -1.220023553697935 -1.220023553698041 1.064e-13 0.8 -0.9650085667061385 -0.9650085667062314 9.281e-14 0.9 -0.7549269499470515 -0.7549269499471327 8.127e-14 1 -0.5772156649015329 -0.5772156649016036 7.072e-14 1.1 -0.4237549404110768 -0.4237549404111393 6.251e-14 1.2 -0.2890398965921883 -0.2890398965922431 5.479e-14 1.3 -0.1691908888667997 -0.1691908888668481 4.841e-14 1.4 -0.06138454458511615 -0.06138454458515841 4.226e-14 1.5 0.03648997397857652 0.03648997397843547 1.411e-13 1.6 0.1260474527734763 0.1260474527733536 1.227e-13 1.7 0.208547874873494 0.2085478748733869 1.071e-13 1.8 0.2849914332938615 0.2849914332937686 9.293e-14 1.9 0.3561841611640597 0.3561841611639783 8.143e-14 2 0.4227843350984671 0.4227843350983961 7.094e-14 dipole_cdf_test(): dipole_cdf() evaluates the Dipole CDF. dipole_cdf_inv() inverts the Dipole CDF. dipole_pdf() evaluates the Dipole PDF. PDF parameter A = 0 PDF parameter B = 1 X PDF CDF CDf_inv -2.26813 -0.020792 0.0146794 -2.27148 -0.3006 0.626523 0.319297 -0.300781 -1.27195 0.159385 0.0574203 -1.27197 -2.87702 -0.00334545 0.00776774 -2.86719 -0.822884 0.412051 0.124653 -0.823242 -1.93886 -0.00848285 0.0218957 -1.94141 0.991511 0.308941 0.907792 0.991211 1.3588 0.122505 0.950007 1.3584 -0.927425 0.347793 0.103276 -0.927734 -0.187832 0.634902 0.383149 -0.187988 PDF parameter A = 0.785398 PDF parameter B = 0.5 X PDF CDF CDf_inv -1.66273 0.0811375 0.151215 -1.66309 0.431053 0.305475 0.56244 0.431152 -7.67 0.00544302 0.0399378 -7.68359 -1.25051 0.142253 0.183674 -1.25 -0.630046 0.276783 0.264075 -0.629883 1.07084 0.179398 0.723814 1.07129 -35.429 0.000256908 0.0089187 -35.5625 0.409778 0.307488 0.555657 0.409668 1.0047 0.194261 0.711145 1.00464 0.589772 0.283929 0.610575 0.589844 PDF parameter A = 1.5708 PDF parameter B = 0 X PDF CDF CDf_inv -3.89136 0.0197185 0.0800666 -3.89258 0.65308 0.223138 0.684154 0.65332 0.768664 0.200089 0.708601 0.768555 0.981222 0.162172 0.746983 0.981445 3.51189 0.0238732 0.911699 3.51172 -1.73529 0.0793543 0.166409 -1.73633 -0.3634 0.281178 0.389049 -0.363281 5.19273 0.0113827 0.939442 5.1875 0.96217 0.16529 0.743864 0.961914 0.949751 0.167353 0.741798 0.949219 dipole_sample_test(): dipole_sample() samples the Dipole distribution. PDF parameter A = 0 PDF parameter B = 1 Sample size = 1000 Sample mean = 0.0320504 Sample variance = 0.980003 Sample maximum = 9.15108 Sample minimum = -4.80672 PDF parameter A = 0.785398 PDF parameter B = 0.5 Sample size = 1000 Sample mean = 2.81749 Sample variance = 10979.1 Sample maximum = 3200.95 Sample minimum = -799.322 PDF parameter A = 1.5708 PDF parameter B = 0 Sample size = 1000 Sample mean = -1.69424 Sample variance = 1017.18 Sample maximum = 434.757 Sample minimum = -744.06 dirichlet_pdf_test(): dirichlet_pdf() evaluates the Dirichlet PDF. Number of components N = 3 PDF parameters A: 0: 0.25 1: 0.5 2: 1.25 PDF arguments X: 0: 0.5 1: 0.125 2: 0.375 PDF value = 0.63907 dirichlet_sample_test(): dirichlet_sample() samples the Dirichlet distribution dirichlet_mean() computes the Dirichlet mean dirichlet_variance() computes the Dirichlet variance. Number of components N = 3 PDF parameters A: 0: 0.25 1: 0.5 2: 1.25 PDF mean, variance: 0 0.125 0.0364583 1 0.25 0.0625 2 0.625 0.078125 Second moment matrix: Col: 0 1 2 Row 0 : 0.0520833 0.0208333 0.0520833 1 : 0.0208333 0.125 0.104167 2 : 0.0520833 0.104167 0.46875 Sample size = 1000 Observed Min, Max, Mean, Variance: 0 1.17746e-13 0.980039 0.115903 0.0337604 1 2.25497e-07 0.974528 0.25084 0.0652326 2 0.00561707 0.999998 0.633258 0.0788432 dirichlet_mix_pdf_test(): dirichlet_mix_pdf() evaluates the Dirichlet Mix PDF. Number of elements ELEM_NUM = 3 Number of components COMP_NUM = 2 PDF parameters A(ELEM,COMP): Col: 0 1 Row 0 : 0.25 1.5 1 : 0.5 0.5 2 : 1.25 2 Component weights: 0: 1 1: 2 PDF value = 2.12288 dirichlet_mix_sample_test(): dirichlet_mix_sample() samples the Dirichlet Mix distribution dirichlet_mix_mean() computes the Dirichlet Mix mean Number of elements ELEM_NUM = 3 Number of components COMP_NUM = 2 PDF parameters A(ELEM,COMP): Col: 0 1 Row 0 : 0.25 1.5 1 : 0.5 0.5 2 : 1.25 2 Component weights: 0: 1 1: 2 PDF mean: 0: 0.291667 1: 0.166667 2: 0.541667 Sample size = 1000 Observed Min, Max, Mean, Variance: 0 3.14711e-13 0.937741 0.290568 0.0571882 1 1.3711e-07 0.976167 0.166736 0.0400228 2 0.0102166 0.999704 0.542696 0.0645849 discrete_cdf_test(): discrete_cdf() evaluates the Discrete CDF discrete_cdf_inv() inverts the Discrete CDF. discrete_pdf() evaluates the Discrete PDF PDF parameter A = 6 PDF parameters B: 0: 1 1: 2 2: 6 3: 2 4: 4 5: 1 X PDF CDF CDf_inv 3 0.375 0.5625 3 3 0.375 0.5625 3 3 0.375 0.5625 3 2 0.125 0.1875 2 5 0.25 0.9375 5 4 0.125 0.6875 4 3 0.375 0.5625 3 2 0.125 0.1875 2 3 0.375 0.5625 3 3 0.375 0.5625 3 discrete_sample_test(): discrete_mean() computes the Discrete mean discrete_sample() samples the Discrete distribution discrete_variance() computes the Discrete variance. PDF parameter A = 6 PDF parameters B: 0: 1 1: 2 2: 6 3: 2 4: 4 5: 1 PDF mean = 3.5625 PDF variance = 1.74609 Sample size = 1000 Sample mean = 3.535 Sample variance = 1.74678 Sample maximum = 6 Sample minimum = 1 disk_sample_test(): disk_mean() returns the Disk mean. disk_sample() samples the Disk distribution. disk_variance() returns the Disk variance. X coordinate of center is A = 10 Y coordinate of center is B = 4 Radius is C = 5 Disk mean = 10 4 Disk variance = 12.5 Sample size = 1000 Sample mean = 9.98451 4.07293 Sample variance = 12.5211 Sample maximum = 14.8677 8.88644 Sample minimum = 5.03501 -0.905225 empirical_discrete_cdf_test(): empirical_discrete_cdf() evaluates the Empirical Discrete CDF empirical_discrete_cdf_inv() inverts the Empirical Discrete CDF. empirical_discrete_pdf() evaluates the Empirical Discrete PDF PDF parameter A = 6 PDF parameter B: 0: 1 1: 1 2: 3 3: 2 4: 1 5: 2 PDF parameter C: 0: 0 1: 1 2: 2 3: 4.5 4: 6 5: 10 X PDF CDF CDf_inv 1 0.1 0.2 1 10 0.2 1 10 2 0.3 0.5 2 2 0.3 0.5 2 2 0.3 0.5 2 2 0.3 0.5 2 2 0.3 0.5 2 10 0.2 1 10 10 0.2 1 10 10 0.2 1 10 empirical_discrete_sample_test(): empirical_discrete_mean() computes the Empirical Discrete mean empirical_discrete_sample() samples the Empirical Discrete distribution empirical_discrete_variance() computes the Empirical Discrete variance. PDF parameter A = 6 PDF parameter B: 0: 1 1: 1 2: 3 3: 2 4: 1 5: 2 PDF parameter C: 0: 0 1: 1 2: 2 3: 4.5 4: 6 5: 10 PDF mean = 4.2 PDF variance = 11.31 Sample size = 1000 Sample mean = 4.2605 Sample variance = 11.4214 Sample maximum = 10 Sample minimum = 0 english_letter_cdf_test(): english_letter_cdf() evaluates the English Letter CDF english_letter_cdf_inv() inverts the English Letter CDF. english_letter_pdf() evaluates the English Letter PDF C PDF CDF CDf_inv 'm' 0.02406 0.54055 'm' 'e' 0.12702 0.29396 'e' 'a' 0.08167 0.08167 'a' 'f' 0.02228 0.31624 'f' 'h' 0.06094 0.39733 'h' 'd' 0.04253 0.16694 'd' 'o' 0.07507 0.68311 'o' 'y' 0.01974 0.99926 'y' 'n' 0.06749 0.60804 'n' 'b' 0.01492 0.09659 'b' english_sentence_length_cdf_test(): english_sentence_length_cdf() evaluates the English Sentence Length CDF english_sentence_length_cdf_inv() inverts the English Sentence Length CDF. english_sentence_length_pdf() evaluates the English Sentence Length PDF X PDF CDF CDf_inv 5 0.0305008 0.0965039 5 24 0.0253187 0.726295 24 13 0.0367753 0.378037 13 24 0.0253187 0.726295 24 5 0.0305008 0.0965039 5 41 0.00625451 0.947468 41 21 0.0287367 0.647141 21 11 0.0357028 0.303294 11 12 0.0379681 0.341262 12 18 0.0331269 0.553935 18 english_sentence_length_sample_test(): english_sentence_length_mean() computes the English Sentence Length mean english_sentence_length_sample() samples the English Sentence Length distribution english_sentence_length_variance() computes the English Sentence Length variance. PDF mean = 19.1147 PDF variance = 147.443 Sample size = 1000 Sample mean = 19.59 Sample variance = 148.326 Sample maximum = 69 Sample minimum = 1 english_word_length_cdf_test(): english_word_length_cdf() evaluates the English Word Length CDF english_word_length_cdf_inv() inverts the English Word Length CDF. english_word_length_pdf() evaluates the English Word Length PDF X PDF CDF CDf_inv 3 0.211926 0.413282 3 6 0.0852426 0.763833 6 6 0.0852426 0.763833 6 5 0.108523 0.67859 5 4 0.156785 0.570067 4 4 0.156785 0.570067 4 3 0.211926 0.413282 3 3 0.211926 0.413282 3 3 0.211926 0.413282 3 4 0.156785 0.570067 4 english_word_length_sample_test(): english_word_length_mean() computes the English Word Length mean english_word_length_sample() samples the English Word Length distribution english_word_length_variance() computes the English Word Length variance. PDF mean = 4.73912 PDF variance = 7.05635 Sample size = 1000 Sample mean = 4.785 Sample variance = 7.48078 Sample maximum = 17 Sample minimum = 1 erlang_cdf_test(): erlang_cdf() evaluates the Erlang CDF. erlang_cdf_inv() inverts the Erlang CDF. erlang_pdf() evaluates the Erlang PDF. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 X PDF CDF CDf_inv 3.33378 0.105981 0.113414 3.33398 5.88579 0.129662 0.441457 5.88574 3.77306 0.120126 0.163256 3.77344 6.64618 0.118396 0.536032 6.64648 13.6837 0.0177069 0.951657 13.6875 7.64295 0.0995643 0.644868 7.64258 8.2496 0.0875544 0.70162 8.25 5.76287 0.131031 0.425433 5.7627 13.4744 0.0190171 0.947814 13.4766 7.3769 0.104803 0.617681 7.37695 erlang_sample_test(): erlang_mean() computes the Erlang mean erlang_sample() samples the Erlang distribution erlang_variance() computes the Erlang variance. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF mean = 7 PDF variance = 12 Sample size = 1000 Sample mean = 6.9918 Sample variance = 11.9041 Sample maximum = 22.0464 Sample minimum = 1.27453 exponential_cdf_test(): exponential_cdf() evaluates the Exponential CDF. exponential_cdf_inv() inverts the Exponential CDF. exponential_pdf() evaluates the Exponential PDF. PDF parameter A = 1 PDF parameter B = 2 X PDF CDF CDf_inv 7.2665 0.0217879 0.956424 7.2665 1.68001 0.355883 0.288234 1.68001 6.11779 0.038695 0.92261 6.11779 2.97851 0.185926 0.628147 2.97851 2.00626 0.302317 0.395366 2.00626 1.55495 0.378846 0.242307 1.55495 3.59261 0.136771 0.726459 3.59261 1.77703 0.339032 0.321936 1.77703 2.86914 0.196378 0.607245 2.86914 2.94604 0.18897 0.62206 2.94604 exponential_sample_test(): exponential_mean() computes the Exponential mean exponential_sample() samples the Exponential distribution exponential_variance() computes the Exponential variance. PDF parameter A = 1 PDF parameter B = 10 PDF mean = 11 PDF variance = 100 Sample size = 1000 Sample mean = 11.627 Sample variance = 110.352 Sample maximum = 93.5184 Sample minimum = 1.03536 exponential_01_cdf_test(): exponential_01_cdf() evaluates the Exponential 01 CDF. exponential_01_cdf_inv() inverts the Exponential 01 CDF. exponential_01_pdf() evaluates the Exponential 01 PDF. X PDF CDF CDf_inv 0.140854 0.868616 0.131384 0.140854 0.102271 0.902785 0.0972149 0.102271 2.66839 0.0693637 0.930636 2.66839 0.830984 0.435621 0.564379 0.830984 0.112361 0.893721 0.106279 0.112361 0.97517 0.377128 0.622872 0.97517 0.495952 0.608991 0.391009 0.495952 1.09272 0.335304 0.664696 1.09272 1.04495 0.351709 0.648291 1.04495 0.114016 0.892243 0.107757 0.114016 exponential_01_sample_test(): exponential_01_mean() computes the Exponential 01 mean exponential_01_sample() samples the Exponential 01 distribution exponential_01_variance() computes the Exponential 01 variance. PDF mean = 1 PDF variance = 1 Sample size = 1000 Sample mean = 1.00495 Sample variance = 0.911533 Sample maximum = 5.68517 Sample minimum = 0.000609758 extreme_values_cdf_test(): extreme_values_cdf() evaluates the Extreme Values CDF extreme_values_cdf_inv() inverts the Extreme Values CDF. extreme_values_pdf() evaluates the Extreme Values PDF PDF parameter A = 2 PDF parameter B = 3 X PDF CDF CDf_inv 3.64835 0.108032 0.56143 3.64835 8.2098 0.0370777 0.881444 8.2098 2.11896 0.122531 0.382464 2.11896 7.73828 0.0424663 0.862714 7.73828 -0.15185 0.0880202 0.128882 -0.15185 10.8583 0.0165135 0.949144 10.8583 2.0582 0.122604 0.375015 2.0582 0.408583 0.103532 0.182732 0.408583 3.09769 0.115546 0.499786 3.09769 2.41862 0.121492 0.419052 2.41862 extreme_values_sample_test(): extreme_values_mean() computes the Extreme Values mean extreme_values_sample() samples the Extreme Values distribution extreme_values_variance() computes the Extreme Values variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 3.73165 PDF variance = 14.8044 Sample size = 1000 Sample mean = 3.86103 Sample variance = 14.8069 Sample maximum = 29.3425 Sample minimum = -4.79169 f_cdf_test(): f_cdf() evaluates the F CDF. f_pdf() evaluates the F PDF. PDF parameter M = 1 PDF parameter N = 1 X M N PDF CDF 3.25561 1 1 0.0414545 0.32218 0.0147142 1 1 2.58606 0.923152 7.01961 1 1 0.014981 0.229759 12.39 1 1 0.00675357 0.176218 0.0724312 1 1 1.10285 0.832632 27.9547 1 1 0.00207923 0.119002 0.959701 1 1 0.165803 0.506546 0.397339 1 1 0.361383 0.641942 0.0157474 1 1 2.49724 0.920527 174.946 1 1 0.000136779 0.04804 f_sample_test(): f_mean() computes the F mean f_sample() samples the F distribution f_variance() computes the F variance. PDF parameter M = 8 PDF parameter N = 6 PDF mean = 1.5 PDF variance = 3.375 Sample size = 1000 Sample mean = 1.49554 Sample variance = 4.08398 Sample maximum = 36.3596 Sample minimum = 0.0684931 fermi_dirac_sample_test(): fermi_dirac_sample() samples the Fermi Dirac distribution. U = 1 V = 1 SAMPLE_NUM = 10000 Sample mean = 0.592572 Sample variance = 0.173831 Maximum value = 2.69086 Minimum value = 0.000195307 U = 2 V = 1 SAMPLE_NUM = 10000 Sample mean = 1.06125 Sample variance = 0.438953 Maximum value = 3.69389 Minimum value = 4.55204e-05 U = 4 V = 1 SAMPLE_NUM = 10000 Sample mean = 2.02179 Sample variance = 1.40192 Maximum value = 5.60152 Minimum value = 1.21596e-05 U = 8 V = 1 SAMPLE_NUM = 10000 Sample mean = 3.99129 Sample variance = 5.43722 Maximum value = 9.31369 Minimum value = 0.00241748 U = 16 V = 1 SAMPLE_NUM = 10000 Sample mean = 7.99045 Sample variance = 21.4382 Maximum value = 17.1659 Minimum value = 0.000756147 U = 32 V = 1 SAMPLE_NUM = 10000 Sample mean = 16.1337 Sample variance = 84.5561 Maximum value = 33.1192 Minimum value = 0.00194394 U = 1 V = 0.25 SAMPLE_NUM = 10000 Sample mean = 0.508274 Sample variance = 0.0884134 Maximum value = 1.37074 Minimum value = 1.904e-05 fisher_pdf_test(): fisher_pdf() evaluates the Fisher PDF. PDF Input: Concentration parameter KAPPA = 0 1 0 0 X PDF 0.108789 0.62696 -0.771418 0.0795775 0.877535 -0.0506549 0.476829 0.0795775 -0.867354 -0.0197346 -0.497301 0.0795775 -0.245018 0.144306 0.958719 0.0795775 0.776817 0.623681 0.08705 0.0795775 -0.31794 0.643904 0.695918 0.0795775 0.991214 0.100396 0.0861128 0.0795775 0.67943 -0.709667 0.186407 0.0795775 -0.826565 0.376404 -0.418461 0.0795775 0.413406 0.889979 0.192438 0.0795775 PDF Input: Concentration parameter KAPPA = 0.5 1 0 0 X PDF -0.750959 0.3069 0.584699 0.0524535 0.971186 0.0901376 0.220619 0.124089 0.381316 0.921868 0.068975 0.0923942 0.0980974 -0.987524 0.123177 0.0801945 -0.0399677 -0.0974652 -0.994436 0.0748452 0.931799 -0.263412 0.24973 0.121669 0.618667 -0.726182 0.299851 0.104036 0.770173 -0.166109 -0.615826 0.112224 -0.827006 -0.393348 -0.401671 0.0504965 0.669679 0.00874967 0.742599 0.106724 PDF Input: Concentration parameter KAPPA = 10 1 0 0 X PDF 0.955629 0.294509 0.00620233 1.02122 0.944009 -0.0559087 0.325149 0.909188 0.910227 0.408625 -0.067174 0.648547 0.824394 -0.088331 -0.559082 0.274898 0.764428 0.612644 0.20079 0.150919 0.856535 -0.00727419 0.516037 0.379106 0.96797 0.191738 0.162085 1.15535 0.998063 -0.0314766 -0.053663 1.56102 0.87657 0.407667 0.255799 0.463201 0.932166 -0.276898 -0.233226 0.807644 fisk_cdf_test(): fisk_cdf() evaluates the Fisk CDF fisk_cdf_inv() inverts the Fisk CDF. fisk_pdf() evaluates the Fisk PDF PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 X PDF CDF CDf_inv 2.92442 0.38843 0.471139 2.92442 3.40592 0.288957 0.635146 3.40592 1.77584 0.201512 0.0551562 1.77584 2.32762 0.395654 0.226308 2.32762 2.82376 0.403464 0.431252 2.82376 2.49923 0.417293 0.296381 2.49923 2.63723 0.419164 0.354248 2.63723 2.81688 0.404349 0.428473 2.81688 2.91174 0.39052 0.466202 2.91174 2.03305 0.309125 0.121116 2.03305 fisk_sample_test(): fisk_mean() computes the Fisk mean fisk_sample() samples the Fisk distribution fisk_variance() computes the Fisk variance. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF mean = 3.4184 PDF variance = 3.82494 Sample size = 1000 Sample mean = 3.45019 Sample variance = 3.6116 Sample maximum = 30.2612 Sample minimum = 1.15328 folded_normal_cdf_test(): folded_normal_cdf() evaluates the Folded Normal CDF. folded_normal_cdf_inv() inverts the Folded Normal CDF. folded_normal_pdf() evaluates the Folded Normal PDF. PDF parameter A = 2 PDF parameter B = 3 X PDF CDF CDf_inv 8.2004 0.0161217 0.980287 8.20463 7.8303 0.0207402 0.973493 7.83104 0.566051 0.210865 0.120153 0.566051 0.285054 0.212432 0.060656 0.285064 1.01654 0.206236 0.214199 1.01656 5.21183 0.0823654 0.849718 5.21202 2.09625 0.185266 0.426737 2.09673 2.49046 0.174594 0.497713 2.49027 2.54713 0.172949 0.50756 2.54768 0.690859 0.209841 0.146408 0.690772 folded_normal_sample_test(): folded_normal_mean() computes the Folded Normal mean folded_normal_sample() samples the Folded Normal distribution folded_normal_variance() computes the Folded Normal variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 2.90672 PDF variance = 4.55099 Sample size = 1000 Sample mean = 2.88644 Sample variance = 4.41072 Sample maximum = 10.1893 Sample minimum = 0.000976889 frechet_cdf_test(): frechet_cdf() evaluates the Frechet CDF frechet_cdf_inv() inverts the Frechet CDF. frechet_pdf() evaluates the Frechet PDF PDF parameter ALPHA = 3 X PDF CDF CDf_inv 0.943732 1.15087 0.304299 0.943732 1.78414 0.248277 0.838549 1.78414 1.29984 0.666518 0.634236 1.29984 0.718072 0.757698 0.0671501 0.718072 2.17865 0.120886 0.907826 2.17865 1.01246 1.08932 0.381539 1.01246 0.759617 0.92026 0.102134 0.759617 1.20244 0.807359 0.562599 1.20244 0.923294 1.15873 0.280687 0.923294 1.03555 1.06012 0.406362 1.03555 frechet_sample_test(): frechet_mean() computes the Frechet mean frechet_sample() samples the Frechet distribution frechet_variance() computes the Frechet variance. PDF parameter ALPHA = 3 PDF mean = 1.35412 PDF variance = 0.845303 Sample size = 1000 Sample mean = 1.34273 Sample variance = 0.73649 Sample maximum = 12.2454 Sample minimum = 0.489784 gamma_cdf_test(): gamma_cdf() evaluates the Gamma CDF. gamma_pdf() evaluates the Gamma PDF. PDF parameter A = 1 PDF parameter B = 1.5 PDF parameter C = 3 X PDF CDF 5.02225 0.164087 0.501833 4.21423 0.179571 0.361917 2.09912 0.0860123 0.0382466 2.70304 0.138059 0.106807 2.56906 0.128141 0.0889623 4.7821 0.170268 0.461665 6.96413 0.098855 0.758371 10.186 0.0273731 0.943341 3.44631 0.173556 0.224649 3.03821 0.158151 0.156638 gamma_sample_test(): gamma_mean() computes the Gamma mean gamma_sample() samples the Gamma distribution gamma_variance() computes the Gamma variance. TEST NUMBER: 0 PDF parameter A = 1 PDF parameter B = 3 PDF parameter C = 2 PDF mean = 7 PDF variance = 18 Sample size = 1000 Sample mean = 7.17728 Sample variance = 19.1559 Sample maximum = 30.0917 Sample minimum = 1.14263 TEST NUMBER: 1 PDF parameter A = 2 PDF parameter B = 0.5 PDF parameter C = 0.5 PDF mean = 2.25 PDF variance = 0.125 Sample size = 1000 Sample mean = 2.23483 Sample variance = 0.111478 Sample maximum = 4.91033 Sample minimum = 2 gamma_inc_values_test(): gamma_inc_values() stores values of the incomplete Gamma function. A X gamma_inc(A,X) 0.100000 0.030000 2.49030283630057 0.100000 0.300000 0.8718369702247978 0.100000 1.500000 0.1079213896175866 0.500000 0.075000 1.238121685818417 0.500000 0.750000 0.3911298052193973 0.500000 3.500000 0.01444722098952533 1.000000 0.100000 0.9048374180359596 1.000000 1.000000 0.3678794411714423 1.000000 5.000000 0.006737946999085467 1.100000 0.100000 0.8827966752611692 1.100000 1.000000 0.3908330082003269 1.100000 5.000000 0.008051456628620992 2.000000 0.150000 0.9898141728888165 2.000000 1.500000 0.5578254003710746 2.000000 7.000000 0.00729505572443613 6.000000 2.500000 114.9574754165633 6.000000 12.000000 2.440923530031405 11.000000 16.000000 280854.6620274718 26.000000 25.000000 8.576480283455533e+24 41.000000 45.000000 2.085031346403364e+47 geometric_cdf_test(): geometric_cdf() evaluates the Geometric CDF geometric_cdf_inv() inverts the Geometric CDF. geometric_pdf() evaluates the Geometric PDF PDF parameter A = 0.25 X PDF CDF CDf_inv 3 0.140625 0.578125 4 1 0.25 0.25 2 1 0.25 0.25 2 5 0.0791016 0.762695 6 2 0.1875 0.4375 3 8 0.033371 0.899887 9 1 0.25 0.25 2 2 0.1875 0.4375 3 4 0.105469 0.683594 5 8 0.033371 0.899887 9 geometric_sample_test(): geometric_mean() computes the Geometric mean geometric_sample() samples the Geometric distribution geometric_variance() computes the Geometric variance. PDF parameter A = 0.25 PDF mean = 4 PDF variance = 12 Sample size = 1000 Sample mean = 4.044 Sample variance = 12.3141 Sample maximum = 26 Sample minimum = 1 gompertz_cdf_test(): gompertz_cdf() evaluates the Gompertz CDF gompertz_cdf_inv() inverts the Gompertz CDF. gompertz_pdf() evaluates the Gompertz PDF PDF parameter A = 2 PDF parameter B = 3 X PDF CDF CDf_inv 1.0281 0.0680781 0.988873 1.0281 0.383602 1.04728 0.732413 0.383602 0.313823 1.30268 0.650661 0.313823 0.367513 1.10259 0.715122 0.367513 0.213352 1.74489 0.498327 0.213352 0.0160798 2.88998 0.0473514 0.0160798 0.0566577 2.62355 0.159161 0.0566577 0.561895 0.563846 0.872681 0.561895 0.517777 0.6629 0.845667 0.517777 0.466119 0.795345 0.808081 0.466119 gompertz_sample_test(): gompertz_sample() samples the Gompertz distribution PDF parameter A = 2 PDF parameter B = 3 Sample size = 1000 Sample mean = 0.267969 Sample variance = 0.056476 Sample maximum = 1.3026 Sample minimum = 0.00023142 gumbel_cdf_test(): gumbel_cdf() evaluates the Gumbel CDF. gumbel_cdf_inv() inverts the Gumbel CDF. gumbel_pdf() evaluates the Gumbel PDF. X PDF CDF CDf_inv -0.615108 0.290908 0.15726 -0.615108 -1.00499 0.177837 0.0650971 -1.00499 -0.500945 0.316847 0.191996 -0.500945 1.2487 0.215331 0.750604 1.2487 1.01295 0.252563 0.695484 1.01295 -0.491884 0.318698 0.194876 -0.491884 0.69596 0.302839 0.607383 0.69596 0.728333 0.297886 0.617107 0.728333 1.68646 0.153873 0.830959 1.68646 0.400324 0.342862 0.511656 0.400324 gumbel_sample_test(): gumbel_mean() computes the Gumbel mean gumbel_sample() samples the Gumbel distribution gumbel_variance() computes the Gumbel variance. PDF mean = 0.577216 PDF variance = 1.64493 Sample size = 1000 Sample mean = 0.611767 Sample variance = 1.75204 Sample maximum = 7.12572 Sample minimum = -1.96731 half_normal_cdf_test(): half_normal_cdf() evaluates the Half Normal CDF. half_normal_cdf_inv() inverts the Half Normal CDF. half_normal_pdf() evaluates the Half Normal PDF. PDF parameter A = 0 PDF parameter B = 2 X PDF CDF CDf_inv 1.37636 0.314826 0.508662 1.37636 0.851428 0.364381 0.329684 0.851428 1.67402 0.281048 0.597412 1.67402 3.3056 0.101793 0.901629 3.3056 2.12985 0.226283 0.71309 2.12985 1.50606 0.300453 0.548568 1.50606 3.83562 0.0634226 0.944865 3.83562 0.566015 0.383282 0.222829 0.566015 2.28598 0.207598 0.746957 2.28598 0.52056 0.385655 0.205352 0.52056 half_normal_sample_test(): half_normal_mean() computes the Half Normal mean half_normal_sample() samples the Half Normal distribution half_normal_variance() computes the Half Normal variance. PDF parameter A = 0 PDF parameter B = 10 PDF mean = 7.97885 PDF variance = 36.338 Sample size = 1000 Sample mean = 7.83659 Sample variance = 37.4879 Sample maximum = 35.6162 Sample minimum = 0.0396396 hypergeometric_cdf_test(): hypergeometric_cdf() evaluates the Hypergeometric CDF. hypergeometric_pdf() evaluates the Hypergeometric PDF. PDF argument X = 7 Total number of balls = 100 Number of white balls = 7 Number of balls taken = 10 PDF value = = 7.49646e-09 CDF value = = 1 hypergeometric_sample_test(): hypergeometric_mean() computes the Hypergeometric mean hypergeometric_sample() samples the Hypergeometric distribution hypergeometric_variance() computes the Hypergeometric variance. PDF parameter N = 10 PDF parameter M = 7 PDF parameter L = 100 PDF mean = 0.7 PDF variance = 0.591818 Sample size = 1000 Sample mean = 0.692 Sample variance = 0.609136 Sample maximum = 4 Sample minimum = 0 i4_choose_test(): i4_choose() evaluates C(N,K). N K CNK 0 0 1 1 0 1 1 1 1 2 0 1 2 1 2 2 2 1 3 0 1 3 1 3 3 2 3 3 3 1 4 0 1 4 1 4 4 2 6 4 3 4 4 4 1 i4_factorial_log_test(): i4_factorial_log() evaluates log(N!). N lfact elfact fact 0 0 1 1 1 0 1 1 2 0.693147 2 2 3 1.79176 6 6 4 3.17805 24 24 5 4.78749 120 120 6 6.57925 720 720 7 8.52516 5040 5040 8 10.6046 40320 40320 9 12.8018 362880 362880 10 15.1044 3.6288e+06 3628800 11 17.5023 3.99168e+07 39916800 12 19.9872 4.79002e+08 479001600 i4_is_power_of_10_test(): i4_is_power_of_10() reports whether an I4 is a power of 10. I i4_is_power_of_10(I) 97 False 98 False 99 False 100 True 101 False 102 False 103 False i4mat_print_test(): i4mat_print() prints an I4MAT. A 5 x 6 integer matrix: Col: 0 1 2 3 4 5 Row 0: 11 12 13 14 15 16 1: 21 22 23 24 25 26 2: 31 32 33 34 35 36 3: 41 42 43 44 45 46 4: 51 52 53 54 55 56 i4mat_print_some_test(): i4mat_print_some() prints some of an I4MAT. Here is I4MAT, rows 0:2, cols 3:5: Col: 3 4 5 Row 0: 14 15 16 1: 24 25 26 2: 34 35 36 i4row_max_test(): i4row_max() computes maximums of an I4ROW. The matrix: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 Row maximums: 0 4 1 8 2 12 i4row_mean_test(): i4row_mean() computes row means of an I4ROW. The matrix: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 The row means: 0: 2.5 1: 6.5 2: 10.5 i4row_min_test(): i4row_min() computes minimums of an I4ROW. The matrix: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 Row minimums: 0 1 1 5 2 9 i4row_variance_test(): i4row_variance() computes variances of an I4ROW. The matrix: Col: 0 1 2 3 Row 0: 1 2 3 4 1: 5 6 7 8 2: 9 10 11 12 The row variances: 0: 1.66667 1: 1.66667 2: 1.66667 i4vec_print_test(): i4vec_print() prints an I4VEC. Here is an I4VEC: 0 91 1 92 2 93 3 94 i4vec_run_count_test(): i4vec_run_count() counts runs in an I4VEC Run Count Sequence 15 1 0 0 1 0 1 0 1 0 1 0 0 1 1 1 0 1 1 0 1 15 1 0 1 0 1 1 0 1 0 0 1 1 1 0 0 1 0 1 0 1 8 1 1 0 1 1 1 1 0 0 1 1 1 0 1 1 1 1 0 0 0 8 1 0 0 1 1 1 1 0 0 0 1 1 0 1 1 0 0 0 0 0 10 0 0 1 0 1 0 0 0 0 0 1 0 1 1 1 0 0 1 1 1 12 1 1 1 0 1 1 0 0 1 0 1 1 1 0 0 1 1 0 1 0 13 1 0 0 1 1 0 0 0 0 1 1 0 1 0 1 1 0 1 0 1 10 0 0 1 0 0 1 1 1 0 0 0 0 0 0 1 0 0 1 0 1 12 1 0 1 0 1 1 0 0 0 1 0 1 1 1 1 0 0 0 1 0 11 0 1 0 1 0 0 0 1 1 0 0 1 0 1 0 0 0 0 0 0 i4vec_unique_count_test(): i4vec_unique_count() counts unique entries in an I4VEC. Input vector: 0 20 1 5 2 12 3 5 4 17 5 10 6 13 7 3 8 6 9 3 10 8 11 15 12 5 13 5 14 16 15 12 16 19 17 8 18 5 19 7 Number of unique entries is 13 inverse_gaussian_cdf_test(): inverse_gaussian_cdf() evaluates the Inverse Gaussian CDF. inverse_gaussian_pdf() evaluates the Inverse Gaussian PDF. PDF parameter A = 5 PDF parameter B = 2 X PDF CDF 0.523386 0.322136 0.0743352 1.64009 0.203963 0.387606 0.808125 0.325444 0.168862 2.51733 0.12808 0.529568 3.13841 0.0970906 0.598828 4.89578 0.052078 0.72378 0.737125 0.332547 0.145486 46.768 0.000396725 0.993867 1.8967 0.176288 0.436276 0.396622 0.266505 0.0364423 inverse_gaussian_sample_test(): inverse_gaussian_mean() computes the Inverse Gaussian mean inverse_gaussian_sample() samples the Inverse Gaussian distribution inverse_gaussian_variance() computes the Inverse Gaussian variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 2 PDF variance = 2.66667 Sample size = 1000 Sample mean = 2.07162 Sample variance = 2.64417 Sample maximum = 11.5493 Sample minimum = 0.259133 laplace_cdf_test(): laplace_cdf() evaluates the Laplace CDF laplace_cdf_inv() inverts the Laplace CDF. laplace_pdf() evaluates the Laplace PDF PDF parameter A = 1 PDF parameter B = 2 X PDF CDF CDf_inv -2.576 0.0418235 0.083647 -2.576 -1.06651 0.0889619 0.177924 -1.06651 4.85179 0.0364363 0.927127 4.85179 0.393742 0.184626 0.369252 0.393742 1.18384 0.228045 0.54391 1.18384 5.39965 0.0277056 0.944589 5.39965 -0.892621 0.0970426 0.194085 -0.892621 0.621058 0.206849 0.413698 0.621058 -0.0815174 0.145577 0.291153 -0.0815174 1.77954 0.169303 0.661395 1.77954 laplace_sample_test(): laplace_mean() computes the Laplace mean laplace_sample() samples the Laplace distribution laplace_variance() computes the Laplace variance. PDF parameter A = 1 PDF parameter B = 2 PDF mean = 1 PDF variance = 8 Sample size = 1000 Sample mean = 1.08243 Sample variance = 7.66657 Sample maximum = 13.1229 Sample minimum = -11.2539 levy_cdf_test(): levy_cdf() evaluates the Levy CDF levy_cdf_inv() inverts the Levy CDF. levy_pdf() evaluates the Levy PDF PDF parameter A = 1 PDF parameter B = 2 X PDF CDF CDf_inv 2.02731 0.204703 0.162929 2.02731 17.0581 0.00823833 0.724153 17.0581 9.22355 0.0211848 0.621902 9.22355 40.4086 0.0022234 0.821763 40.4086 2.25117 0.181275 0.206116 2.25117 3.55839 0.0932657 0.376609 3.55839 2.86096 0.129852 0.299883 2.86096 10.5591 0.0171935 0.647377 10.5591 2.02375 0.205076 0.1622 2.02375 54.1307 0.00142966 0.846162 54.1307 log_normal_cdf_test(): log_normal_cdf() evaluates the Log Normal CDF log_normal_cdf_inv() inverts the Log Normal CDF. log_normal_pdf() evaluates the Log Normal PDF PDF parameter A = 10 PDF parameter B = 2.25 X PDF CDF CDf_inv 139314 9.09504e-07 0.793827 139314 2040.26 4.96875e-05 0.145163 2040.26 2995.94 3.99466e-05 0.187631 2995.94 65280.6 2.4172e-06 0.685405 65280.6 63278.7 2.51016e-06 0.680474 63278.7 9231.92 1.78238e-05 0.349571 9231.92 6266.85 2.42052e-05 0.288199 6266.85 16507 1.06535e-05 0.448993 16507 128446 1.01543e-06 0.783384 128446 34652.5 5.01403e-06 0.579803 34652.5 log_normal_sample_test(): log_normal_mean() computes the Log Normal mean log_normal_sample() samples the Log Normal distribution log_normal_variance() computes the Log Normal variance. PDF parameter A = 1 PDF parameter B = 2 PDF mean = 20.0855 PDF variance = 21623 Sample size = 1000 Sample mean = 18.4705 Sample variance = 5168.86 Sample maximum = 1180.72 Sample minimum = 0.00317592 log_series_cdf_test(): log_series_cdf() evaluates the Log Series CDF log_series_cdf_inv() inverts the Log Series CDF. log_series_pdf() evaluates the Log Series PDF PDF parameter A = 0.25 X PDF CDF CDf_inv 3 0.0181045 0.995746 4 2 0.108627 0.977642 3 1 0.869015 0.869015 2 1 0.869015 0.869015 2 1 0.869015 0.869015 2 1 0.869015 0.869015 2 2 0.108627 0.977642 3 1 0.869015 0.869015 2 1 0.869015 0.869015 2 1 0.869015 0.869015 2 log_series_sample_test(): log_series_mean() computes the Log Series mean log_series_variance() computes the Log Series variance log_series_sample() samples the Log Series distribution. PDF parameter A = 0.25 PDF mean = 1.15869 PDF variance = 0.202361 Sample size = 1000 Sample mean = 1.178 Sample variance = 0.262316 Sample maximum = 8 Sample minimum = 1 log_uniform_cdf_test(): log_uniform_cdf() evaluates the Log Uniform CDF log_uniform_cdf_inv() inverts the Log Uniform CDF. log_uniform_pdf() evaluates the Log Uniform PDF PDF parameter A = 2 PDF parameter B = 20 X PDF CDF CDf_inv 4.78621 0.0907388 0.378961 4.78621 3.36298 0.12914 0.225694 3.36298 2.05816 0.211011 0.0124493 2.05816 2.58935 0.167723 0.112162 2.58935 4.74037 0.0916162 0.374782 4.74037 19.1776 0.0226459 0.981765 19.1776 13.8051 0.031459 0.839009 13.8051 4.92686 0.0881483 0.391541 4.92686 3.20856 0.135355 0.20528 3.20856 2.30973 0.188028 0.0625307 2.30973 log_uniform_sample_test(): log_uniform_mean() computes the Log Uniform mean log_uniform_sample() samples the Log Uniform distribution log_uniform_variance() computes the Log Uniform variance PDF parameter A = 2 PDF parameter B = 20 PDF mean = 7.8173 PDF variance = 24.8801 Sample size = 1000 Sample mean = 7.8688 Sample variance = 25.4987 Sample maximum = 19.9435 Sample minimum = 2.00509 logistic_cdf_test(): logistic_cdf() evaluates the Logistic CDF logistic_cdf_inv() inverts the Logistic CDF. logistic_pdf() evaluates the Logistic PDF PDF parameter A = 1 PDF parameter B = 2 X PDF CDF CDf_inv 0.415404 0.122368 0.427441 0.415404 2.54455 0.10807 0.684012 2.54455 0.984684 0.124998 0.498086 0.984684 3.73468 0.0809103 0.79695 3.73468 6.06192 0.034141 0.926284 6.06192 5.42877 0.0443852 0.901534 5.42877 -5.4352 0.0185129 0.0385087 -5.4352 3.49982 0.0865568 0.777284 3.49982 -4.61022 0.026897 0.0570486 -4.61022 -5.82907 0.0154151 0.0318442 -5.82907 logistic_sample_test(): logistic_mean() computes the Logistic mean logistic_sample() samples the Logistic distribution logistic_variance() computes the Logistic variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 2 PDF variance = 29.6088 Sample size = 1000 Sample mean = 1.53439 Sample variance = 29.7635 Sample maximum = 22.1031 Sample minimum = -18.9141 lorentz_cdf_test(): lorentz_cdf() evaluates the Lorentz CDF lorentz_cdf_inv() inverts the Lorentz CDF. lorentz_pdf() evaluates the Lorentz PDF X PDF CDF CDf_inv 0.293155 0.293119 0.590771 0.293155 -4.22693 0.0168713 0.0739457 -4.22693 0.667898 0.220118 0.687438 0.667898 -2.7302 0.037652 0.111758 -2.7302 -1.71151 0.0810098 0.168316 -1.71151 -25.8743 0.000474748 0.012296 -25.8743 0.104757 0.314855 0.533224 0.104757 -2.23037 0.0532775 0.134163 -2.23037 -3.72957 0.0213492 0.0833862 -3.72957 -0.614922 0.230973 0.32451 -0.614922 lorentz_sample_test(): lorentz_mean() computes the Lorentz mean lorentz_variance() computes the Lorentz variance lorentz_sample() samples the Lorentz distribution. PDF mean = 0 PDF variance = 1.79769e+308 Sample size = 1000 Sample mean = -1.72588 Sample variance = 3711.49 Sample maximum = 190.989 Sample minimum = -1846.97 maxwell_cdf_test(): maxwell_cdf() evaluates the Maxwell CDF. maxwell_cdf_inv() inverts the Maxwell CDF. maxwell_pdf() evaluates the Maxwell PDF. PDF parameter A = 2 X PDF CDF CDf_inv 2.03418 0.246035 0.434658 2.03418 2.65605 0.291307 0.552252 2.65625 2.67743 0.291826 0.555923 2.67773 3.23424 0.282187 0.6431 3.23438 3.11023 0.287932 0.625059 3.11035 3.21849 0.283019 0.64085 3.21875 2.33992 0.275435 0.495083 2.33984 2.25137 0.268276 0.478098 2.25098 5.52162 0.0672774 0.862649 5.52148 3.94644 0.221706 0.732694 3.94727 maxwell_sample_test(): maxwell_mean() computes the Maxwell mean maxwell_variance() computes the Maxwell variance maxwell_sample() samples the Maxwell distribution. PDF parameter A = 2 PDF mean = 3.19154 PDF mean = 1.81408 Sample size = 1000 Sample mean = 3.14091 Sample variance = 1.84345 Sample maximum = 7.99228 Sample minimum = 0.178666 multinomial_coef_test(): multinomial_coef1 computes multinomial coefficients using the Gamma function multinomial_coef2 computes multinomial coefficients directly. Line 10 of the BINOMIAL table: 0 10 1 1 1 9 10 10 2 8 45 45 3 7 120 120 4 6 210 210 5 5 252 252 6 4 210 210 7 3 120 120 8 2 45 45 9 1 10 10 10 0 1 1 Level 5 of the TRINOMIAL coefficients: 0 0 5 1 1 0 1 4 5 5 0 2 3 10 10 0 3 2 10 10 0 4 1 5 5 0 5 0 1 1 1 0 4 5 5 1 1 3 20 20 1 2 2 30 30 1 3 1 20 20 1 4 0 5 5 2 0 3 10 10 2 1 2 30 30 2 2 1 30 30 2 3 0 10 10 3 0 2 10 10 3 1 1 20 20 3 2 0 10 10 4 0 1 5 5 4 1 0 5 5 5 0 0 1 1 multinomial_pdf_test(): multinomial_pdf() evaluates the Multinomial PDF. PDF argument X: 0 0 1 2 2 3 PDF parameter A = 5 PDF parameter B = 3 PDF parameter C: 0: 0.1 1: 0.5 2: 0.4 PDF value = 0.16 multinomial_sample_test(): multinomial_mean() computes the Multinomial mean multinomial_sample() samples the Multinomial distribution multinomial_variance() computes the Multinomial variance PDF parameter A = 5 PDF parameter B = 3 PDF parameter C: 0: 0.125 1: 0.5 2: 0.375 PDF means and variances: 0.625 0.546875 2.5 1.25 1.875 1.17188 Sample size = 1000 Component Min, Max, Mean, Variance: 0 0 4 0.681 0.597839 1 0 5 2.451 1.25886 2 0 5 1.868 1.15573 multinoulli_pdf_test(): multinoulli_pdf() evaluates the Multinoulli PDF. X pdf(X) -1 0 0 0.411237 1 0.0722834 2 0.0343151 3 0.12871 4 0.353455 5 0 nakagami_cdf_test(): nakagami_cdf() evaluates the Nakagami CDF nakagami_pdf() evaluates the Nakagami PDF X PDF CDF CDf_inv 3.18257 0.586699 0.692394 3.18258 3.2582 0.540746 0.735053 3.2582 3.31623 0.503068 0.765346 3.31623 3.36515 0.470346 0.789159 3.36515 3.40825 0.441184 0.808803 3.40825 3.44721 0.414809 0.82548 3.44721 3.48305 0.390724 0.839911 3.48305 3.5164 0.36858 0.852572 3.5164 3.54772 0.348118 0.863796 3.54772 3.57735 0.329135 0.873828 3.57735 nakagami_sample_test(): nakagami_mean() computes the Nakagami mean nakagami_variance() computes the Nakagami variance. PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF mean = 2.91874 PDF variance = 0.318446 negative_binomial_cdf_test(): negative_binomial_cdf() evaluates the Negative Binomial CDF. negative_binomial_cdf_inv() inverts the Negative Binomial CDF. negative_binomial_pdf() evaluates the Negative Binomial PDF. PDF parameter A = 2 PDF parameter B = 0.25 X PDF CDF CDf_inv 9 0.0667419 0.699661 9 11 0.0469279 0.802903 11 18 0.010649 0.960536 18 7 0.0889893 0.555054 7 3 0.09375 0.15625 3 9 0.0667419 0.699661 9 9 0.0667419 0.699661 9 7 0.0889893 0.555054 7 9 0.0667419 0.699661 9 5 0.105469 0.367188 5 negative_binomial_sample_test(): negative_binomial_mean() computes the Negative Binomial mean negative_binomial_sample() samples the Negative Binomial distribution negative_binomial_variance() computes the Negative Binomial variance. PDF parameter A = 2 PDF parameter B = 0.75 PDF mean = 2.66667 PDF variance = 0.888889 Sample size = 1000 Sample mean = 2.606 Sample variance = 0.802764 Sample maximum = 9 Sample minimum = 2 normal_01_cdf_test(): normal_01_cdf() evaluates the Normal 01 CDF normal_01_cdf_inv() inverts the Normal 01 CDF. normal_01_pdf() evaluates the Normal 01 PDF X PDF CDF CDf_inv 0.731905 0.305202 0.767887 0.731905 0.927627 0.259452 0.823199 0.927627 -0.66087 0.32068 0.254348 -0.66087 -2.05777 0.0480191 0.0198059 -2.05777 -1.30137 0.171063 0.0965655 -1.30137 -2.66742 0.0113731 0.00382178 -2.66742 -0.435147 0.362905 0.331728 -0.435147 -0.306854 0.380595 0.379477 -0.306854 0.338626 0.376713 0.632554 0.338626 0.111619 0.396465 0.544437 0.111619 normal_01_cdf_values_test(): normal_01_cdf_values() stores values of the unit normal CDF. X normal_01_cdf(X) 0.000000 0.5000000000000000 0.100000 0.5398278372770290 0.200000 0.5792597094391030 0.300000 0.6179114221889526 0.400000 0.6554217416103242 0.500000 0.6914624612740131 0.600000 0.7257468822499270 0.700000 0.7580363477769270 0.800000 0.7881446014166033 0.900000 0.8159398746532405 1.000000 0.8413447460685429 1.500000 0.9331927987311419 2.000000 0.9772498680518208 2.500000 0.9937903346742240 3.000000 0.9986501019683699 3.500000 0.9997673709209645 4.000000 0.9999683287581669 normal_01_sample_test(): normal_01_mean() computes the Normal 01 mean normal_01_sample() samples the Normal 01 distribution normal_01_variance() returns the Normal 01 variance. PDF mean = 0 PDF variance = 1 Sample size = 1000 Sample mean = -0.00479616 Sample variance = 1.01771 Sample maximum = 3.5583 Sample minimum = -3.07723 normal_cdf_test(): normal_cdf() evaluates the Normal CDF normal_cdf_inv() inverts the Normal CDF. normal_pdf() evaluates the Normal PDF PDF parameter MU = 100 PDF parameter SIGMA = 15 X PDF CDF CDf_inv 103.96 0.0256853 0.604111 103.96 101.664 0.0264331 0.544153 101.664 92.411 0.0234011 0.306451 92.411 75.9775 0.00737695 0.0546328 75.9775 132.432 0.00256863 0.984696 132.432 98.5354 0.0264697 0.46111 98.5354 82.5667 0.0135365 0.122573 82.5667 101.732 0.0264195 0.545952 101.732 91.0828 0.0222883 0.276096 91.0828 74.8551 0.00652562 0.0468374 74.8551 normal_sample_test(): normal_mean() computes the Normal mean normal_sample samples() the Normal distribution normal_variance() returns the Normal variance. PDF parameter MU = 100 PDF parameter SIGMA = 15 PDF mean = 100 PDF variance = 225 Sample size = 1000 Sample mean = 100.28 Sample variance = 220.49 Sample maximum = 152.883 Sample minimum = 56.0497 normal_truncated_ab_cdf_test(): normal_truncated_ab_cdf() evaluates the Normal Truncated AB CDF. normal_truncated_ab_cdf_inv() inverts the Normal Truncated AB CDF. normal_truncated_ab_pdf() evaluates the Normal Truncated AB PDF. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [50,150] X PDF CDF CDf_inv 120.264 0.0120372 0.805073 120.264 111.193 0.015124 0.681059 111.193 100.303 0.0167172 0.505074 100.303 122.218 0.0112637 0.827843 122.218 95.4756 0.0164468 0.42477 95.4756 98.0584 0.016668 0.467572 98.0584 90.0357 0.0154418 0.33772 90.0357 95.2406 0.0164181 0.420908 95.2406 77.245 0.0110483 0.16617 77.245 71.9604 0.00891321 0.113431 71.9604 normal_truncated_ab_sample_test(): normal_truncated_ab_mean() computes the Normal Truncated AB mean normal_truncated_ab_sample() samples the Normal Truncated AB distribution normal_truncated_ab_variance() computes the Normal Truncated AB variance. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [50,150] PDF mean = 100 PDF variance = 483.588 Sample size = 1000 Sample mean = 99.7945 Sample variance = 484.216 Sample maximum = 149.986 Sample minimum = 50.2396 normal_truncated_a_cdf_test(): normal_truncated_a_cdf() evaluates the Normal Truncated A CDF. normal_truncated_a_cdf_inv() inverts the Normal Truncated A CDF. normal_truncated_a_pdf() evaluates the Normal Truncated A PDF. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [50,+oo) X PDF CDF CDf_inv 99.8934 0.016329 0.486619 99.8934 116.591 0.0131017 0.740641 116.591 86.4801 0.0141077 0.277896 86.4801 105.24 0.0159744 0.573301 105.24 121.438 0.0113053 0.799871 121.438 86.1441 0.0140043 0.273173 86.1441 70.9834 0.00832602 0.102471 70.9834 122.095 0.0110496 0.807216 122.095 112.367 0.0144486 0.68236 112.367 94.6398 0.0159581 0.401498 94.6398 normal_truncated_a_sample_test(): normal_truncated_a_mean() computes the Normal Truncated A mean normal_truncated_a_sample() samples the Normal Truncated A distribution normal_truncated_a_variance() computes the Normal Truncated A variance. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [50,+oo] PDF mean = 101.381 PDF variance = 554.032 Sample size = 1000 Sample mean = 102.039 Sample variance = 571.869 Sample maximum = 186.228 Sample minimum = 50.0427 normal_truncated_b_cdf_test(): normal_truncated_b_cdf() evaluates the Normal Truncated B CDF. normal_truncated_b_cdf_inv() inverts the Normal Truncated B CDF. normal_truncated_b_pdf() evaluates the Normal Truncated B PDF. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [-oo,150] X PDF CDF CDf_inv 108.809 0.0153464 0.652556 108.809 89.5094 0.014953 0.345234 89.5094 94.7744 0.0159763 0.426927 94.7744 93.6908 0.0158174 0.409699 93.6908 70.4872 0.00813479 0.121666 70.4872 89.919 0.0150541 0.35138 89.919 105.38 0.0159554 0.598814 105.38 112.096 0.0145254 0.70172 112.096 64.0947 0.00582177 0.0772282 64.0947 78.9037 0.0114376 0.204017 78.9037 normal_truncated_b_sample_test(): normal_truncated_b_mean() computes the Normal Truncated B mean normal_truncated_b_sample() samples the Normal Truncated B distribution normal_truncated_b_variance() computes the Normal Truncated B variance. The "parent" normal distribution has mean = 100 standard deviation = 25 The parent distribution is truncated to the interval [-oo,150] PDF mean = 98.6188 PDF variance = 554.032 Sample size = 1000 Sample mean = 97.2759 Sample variance = 546.985 Sample maximum = 149.781 Sample minimum = 21.4415 owen_values_test(): owen_values() stores values of the OWEN function. H A T 0.062500 0.250000 0.038912 6.500000 0.437500 0.000000 7.000000 0.968750 0.000000 4.781250 0.062500 0.000000 2.000000 0.500000 0.008625 1.000000 0.999997 0.066742 1.000000 0.500000 0.043065 1.000000 1.000000 0.066742 1.000000 2.000000 0.078468 1.000000 3.000000 0.079300 0.500000 0.500000 0.064489 0.500000 1.000000 0.106671 0.500000 2.000000 0.141581 0.500000 3.000000 0.151084 0.250000 0.500000 0.071347 0.250000 1.000000 0.120129 0.250000 2.000000 0.166613 0.250000 3.000000 0.184750 0.125000 0.500000 0.073173 0.125000 1.000000 0.123763 0.125000 2.000000 0.173744 0.125000 3.000000 0.195119 0.007812 0.500000 0.073789 0.007812 1.000000 0.124995 0.007812 2.000000 0.176198 0.007812 3.000000 0.198777 0.007812 10.000000 0.234089 0.007812 100.000000 0.247946 pareto_cdf_test(): pareto_cdf() evaluates the Pareto CDF pareto_cdf_inv() inverts the Pareto CDF. pareto_pdf() evaluates the Pareto PDF PDF parameter A = 0.5 PDF parameter B = 5 X PDF CDF CDf_inv 0.517829 8.10405 0.160697 0.517829 0.874195 0.350083 0.938792 0.874195 0.659742 1.89485 0.749978 0.659742 0.780058 0.693517 0.891803 0.780058 0.550694 5.60218 0.382982 0.550694 0.737214 0.973329 0.85649 0.737214 0.500328 9.96069 0.00327658 0.500328 0.578963 4.14874 0.519607 0.578963 0.526532 7.33285 0.227805 0.526532 0.580768 4.07197 0.527026 0.580768 pareto_sample_test(): pareto_mean() computes the Pareto mean pareto_sample() samples the Pareto distribution pareto_variance() computes the Pareto variance. PDF parameter A = 0.5 PDF parameter B = 5 PDF mean = 0.625 PDF variance = 0.0260417 Sample size = 1000 Sample mean = 0.629221 Sample variance = 0.0242369 Sample maximum = 2.1028 Sample minimum = 0.500162 pearson_05_pdf_test(): pearson_05_pdf() evaluates the Pearson 05 PDF. PDF argument X = 5 PDF parameter A = 1 PDF parameter B = 2 PDF parameter C = 3 PDF value = 0.0758163 planck_pdf_test(): planck_pdf() evaluates the Planck PDF. PDF parameter A = 2 PDF parameter B = 3 X PDF 2.15477 0.3358 3.14647 0.142199 1.84301 0.396647 1.01347 0.389138 2.06509 0.354613 0.9367 0.367477 1.45399 0.437275 0.987541 0.382279 0.299396 0.0806452 1.60952 0.427972 planck_sample_test(): planck_mean() returns the mean of the Planck distribution. planck_sample() samples the Planck distribution. planck_variance() returns the variance of the Planck distribution. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 3.83223 PDF variance = 4.11319 Sample size = 1000 Sample mean = 1.90693 Sample variance = 1.06136 Sample maximum = 8.31417 Sample minimum = 0.147987 poisson_cdf_test(): poisson_cdf() evaluates the Poisson CDF, poisson_cdf_inv() inverts the Poisson CDF. poisson_pdf() evaluates the Poisson PDF. PDF parameter A = 10 X PDF CDF CDf_inv 12 0.0947803 0.791556 12 11 0.113736 0.696776 11 14 0.0520771 0.916542 14 6 0.0630555 0.130141 6 11 0.113736 0.696776 11 14 0.0520771 0.916542 14 9 0.12511 0.45793 9 10 0.12511 0.58304 10 12 0.0947803 0.791556 12 7 0.0900792 0.220221 7 poisson_sample_test(): poisson_mean() computes the Poisson mean poisson_sample() samples the Poisson distribution poisson_variance() computes the Poisson variance. PDF parameter A = 10 PDF mean = 10 PDF variance = 10 Sample size = 1000 Sample mean = 10.019 Sample variance = 10.4066 Sample maximum = 21 Sample minimum = 2 power_cdf_test(): power_cdf() evaluates the Power CDF power_cdf_inv() inverts the Power CDF. power_pdf() evaluates the Power PDF PDF parameter A = 2 PDF parameter B = 3 X PDF CDF CDf_inv 1.67021 0.371158 0.309956 1.67021 1.96945 0.437656 0.430971 1.96945 1.6985 0.377445 0.320546 1.6985 2.51695 0.559322 0.703893 2.51695 2.8617 0.635933 0.909925 2.8617 1.81592 0.403537 0.366396 1.81592 2.53173 0.562606 0.712182 2.53173 2.97125 0.660277 0.980922 2.97125 2.86342 0.636315 0.911019 2.86342 0.953389 0.211864 0.100994 0.953389 power_sample_test(): power_mean() computes the Power mean power_sample() samples the Power distribution power_variance() computes the Power variance. PDF parameter A = 2 PDF parameter B = 3 PDF mean = 2 PDF variance = 0.5 Sample size = 1000 Sample mean = 2.03268 Sample variance = 0.476582 Sample maximum = 2.99949 Sample minimum = 0.121254 psi_values_test(): psi_values() stores values of the PSI function. X PSI(X) 0.100000 -10.4237549404110794 0.200000 -5.2890398965921879 0.300000 -3.5025242222001332 0.400000 -2.5613845445851160 0.500000 -1.9635100260214231 0.600000 -1.5406192138931900 0.700000 -1.2200235536979349 0.800000 -0.9650085667061385 0.900000 -0.7549269499470515 1.000000 -0.5772156649015329 1.100000 -0.4237549404110768 1.200000 -0.2890398965921883 1.300000 -0.1691908888667997 1.400000 -0.0613845445851161 1.500000 0.0364899739785765 1.600000 0.1260474527734763 1.700000 0.2085478748734940 1.800000 0.2849914332938615 1.900000 0.3561841611640597 2.000000 0.4227843350984671 quasigeometric_cdf_test(): quasigeometric_cdf() evaluates the Quasigeometric CDF quasigeometric_cdf_inv() inverts the Quasigeometric CDF. quasigeometric_pdf() evaluates the Quasigeometric PDF PDF parameter A = 0.4825 PDF parameter B = 0.5893 X PDF CDF CDf_inv 1 0.212537 0.695037 2 0 0.4825 0.4825 1 6 0.0151049 0.978327 7 1 0.212537 0.695037 2 0 0.4825 0.4825 1 0 0.4825 0.4825 1 0 0.4825 0.4825 1 0 0.4825 0.4825 1 0 0.4825 0.4825 1 1 0.212537 0.695037 2 quasigeometric_sample_test(): quasigeometric_mean() computes the Quasigeometric mean quasigeometric_sample() samples the Quasigeometric distribution quasigeometric_variance() computes the Quasigeometric variance. PDF parameter A = 0.4825 PDF parameter B = 0.5893 PDF parameter A = 0.4825 PDF mean = 1.26004 PDF variance = 3.28832 Sample size = 1000 Sample mean = 1.241 Sample variance = 3.37492 Sample maximum = 19 Sample minimum = 0 r8_beta_test(): r8_beta() evaluates the Beta function. X Y BETA(X,Y) r8_beta(X,Y) tabulated computed. 0.2 1 5 5 0.4 1 2.5 2.5 0.6 1 1.666666666666667 1.666666666666667 0.8 1 1.25 1.25 1 0.2 5 5 1 0.4 2.5 2.5 1 1 1 1 2 2 0.1666666666666667 0.1666666666666667 3 3 0.03333333333333333 0.03333333333333333 4 4 0.007142857142857143 0.007142857142857143 5 5 0.001587301587301587 0.001587301587301587 6 2 0.02380952380952381 0.02380952380952381 6 3 0.005952380952380952 0.005952380952380952 6 4 0.001984126984126984 0.001984126984126984 6 5 0.0007936507936507937 0.0007936507936507937 6 6 0.0003607503607503608 0.0003607503607503608 7 7 8.325008325008325e-05 8.325008325008325e-05 r8_csc_test(): r8_csc() computes the cosecant of an angle. ANGLE r8_csc(ANGLE) 0.00 Undefined 15.00 1.46037 30.00 -1.00194 45.00 1.29876 60.00 -8.05172 75.00 -1.69796 90.00 1.0177 105.00 -1.18588 120.00 4.05727 135.00 2.06607 150.00 -1.05049 165.00 1.10683 180.00 Undefined 195.00 -2.69119 210.00 1.10311 225.00 -1.05294 240.00 2.09321 255.00 3.94632 270.00 -1.18053 285.00 1.0191 300.00 -1.71484 315.00 -7.61642 330.00 1.29113 345.00 -1.00241 360.00 Undefined r8_erf_test(): r8_erf() evaluates the error function. X ERF(X) r8_erf(X) 0 0 0 0.1 0.1124629160182849 0.1124629160182849 0.2 0.2227025892104785 0.2227025892104785 0.3 0.3286267594591274 0.3286267594591273 0.4 0.4283923550466685 0.4283923550466684 0.5 0.5204998778130465 0.5204998778130465 0.6 0.6038560908479259 0.6038560908479259 0.7 0.6778011938374185 0.6778011938374184 0.8 0.7421009647076605 0.7421009647076605 0.9 0.7969082124228321 0.7969082124228322 1 0.8427007929497149 0.8427007929497148 1.1 0.8802050695740817 0.8802050695740817 1.2 0.9103139782296354 0.9103139782296354 1.3 0.9340079449406524 0.9340079449406524 1.4 0.9522851197626488 0.9522851197626487 1.5 0.9661051464753106 0.9661051464753108 1.6 0.976348383344644 0.976348383344644 1.7 0.9837904585907746 0.9837904585907746 1.8 0.9890905016357306 0.9890905016357308 1.9 0.9927904292352575 0.9927904292352574 2 0.9953222650189527 0.9953222650189527 r8_gamma_inc_test(): r8_gamma_inc() evaluates the normalized incomplete Gamma function P(A,X). A X Exact F r8_gamma_inc(A,X) 0.1 0.03 2.4903 0.738235 0.1 0.3 0.871837 0.908358 0.1 1.5 0.107921 0.988656 0.5 0.075 1.23812 0.301465 0.5 0.75 0.39113 0.779329 0.5 3.5 0.0144472 0.991849 1 0.1 0.904837 0.0951626 1 1 0.367879 0.632121 1 5 0.00673795 0.993262 1.1 0.1 0.882797 0.0720597 1.1 1 0.390833 0.589181 1.1 5 0.00805146 0.991537 2 0.15 0.989814 0.0101858 2 1.5 0.557825 0.442175 2 7 0.00729506 0.992705 6 2.5 114.957 0.042021 6 12 2.44092 0.979659 11 16 280855 0.922604 26 25 8.57648e+24 0.447079 41 45 2.08503e+47 0.744455 r8_zeta_test(): r8_zeta() estimates the Zeta function. P r8_zeta(P) 1 1.79769e+308 2 1.64493 3 1.20206 4 1.08232 5 1.03693 6 1.01734 7 1.00835 8 1.00408 9 1.00201 10 1.00099 11 1.00049 12 1.00025 13 1.00012 14 1.00006 15 1.00003 16 1.00002 17 1.00001 18 1 19 1 20 1 21 1 22 1 23 1 24 1 25 1 3 1.20206 3.125 1.17905 3.25 1.15915 3.375 1.14185 3.5 1.12673 3.625 1.11347 3.75 1.10179 3.875 1.09147 4 1.08232 r8poly_print_test(): r8poly_print() prints an R8POLY. The R8POLY: p(x) = 9 * x^5 + 0.78 * x^4 + 56 * x^2 - 3.4 * x + 12 r8poly_value_horner_test(): r8poly_value_horner() evaluates a polynomial at a point using Horner's method. The polynomial coefficients: p(x) = 1 * x^4 - 10 * x^3 + 35 * x^2 - 50 * x + 24 I X P(X) 0 0.0000 24 1 0.3333 10.8642 2 0.6667 3.45679 3 1.0000 0 4 1.3333 -0.987654 5 1.6667 -0.691358 6 2.0000 0 7 2.3333 0.493827 8 2.6667 0.493827 9 3.0000 0 10 3.3333 -0.691358 11 3.6667 -0.987654 12 4.0000 0 13 4.3333 3.45679 14 4.6667 10.8642 15 5.0000 24 r8row_max_test(): r8row_max() computes maximums of an R8ROW. The matrix: Col: 0 1 2 3 Row 0 : 1 2 3 4 1 : 5 6 7 8 2 : 9 10 11 12 Row maximums: 0: 4 1: 8 2: 12 r8row_mean_test(): r8row_mean() computes row means of an R8ROW. The matrix: Col: 0 1 2 3 Row 0 : 1 2 3 4 1 : 5 6 7 8 2 : 9 10 11 12 The row means: 0: 2.5 1: 6.5 2: 10.5 r8row_min_test(): r8row_min() computes minimums of an R8ROW. The matrix: Col: 0 1 2 3 Row 0 : 1 2 3 4 1 : 5 6 7 8 2 : 9 10 11 12 Row minimums: 0: 1 1: 5 2: 9 r8row_variance_test(): r8row_variance() computes variances of an R8ROW. The matrix: Col: 0 1 2 3 Row 0 : 1 2 3 4 1 : 5 6 7 8 2 : 9 10 11 12 The row variances: 0: 1.66667 1: 1.66667 2: 1.66667 r8vec_dot_product_test(): r8vec_dot_product() computes the dot product of two R8VEC's. V1 and V2: 0: 0.981491 0.457154 1: 0.817016 0.514551 2: 0.467281 0.216354 3: 0.320782 0.164262 4: 0.137316 0.230044 5: 0.412487 0.304326 6: 0.61192 0.129119 7: 0.154362 0.825986 8: 0.880713 0.623751 9: 0.218678 0.970163 V1 dot V2 = 2.14801 r8vec_transpose_print_test(): r8vec_transpose_print() prints an R8VEC "tranposed", that is, placing multiple entries on a line. The vector X: 1 2.1 3.2 4.3 5.4 6.5 7.6 8.7 9.8 10.9 11 r8vec2_print_test(): r8vec2_print() prints a pair of R8VEC's. Print a pair of R8VEC's: 0: 0 0 1: 0.2 0.04 2: 0.4 0.16 3: 0.6 0.36 4: 0.8 0.64 5: 1 1 rayleigh_cdf_test(): rayleigh_cdf() evaluates the Rayleigh CDF rayleigh_cdf_inv() inverts the Rayleigh CDF. rayleigh_pdf() evaluates the Rayleigh PDF PDF parameter A = 2 X PDF CDF CDf_inv 1.7593 0.298714 0.320836 1.7593 2.40392 0.291841 0.514392 2.40392 4.26433 0.109802 0.897004 4.26433 1.02703 0.225042 0.123528 1.02703 2.13505 0.301915 0.434365 2.13505 2.8584 0.257345 0.639876 2.8584 0.670533 0.158472 0.0546517 0.670533 2.40435 0.291818 0.514517 2.40435 4.57177 0.0838249 0.926659 4.57177 0.960424 0.213958 0.108903 0.960424 rayleigh_sample_test(): rayleigh_mean() computes the Rayleigh mean rayleigh_sample() samples the Rayleigh distribution rayleigh_variance() computes the Rayleigh variance. PDF parameter A = 2 PDF mean = 2.50663 PDF variance = 1.71681 Sample size = 1000 Sample mean = 2.55451 Sample variance = 1.79929 Sample maximum = 7.21056 Sample minimum = 0.18761 reciprocal_cdf_test(): reciprocal_cdf() evaluates the Reciprocal CDF. reciprocal_cdf_inv() inverts the Reciprocal CDF. reciprocal_pdf() evaluates the Reciprocal PDF. PDF parameter A = 1 PDF parameter B = 3 X PDF CDF CDf_inv 2.01429 0.451891 0.63741 2.01429 1.44494 0.629948 0.335031 1.44494 2.02845 0.448737 0.643785 2.02845 2.44622 0.3721 0.81425 2.44622 1.15282 0.789574 0.129449 1.15282 2.89565 0.314347 0.967775 2.89565 1.70661 0.53336 0.486533 1.70661 2.6004 0.350039 0.869883 2.6004 1.76898 0.514557 0.519202 1.76898 1.59175 0.571847 0.423112 1.59175 reciprocal_sample_test(): reciprocal_mean() computes the Reciprocal mean reciprocal_sample() samples the Reciprocal distribution reciprocal_variance() computes the Reciprocal variance. PDF parameter A = 1 PDF parameter B = 3 PDF mean = 1.82048 PDF variance = 0.326815 Sample size = 1000 Sample mean = 1.81954 Sample variance = 0.348448 Sample maximum = 2.99276 Sample minimum = 1.00061 sech_cdf_test(): sech_cdf() evaluates the Sech CDF. sech_cdf_inv() inverts the Sech CDF. sech_pdf() evaluates the Sech PDF. PDF parameter A = 3 PDF parameter B = 2 X PDF CDF CDf_inv 4.90852 0.106749 0.765983 4.90852 2.91907 0.159025 0.487123 2.91907 -5.48413 0.00457568 0.00915262 -5.48413 2.86971 0.158818 0.479278 2.86971 1.09425 0.106859 0.234313 1.09425 2.27907 0.149347 0.387668 2.27907 2.93454 0.15907 0.489584 2.93454 9.07687 0.0152152 0.969523 9.07687 2.46029 0.153531 0.415126 2.46029 4.70575 0.114807 0.74352 4.70575 sech_sample_test(): sech_mean() computes the Sech mean sech_sample() samples the Sech distribution sech_variance() computes the Sech variance. PDF parameter A = 3 PDF parameter B = 2 PDF mean = 3 PDF variance = 9.8696 Sample size = 1000 Sample mean = 3.06052 Sample variance = 9.57945 Sample maximum = 18.1325 Sample minimum = -10.4983 semicircular_cdf_test(): semicircular_cdf() evaluates the Semicircular CDF. semicircular_cdf_inv() inverts the Semicircular CDF. semicircular_pdf() evaluates the Semicircular PDF. PDF parameter A = 3 PDF parameter B = 2 X PDF CDF CDf_inv 3.90775 0.283635 0.778693 3.90771 2.29275 0.297743 0.279658 2.29297 4.90892 0.0949666 0.994207 4.9082 2.80739 0.31683 0.438786 2.80762 2.58979 0.311542 0.370346 2.58984 3.52745 0.307041 0.665926 3.52734 1.74414 0.247731 0.128345 1.74414 4.29198 0.242981 0.880518 4.29199 3.14053 0.317523 0.544695 3.14062 3.59507 0.303894 0.686582 3.59521 semicircular_sample_test(): semicircular_mean() computes the Semicircular mean semicircular_sample() samples the Semicircular distribution semicircular_variance() computes the Semicircular variance. PDF parameter A = 3 PDF parameter B = 2 PDF mean = 3 PDF variance = 1 Sample size = 1000 Sample mean = 3.01667 Sample variance = 0.980586 Sample maximum = 4.87918 Sample minimum = 1.0144 sin_power_int_test(): sin_power_int() returns values of the integral of SIN(X)^N from A to B. A B N Exact Computed 10.000000 20.000000 0 1.000000e+01 1.000000e+01 0.000000 1.000000 1 4.596977e-01 4.596977e-01 0.000000 1.000000 2 2.726756e-01 2.726756e-01 0.000000 1.000000 3 1.789406e-01 1.789406e-01 0.000000 1.000000 4 1.240256e-01 1.240256e-01 0.000000 1.000000 5 8.897440e-02 8.897440e-02 0.000000 2.000000 5 9.039312e-01 9.039312e-01 1.000000 2.000000 5 8.149568e-01 8.149568e-01 0.000000 1.000000 10 2.188752e-02 2.188752e-02 0.000000 1.000000 11 1.702344e-02 1.702344e-02 sin_power_int_values_test(): sin_power_int_values() stores values of the sine power integral. A B N F 10.000000 20.000000 0 10 0.000000 1.000000 1 0.4596976941318603 0.000000 1.000000 2 0.2726756432935796 0.000000 1.000000 3 0.1789405625488581 0.000000 1.000000 4 0.1240255653152068 0.000000 1.000000 5 0.08897439645157594 0.000000 2.000000 5 0.9039312384814995 1.000000 2.000000 5 0.8149568420299235 0.000000 1.000000 10 0.02188752242172985 0.000000 1.000000 11 0.01702343937406933 stirling2_number_test(): stirling2_number() calculates a Stirling number S2(n,k) of the second kind. 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 3 1 0 0 0 0 0 0 1 7 6 1 0 0 0 0 0 1 15 25 10 1 0 0 0 0 1 31 90 65 15 1 0 0 0 1 63 301 350 140 21 1 0 0 1 127 966 1701 1050 266 28 1 student_cdf_test(): student_cdf() evaluates the Student CDF. student_pdf() evaluates the Student PDF. PDF argument X = 2.447 PDF parameter A = 0.5 PDF parameter B = 2 PDF parameter C = 6 PDF value = 0.14754 CDF value = 0.816049 student_sample_test(): student_mean() computes the Student mean student_sample() samples the Student distribution student_variance() computes the Student variance. PDF parameter A = 0.5 PDF parameter B = 2 PDF parameter C = 6 PDF mean = 0.5 PDF variance = 6 Sample size = 1000 Sample mean = 0.485189 Sample variance = 2.11725 Sample maximum = 14.9243 Sample minimum = -9.02625 student_noncentral_cdf_test(): student_noncentral_cdf() evaluates the Student Noncentral CDF PDF argument X = 0.5 PDF parameter IDF = 10 PDF parameter B = 1 CDF value = 0.30528 tfn_test(): tfn() evaluates Owen's T function. H A T(H,A) Exact 0.0625 0.25 0.0389119 0.0389119 6.5 0.4375 2.00058e-11 2.00058e-11 7 0.96875 6.39906e-13 6.39906e-13 4.78125 0.0625 1.0633e-07 1.0633e-07 2 0.5 0.00862508 0.00862508 1 0.999997 0.0667418 0.0667418 1 0.5 0.0430647 0.0430647 1 1 0.0667419 0.0667419 1 2 0.0784682 0.0784682 1 3 0.0792995 0.0792995 0.5 0.5 0.0644886 0.0644886 0.5 1 0.106671 0.106671 0.5 2 0.141581 0.141581 0.5 3 0.151084 0.151084 0.25 0.5 0.0713466 0.0713466 0.25 1 0.120129 0.120129 0.25 2 0.166613 0.166613 0.25 3 0.18475 0.18475 0.125 0.5 0.0731727 0.0731727 0.125 1 0.123763 0.123763 0.125 2 0.173744 0.173744 0.125 3 0.195119 0.195119 0.0078125 0.5 0.0737894 0.0737894 0.0078125 1 0.124995 0.124995 0.0078125 2 0.176198 0.176198 0.0078125 3 0.198777 0.198777 0.0078125 10 0.234074 0.234089 0.0078125 100 0.233737 0.247946 triangle_cdf_test(): triangle_cdf() evaluates the Triangle CDF triangle_cdf_inv() inverts the Triangle CDF. triangle_pdf() evaluates the Triangle PDF PDF parameter A = 1 PDF parameter B = 3 PDF parameter C = 10 X PDF CDF CDf_inv 4.41727 0.177229 0.505288 4.41727 4.2626 0.18214 0.477496 4.2626 3.59511 0.20333 0.348847 3.59511 4.64378 0.170039 0.544617 4.64378 2.55713 0.173015 0.134704 2.55713 7.29103 0.085999 0.883516 7.29103 2.36015 0.151128 0.102779 2.36015 3.04705 0.220729 0.232642 3.04705 5.04919 0.157169 0.610944 5.04919 8.10379 0.0601971 0.942927 8.10379 triangle_sample_test(): triangle_mean() returns the Triangle mean triangle_sample samples the Triangle distribution triangle_variance returns the Triangle variance PDF parameter A = 1 PDF parameter B = 3 PDF parameter C = 10 PDF parameter MEAN = 4.66667 PDF parameter VARIANCE = 3.72222 Sample size = 1000 Sample mean = 4.70524 Sample variance = 3.87273 Sample maximum = 9.83817 Sample minimum = 1.05426 triangular_cdf_test(): triangular_cdf() evaluates the Triangular CDF triangular_cdf_inv() inverts the Triangular CDF. triangular_pdf() evaluates the Triangular PDF PDF parameter A = 1 PDF parameter B = 10 X PDF CDF CDf_inv 2.1247 0.0555407 0.0312332 2.1247 3.55199 0.126024 0.160807 3.55199 7.13908 0.14128 0.797905 7.13908 7.80497 0.108396 0.881034 7.80497 6.59894 0.167954 0.714389 6.59894 8.29929 0.0839859 0.928582 8.29929 2.99857 0.0986946 0.0986239 2.99857 5.0976 0.202351 0.414576 5.0976 6.50571 0.172558 0.698517 6.50571 7.10991 0.14272 0.793763 7.10991 triangular_sample_test(): triangular_mean() computes the Triangular mean triangular_sample() samples the Triangular distribution triangular_variance() computes the Triangular variance. PDF parameter A = 1 PDF parameter B = 10 PDF mean = 5.5 PDF variance = 3.375 Sample size = 1000 Sample mean = 5.40015 Sample variance = 3.55048 Sample maximum = 9.67196 Sample minimum = 1.0477 trigamma_test(): trigamma() evaluates the trigamma function. X FX FX Tabulated Computed 1 1.644934066848226 1.644934065473016 1.1 1.433299150792759 1.43329914968199 1.2 1.267377205423779 1.267377204522996 1.3 1.134253434996619 1.134253434263296 1.4 1.025356590529597 1.025356589930374 1.5 0.9348022005446793 0.9348022000532704 1.6 0.8584318931245799 0.8584318927201864 1.7 0.7932328301639984 0.793232829830095 1.8 0.7369741375017002 0.7369741372251055 1.9 0.6879720582426356 0.6879720580127948 2 0.6449340668482264 0.6449340654730159 trigamma_values_test(): trigamma_values() stores values of the trigamma function. X trigamma(X) 1.000000 1.6449340668482260 1.100000 1.4332991507927590 1.200000 1.2673772054237791 1.300000 1.1342534349966189 1.400000 1.0253565905295969 1.500000 0.9348022005446793 1.600000 0.8584318931245799 1.700000 0.7932328301639984 1.800000 0.7369741375017002 1.900000 0.6879720582426356 2.000000 0.6449340668482264 uniform_01_cdf_test(): uniform_01_cdf() evaluates the Uniform 01 CDF uniform_01_cdf_inv() inverts the Uniform 01 CDF. uniform_01_pdf() evaluates the Uniform 01 PDF X PDF CDF CDf_inv 0.854758 1 0.854758 0.854758 0.778721 1 0.778721 0.778721 0.183095 1 0.183095 0.183095 0.782964 1 0.782964 0.782964 0.790172 1 0.790172 0.790172 0.132199 1 0.132199 0.132199 0.45208 1 0.45208 0.45208 0.824666 1 0.824666 0.824666 0.575758 1 0.575758 0.575758 0.229895 1 0.229895 0.229895 uniform_01_sample_test(): uniform_01_mean() computes the Uniform 01 mean uniform_01_sample() samples the Uniform 01 distribution uniform_01_variance() computes the Uniform 01 variance. PDF mean = 0.5 PDF variance = 0.0833333 Sample size = 1000 Sample mean = 0.500878 Sample variance = 0.0822088 Sample maximum = 0.999074 Sample minimum = 0.000431565 uniform_01_order_sample_test(): uniform_order_sample() samples the Uniform 01 Order distribution. Ordered sample: 0: 0.0534981 1: 0.40119 2: 0.429056 3: 0.469158 4: 0.572025 5: 0.650874 6: 0.750412 7: 0.890893 8: 0.894383 9: 0.987504 uniform_cdf_test(): uniform_cdf() evaluates the Uniform CDF uniform_cdf_inv() inverts the Uniform CDF. uniform_pdf() evaluates the Uniform PDF PDF parameter A = 1 PDF parameter B = 10 X PDF CDF CDf_inv 1.53474 0.111111 0.0594158 1.53474 7.23591 0.111111 0.692879 7.23591 1.23307 0.111111 0.0258966 1.23307 3.77055 0.111111 0.307839 3.77055 1.57354 0.111111 0.0637267 1.57354 6.97557 0.111111 0.663953 6.97557 4.42369 0.111111 0.38041 4.42369 1.82912 0.111111 0.0921244 1.82912 2.61971 0.111111 0.179968 2.61971 2.53865 0.111111 0.170961 2.53865 uniform_sample_test(): uniform_mean() computes the Uniform mean uniform_sample() samples the Uniform distribution uniform_variance() computes the Uniform variance. PDF parameter A = 1 PDF parameter B = 10 PDF mean = 5.5 PDF variance = 6.75 Sample size = 1000 Sample mean = 5.4531 Sample variance = 6.87191 Sample maximum = 9.98592 Sample minimum = 1.00936 uniform_discrete_cdf_test(): uniform_discrete_cdf() evaluates the Uniform Discrete CDF uniform_discrete_cdf_inv() inverts the Uniform Discrete CDF. uniform_discrete_pdf() evaluates the Uniform Discrete PDF PDF parameter A = 1 PDF parameter B = 6 X PDF CDF CDf_inv 3 0.166667 0.5 4 5 0.166667 0.833333 6 3 0.166667 0.5 4 4 0.166667 0.666667 5 2 0.166667 0.333333 3 5 0.166667 0.833333 6 3 0.166667 0.5 4 3 0.166667 0.5 4 2 0.166667 0.333333 3 1 0.166667 0.166667 2 uniform_discrete_sample_test(): uniform_discrete_mean() computes the Uniform Discrete mean uniform_discrete_sample() samples the Uniform Discrete distribution uniform_discrete_variance() computes the Uniform Discrete variance. PDF parameter A = 1 PDF parameter B = 6 PDF mean = 3.5 PDF variance = 2.91667 Sample size = 1000 Sample mean = 3.532 Sample variance = 2.94098 Sample maximum = 6 Sample minimum = 1 uniform_nsphere_sample_test(): uniform_nsphere_sample() samples the Uniform Nsphere distribution. Dimension N of sphere = 3 Points on the sphere: -0.6738 -0.736473 -0.0600108 0.816109 -0.507201 -0.276972 0.517105 0.850045 0.100133 -0.390761 0.250376 -0.885787 -0.407681 0.904647 -0.12414 0.273119 -0.0989639 -0.956876 0.0779302 0.687027 -0.722441 -0.956166 -0.193383 -0.219885 0.327232 -0.371539 0.868837 -0.577268 -0.814865 -0.0525102 von_mises_cdf_test(): von_mises_cdf() evaluates the Von Mises CDF. von_mises_cdf_inv() inverts the Von Mises CDF. von_mises_pdf() evaluates the Von Mises PDF. PDF parameter A = 1 PDF parameter B = 2 X PDF CDF CDf_inv 0.189331 0.276971 0.155949 0.189291 1.42213 0.432818 0.705601 1.42223 0.27751 0.313001 0.181953 0.277495 0.684747 0.467463 0.342572 0.684767 -0.576173 0.0690708 0.0371501 -0.576932 1.37682 0.448346 0.685632 1.37698 1.37903 0.447616 0.686621 1.37889 2.0832 0.178198 0.90553 2.08299 1.32679 0.464073 0.662801 1.32674 1.28078 0.477024 0.641144 1.28072 von_mises_sample_test(): von_mises_mean() computes the Von Mises mean von_mises_sample() samples the Von Mises distribution. von_mises_circular_variance() computes the Von Mises circular variance PDF parameter A = 1 PDF parameter B = 2 PDF mean = 1 PDF circular variance = 0.302225 Sample size = 1000 Sample mean = 0.983071 Sample circular variance = 0.294105 Sample maximum = 4.12839 Sample minimum = -2.1298 weibull_cdf_test(): weibull_cdf() evaluates the Weibull CDF weibull_cdf_inv() inverts the Weibull CDF. weibull_pdf() evaluates the Weibull PDF PDF parameter A = 2 PDF parameter B = 3 PDF parameter C = 4 X PDF CDF CDf_inv 4.05575 0.344133 0.197877 4.05575 4.55178 0.486146 0.407536 4.55178 4.828 0.507077 0.545993 4.828 4.85037 0.506243 0.55733 4.85037 3.09157 0.0631124 0.0173747 3.09157 4.12048 0.366841 0.220892 4.12048 3.55535 0.172855 0.0697001 3.55535 4.53979 0.484031 0.40172 4.53979 4.59243 0.492631 0.427432 4.59243 3.93919 0.302423 0.16019 3.93919 weibull_sample_test(): weibull_mean() computes the Weibull mean weibull_sample() samples the Weibull distribution weibull_variance() computes the Weibull variance. PDF parameter A = 2 PDF parameter B = 3 PDF parameter C = 4 PDF mean = 4.71921 PDF variance = 0.581953 Sample size = 1000 Sample mean = 4.70408 Sample variance = 0.568552 Sample maximum = 7.2369 Sample minimum = 2.45623 weibull_discrete_cdf_test(): weibull_discrete_cdf() evaluates the Weibull Discrete CDF weibull_discrete_cdf_inv() inverts the Weibull Discrete CDF. weibull_discrete_pdf() evaluates the Weibull Discrete PDF PDF parameter A = 0.5 PDF parameter B = 1.5 X PDF CDF CDf_inv 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 1 0.359214 0.859214 1 2 0.113508 0.972723 3 1 0.359214 0.859214 1 weibull_discrete_sample_test(): weibull_discrete_sample() samples the Weibull Discrete distribution PDF parameter A = 0.5 PDF parameter B = 1.5 Sample size = 1000 Sample mean = 1.15 Sample variance = 0.1695 Sample maximum = 4 Sample minimum = 1 zipf_cdf_test(): zipf_pdf() evaluates the Zipf PDF. zipf_cdf() evaluates the Zipf CDF. zipf_cdf_inv() inverts the Zipf CDF. PDF parameter A = 2 X PDF(X) CDF(X) CDf_inv(CDF) 1 0.607927 0.607927 1 2 0.151982 0.759909 2 3 0.0675475 0.827456 3 4 0.0379954 0.865452 4 5 0.0243171 0.889769 5 6 0.0168869 0.906656 6 7 0.0124067 0.919062 7 8 0.00949886 0.928561 8 9 0.00750527 0.936067 9 10 0.00607927 0.942146 10 11 0.00502419 0.94717 11 12 0.00422172 0.951392 12 13 0.0035972 0.954989 13 14 0.00310167 0.958091 14 15 0.0027019 0.960792 15 16 0.00237472 0.963167 16 17 0.00210355 0.965271 17 18 0.00187632 0.967147 18 19 0.00168401 0.968831 19 20 0.00151982 0.970351 20 zipf_sample_test(): zipf_mean() returns the mean of the Zipf distribution. zipf_sample() samples the Zipf distribution. zipf_variance() returns the variance of the Zipf distribution. PDF parameter A = 4 PDF mean = 1.11063 PDF variance = 0.286326 Sample size = 1000 Sample mean = 1.118 Sample variance = 0.220076 Sample maximum = 7 Sample minimum = 1 prob_test(): Normal end of execution. Mon Nov 14 18:19:08 2022