      SUBROUTINE BSSLZ1(BES,N)
C
      DIMENSION BES(N)
      DIMENSION BZ(50)
C
      DATA PI/3.14159265358979E0/
      DATA BZ         / 2.4048255577E0, 5.5200781103E0,
     $  8.6537279129E0,11.7915344391E0,14.9309177086E0,18.0710639679E0,
     $ 21.2116366299E0,24.3524715308E0,27.4934791320E0,30.6346064684E0,
     $ 33.7758202136E0,36.9170983537E0,40.0584257646E0,43.1997917132E0,
     $ 46.3411883717E0,49.4826098974E0,52.6240518411E0,55.7655107550E0,
     $ 58.9069839261E0,62.0484691902E0,65.1899648002E0,68.3314693299E0,
     $ 71.4729816036E0,74.6145006437E0,77.7560256304E0,80.8975558711E0,
     $ 84.0390907769E0,87.1806298436E0,90.3221726372E0,93.4637187819E0,
     $ 96.6052679510E0,99.7468198587E0,102.888374254E0,106.029930916E0,
     $ 109.171489649E0,112.313050280E0,115.454612653E0,118.596176630E0,
     $ 121.737742088E0,124.879308913E0,128.020877005E0,131.162446275E0,
     $ 134.304016638E0,137.445588020E0,140.587160352E0,143.728733573E0,
     $ 146.870307625E0,150.011882457E0,153.153458019E0,156.295034268E0/
      NN=N
      SELECT CASE (N)
       CASE(:50)
         BES(1:NN)=BZ(1:N)
       CASE(51:)
         BES(50)=BZ(50)
         BES(51:N)=BES(50:N-1)+PI
         BES(1:49)=BZ(1:49)
      END SELECT
      RETURN
      END
