 Then $$\label{eqn:JSbisection}
    W\left(\epsilon_{\rm S*}-\epsilon\right)=\frac{p_{\rm *}^{2}-p^{2}}{2},$$ where $\epsilon_{\rm S*}=\epsilon_{\rm *}(W)$ is implicitly a function of $W$ (via Equations \[eqn:eoscopy\]-\[eqn:JScopy\]). The right hand side of Equation \[eqn:JSbisection\] is a positive constant in the iteration of $J_{\rm S}$. The derivative of the left hand side is $$\frac{d\text{LHS}}{dW}=\epsilon_{\rm *}(W)-\epsilon+W\left(\frac{\partial\epsilon_{\rm S*}}{\partial\rho_{\rm S*}}\right)_{p_{\rm *}}\frac{d\rho_{\rm S*}}{dW},$$ from Equation \[eqn:rho\*copy\] $$\begin{aligned}
    \frac{d\rho_{\rm S*}}{dW}=-\rho_{\rm S*}^2\frac{p_{\rm *}-p}{W^2}<0,\end{aligned}$$ and from Equation \[eqn:unique1\] and \[eqn:unique2\], $$\begin{aligned}
    \left(\frac{\partial\epsilon}{\partial\rho}\right)_{p}=-\left(\frac{\partial\epsilon}{\partial p}\right)_{\rho}\bigg/\left(\frac{\partial\rho}{\partial p}\right)_{\epsilon}=\epsilon_{\rho}\equiv\epsilon_{\rho_{\rm S*}}<0.\end{aligned}$$ Thus, $d\text{LHS}/dW\ge0$ and only equals 0 when $W=J_{\rm S}^{2}=0$, $\text{LHS}=0$ when $W=J_{\rm S}^{2}=0$, so there is exactly one solution of $W$ to Equation \[eqn:JSbisection\]. In Section \[sec:JS\], we have chosen the $J_{\rm S,max}$ to be large enough to ensure the solution $J_{\rm S}\le J_{\rm S,max}$.

