Add regularized incomplete beta and gamma functions. Use them and their auxiliary functions as much possible.