fn gaussian_g(kappa: f64, tau: f64) -> f64
(1 − e^{−κτ}) / κ, τ itself at zero mean reversion.
(1 − e^{−κτ}) / κ
τ