code 35164;
    procedure BESS Y(X, N, Y); value X, N;
    real X; integer N; array Y;
    begin integer I; real Y0, Y1, Y2;
      BESS Y01(X, Y0, Y1); Y[0]:= Y0;
      if N > 0 then Y[1]:= Y1 ;
      X:= 2/X;
      for I:=2 step 1 until N do 
      begin Y[I]:= Y2:= (I-1)*X*Y1 - Y0;
        Y0:= Y1; Y1:= Y2
      end 
    end BESS Y

        eop