.


           


6.1.

,

. , . .

k = 0,

.

k = 0, c

n, n.

.

Function

Extent(x : real; n : integer) : real;

var

i : integer;

E : real;

begin

E := 1;

if n = 0 then E := 1

else for

i := 1 to n do E := E*x/2;

Extent := E

end;

.

Function Fakt(n : integer) : longint;

var

i : integer;

F : longint;

begin

F := 1;

if n = 0 then F := 1

else for

i := 1 to n do F := F*i;

Fakt := F

end;

Function J(x, eps : real; n : integer) : real;

var

y, jj : real;

k : integer;

begin

k := 0; y := Extent(x, n)/Fakt(n); jj := 0;

repeat

jj := jj + y;

k := k + 1;

y := -y*x*x/(4*k*(n + k))

until abs(y) < eps;

J := jj

end;

Program Bessel;

uses WinCrt;

var

n : integer;

x, eps : real;

{----------------------------------------------------------------------------------------}

Function

t(eps : real) : integer;

var

k : integer;

begin

k := -1;

repeat