数学算法

				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;