function q=find_q(delta)
% q=find_q(delta)
% Returns a value q such that equation for epsilon(q,q)<delta 
% on page  531 in GVL is satisfied. 
%

for q=1:100
    p=q;
    e=2^(3-(p+q))*(factorial(p)*factorial(q))/(factorial(p+q)*factorial(p+q+1));

    if (e<delta)
        return;
    end
end