pub fn prs<F>( f: F, n: usize, (a, b): (f64, f64), m: usize, eps: f64 ) -> Vec<f64>where F: Fn(f64) -> f64 + Copy,
Piecewise Rejection Sampling