From c15dafb8da632cf6e60f548374bfbec754232ec4 Mon Sep 17 00:00:00 2001
From: eduard-kuhn <164213309+eduard-kuhn@users.noreply.github.com>
Date: Fri, 31 Jan 2025 15:16:21 +0100
Subject: [PATCH] Update ct_distributions.jl

---
 src/ct_distributions.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/ct_distributions.jl b/src/ct_distributions.jl
index 23f514e3..1c283e13 100644
--- a/src/ct_distributions.jl
+++ b/src/ct_distributions.jl
@@ -54,7 +54,7 @@ The incomplete Fermi-Dirac integral of order 1/2, implemented according to [Bedn
 function FermiDiracOneHalfBednarczyk(x::Real)
 
     a = x^4 + 33.6*x * (1.0 - 0.68*exp(-0.17*(x+1)^2)) + 50
-    sqrt(pi) / ( 2 * (3/4*sqrt(pi) * a^(-3/8) + exp(-x) ) )
+    1.0 / (3/4*sqrt(pi) * a^(-3/8) + exp(-x))
 
 end