Simplified initial C

1 files changed, 1 insertions, 1 deletions

diff --git a/log-fourier_integrator.cpp b/log-fourier_integrator.cpp
index f553569..30354a6 100644
--- a/log-fourier_integrator.cpp
+++ b/log-fourier_integrator.cpp
@@ -114,7 +114,7 @@ int main(int argc, char* argv[]) {
         if (τ₀ == 2) {
           Cₜ₋₁[n] = Γ₀ * std::exp(-fft.t(n) / 2) * (1 + fft.t(n) / 2);
         } else {
-          Cₜ₋₁[n] = Γ₀ * (std::exp(-μ₀ * fft.t(n)) - μ₀ * τ₀ * std::exp(-fft.t(n) / τ₀)) / (μ₀ - pow(μ₀, 3) * pow(τ₀, 2));
+          Cₜ₋₁[n] = Γ₀ * (std::exp(-μ₀ * fft.t(n)) / μ₀ - τ₀ * std::exp(-fft.t(n) / τ₀)) / (1 - pow(μ₀ * τ₀, 2));
         }
       } else {
         Cₜ₋₁[n] = Γ₀ * std::exp(-μ₀ * fft.t(n)) / μ₀;