diff --git a/nucleons.jl b/nucleons.jl
index 372c717..35f61df 100644
--- a/nucleons.jl
+++ b/nucleons.jl
@@ -96,7 +96,7 @@ end
 "Find all bound energies between E_min (=850.0) and E_max (=938.0) where
     tol_digits determines the precision for root finding and the threshold for identifying duplicates,
     the other parameters are the same from dirac!(...)."
-function findEs(κ, p::Bool, s::system, E_min=850.0, E_max=938.0, tol_digits=5)
+function findEs(κ, p::Bool, s::system, E_min=850.0, E_max=938.0; tol_digits=8)
     func = determinantFunc(κ, p, s)
     Es = find_all_zeros(func, E_min, E_max; partitions=20, tol=1/10^tol_digits)
     return unique(E -> round(E; digits=tol_digits), Es)