code 35179;
    procedure NONEXP BESS K(X, N, K); value X, N;
    real X; integer N; array K;
    begin  integer I; real K0, K1, K2;
        NONEXP BESS K01(X, K0, K1);
        K[0]:= K0; if N> 0 then K[1]:= K1; X:= 2 / X;
        for I:= 2 step 1 until N do 
        begin K[I]:= K2:= K0 + X * (I-1)* K1;
            K0:= K1; K1:= K2
        end 
    end NONEXP BESS K;
        eop