diff --git a/system.jl b/system.jl
index ac61ab8..89d23ec 100644
--- a/system.jl
+++ b/system.jl
@@ -46,7 +46,7 @@ function solve_system(s::system, initial_dens=nothing, initial_flds=(zeros(1 + s
     E_total_previous = NaN
 
     while true
-        (Φ0s, W0s, B0s, A0s) = solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, s.r_max, s.divs, 50; initial_sol = (Φ0s, W0s, B0s, A0s))
+        (Φ0s, W0s, B0s, A0s) = solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, s.r_max, s.divs, isnan(E_total_previous) ? 50 : 5; initial_sol = (Φ0s, W0s, B0s, A0s))
 
         S_interp = linear_interpolation(rs(s), Φ0s)
         V_interp = linear_interpolation(rs(s), W0s)