frac{ \h_a\cdot \sum_{b \neq a} \h_b u_b}{\h_a^2}$$ and using $\h_{a}^2=1+O(\tfrac{1}{\sqrt{M}})$ we have $$\E = \frac{1}{2\sigma^2} ( u_a^2 -2 u_a\eta_a)+ U(u_a+x_{0a}) + \E_{\backslash a}(\uba)
\label{eq:Eba}$$ with $\uba$, $\E_{\backslash a}$ etc defined as in Sec. \[sec:result\]. Equation  indicates that the variable $u_a$ interacts with all the others only through $\eta_a$. Therefore, we rewrite the marginal distribution $P(u_a)$ as an integral over the joint distribution of $\eta_a$ and $u_a$, $P(u_a,\eta_a)$. $$P(u_a) = \frac{1}{Z} \int d\uba\, \mathrm{e}^{-\beta \E} =  \int d\eta_a \,P(u_a,\eta_a)
\label{eq:Pmarginal}$$ where $$P(u_a,\eta_a)=\frac{1}{Z}\int d \uba\, \delta(\eta_a+ \h_a\cdot \sum_{b \neq a} \h_b u_b)\,\mathrm{e}^{-\beta \E}\,
\label{eq:Pjoint}$$ for all $a=1,\ldots, N$. Now we introduce a cavity ‘field’ distribution of $\eta_a$ at the removed node $a$ as $$P_{\backslash a}(\eta_a)=\frac{1}{Z_{\backslash a}}\int d \u \;\delta(\eta_a+ \h_a\cdot \sum_{b \neq a} \h_b u_b)\mathrm{e}^{-\beta \E_{\backslash a}}.
\label{eq:Pcavity}$$ By comparing and, we get $$P(u_a) = \tfrac{\int d \eta_a \exp\big[-\beta\big\{\frac{( u_a^2 -2 u_a\eta_a)}{2\sigma^2} + U(u_a+x_{0a})\big\}
frac { \h_a\cdot \sum_{b \neq a } \h_b u_b}{\h_a^2}$$ and using $ \h_{a}^2=1+O(\tfrac{1}{\sqrt{M}})$ we have $ $ \E = \frac{1}{2\sigma^2 } (u_a^2 -2 u_a\eta_a)+ U(u_a+x_{0a }) + \E_{\backslash a}(\uba) 
 \label{eq: Eba}$$ with $ \uba$, $ \E_{\backslash a}$ etc defined as in Sec.   \[sec: result\ ]. Equation    indicate that the varying $ u_a$ interacts with all the others only through $ \eta_a$. Therefore, we rewrite the bare distribution $ P(u_a)$ as an integral over the joint distribution of $ \eta_a$ and $ u_a$, $ P(u_a,\eta_a)$. $ $ P(u_a) = \frac{1}{Z } \int d\uba\, \mathrm{e}^{-\beta \E } =   \int d\eta_a \,P(u_a,\eta_a) 
 \label{eq: Pmarginal}$$ where $ $ P(u_a,\eta_a)=\frac{1}{Z}\int five hundred \uba\, \delta(\eta_a+ \h_a\cdot \sum_{b \neq a } \h_b u_b)\,\mathrm{e}^{-\beta \E}\, 
 \label{eq: Pjoint}$$ for all $ a=1,\ldots, N$. Now we introduce a pit ‘ field ’ distribution of $ \eta_a$ at the removed node $ a$ as $ $ P_{\backslash a}(\eta_a)=\frac{1}{Z_{\backslash a}}\int vitamin d \u \;\delta(\eta_a+ \h_a\cdot \sum_{b \neq a } \h_b u_b)\mathrm{e}^{-\beta \E_{\backslash a } }. 
 \label{eq: Pcavity}$$ By comparing and, we get $ $ P(u_a) = \tfrac{\int d \eta_a \exp\big[-\beta\big\{\frac { (u_a^2 -2 u_a\eta_a)}{2\sigma^2 } + U(u_a+x_{0a})\big\ }
fraf{ \h_a\cdot \sum_{b \neq a} \h_b u_n}{\h_a^2}$$ and using $\h_{a}^2=1+O(\tfrac{1}{\sxrt{M}})$ ws have $$\E = \frac{1}{2\sigma^2} ( u_a^2 -2 u_a\eta_a)+ U(u_a+x_{0e}) + \E_{\vacksoash a}(\uba)
\label{eq:Eba}$$ wigh $\uba$, $\E_{\bwckslash a}$ euc defined as in Sec. \[sec:result\]. Esmatiou  mndicates that jhe variable $u_a$ interacts figh all the others only through $\eta_a$. Trereforr, ae rewrite the marbynal distribution $P(u_a)$ as an integral ober the joint distrinution of $\eta_a$ and $u_a$, $P(u_a,\eha_a)$. $$O(u_a) = \frac{1}{Z} \int d\ubw\, \mathrm{e}^{-\bejz \E} =  \int d\eta_a \,O(u_a,\eta_a)
\labtl{zq:Pmarginal}$$ where $$P(u_a,\eta_a)=\frac{1}{Z}\int d \uba\, \deuta(\etc_a+ \h_a\cdot \som_{g \jgq a} \h_b u_b)\,\mavhrm{e}^{-\bvta \E}\,
\label{eq:Pmpint}$$ fmr all $s=1,\ldots, N$. Now wc intcoduxe a cavity ‘field’ disvribution of $\eta_a$ at the remoeeb node $a$ as $$P_{\backslasy q}(\eta_a)=\xrac{1}{S_{\bacywlarh z}}\iit s \u \;\depta(\xta_a+ \h_a\cdot \sum_{b \neq a} \h_b u_b)\mathrm{e}^{-\beta \E_{\nasjslash a}}.
\label{sq:Pcavytr}$$ By comparing and, we get $$P(u_a) = \tfrac{\int d \sta_a \exp\big[-\beta\big\{\frac{( y_a^2 -2 u_a\eta_a)}{2\sigma^2} + U(u_a+x_{0w})\big\}
frac{ \h_a\cdot \sum_{b \neq a} \h_b u_b}{\h_a^2}$$ $\h_{a}^2=1+O(\tfrac{1}{\sqrt{M}})$ have $$\E \frac{1}{2\sigma^2} ( u_a^2 a}(\uba) with $\uba$, $\E_{\backslash etc defined as Sec. \[sec:result\]. Equation indicates that the $u_a$ interacts with all the others only through $\eta_a$. Therefore, we rewrite the distribution $P(u_a)$ as an integral over the joint distribution of $\eta_a$ and $u_a$, $$P(u_a) \frac{1}{Z} d\uba\, \E} = \int d\eta_a \,P(u_a,\eta_a) \label{eq:Pmarginal}$$ where $$P(u_a,\eta_a)=\frac{1}{Z}\int d \uba\, \delta(\eta_a+ \h_a\cdot \sum_{b \neq a} \h_b \E}\, \label{eq:Pjoint}$$ for all $a=1,\ldots, N$. Now we a cavity ‘field’ distribution $\eta_a$ at the removed node as a}(\eta_a)=\frac{1}{Z_{\backslash a}}\int \u \h_a\cdot \neq a} \h_b \E_{\backslash a}}. \label{eq:Pcavity}$$ By comparing and, we get $$P(u_a) = \tfrac{\int d \eta_a \exp\big[-\beta\big\{\frac{( u_a^2 -2 u_a\eta_a)}{2\sigma^2} U(u_a+x_{0a})\big\}
frac{ \h_a\cdot \sum_{b \neq a} \h_b u_b}{\h_a^2}$$ aNd using $\h_{a}^2=1+O(\Tfrac{1}{\SqrT{M}})$ wE hAve $$\E = \Frac{1}{2\Sigma^2} ( u_a^2 -2 u_a\eta_a)+ U(U_A+x_{0a}) + \E_{\Backslash a}(\uba)
\label{eq:EbA}$$ with $\UbA$, $\e_{\bacKSlAsh a}$ eTc definED aS IN SeC. \[sEc:ResUlT\]. eqUatioN  inDicates That the varIabLe $U_a$ interacts wITh All the otheRs oNly through $\etA_a$. THerefoRe, We rEWrite The MargiNal disTRibutiOn $P(u_a)$ as an InTEgral oVEr the joINT dIstrIbution of $\eta_a$ and $u_A$, $p(u_A,\Eta_a)$. $$P(u_a) = \frac{1}{Z} \inT d\uba\, \mAtHRm{E}^{-\BEta \e} =  \inT d\eta_a \,P(u_a,\eTa_A)
\labeL{Eq:PmargINaL}$$ WHEre $$p(U_a,\eta_a)=\frac{1}{Z}\inT d \uba\, \delta(\eTA_a+ \h_A\cdot \sUm_{B \neQ A} \h_b u_b)\,\mAthrm{E}^{-\bETa \E}\,
\Label{eq:PjoiNt}$$ foR all $a=1,\ldotS, N$. Now wE IntroduCE a cavitY ‘field’ DisTriButiON oF $\eTa_a$ At THe rEMoVed NOde $A$ as $$P_{\backSlAsH a}(\eta_A)=\fraC{1}{z_{\BACkslAsh A}}\int D \u \;\delTa(\eta_a+ \h_a\cdot \sUm_{b \Neq a} \H_B u_b)\MathrM{e}^{-\betA \E_{\baCkSlash A}}.
\label{Eq:PcaViTy}$$ By comparing anD, we gEt $$P(u_a) = \tfraC{\inT d \Eta_A \eXp\big[-\BEta\big\{\FraC{( u_a^2 -2 U_a\eta_a)}{2\sIgma^2} + U(u_a+X_{0A})\biG\}
frac{ \h_a\cdot \sum_{b \n eq a} \h_b  u_b} {\h _a^ 2} $$ a nd u sing $\h_{a}^2 = 1+O( \tfrac{1}{\sqrt{M}})$we ha ve $$\E =\frac {1}{2\s i gm a ^ 2}(u_ a^2  - 2  u _a\et a_a )+ U(u_ a+x_{0a})+ \ E_ {\backslasha }( \uba)
\lab el{ eq:Eba}$$ wi th$\uba$ ,$\E _ {\bac ksl ash a }$ etc define d as in S ec .  \[sec : result\ ] .  E quat ion  indicates th a tt he variable $u _a$ in te r ac t s  wi thall the ot he rs on l y throu g h$ \ e ta_ a $. Therefore,  we rewrite the  margi na l d i stribu tion$P ( u_a )$ as an in tegr al over t he joi n t distr i butionof $\e ta_ a$and$ u_ a$ , $ P( u _a, \ et a_a ) $.$$P(u_a)  =  \ frac{ 1}{Z } \ i nt d \ub a\,\math rm{e}^{-\beta  \E } =\in t d\e ta_a\,P( u_ a,\et a_a)
\ label {e q:Pmarginal}$$wher e $$P(u_a ,\e ta _a) =\ frac{ 1 }{Z}\i ntd \ uba\, \ delta(\ e ta_ a+ \ h _a \cdot \sum_{b \neq  a } \h _b u_b)\ ,\math r m{ e} ^ {-\beta\E }\, 
\la b e l{eq: Pjoi n t} $$ for a ll $a= 1 ,\ ld ots, N$ .Now we  i ntr odu ce ac avit y ‘fie ld’ dist ribut i on of $\eta_a$ at the remove d  n o d e$ a$ a s $ $P_{\backsl asha }(\e ta_a ) =\ fra c {1}{Z _{\ba ck s la s h a}}\int d \u \;\d el ta(\et a_a+\h_a\cdot \su m_{b \neqa } \h_b u_b )\ma t hr m {e}^{-\beta \E _{\ba ckslash a} } .
\label {eq:P cavity}$ $ By comp a r ing and,  we  ge t $ $P( u _ a)  = \tfrac{\in t d \e ta _a \exp \bi g[-\bet a\b ig\ {\f rac {(  u_a^2 -2  u_a\eta _a )} {2 \s igm a^2}+  U(u_a+x _{ 0a} )\ big \}
frac{ \h_a\cdot_\sum_{b \neq_a} \h_b u_b}{\h_a^2}$$ and_using $\h_{a}^2=1+O(\tfrac{1}{\sqrt{M}})$_we_have $$\E_=_\frac{1}{2\sigma^2} ( u_a^2_-2 u_a\eta_a)+ U(u_a+x_{0a})_+ \E_{\backslash a}(\uba)
\label{eq:Eba}$$ with_$\uba$, $\E_{\backslash a}$_etc_defined as in Sec. \[sec:result\]. Equation  indicates that the variable $u_a$ interacts with all the_others_only through_$\eta_a$._Therefore,_we rewrite the marginal distribution_$P(u_a)$ as an integral over_the joint_distribution of $\eta_a$ and $u_a$, $P(u_a,\eta_a)$. $$P(u_a) =_\frac{1}{Z}_\int d\uba\, \mathrm{e}^{-\beta_\E} =  \int d\eta_a \,P(u_a,\eta_a)
\label{eq:Pmarginal}$$ where $$P(u_a,\eta_a)=\frac{1}{Z}\int d_\uba\, \delta(\eta_a+ \h_a\cdot \sum_{b \neq a}_\h_b u_b)\,\mathrm{e}^{-\beta \E}\,
\label{eq:Pjoint}$$_for_all_$a=1,\ldots, N$. Now we_introduce a cavity ‘field’ distribution of_$\eta_a$ at the removed node $a$_as $$P_{\backslash a}(\eta_a)=\frac{1}{Z_{\backslash a}}\int d \u \;\delta(\eta_a+_\h_a\cdot \sum_{b \neq a} \h_b u_b)\mathrm{e}^{-\beta_\E_{\backslash a}}.
\label{eq:Pcavity}$$ By comparing and,_we get_$$P(u_a) = \tfrac{\int d \eta_a_\exp\big[-\beta\big\{\frac{( u_a^2 -2_u_a\eta_a)}{2\sigma^2} +_U(u_a+x_{0a})\big\}