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