yapanuwan
/
NuclearRMF
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Close solution found by trial and error

This commit is contained in:
Nuwan Yapa 2025-01-30 19:15:21 -05:00
parent c9ec9c6203
commit 8a42579594
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
mesons.jl
View File

@ -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