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