function g = gamma_m(m,tol)

% Find n, such that m*(m+1)*...*(n-1)>1/tol
alpha = m; k = 1;
while alpha<1/tol
alpha = alpha*(m+k);
k = k+1;
end
n = m+k+1


% Backward-iteration
a = 0;
for k=n-1:-1:m
    a=(a+1)/k;
end
g=a;