JavaScript实现伽马函数（gamma()）

前言

伽玛函数（Gamma函数），也叫欧拉第二积分，是阶乘函数在实数与复数上扩展的一类函数。该函数在分析学、概率论、偏微分方程和组合数学中有重要的应用。与之有密切联系的函数是贝塔函数，也叫第一类欧拉积分。可以用来快速计算同伽马函数形式相类似的积分。

伽玛函数不是初等函数，而是用积分形式定义的超越函数，怎么看都让人觉得不如初等函数自然亲切。然而伽玛函数也被称为阶乘函数，高等数学会告诉我们一个基本结论：伽玛函数是阶乘的推广。

伽马函数（Gamma）

由于伽玛函数在整个实数轴上都有定义，于是可以看做阶乘概念在实数集上的延拓。

$$gamma(x)={\int }_0^{\infty }{t}^{x-1}{e}^{-t}dt$$

gamma(x) = integral from 0 to inf of t^(x-1) exp(-t) dt.

• 当x为int时
  gamma(x)=(x-1)!
• 当x为float或double时
  gamma(x)是积分（面积）

令$t=\frac{t}{1-t}$则
$${\int }_0^{\infty }{t}^{x-1}{e}^{-t}dt={\int }_0^1(\frac{t}{1-t}{)}^{x-1}{e}^{-\frac{t}{1-t}}\frac{dt}{(1-t{)}^2}$$

实现

function gamma(x) {
  if (x < 0.5) {
    // 使用反正切公式计算伽马函数
    return Math.PI / (Math.sin(Math.PI * x) * gamma(1 - x));
  } else {
    // 使用连乘积公式计算伽马函数
    let y = 1;
    while (x > 1.5) {
      x--;
      y *= x;
    }
    return y * Math.sqrt(Math.PI) * Math.pow(2, x - 1) * gamma(x - 0.5);
  }
}

在这个实现中，如果x小于0.5，则使用反正切公式计算伽马函数；否则，使用连乘积公式计算伽马函数。其中，反正切公式和连乘积公式是两种常用的伽马函数计算方法。

使用这个函数，您可以计算任意x值的伽马函数，例如：

console.log(gamma(0.5)); // 1.77245385091
console.log(gamma(1)); // 1
console.log(gamma(2)); // 1
console.log(gamma(3)); // 2
console.log(gamma(4)); // 6
console.log(gamma(5)); // 24

需要注意的是，伽马函数在计算时可能会出现数值不稳定的情况，特别是在x较小或较大时。因此，在实际应用中，需要对计算结果进行检验和调整，以确保计算的正确性和稳定性。