213 attempted solution

I tried to solve #213, the result is compatible with the mathematica result even if some ratios of gamma functions are not simplified properly
Why do we need a seperate entry for `IncompleteGamma`? Isn't this already handled by `Gamma`. Something to do with `mpmath` evaluation? Might be worth adding a comment here to explain.