code 35051;
procedure IBPPLUSN(X,P,Q,NMAX,EPS,I); value X,P,Q,NMAX,EPS;
integer NMAX; real X,P,Q,EPS; array I;
begin integer N;

    if X=0 or X=1 then 
    begin for N:= 0 step 1 until NMAX do I[N]:= X end 
    else 
    begin if X <=.5 then IXQFIX(X,P,Q,NMAX,EPS,I) else 
        begin IXPFIX(1-X,Q,P,NMAX,EPS,I);
            for N:= 0 step 1 until NMAX do I[N]:= 1-I[N]
        end 
    end 
end IBPPLUSN

        eop