From db3a9304e795341adc8eedc6afa782c1dc61a344 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Fri, 17 Jan 2025 12:56:01 -0500
Subject: [PATCH] Refactor

---
 mesons.jl | 5 +----
 1 file changed, 1 insertion(+), 4 deletions(-)

diff --git a/mesons.jl b/mesons.jl
index 9c73d44..675a34a 100644
--- a/mesons.jl
+++ b/mesons.jl
@@ -37,10 +37,7 @@ function mesons!(ddu, du, u, (ρ_sp, ρ_vp, ρ_sn, ρ_vn), r)
     ddB0 = -(2/(r + r_reg)) * dB0 + m_ρ^2 * B0 / ħc + g2_ρ * (2 * Λv * W0^2 * B0 - (ρ_vp(r) - ρ_vn(r)) / 2) / ħc
     ddA0 = -(2/(r + r_reg)) * dA0 - ρ_vp(r) / ħc # e goes here?
 
-    ddu[1] = ddΦ0
-    ddu[2] = ddW0
-    ddu[3] = ddB0 
-    ddu[4] = ddA0
+    ddu .= [ddΦ0, ddW0, ddB0, ddA0]
 end
 
 "Solve meson equations and return the wave functions u(r)=[Φ0(r), W0(r), B0(r), A0(r)] where