code 35081;
    procedure EIALPHA(X,N,ALPHA);
    value X,N;real X;integer N;array ALPHA;
    begin real A,B,C;integer K;
        C:=1/X;A:=EXP(-X);
        B:=ALPHA[0]:=A*C;
        for K:=1 step 1 until N do 
        ALPHA[K]:=B:=(A+K*B)*C
    end EIALPHA;
        eop