From 8a425795941d5cdf0cac96742e55567fec70c96a Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Thu, 30 Jan 2025 19:15:21 -0500 Subject: [PATCH] Close solution found by trial and error --- mesons.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/mesons.jl b/mesons.jl index 0a7ec90..f0a4602 100644 --- a/mesons.jl +++ b/mesons.jl @@ -56,13 +56,13 @@ function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3; in A0 .= solveKG(m_γ, src_A0, r_max) for _ in 1:iterations - @. src_Φ0 = g2_s * ((κ/2) * (Φ0/ħc)^2 + (λ/6) * (Φ0/ħc)^3) - g2_s * (ρ_sp + ρ_sn) * ħc + @. src_Φ0 = g2_s * ((κ/2) * (Φ0/ħc)^2 + (λ/6) * (Φ0/ħc)^3) - g2_s * ħc * (ρ_sp + ρ_sn) Φ0 .= solveKG(m_s, src_Φ0, r_max) - @. src_W0 = g2_v * ((ζ/6) * (W0/ħc)^3 + 2 * Λv * (B0/ħc)^2 * (W0/ħc)) - g2_v * (ρ_vp + ρ_vn) * ħc + @. src_W0 = g2_v * ħc * ((ζ/6) * (W0/ħc)^3 + 2Λv * (B0/ħc)^2 * (W0/ħc)) - g2_v * (ρ_vp + ρ_vn) * ħc W0 .= solveKG(m_ω, src_W0, r_max) - @. src_B0 = g2_ρ * (2 * Λv * (W0/ħc)^2 * (B0/ħc)) - (g2_ρ / 2) * (ρ_vp - ρ_vn) * ħc + @. src_B0 = 2Λv * g2_ρ * W0^2 * B0 / ħc^2 - (g2_ρ / 4) * (ρ_vp - ρ_vn) * ħc B0 .= solveKG(m_ρ, src_B0, r_max) end