2}}}\,({{\rho_{B}}})^{2}
-{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{d,s}{\frac{|Q_{f}|B}{2\pi^{2}}}  \sum_{n=0}^{n^{f}_{max}} 3(2-\delta_{n0})
\int_{0}^{k^{f}_{z,F}} dk_{z} \,{\frac{{k_{z}}^{2}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}
\label{parallelpressuremagon}$$ and $$p_{\perp}={\frac{27{g_{h}}^{2}}{16{m_{G}}^{2}}}\,({{\rho_{B}}})^{2}
-{\mathcal{B}}_{QCD}+{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{d,s}{\frac{|Q_{f}|^{2}B^{2}}{2\pi^{2}}}  \sum_{n=0}^{n^{f}_{max}} 3(2-\delta_{n0})
n \int_{0}^{k^{f}_{z,F}}   {\frac{dk_{z}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}  
\label{perppressuremagon}$$ The baryon density ($\rho_{B}$) is given by [@we16; @we17]: $$\rho_{B}=\sum_{f=u}^{d,s}\,{\frac{|Q_{f}|B}{2\pi^{2}}} \,
\sum_{n=0}^{n^{f}_{max}}(2-\delta_{n0}) \, \sqrt{{\nu_{f}}^{2}-m_{f}^{2}-2n|Q_{f}|B}
\hspace{0.7cm} \textrm{with} \hspace{0.5cm}
n\leq n^{f}_{max}=int\Bigg[ {\frac{{\nu_{f}}^{2}-m_{f}^{2}}{2|Q_{f}|B}} \Bigg]
\label{magbaryondens}$$ where $ {\it{int}}[a]$ denotes the integer part of $a$ and ${\nu_{f}}$ is the chemical potential for the quark $f$. As in [@we16] we
2}}}\,({{\rho_{B}}})^{2 } 
 -{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi } } 
 + \sum_{f = u}^{d, s}{\frac{|Q_{f}|B}{2\pi^{2 } } }   \sum_{n=0}^{n^{f}_{max } } 3(2-\delta_{n0 }) 
 \int_{0}^{k^{f}_{z, F } } dk_{z } \,{\frac{{k_{z}}^{2}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B } } } 
 \label{parallelpressuremagon}$$ and $ $ p_{\perp}={\frac{27{g_{h}}^{2}}{16{m_{G}}^{2}}}\,({{\rho_{B}}})^{2 } 
 -{\mathcal{B}}_{QCD}+{\frac{B^{2}}{8\pi } } 
 + \sum_{f = u}^{d, s}{\frac{|Q_{f}|^{2}B^{2}}{2\pi^{2 } } }   \sum_{n=0}^{n^{f}_{max } } 3(2-\delta_{n0 }) 
 n \int_{0}^{k^{f}_{z, F } }    { \frac{dk_{z}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B } } }  
 \label{perppressuremagon}$$ The baryon density ($ \rho_{B}$) is given by [ @we16; @we17 ]: $ $ \rho_{B}=\sum_{f = u}^{d, s}\,{\frac{|Q_{f}|B}{2\pi^{2 } } } \, 
 \sum_{n=0}^{n^{f}_{max}}(2-\delta_{n0 }) \, \sqrt{{\nu_{f}}^{2}-m_{f}^{2}-2n|Q_{f}|B } 
 \hspace{0.7 cm } \textrm{with } \hspace{0.5 cm } 
 n\leq n^{f}_{max}=int\Bigg [ { \frac{{\nu_{f}}^{2}-m_{f}^{2}}{2|Q_{f}|B } } \Bigg ] 
 \label{magbaryondens}$$ where $ { \it{int}}[a]$ denotes the integer character of $ a$ and $ { \nu_{f}}$ is the chemical electric potential for the quark $ f$. As in [ @we16 ] we
