SUBROUTINE BSSLZ1(BES,N) C IMPLICIT REAL*8(A-H,O-Z) DIMENSION BES(N) DIMENSION BZ(50) C DATA BZ / 2.40482555770, 5.52007811030, $ 8.65372791290,11.79153443910,14.93091770860,18.07106396790, $ 21.21163662990,24.35247153080,27.49347913200,30.63460646840, $ 33.77582021360,36.91709835370,40.05842576460,43.19979171320, $ 46.34118837170,49.48260989740,52.62405184110,55.76551075500, $ 58.90698392610,62.04846919020,65.18996480020,68.33146932990, $ 71.47298160360,74.61450064370,77.75602563040,80.89755587110, $ 84.03909077690,87.18062984360,90.32217263720,93.46371878190, $ 96.60526795100,99.74681985870,102.8883742540,106.0299309160, $ 109.1714896490,112.3130502800,115.4546126530,118.5961766300, $ 121.7377420880,124.8793089130,128.0208770050,131.1624462750, $ 134.3040166380,137.4455880200,140.5871603520,143.7287335730, & 146.8703076250,150.0118824570,153.1534580190,156.2950342680/ NN=N IF(N.LE.50) GO TO 12 BES(50)=BZ(50) DO 5 J=51,N 5 BES(J)=BES(J-1)+PI NN=49 12 DO 15 J=1,NN 15 BES(J)=BZ(J) RETURN END SUBROUTINE BSSLZ1