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