To apply the bisection root finding algorithm to Equation \[eqn:JSbisection\], we pick a $J_{\rm S}$ in the range of $[J_{\rm S,min},J_{\rm S,max}]$. We then calculate $\rho_{\rm *}$ using Equation \[eqn:rho
Then $ $ \label{eqn: JSbisection } 
     W\left(\epsilon_{\rm S*}-\epsilon\right)=\frac{p_{\rm * } ^{2}-p^{2}}{2},$$ where $ \epsilon_{\rm S*}=\epsilon_{\rm * } (W)$ is implicitly a function of $ W$ (via Equations   \[eqn: eoscopy\]-\[eqn: JScopy\ ]). The right hand side of Equation \[eqn: JSbisection\ ] is a incontrovertible constant in the iteration of $ J_{\rm S}$. The derived function of the left hand side is $ $ \frac{d\text{LHS}}{dW}=\epsilon_{\rm * } (W)-\epsilon+W\left(\frac{\partial\epsilon_{\rm S*}}{\partial\rho_{\rm S*}}\right)_{p_{\rm * } } \frac{d\rho_{\rm S*}}{dW},$$ from Equation \[eqn: rho\*copy\ ] $ $ \begin{aligned } 
     \frac{d\rho_{\rm S*}}{dW}=-\rho_{\rm S*}^2\frac{p_{\rm * } -p}{W^2}<0,\end{aligned}$$ and from Equation \[eqn: unique1\ ] and \[eqn: unique2\ ], $ $ \begin{aligned } 
     \left(\frac{\partial\epsilon}{\partial\rho}\right)_{p}=-\left(\frac{\partial\epsilon}{\partial p}\right)_{\rho}\bigg/\left(\frac{\partial\rho}{\partial p}\right)_{\epsilon}=\epsilon_{\rho}\equiv\epsilon_{\rho_{\rm S*}}<0.\end{aligned}$$ therefore, $ d\text{LHS}/dW\ge0 $ and only equals 0 when $ W = J_{\rm S}^{2}=0 $, $ \text{LHS}=0 $ when $ W = J_{\rm S}^{2}=0 $, so there is precisely one solution of $ W$ to Equation \[eqn: JSbisection\ ]. In Section \[sec: JS\ ], we have chosen the $ J_{\rm S, max}$ to be large enough to guarantee the solution $ J_{\rm S}\le J_{\rm S, max}$. 

 To apply the bisection root receive algorithm to Equation \[eqn: JSbisection\ ], we pick a $ J_{\rm S}$ in the range of $ [ J_{\rm S, min},J_{\rm S, max}]$. We then forecast $ \rho_{\rm * } $ use Equation \[eqn: rho
 Thfn $$\label{eqn:JSbisection}
    W\left(\epsilon_{\rm S*}-\epsinon\riggt)=\frac{p_{\ro *}^{2}-p^{2}}{2},$$ where $\epsilon_{\rm S*}=\epsilon_{\cm *}(W)$ is inplicitly a function ow $W$ (via Eeuations \[wqn:eiwcopy\]-\[eqn:JSrkpy\]). The right mand vmde of Equation \[eqn:JSbisewtion\] is a poshtkvz constant in the iteration of $J_{\rm S}$. The detigative of the jeft rand side is $$\frac{d\text{LHS}}{dW}=\epsilon_{\rm *}(W)-\spsilon+X\left(\frac{\partiak\epsilon_{\rm S*}}{\partial\rho_{\rm S*}}\gighh)_{p_{\rm *}}\frac{d\rho_{\rm S*}}{dA},$$ from Equajjon \[wqn:rho\*copy\] $$\bdgin{aligned}
    \frac{d\rhk_{\rm S*}}{dW}=-\rho_{\rm S*}^2\frac{p_{\rm *}-p}{W^2}<0,\end{aligved}$$ aud from Equqtuon \[gqn:unique1\] anv \[eqn:ugique2\], $$\begin{aligned}
    \left(\ftac{\partial\epsikon}{\'artual\rho}\right)_{p}=-\left(\frac{\pertial\epsilon}{\partial p}\right)_{\rhm}\bngg/\left(\frac{\partial\rho}{\paetial p}\richt)_{\eowilun}=\eksikoh_{\rho}\eqkiv\xpsilon_{\rho_{\rj S*}}<0.\end{alignwd}$$ Thus, $d\text{LHS}/dW\gt0$ age only equals 0 when $W=T_{\rm S}^{2}=0$, $\text{LHS}=0$ when $W=J_{\rm S}^{2}=0$, so there is eqactmy one solution of $W$ to Equation \[eqn:JSbisectlon\]. In Sestion \[sec:JS\], we have chosen the $J_{\rm S,max}$ to be larce enkjgh tj wndure the solution $J_{\rm S}\le J_{\rm S,max}$.

To apply thq bosvction root findikg algorithm to Eqiahipg \[eqn:JSbisectkon\], we pidk a $J_{\rm S}$ in the gange os $[J_{\rm S,min},J_{\rm F,max}]$. We then calculate $\rho_{\rm *}$ uwing Equatiou \[ewn:rho
Then $$\label{eqn:JSbisection} W\left(\epsilon_{\rm S*}-\epsilon\right)=\frac{p_{\rm *}^{2}-p^{2}}{2},$$ where $\epsilon_{\rm is a function $W$ (via Equations of \[eqn:JSbisection\] is a constant in the of $J_{\rm S}$. The derivative of left hand side is $$\frac{d\text{LHS}}{dW}=\epsilon_{\rm *}(W)-\epsilon+W\left(\frac{\partial\epsilon_{\rm S*}}{\partial\rho_{\rm S*}}\right)_{p_{\rm *}}\frac{d\rho_{\rm S*}}{dW},$$ from Equation \[eqn:rho\*copy\] \frac{d\rho_{\rm S*}}{dW}=-\rho_{\rm S*}^2\frac{p_{\rm *}-p}{W^2}<0,\end{aligned}$$ and from Equation \[eqn:unique1\] and \[eqn:unique2\], $$\begin{aligned} \left(\frac{\partial\epsilon}{\partial\rho}\right)_{p}=-\left(\frac{\partial\epsilon}{\partial p}\right)_{\rho}\bigg/\left(\frac{\partial\rho}{\partial S*}}<0.\end{aligned}$$ $d\text{LHS}/dW\ge0$ only 0 when $W=J_{\rm S}^{2}=0$, $\text{LHS}=0$ when $W=J_{\rm S}^{2}=0$, so there is exactly one solution of $W$ Equation \[eqn:JSbisection\]. In Section \[sec:JS\], we have chosen $J_{\rm S,max}$ to be enough to ensure the solution S}\le S,max}$. To the root algorithm to Equation we pick a $J_{\rm S}$ in the range of $[J_{\rm S,min},J_{\rm S,max}]$. We then calculate $\rho_{\rm *}$ Equation \[eqn:rho
 Then $$\label{eqn:JSbisection}
    W\lEft(\epsilon_{\Rm S*}-\epSilOn\rIgHt)=\frAc{p_{\rM *}^{2}-p^{2}}{2},$$ where $\epsilon_{\RM S*}=\epSilon_{\rm *}(W)$ is implicitly a fUnctiOn OF $W$ (viA eqUatioNs \[eqn:eoSCoPY\]-\[Eqn:jSCoPy\]). THe RIgHt hanD siDe of EquAtion \[eqn:JSBisEcTion\] is a positIVe Constant in The Iteration of $J_{\Rm S}$. the derIvAtiVE of thE leFt hanD side iS $$\Frac{d\tExt{LHS}}{dW}=\ePsILon_{\rm *}(W)-\EPsilon+W\LEFt(\Frac{\Partial\epsilon_{\rm S*}}{\PArTIal\rho_{\rm S*}}\right)_{P_{\rm *}}\fraC{d\RHo_{\RM s*}}{dW},$$ FroM Equation \[eQn:Rho\*coPY\] $$\begin{aLIgNED}
    \FraC{D\rho_{\rm S*}}{dW}=-\rho_{\rM S*}^2\frac{p_{\rm *}-p}{W^2}<0,\ENd{aLigned}$$ AnD frOM EquatIon \[eqN:uNIquE1\] and \[eqn:uniqUe2\], $$\beGin{aligneD}
    \left(\fRAc{\partiAL\epsiloN}{\partiAl\rHo}\rIght)_{P}=-\LeFt(\FraC{\pARtiAL\ePsiLOn}{\pArtial p}\rIgHt)_{\Rho}\biGg/\leFT(\FRAc{\paRtiAl\rhO}{\partIal p}\right)_{\epsiLon}=\EpsiLOn_{\rHo}\equIv\epsIlon_{\RhO_{\rm S*}}<0.\eNd{aligNed}$$ ThUs, $D\text{LHS}/dW\ge0$ and Only Equals 0 wheN $W=J_{\Rm s}^{2}=0$, $\teXt{lHS}=0$ whEN $W=J_{\rm S}^{2}=0$, So tHerE is exacTly one sOLutIoN OF $w$ tO Equation \[eqn:JSbiseCtION\]. IN Section \[Sec:JS\], wE HaVe CHosen the $j_{\rM S,mAx}$ to BE Large EnouGH tO ensure tHe soluTIoN $J_{\Rm S}\le J_{\rM S,Max}$.

To aPpLy tHe bIsectIOn roOt findIng algorIthm tO equation \[eqn:JSbISection\], we pick A $j_{\rM s}$ In THe raNge Of $[J_{\rm S,min},J_{\rM S,maX}]$. we thEn caLCuLatE $\Rho_{\rm *}$ Using eqUAtIOn \[eqn:rho
 Then $$\label{eqn:JSbisec tion}
   W\lef t(\ eps il on_{ \rmS*}-\epsilon\r i ght) =\frac{p_{\rm *}^{2}-p ^{2}} {2 } ,$$w he re $\ epsilon _ {\ r m  S* }= \e psi lo n _{ \rm * }(W )$ is i mplicitlya f un ction of $W$ (v ia Equatio ns\[eqn:eoscop y\] -\[eqn :J Sco p y\]).  Th e rig ht han d  sideof Equati on \[eqn: J Sbisect i o n\ ] is  a positive const a nt in the iterati on of$J _ {\ r m  S} $.The deriva ti ve of the lef t  h a n d  si d e is $$\frac{ d\text{LHS} } {dW }=\eps il on_ { \rm *} (W)-\ ep s ilo n+W\left(\f rac{ \partial\ epsilo n _{\rm S * }}{\par tial\r ho_ {\r m S* } }\ ri ght )_ { p_{ \ rm  *} } \fr ac{d\rho _{ \r m S*} }{dW } , $ $  fro m E quat ion \ [eqn:rho\*cop y\]  $$\ b egi n{ali gned} 
   \ frac{ d\rho_ {\rmS* }}{dW}=-\rho_{\ rm S *}^2\frac {p_ {\ rm*} -p}{W ^ 2}<0,\ end {al igned}$ $ and f r omEq u a t io n \[eqn:unique1\]an d \[ eqn:uniq ue2\], $$ \b e gin{alig ne d}   \ l eft(\ frac { \p artial\e psilon } {\ pa rtial\r ho }\righ t) _{p }=- \left ( \fra c{\par tial\eps ilon} { \partial p}\ri g ht)_{\rho}\bi g g/ \ l ef t (\fr ac{ \partial\rh o}{\ p arti al p } \r igh t )_{\e psilo n} = \e p silon_{\rho}\equiv\ ep silon_ {\rho _{\rm S*}}<0. \end{align e d } $$ Thus,  $d\ t ex t {LHS}/dW\ge0$and o nly equals 0 when $ W=J_{ \rm S}^{ 2}=0$, $\ t e xt{LHS}= 0$whe n $ W=J _ { \r m S}^{2}=0$,s o  the re  is exa ctl y one s olu tio n o f $ W$  to Equat ion \[eq n: JS bi se cti on\]. In Secti on  \[ se c:J S\],w e have  chos en t he  $ J _{\ rm S,ma x }$ t o be  l ar ge e nou gh  to e nsur e  th e solut ion $J_{\ rmS }\le  J _{ \rm S,m ax}$.

To app ly  the bisec ti onroot f i n ding alg orithm to Equation \[eq n :JSbise cti on\],  wepick a $J _{\ rm S}$  in the ra nge of  $[J_ {\ rmS , min}, J _ {\ rmS, max}]$. We t hen  calc ul ate$\rho_{ \rm *}$ using Equa t ion  \[eqn:rho
 Then_$$\label{eqn:JSbisection}
 _  W\left(\epsilon_{\rm S*}-\epsilon\right)=\frac{p_{\rm_*}^{2}-p^{2}}{2},$$ where_$\epsilon_{\rm_S*}=\epsilon_{\rm *}(W)$_is_implicitly a function_of $W$ (via_Equations \[eqn:eoscopy\]-\[eqn:JScopy\]). The right hand_side of Equation_\[eqn:JSbisection\]_is a positive constant in the iteration of $J_{\rm S}$. The derivative of the_left_hand side_is_$$\frac{d\text{LHS}}{dW}=\epsilon_{\rm_*}(W)-\epsilon+W\left(\frac{\partial\epsilon_{\rm S*}}{\partial\rho_{\rm S*}}\right)_{p_{\rm *}}\frac{d\rho_{\rm S*}}{dW},$$_from Equation \[eqn:rho\*copy\] $$\begin{aligned}
 _ _\frac{d\rho_{\rm S*}}{dW}=-\rho_{\rm S*}^2\frac{p_{\rm *}-p}{W^2}<0,\end{aligned}$$ and from Equation \[eqn:unique1\]_and_\[eqn:unique2\], $$\begin{aligned}
 _  \left(\frac{\partial\epsilon}{\partial\rho}\right)_{p}=-\left(\frac{\partial\epsilon}{\partial p}\right)_{\rho}\bigg/\left(\frac{\partial\rho}{\partial p}\right)_{\epsilon}=\epsilon_{\rho}\equiv\epsilon_{\rho_{\rm S*}}<0.\end{aligned}$$ Thus, $d\text{LHS}/dW\ge0$ and_only equals 0 when $W=J_{\rm S}^{2}=0$,_$\text{LHS}=0$ when $W=J_{\rm_S}^{2}=0$,_so_there is exactly one_solution of $W$ to Equation \[eqn:JSbisection\]._In Section \[sec:JS\], we have chosen_the $J_{\rm S,max}$ to be large enough_to ensure the solution $J_{\rm S}\le_J_{\rm S,max}$.

To apply the bisection_root finding_algorithm to Equation \[eqn:JSbisection\], we_pick a $J_{\rm_S}$ in_the range of_$[J_{\rm S,min},J_{\rm S,max}]$. We then calculate_$\rho_{\rm *}$ using_Equation \[eqn:rho