Need a robust implementation of the 2F1 hypergeometric function.
Currently, the t distribution CDF is calculated with a numerical approximation of the integral of the PDF. This is extremely slow, in particular when calculating the PPF, which executes CDF many times.