Increase divs
This commit is contained in:
parent
9e24943a2d
commit
e389adf885
|
|
@ -21,7 +21,7 @@ N = 82
|
||||||
r_max = maximum(xs)
|
r_max = maximum(xs)
|
||||||
E_min = 880
|
E_min = 880
|
||||||
E_max = 939
|
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)
|
(ρ_s, ρ_v) = calculateNucleonDensity(N, p, S_interp, V_interp, R_interp, A_interp, r_max, divs, E_min, E_max)
|
||||||
|
|
||||||
|
|
|
||||||
|
|
@ -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
|
groundE = findEs(κ, p, S_interp, V_interp, R_interp, A_interp, r_max, E_min, E_max) |> minimum
|
||||||
println("ground state E = $groundE")
|
println("ground state E = $groundE")
|
||||||
|
|
||||||
divs = 50
|
divs = 400
|
||||||
wf = solveWf(κ, p, groundE, S_interp, V_interp, R_interp, A_interp, r_max, divs)
|
wf = solveWf(κ, p, groundE, S_interp, V_interp, R_interp, A_interp, r_max, divs)
|
||||||
rs = range(0, r_max, length=divs+1)
|
rs = range(0, r_max, length=divs+1)
|
||||||
gs = wf[1, :]
|
gs = wf[1, :]
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue