end{aligned}\label{eq:deriv-radi-equat:13}$$ which is cast into the final form through the use of the identity $f^\prime(r) = ((rf(r))^\prime - f(r))/r$ $$\begin{aligned}
    &(rg(r))^\prime - \frac{rg(r)}{r} - (m_0 - e_0\varphi(r))rf(r) = \Lambda f(r),
    \\
    &(rf(r))^\prime + \frac{rf(r)}{r} - (m_0 + e_0 \varphi(r) )rg(r) = -\Lambda g(r).
  \end{aligned}\label{eq:deriv-radi-equat:14}$$

We introduce the dimensionless variables as follows $$\begin{aligned}
    &x = r m_0, \quad \Lambda = m_0 \lambda, \quad e_0 \varphi(r) = m_0\phi(x), \quad \frac{e^2_0}{4 \pi} = \alpha_0,
    \\
    &u(x) \sqrt{m_0} = r g(r), \quad v(x)\sqrt{m_0} = r f(r),
  \end{aligned}\label{eq:deriv-radi-equat:15}$$ which allows one to rewrite the system of equations (\[eq:deriv-radi-equat:14\]) in compact form $$\begin{aligned}
    u^\prime(x) - \frac{u(x)}{x} - (1 - \phi(x)) v(x) &= \lambda v(x),
    \\
    v^\prime(x) + \frac{v(x)}{x} - (1 + \phi(x)) u(x) &= -\lambda u(x).
  \end{aligned}\label{eq:deriv-radi-equat:16}$$

