bessk.txt

来自「《Delphi常用数值算法集》的配书源码」· 文本 代码 · 共 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;