factrl.txt

来自「1.г函数、贝塔函数、阶乘及二项式系数 2.不完全г函数、误差函数 3.不完」· 文本 代码 · 共 21 行

Function FACTRL(N:integer):REAL;
var
  GLNTOP:integer;
  GLA:ARRAY[1..33] OF real;
  J:integer;
BEGIN
  GLNTOP:=0;
  GLA[1]:= 1;
  If N < 0 Then  ShowMessage('警告：现在N<0，请修改程序');
  If N = 0 Then
    FACTRL:= GLA[N+1]
  else If N <= 32 Then
  begin
      For J:= GLNTOP + 1 To N do
        GLA[J + 1]:= J * GLA[J];
      GLNTOP:= N;
      FACTRL:= GLA[N + 1];
  end
  else If N>32 then
    FACTRL:= Exp(GAMMLN(N + 1));
END;