2}}}\,({{\rho_{H}}})^{2}
-{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{a,s}{\frac{|Q_{f}|B}{2\pi^{2}}}  \sum_{n=0}^{n^{f}_{max}} 3(2-\delta_{h0})
\int_{0}^{k^{f}_{z,F}} dk_{z} \,{\frac{{k_{z}}^{2}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}
\label{'araolelpeessuremagon}$$ and $$p_{\perp}={\wrac{27{g_{h}}^{2}}{16{m_{G}}^{2}}}\,({{\rjo_{B}}})^{2}
-{\mathcql{B}}_{QRD}+{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{d,s}{\hdac{|Q_{f}|^{2}B^{2}}{2\pl^{2}}}  \suj_{k=0}^{n^{f}_{mar}} 3(2-\velta_{n0})
n \int_{0}^{k^{f}_{z,F}}   {\frac{dk_{z}}{\vqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}  
\nacep{perppressuremagon}$$ The baryon densiey ($\rho_{B}$) id given by [@we16; @re17]: $$\rnj_{B}=\suj_{f=u}^{d,s}\,{\frac{|Q_{f}|B}{2\pi^{2}}} \,
\sum_{n=0}^{n^{f}_{max}}(2-\delta_{n0}) \, \sqdt{{\nu_{f}}^{2}-m_{f}^{2}-2i|Q_{f}|B}
\hspace{0.7cm} \tectrm{with} \hspace{0.5cm}
n\leq n^{f}_{mad}=int\Higg[ {\frac{{\nu_{f}}^{2}-m_{f}^{2}}{2|Q_{f}|B}} \Higg]
\label{matbarrindens}$$ where $ {\it{int}}[a]$ dtnmtes the ihteger part of $a$ and ${\nu_{f}}$ is the chemncal potentuao flt the quark $h$. As ig [@we16] we
2}}}\,({{\rho_{B}}})^{2} -{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi}} +\sum_{f=u}^{d,s}{\frac{|Q_{f}|B}{2\pi^{2}}} \sum_{n=0}^{n^{f}_{max}} 3(2-\delta_{n0}) \int_{0}^{k^{f}_{z,F}} dk_{z} and -{\mathcal{B}}_{QCD}+{\frac{B^{2}}{8\pi}} +\sum_{f=u}^{d,s}{\frac{|Q_{f}|^{2}B^{2}}{2\pi^{2}}} 3(2-\delta_{n0}) n \int_{0}^{k^{f}_{z,F}} ($\rho_{B}$) given by [@we16; $$\rho_{B}=\sum_{f=u}^{d,s}\,{\frac{|Q_{f}|B}{2\pi^{2}}} \, \sum_{n=0}^{n^{f}_{max}}(2-\delta_{n0}) \sqrt{{\nu_{f}}^{2}-m_{f}^{2}-2n|Q_{f}|B} \hspace{0.7cm} \textrm{with} \hspace{0.5cm} n\leq n^{f}_{max}=int\Bigg[ \Bigg] \label{magbaryondens}$$ where $ {\it{int}}[a]$ denotes the integer part of $a$ and ${\nu_{f}}$ the chemical potential for the quark $f$. As in [@we16] we
2}}}\,({{\rho_{B}}})^{2}
-{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi}}
+\sum_{F=u}^{d,s}{\frac{|Q_{f}|b}{2\pi^{2}}}  \suM_{n=0}^{n^{F}_{maX}} 3(2-\dElta_{N0})
\int_{0}^{K^{f}_{z,F}} dk_{z} \,{\frac{{k_{z}}^{2}}{\sQRt{m_{f}^{2}+K_{z}^{2}+2n|Q_{f}|B}}}
\label{parallelpreSsureMaGOn}$$ anD $$P_{\pErp}={\frAc{27{g_{h}}^{2}}{16{m_{G}}^{2}}}\,({{\rHO_{B}}})^{2}
-{\MAThcAl{b}}_{QcD}+{\fRaC{b^{2}}{8\pI}}
+\sum_{f=U}^{d,s}{\Frac{|Q_{f}|^{2}B^{2}}{2\Pi^{2}}}  \sum_{n=0}^{n^{f}_{maX}} 3(2-\deLtA_{n0})
n \int_{0}^{k^{f}_{z,F}}   {\frAC{dK_{z}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2N|Q_{f}|b}}}  
\label{perpprEssUremagOn}$$ the BAryon DenSity ($\rHo_{B}$) is gIVen by [@wE16; @we17]: $$\rho_{B}=\suM_{f=U}^{D,s}\,{\frac{|q_{F}|B}{2\pi^{2}}} \,
\sum_{N=0}^{N^{F}_{mAx}}(2-\deLta_{n0}) \, \sqrt{{\nu_{f}}^{2}-m_{f}^{2}-2n|Q_{f}|B}
\HSpACe{0.7cm} \textrm{with} \Hspace{0.5Cm}
N\LeQ N^{F}_{maX}=inT\Bigg[ {\frac{{\nU_{f}}^{2}-M_{f}^{2}}{2|Q_{f}|B}} \bIgg]
\labeL{MaGBARyoNDens}$$ where $ {\it{inT}}[a]$ denotes thE IntEger paRt Of $a$ ANd ${\nu_{f}}$ iS the cHeMIcaL potential fOr thE quark $f$. As In [@we16] we
2}}}\,({{\rho_{B}}})^{2}
- {\mathcal{ B}}_{ QCD }-{ \f rac{ B^{2 }}{8\pi}}
+\su m _{f= u}^{d,s}{\frac{|Q_{f}| B}{2\ pi ^ {2}} }  \sum_ {n=0}^{ n ^{ f } _{m ax }}  3( 2- \ de lta_{ n0} )
\int_ {0}^{k^{f} _{z ,F }} dk_{z} \, { \f rac{{k_{z} }^{ 2}}{\sqrt{m_ {f} ^{2}+k _{ z}^ { 2}+2n |Q_ {f}|B }}}
\l a bel{pa rallelpre ss u remago n }$$ and $ $p _{\p erp}={\frac{27{g_ { h} } ^{2}}{16{m_{G} }^{2}} }\ , ({ { \ rho _{B }}})^{2}
- {\ mathc a l{B}}_{ Q CD } + { \fr a c{B^{2}}{8\pi }}
+\sum_{f = u}^ {d,s}{ \f rac { |Q_{f} |^{2} B^ { 2}} {2\pi^{2}}}   \s um_{n=0}^ {n^{f} _ {max}}3 (2-\del ta_{n0 })n \ int_ { 0} ^{ k^{ f} _ {z, F }}   { \fr ac{dk_{z }} {\ sqrt{ m_{f } ^ { 2 }+k_ {z} ^{2} +2n|Q _{f}|B}}}  
\ lab el{p e rpp ressu remag on}$ $The b aryondensi ty  ($\rho_{B}$) i s gi ven by [@ we1 6;  @w e1 7]: $ $ \rho_{ B}= \su m_{f=u} ^{d,s}\ , {\f ra c { | Q_ {f}|B}{2\pi^{2}}}\, \ su m_{n=0}^ {n^{f} _ {m ax } }(2-\del ta _{n 0})\ ,  \sqr t{{\ n u_ {f}}^{2} -m_{f} ^ {2 }- 2n|Q_{f }| B}
\hs pa ce{ 0.7 cm} \ t extr m{with } \hspac e{0.5 c m}
n\leq n^{f} _ {max}=int\Big g [{ \ fr a c{{\ nu_ {f}}^{2}-m_ {f}^ { 2}}{ 2|Q_ { f} |B} }  \Big g]
\l ab e l{ m agbaryondens}$$ whe re  $ {\i t{int }}[a]$ denote s the inte g e r  part of  $a$ an d  ${\nu_{f}}$ i s the  chemicalp otential  forthe quar k $f$. As i n [@we16 ] w e
2}}}\,({{\rho_{B}}})^{2}
-{\mathcal{B}}_{QCD}-{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{d,s}{\frac{|Q_{f}|B}{2\pi^{2}}} _\sum_{n=0}^{n^{f}_{max}} 3(2-\delta_{n0})
\int_{0}^{k^{f}_{z,F}}_dk_{z} \,{\frac{{k_{z}}^{2}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}
\label{parallelpressuremagon}$$ and $$p_{\perp}={\frac{27{g_{h}}^{2}}{16{m_{G}}^{2}}}\,({{\rho_{B}}})^{2}
-{\mathcal{B}}_{QCD}+{\frac{B^{2}}{8\pi}}
+\sum_{f=u}^{d,s}{\frac{|Q_{f}|^{2}B^{2}}{2\pi^{2}}}_ \sum_{n=0}^{n^{f}_{max}}_3(2-\delta_{n0})
n_\int_{0}^{k^{f}_{z,F}} __{\frac{dk_{z}}{\sqrt{m_{f}^{2}+k_{z}^{2}+2n|Q_{f}|B}}}  
\label{perppressuremagon}$$_The baryon density_($\rho_{B}$) is given by_[@we16; @we17]: $$\rho_{B}=\sum_{f=u}^{d,s}\,{\frac{|Q_{f}|B}{2\pi^{2}}}_\,
\sum_{n=0}^{n^{f}_{max}}(2-\delta_{n0})_\, \sqrt{{\nu_{f}}^{2}-m_{f}^{2}-2n|Q_{f}|B}
\hspace{0.7cm} \textrm{with} \hspace{0.5cm}
n\leq n^{f}_{max}=int\Bigg[ {\frac{{\nu_{f}}^{2}-m_{f}^{2}}{2|Q_{f}|B}} \Bigg]
\label{magbaryondens}$$ where $ {\it{int}}[a]$ denotes the integer part_of_$a$ and_${\nu_{f}}$_is_the chemical potential for the_quark $f$. As in [@we16]_we