Add regularized incomplete beta and gamma functions.

Try to use them and their auxiliary functions as much possible. It is possible that we can exploit the auxiliary functions more. I have also added discrete and continuous quantile fallbacks as well as a rand by inversion fallback.