use super::variance::expect_psi_prime;
use crate::rho::RhoFunction;
pub fn influence_function<'a>(
rho: &'a dyn RhoFunction,
quad_points: usize,
) -> impl Fn(f64) -> f64 + 'a {
let c = expect_psi_prime(rho, quad_points);
move |x: f64| rho.psi(x) / c
}