plot (Γ(1 + b) - Γ(1 + b, -log(b))) (-log(b))^(-1 - b) from b=0 to 10