code 35053;
procedure IXQFIX(X,P,Q,NMAX,EPS,I); value X,P,Q,NMAX,EPS;
real X,P,Q,EPS; integer NMAX; array I;
begin integer M,MMAX; real S,IQ0,IQ1,Q0;

    M:= ENTIER(Q); S:= Q-M; Q0:= if S>0 then S else S+1;
    MMAX:= if S>0 then M else M-1;
    IQ0:= INCBETA(X,P,Q0,EPS);
    if MMAX>0 then IQ1:= INCBETA(X,P,Q0+1,EPS);
    begin array IQ[0:MMAX];
        FORWARD(X,P,Q0,IQ0,IQ1,MMAX,IQ);
        BACKWARD(X,P,Q,IQ[MMAX],NMAX,EPS,I)
    end 
end IXQFIX

        eop