pub fn exponential_decay(x: f64, a: f64) -> f64 {
    f64::exp(-a * x)
}