BSSLZ1.inc

1       SUBROUTINE BSSLZ1(BES,N)
2 C
3       IMPLICIT REAL*8(A-H,O-Z)
4       DIMENSION BES(N)
5       DIMENSION BZ(50)
6 C
7       DATA BZ         / 2.40482555770, 5.52007811030,
8      $  8.65372791290,11.79153443910,14.93091770860,18.07106396790,
9      $ 21.21163662990,24.35247153080,27.49347913200,30.63460646840,
10      $ 33.77582021360,36.91709835370,40.05842576460,43.19979171320,
11      $ 46.34118837170,49.48260989740,52.62405184110,55.76551075500,
12      $ 58.90698392610,62.04846919020,65.18996480020,68.33146932990,
13      $ 71.47298160360,74.61450064370,77.75602563040,80.89755587110,
14      $ 84.03909077690,87.18062984360,90.32217263720,93.46371878190,
15      $ 96.60526795100,99.74681985870,102.8883742540,106.0299309160,
16      $ 109.1714896490,112.3130502800,115.4546126530,118.5961766300,
17      $ 121.7377420880,124.8793089130,128.0208770050,131.1624462750,
18      $ 134.3040166380,137.4455880200,140.5871603520,143.7287335730,
19      & 146.8703076250,150.0118824570,153.1534580190,156.2950342680/
20       NN=N
21       IF(N.LE.50) GO TO 12
22       BES(50)=BZ(50)
23       DO 5 J=51,N
24     5 BES(J)=BES(J-1)+PI
25       NN=49
26    12 DO 15 J=1,NN
27    15 BES(J)=BZ(J)
28       RETURN
29       END SUBROUTINE BSSLZ1