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);
procedure GCF(var GAMMCF:real; A:real; X:real; var GLN:real);
label 1;
const
ITMAX=100; EPS=0.0000003;
var
N:integer; GOLD,G,FAC,B1,B0,ANF,ANA,AN,A1,A0:real;
begin
GLN:=GAMMLN(A);
G
Function BESSI(N:integer; X:real):real;
const
IACC = 40; BIGNO = 1.e10;
BIGNI = 1.e-10;
var
TOX,BIP,BI,BIM,ANS:real; M,J:integer;
begin
If N < 2 Then
ShowMessage('ba
