diff --git a/test/Pb208_density.jl b/test/Pb208_density.jl index bb33597..faff79c 100644 --- a/test/Pb208_density.jl +++ b/test/Pb208_density.jl @@ -21,7 +21,7 @@ N = 82 r_max = maximum(xs) E_min = 880 E_max = 939 -divs = 50 +divs = 400 (ρ_s, ρ_v) = calculateNucleonDensity(N, p, S_interp, V_interp, R_interp, A_interp, r_max, divs, E_min, E_max) diff --git a/test/Pb208_wf.jl b/test/Pb208_wf.jl index cebc05f..1a8a5ee 100644 --- a/test/Pb208_wf.jl +++ b/test/Pb208_wf.jl @@ -24,7 +24,7 @@ E_max = 939 groundE = findEs(κ, p, S_interp, V_interp, R_interp, A_interp, r_max, E_min, E_max) |> minimum println("ground state E = $groundE") -divs = 50 +divs = 400 wf = solveWf(κ, p, groundE, S_interp, V_interp, R_interp, A_interp, r_max, divs) rs = range(0, r_max, length=divs+1) gs = wf[1, :]