1 | "use strict";
|
2 |
|
3 | Object.defineProperty(exports, "__esModule", {
|
4 | value: true
|
5 | });
|
6 | exports.gammaNumber = gammaNumber;
|
7 | exports.gammaP = exports.gammaG = void 0;
|
8 |
|
9 | var _number = require("../../utils/number.js");
|
10 |
|
11 | var _product = require("../../utils/product.js");
|
12 |
|
13 |
|
14 | function gammaNumber(n) {
|
15 | var x;
|
16 |
|
17 | if ((0, _number.isInteger)(n)) {
|
18 | if (n <= 0) {
|
19 | return isFinite(n) ? Infinity : NaN;
|
20 | }
|
21 |
|
22 | if (n > 171) {
|
23 | return Infinity;
|
24 | }
|
25 |
|
26 | return (0, _product.product)(1, n - 1);
|
27 | }
|
28 |
|
29 | if (n < 0.5) {
|
30 | return Math.PI / (Math.sin(Math.PI * n) * gammaNumber(1 - n));
|
31 | }
|
32 |
|
33 | if (n >= 171.35) {
|
34 | return Infinity;
|
35 | }
|
36 |
|
37 | if (n > 85.0) {
|
38 |
|
39 | var twoN = n * n;
|
40 | var threeN = twoN * n;
|
41 | var fourN = threeN * n;
|
42 | var fiveN = fourN * n;
|
43 | return Math.sqrt(2 * Math.PI / n) * Math.pow(n / Math.E, n) * (1 + 1 / (12 * n) + 1 / (288 * twoN) - 139 / (51840 * threeN) - 571 / (2488320 * fourN) + 163879 / (209018880 * fiveN) + 5246819 / (75246796800 * fiveN * n));
|
44 | }
|
45 |
|
46 | --n;
|
47 | x = gammaP[0];
|
48 |
|
49 | for (var i = 1; i < gammaP.length; ++i) {
|
50 | x += gammaP[i] / (n + i);
|
51 | }
|
52 |
|
53 | var t = n + gammaG + 0.5;
|
54 | return Math.sqrt(2 * Math.PI) * Math.pow(t, n + 0.5) * Math.exp(-t) * x;
|
55 | }
|
56 |
|
57 | gammaNumber.signature = 'number';
|
58 |
|
59 | var gammaG = 4.7421875;
|
60 | exports.gammaG = gammaG;
|
61 | var gammaP = [0.99999999999999709182, 57.156235665862923517, -59.597960355475491248, 14.136097974741747174, -0.49191381609762019978, 0.33994649984811888699e-4, 0.46523628927048575665e-4, -0.98374475304879564677e-4, 0.15808870322491248884e-3, -0.21026444172410488319e-3, 0.21743961811521264320e-3, -0.16431810653676389022e-3, 0.84418223983852743293e-4, -0.26190838401581408670e-4, 0.36899182659531622704e-5];
|
62 | exports.gammaP = gammaP; |
\ | No newline at end of file |