Function cdf 

pub fn cdf(f: f64, d1: f64, d2: f64) -> f64

cdf P(X ≤ x) = I_x(d₁/2, d₂/2), x = d₁f/(d₁f + d₂).