pub fn gaussian_int(n: f64, alpha: f64) -> f64
Computes the Gaussian integral:
$$ \int_0^{+\inf} x^n e^{-\alpha x^2} dx = \frac{n!}{2\alpha^{(n+1)/2}} $$