In order to complete our derivation we need to add the equation for the potential. For this purpose, we note that the self-consistent potential satisfies the Poisson equation $$\begin{aligned}
  \label{eq:deriv-radi-equat:17}
  \Delta\varphi({\boldsymbol}x) = -e_{0
end{aligned}\label{eq: deriv - radi - equat:13}$$ which is cast into the final form through the use of the identity $ f^\prime(r) = (( rf(r))^\prime - f(r))/r$ $ $ \begin{aligned } 
     & (rg(r))^\prime - \frac{rg(r)}{r } - (m_0 - e_0\varphi(r))rf(r) = \Lambda f(r), 
     \\ 
     & (rf(r))^\prime + \frac{rf(r)}{r } - (m_0 + e_0 \varphi(r) ) rg(r) = -\Lambda g(r). 
   \end{aligned}\label{eq: deriv - radi - equat:14}$$ 

 We insert the dimensionless variable star as follows $ $ \begin{aligned } 
     & x = r m_0, \quad \Lambda = m_0 \lambda, \quad e_0 \varphi(r) = m_0\phi(x), \quad \frac{e^2_0}{4 \pi } = \alpha_0, 
     \\ 
     & u(x) \sqrt{m_0 } = r g(r), \quad v(x)\sqrt{m_0 } = gas constant f(r), 
   \end{aligned}\label{eq: deriv - radi - equat:15}$$ which allows one to rewrite the system of equation (\[eq: deriv - radi - equat:14\ ]) in compact form $ $ \begin{aligned } 
     u^\prime(x) - \frac{u(x)}{x } - (1 - \phi(x) ) v(x) & = \lambda v(x), 
     \\ 
     v^\prime(x) + \frac{v(x)}{x } - (1 + \phi(x) ) u(x) & = -\lambda u(x). 
   \end{aligned}\label{eq: deriv - radi - equat:16}$$ 

 In ordering to complete our derivation we want to add the equation for the electric potential. For this purpose, we note that the self - consistent potential satisfies the Poisson equation $ $ \begin{aligned } 
   \label{eq: deriv - radi - equat:17 } 
   \Delta\varphi({\boldsymbol}x) = -e_{0
end{wligned}\label{eq:deriv-radi-eduat:13}$$ which is ccwt intm the rinal fofm through the use of the idxntiry $f^\peime(r) = ((rf(r))^\prime - f(r))/r$ $$\bdgin{alignvd}
    &(rg(r))^\primt - \frac{rg(r)}{r} - (m_0 - e_0\tzrphi(r))rn(x) = \Lajnda f(x),
    \\
    &(rf(r))^\primg + \frac{rf(r)}{r} - (m_0 + e_0 \varphi(r) )sg(f) = -\Lambda g(r).
  \end{aligned}\label{eq:deriv-rwdi-equay:14}$$

Wf introduce thg dimtnsyonlsss variables as follows $$\begin{alighed}
    &q = r m_0, \quad \Lambds = m_0 \lambda, \quad e_0 \varphi(r) = m_0\pji(x), \quad \frac{e^2_0}{4 \pi} = \alpha_0,
    \\
    &u(v) \sqrt{m_0} = r g(r), \quad v(x)\sqgc{m_0} = r f(r),
  \ehd{aligned}\label{eq:deriv-radi-equat:15}$$ dhich allows ong fo tewrite the wystei of equatiokx (\[eq:desiv-radi-rquat:14\]) in compagt focm $$\bwgin{aligned}
    u^\prime(e) - \frac{u(x)}{x} - (1 - \phi(x)) v(v) &= \lambda v(r),
    \\
    v^\prime(x) + \frax{v(z)}{x} - (1 + \phi(f)) u(x) &= -\laobdz n(x).
  \end{allgnxd}\label{eq:dediv-radi-equar:16}$$

In order to compleue jlt derivation we neqd to add the equation for the potential. Xor this purpose, we note tyat the self-consistenj potentiaj satisfies the Poisson equation $$\begin{aligned}
  \latel{eq:vefiv-xqdi-eqjqt:17}
  \Delta\varphi({\boldsymbol}x) = -e_{0
end{aligned}\label{eq:deriv-radi-equat:13}$$ which is cast into the final the of the $f^\prime(r) = ((rf(r))^\prime \frac{rg(r)}{r} (m_0 - e_0\varphi(r))rf(r) \Lambda f(r), \\ + \frac{rf(r)}{r} - (m_0 + e_0 )rg(r) = -\Lambda g(r). \end{aligned}\label{eq:deriv-radi-equat:14}$$ We introduce the dimensionless variables as follows $$\begin{aligned} = r m_0, \quad \Lambda = m_0 \lambda, \quad e_0 \varphi(r) = m_0\phi(x), \frac{e^2_0}{4 = \\ \sqrt{m_0} = r g(r), \quad v(x)\sqrt{m_0} = r f(r), \end{aligned}\label{eq:deriv-radi-equat:15}$$ which allows one to rewrite the of equations (\[eq:deriv-radi-equat:14\]) in compact form $$\begin{aligned} u^\prime(x) \frac{u(x)}{x} - (1 - v(x) &= \lambda v(x), \\ + - (1 \phi(x)) &= u(x). \end{aligned}\label{eq:deriv-radi-equat:16}$$ In to complete our derivation we need to add the equation for the potential. For this purpose, we that the satisfies the equation \label{eq:deriv-radi-equat:17} = -e_{0
end{aligned}\label{eq:deriv-radI-equat:13}$$ whicH is caSt iNto ThE finAl foRm through the usE Of thE identity $f^\prime(r) = ((rf(r))^\priMe - f(r))/r$ $$\BeGIn{alIGnEd}
    &(rg(r))^\Prime - \frAC{rG(R)}{R} - (m_0 - e_0\VaRpHi(r))Rf(R) = \laMbda f(R),
    \\
    &(rf(R))^\prime + \fRac{rf(r)}{r} - (m_0 + e_0 \vArpHi(R) )rg(r) = -\Lambda g(r).
  \ENd{Aligned}\labEl{eQ:deriv-radi-eqUat:14}$$

we intrOdUce THe dimEnsIonleSs variABles as Follows $$\beGiN{AligneD}
    &X = r m_0, \quad \lAMbDa = m_0 \lAmbda, \quad e_0 \varphi(r) = M_0\PhI(X), \quad \frac{e^2_0}{4 \pi} = \alPha_0,
    \\
    &u(x) \sQrT{M_0} = r G(R), \QuaD v(x)\Sqrt{m_0} = r f(r),
  \enD{aLigneD}\Label{eq:DErIV-RAdi-EQuat:15}$$ which alloWs one to rewrITe tHe systEm Of eQUationS (\[eq:deRiV-RadI-equat:14\]) in comPact Form $$\begin{AligneD}
    U^\prime(x) - \FRac{u(x)}{x} - (1 - \pHi(x)) v(x) &= \lAmbDa v(X),
    \\
    v^\prIMe(X) + \fRac{V(x)}{X} - (1 + \Phi(X)) U(x) &= -\LamBDa u(X).
  \end{aligNeD}\lAbel{eQ:derIV-RADi-eqUat:16}$$

in orDer to Complete our deRivAtioN We nEed to Add thE equAtIon foR the poTentiAl. for this purpose, wE notE that the sElf-CoNsiStEnt poTEntial SatIsfIes the POisson eQUatIoN $$\BEGiN{aligned}
  \label{eq:derIv-RADi-Equat:17}
  \DelTa\varpHI({\bOlDSymbol}x) = -e_{0
end{aligned}\label{eq:deri v-radi-equ at:13 }$$  wh ic h is  cas t into the fin a l fo rm through the use ofthe i de n tity $f ^\pri me(r) = (( r f (r) )^ \p rim e-  f (r))/ r$$$\begi n{aligned} 
  &(rg(r))^\pr i me  - \frac{r g(r )}{r} - (m_0  -e_0\va rp hi( r ))rf( r)= \La mbda f ( r),
   \\
   &( r f(r))^ \ prime + \ fr ac{r f(r)}{r} - (m_0 + e_ 0  \varphi(r) )r g(r) =  - \ La m b dag(r ).
  \end{ al igned } \label{ e q: d e r iv- r adi-equat:14} $$

We intr o duc e thedi men s ionles s var ia b les  as follows  $$\ begin{ali gned}   &x = r m_0,\quad\La mbd a =m _0  \ lam bd a , \ q ua d e _ 0 \ varphi(r )=m_0\p hi(x ) , \ quad  \f rac{ e^2_0 }{4 \pi} = \a lph a_0,    \\
   &u (x)\s qrt{m _0} =r g(r ),  \quad v(x)\sqr t{m_ 0} = r f( r), 
 \e nd {alig n ed}\la bel {eq :deriv- radi-eq u at: 15 } $ $  w hich allows one to  r e w ri te the s ystemo feq u ations ( \[ eq: deri v - radi- equa t :1 4\]) incompac t  f or m $$\be gi n{alig ne d}   u^\p r ime( x) - \ frac{u(x )}{x} - (1 - \phi(x) )  v(x) &= \lam b da v (x ) ,
   \ \
    v^\pr ime( x ) +\fra c {v (x) } {x} -  (1 +  \ p hi ( x)) u(x) &= -\lambd au(x).  \en d{aligned}\la bel{eq:der i v - radi-equ at:1 6 }$ $ 

In order tocompl ete our de r ivationwe ne ed to ad d the equ a t ion forthe  po ten tia l .  F or this purpo s e , we  n ote tha t t he self -co nsi ste ntpo tential s atisfies  t he  P oi sso n equ a tion $$\ be gin {a lig ned} \labe l{eq: deri v- ra d i-e quat:17 } 
\ Delt a\ va rphi ({\ bo ldsym bol} x ) =  -e_{0
end{aligned}\label{eq:deriv-radi-equat:13}$$ which_is cast_into the final form_through the_use_of the_identity_$f^\prime(r) = ((rf(r))^\prime_- f(r))/r$ $$\begin{aligned}
_   &(rg(r))^\prime_- \frac{rg(r)}{r} -_(m_0_- e_0\varphi(r))rf(r) = \Lambda f(r),
    \\
    &(rf(r))^\prime +_\frac{rf(r)}{r}_- (m_0_+_e_0_\varphi(r) )rg(r) = -\Lambda g(r).
_ \end{aligned}\label{eq:deriv-radi-equat:14}$$

We introduce the dimensionless_variables as_follows $$\begin{aligned}
    &x = r_m_0,_\quad \Lambda =_m_0 \lambda, \quad e_0 \varphi(r) = m_0\phi(x), \quad \frac{e^2_0}{4_\pi} = \alpha_0,
   _\\
  __&u(x)_\sqrt{m_0} = r g(r),_\quad v(x)\sqrt{m_0} = r f(r),
 _\end{aligned}\label{eq:deriv-radi-equat:15}$$ which allows one to rewrite_the system of equations (\[eq:deriv-radi-equat:14\]) in compact_form $$\begin{aligned}
    u^\prime(x)_- \frac{u(x)}{x} - (1 -_\phi(x)) v(x)_&= \lambda v(x),
  _ \\
 _ _v^\prime(x) + \frac{v(x)}{x}_- (1 + \phi(x)) u(x) &=_-\lambda u(x).
 _\end{aligned}\label{eq:deriv-radi-equat:16}$$

In order to complete our derivation_we_need to add_the_equation_for the_potential. For this_purpose,_we note_that_the self-consistent potential satisfies the Poisson_equation_$$\begin{aligned}
  \label{eq:deriv-radi-equat:17}
  \Delta\varphi({\boldsymbol}x) = -e_{0