pub fn number_state_probability(coherent_alpha: Complex, n: u32) -> f64
Probability of finding n photons in a coherent state |alpha>: P(n) = |alpha|^{2n} e^{-|alpha|^2} / n!