bessk.txt

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

Function BESSK(N:integer; X:real):real;
var
    TOX,BKM,BK,BKP:real; J:integer;
begin
    If N < 2 Then
        ShowMessage('bad argument N in BESSK');
    TOX:=2 / X;
    BKM:=BESSK0(X);
    BK:=BESSK1(X);
    For J:=1 To N - 1 do
    begin
        BKP:=BKM + J * TOX * BK;
        BKM:=BK;
        BK:=BKP;
    end;
    BESSK:=BK;
end;