code 35182;
  procedure BESS YAPLUSN(A, X, NMAX, YAN); value A, X, NMAX;
  real A, X; integer NMAX; array YAN;
  begin integer N; real Y1;
    BESS YA01(A, X, YAN[0], Y1); A:= A-1; X:= 2/X;
    if NMAX > 0 then YAN[1]:= Y1;
    for N:= 2 step 1 until NMAX do 
    YAN[N]:= -YAN[N-2] + (A+N)*X*YAN[N-1]
  end BESS YAPLUSN;
        eop