Several formulas can approximate the gamma function on the complex domain, but all of them have issues. Please, look at its plot before starting. Nonetheless, there is a trick that makes them all work; and you must find it (I've lost some time trying different formulas, adding more terms, but it is not the way to go).
I've found a 1968-paper by Wrench, Concerning Two Series for the Gamma Function, which gives us the 20th first terms of one formula for instance. https://www.ams.org/journals/mcom/1968-22-103/S0025-5718-1968-0237078-4/ (It would be more than enough if there weren't issues with all approximations.)
A clever and simple approach to deal with the pitfalls in the test suite. It can be easily extended to deal with other problematic values of the argument.
It took me a while to realize the same thing...
815 Solvers
Determine Whether an array is empty
653 Solvers
220 Solvers
Side of an equilateral triangle
2884 Solvers
228 Solvers
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!