code 35153;
  comment MODIFIED SPHERICAL BESSEL FUNCTIONS
      K[.5](X),   , K[N+.5](X);
  procedure SPHER BESS K(X, N, K); value X, N;
  integer N; real X; array K;
  begin real EXPX;
    EXPX:= EXP(-X);
    NONEXP SPHER BESS K(X, N, K);
    for N:= N step -1 until 0 do K[N]:= K[N] * EXPX
  end 

        eop