bessy.txt

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

Function BESSY(var N:integer; X:real):real;
var
    TOX,BY,BYM,BYP:real;j:integer;
begin
    If N < 2 Then
        ShowMessage('bad argument N in BESSY');
    TOX:= 2 / X;
    By:= BESSY1(X);
    BYM:= BESSY0(X);
    For J:= 1 To N - 1 do
    begin
        BYP:= J * TOX * By - BYM;
        BYM:= By;
        By:= BYP;
    end;
    BESSY:= By;
end;