Done! No more epsilon! The result is as accurate as possible given a double representation.

*BTW Tomek Czajka pointed out that my final version is buggy as well. It doesn't work for c = 1e60. The problem is the number of iterations is too small, and for it to work correctly it needs about 1100 steps.*

*BTW Tomek Czajka(of topcoder fame) pointed out that my final version is buggy as well. It doesn't work for c = 1e60. The problem is the number of iterations is too small, and for it to work correctly it needs about 1100 steps.*

We’ve addressed negative numbers, numbers in (0, 1), overflow problems, absolute and relative error problems.