Function could return answers which are orders of magnitude off. It could be due to the fact that we revert to bisection which is linear instead of cubic Halley.
>>> (\a b -> incompleteBeta a b . invIncompleteBeta a b) 90 90 1e-200
9.866213992257607e-197
Function could return answers which are orders of magnitude off. It could be due to the fact that we revert to bisection which is linear instead of cubic Halley.