Better error reporting for test case
Use incomplete beta as approximation for CDF of binomial distribution

It offer lower precision for some inputs but much faster for large N
and works for very large N

Fixes #53