Function rgsl::randist::f_distribution::fdist

pub fn fdist(r: &Rng, nu1: f64, nu2: f64) -> f64

This function returns a random variate from the F-distribution with degrees of freedom nu1 and nu2. The distribution function is,

p(x) dx =
{ Gamma((\nu_1 + \nu_2)/2)

        over Gamma(nu_1/2) Gamma(nu_2/2) }

\nu_1^{\nu_1/2} \nu_2^{\nu_2/2}

   x^{\nu_1/2 - 1} (\nu_2 + \nu_1 x)^{-\nu_1/2 -\nu_2/2}

for x >= 0.