An Interpolatory Method of Approximating Distribution Functions with Applications to the Gamma and Related Variates