diff --git a/nucleons.jl b/nucleons.jl
index d465af4..4796bbf 100644
--- a/nucleons.jl
+++ b/nucleons.jl
@@ -86,7 +86,7 @@ function solveNucleonWf(κ, p::Bool, E, s::system; shooting=true, normalize=true
     end
 
     if normalize
-        wf ./= norm(wf) * Δr(s) # integration by Reimann sum
+        wf ./= norm(wf) * sqrt(Δr(s)) # integration by Reimann sum
     end
 
     return wf