code 35155;
  procedure NONEXP SPHER BESS K(X, N, K); value X, N;
  real X; integer N; array K;
  begin integer I; real KI, KI1, KI2;
    X:= 1/X; K[0]:= KI2:= X*1.5707963267949;
    if N=0 then "GO TO" EXIT;
    K[1]:= KI1:= KI2 * (1+X);
    for I:= 2 step 1 until N do 
    begin K[I]:= KI:= KI2 + (I+I-1) * X * KI1;
      KI2:= KI1; KI1:= KI end;
    EXIT:
  end;
        eop