Function sample_truncated_normal 

pub fn sample_truncated_normal(
    rng: &mut Xoshiro256,
    mu: f64,
    sigma: f64,
    a: f64,
    b: f64,
) -> f64

Sample from a truncated normal distribution N(mu, sigma²) on [a, b] using rejection sampling.