From 2ee76f2c1bf5416deb2d2a5f2d8dbd0698a79e4c Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Tue, 4 Feb 2025 14:08:30 -0500
Subject: [PATCH] Normalization fix

---
 nucleons.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/nucleons.jl b/nucleons.jl
index 0601768..18d791d 100644
--- a/nucleons.jl
+++ b/nucleons.jl
@@ -53,7 +53,7 @@ function solveNucleonWf(κ, p, E, Φ0, W0, B0, A0, r_max, divs; shooting=true, n
 
     if normalize
         norm = sum(wf .* wf) * r_max / divs # integration by Reimann sum
-        wf = wf ./ sqrt(norm)
+        wf = wf ./ sqrt(norm * 2) # WHY FACTOR OF 2?
     end
 
     return wf