pub fn exponential_Pinv(P: f64, mu: f64) -> f64