code 35151;
  comment SPHERICAL BESSEL FUNCTIONS Y[.5](X),   , Y[N+.5](X);
  procedure SPHER BESS Y(X, N, Y); value X, N;
  integer N; real X; array Y;
  if N=0 then  Y[0]:= - COS(X)/X else 
  begin real YI, YI1, YI2; integer I;
    YI2:= Y[0]:= -COS(X)/X; YI1:= Y[1]:= (YI2 - SIN(X))/X;
    for I:= 2 step 1 until N do 
    begin Y[I]:= YI:= -YI2 + (I+I-1) * YI1/X;
      YI2:= YI1; YI1:= YI
    end 
  end 

        eop