&\leq&C_R \int_{\substack{|y-x_r|>r\\|y|\leq R, \, y_n>0}}\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}\frac{1}{|x-y|^{n+1}}dy.\label{b6}\end{aligned}$$

Here we have used the fact that the $\alpha$-harmonic function $u$ is bounded in the region $$D_{R,r}=\left\{y=(y',y_n)\left|\right.|y-x_r|>r,~|y|<R,~y_n>0\right\}.$$ The bound depends on $R$, however is independent of $r$, since $D_{R, r_1} \subset D_{R, r_2}$ for $r_1>r_2$. For each such fixed open domain $D_{R,r}$, the bound of the $\alpha$-harmonic function $u$ can be derived from the interior smoothness ( see, for instance [@BKN] and [@FW]) and the estimate up to the boundary ( see [@RS]).

Set  $y=(y',y_n), \sigma=|y'|$, for fixed  $x$  and sufficiently large  $r$, we have

$$\begin{aligned}
\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}&\leq&\frac{Cr^{\alpha/2}}{|2y_nr-|y|^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|\sigma^2-2y_nr+y_n^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|(y_n-r)^2+\sigma^2-r^2|^{\alpha/2}}.\label{b7}\end{aligned}$$

and $$\frac{1}{|x-y|^{n+1}}\leq\frac{C}{(1+|y|)^{n+1}}\leq\frac{C}{(1+|y'|)^{n+1}}=\frac{C}{(1+\sigma)^{n+1}}.\label{b8}$$

For convenience of estimate, we amplify the domain $D_{
& \leq&C_R \int_{\substack{|y - x_r|>r\\|y|\leq R, \, y_n>0}}\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}\frac{1}{|x - y|^{n+1}}dy.\label{b6}\end{aligned}$$ 

 Here we have used the fact that the $ \alpha$-harmonic function $ u$ is bounded in the region $ $ D_{R, r}=\left\{y=(y',y_n)\left|\right.|y - x_r|>r,~|y|<R,~y_n>0\right\}.$$ The bounce depend on $ R$, however is independent of $ r$, since $ D_{R, r_1 } \subset D_{R, r_2}$ for $ r_1 > r_2$. For each such fixed open world $ D_{R, r}$, the bound of the $ \alpha$-harmonic function $ u$ can be derive from the inner smoothness (attend, for instance [ @BKN ] and [ @FW ]) and the estimate up to the limit (see [ @RS ]). 

 Set   $ y=(y',y_n), \sigma=|y'|$, for fixed   $ x$    and sufficiently large   $ r$, we accept 

 $ $ \begin{aligned } 
 \left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}&\leq&\frac{Cr^{\alpha/2}}{|2y_nr-|y|^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|\sigma^2 - 2y_nr+y_n^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|(y_n - r)^2+\sigma^2 - r^2|^{\alpha/2}}.\label{b7}\end{aligned}$$ 

 and $ $ \frac{1}{|x - y|^{n+1}}\leq\frac{C}{(1+|y|)^{n+1}}\leq\frac{C}{(1+|y'|)^{n+1}}=\frac{C}{(1+\sigma)^{n+1}}.\label{b8}$$ 

 For convenience of estimate, we amplify the domain $ D _ {
&\leq&F_R \int_{\substack{|y-x_r|>r\\|y|\leq R, \, y_n>0}}\left[\frac{2x_nr-|r|^2}{2t_nr-|y|^2}\rijht]^{\alphz/2}\frac{1}{|x-y|^{n+1}}ay.\label{b6}\end{aligned}$$

Here we hate uwed tye fact that the $\alpha$-farmonic vunction $u$ iw bounded ii the rennon $$D_{D,v}=\left\{v=(y',b_n)\left|\right.|y-x_r|>r,~|i|<R,~y_n>0\right\}.$$ Tha bound dependv un $R$, however is independent of $r$, sincq $D_{R, r_1} \xuhset D_{R, r_2}$ for $t_1>r_2$. Fog qach such fixed open domain $D_{R,r}$, the bohnd of uhe $\alpha$-harmonic gunction $u$ can be derived vrom the interior smoohhness ( see, for unstance [@BKN] and [@FW]) and the estimaje up to the boundary ( see [@RS]).

Set  $y=(y',y_n), \sigma=|y'|$, fot ridgd  $x$  and sufhiciennly large  $r$, wc have

$$\bagin{alibned}
\left[\frac{2x_nv-|x|^2}{2y_nr-|b|^2}\rigyt]^{\alpha/2}&\leq&\frac{Cr^{\alpha/2}}{|2b_nr-|y|^2|^{\alpha/2}}\nonumber\\&=&\fras{Cr^{\alpha/2}}{|\shgja^2-2y_nr+y_n^2|^{\alpha/2}}\nonumvee\\&=&\frac{Wr^{\al[ha/2}}{|(y_v-e)^2+\siema^2-d^2|^{\akpga/2}}.\labep{b7}\eid{aligned}$$

ans $$\frac{1}{|x-y|^{n+1}}\lew\frac{C}{(1+|y|)^{n+1}}\leq\frac{C}{(1+|y'|)^{n+1}}=\ftas{B}{(1+\xigma)^{n+1}}.\label{b8}$$

Ror cogvqnience of estimate, we amplify the domapn $D_{
&\leq&C_R \int_{\substack{|y-x_r|>r\\|y|\leq R, \, y_n>0}}\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}\frac{1}{|x-y|^{n+1}}dy.\label{b6}\end{aligned}$$ Here we the that the function $u$ is The depends on $R$, is independent of since $D_{R, r_1} \subset D_{R, r_2}$ $r_1>r_2$. For each such fixed open domain $D_{R,r}$, the bound of the $\alpha$-harmonic $u$ can be derived from the interior smoothness ( see, for instance [@BKN] [@FW]) the up the boundary ( see [@RS]). Set $y=(y',y_n), \sigma=|y'|$, for fixed $x$ and sufficiently large $r$, we $$\begin{aligned} \left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}&\leq&\frac{Cr^{\alpha/2}}{|2y_nr-|y|^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|\sigma^2-2y_nr+y_n^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|(y_n-r)^2+\sigma^2-r^2|^{\alpha/2}}.\label{b7}\end{aligned}$$ and $$\frac{1}{|x-y|^{n+1}}\leq\frac{C}{(1+|y|)^{n+1}}\leq\frac{C}{(1+|y'|)^{n+1}}=\frac{C}{(1+\sigma)^{n+1}}.\label{b8}$$ For convenience of estimate, amplify the domain $D_{
&\leq&C_R \int_{\substack{|y-x_r|>r\\|y|\leq R, \, Y_n>0}}\left[\frac{2X_nr-|x|^2}{2y_Nr-|y|^2}\RigHt]^{\AlphA/2}\fraC{1}{|x-y|^{n+1}}dy.\label{b6}\enD{AligNed}$$

Here we have used the faCt thaT tHE $\alpHA$-hArmonIc functIOn $U$ IS boUnDeD in ThE ReGion $$D_{r,r}=\lEft\{y=(y',y_n)\Left|\right.|y-X_r|>r,~|Y|<R,~Y_n>0\right\}.$$ The boUNd Depends on $R$, HowEver is indepeNdeNt of $r$, sInCe $D_{r, R_1} \subsEt D_{r, r_2}$ for $R_1>r_2$. For eACh such Fixed open DoMAin $D_{R,r}$, THe bound OF ThE $\alpHa$-harmonic functioN $U$ cAN be derived from The intErIOr SMOotHneSs ( see, for inStAnce [@Bkn] and [@FW]) aND tHE EStiMAte up to the bouNdary ( see [@RS]).

SET  $y=(y',Y_n), \sigmA=|y'|$, For FIxed  $x$  aNd sufFiCIenTly large  $r$, we Have

$$\Begin{aligNed}
\lefT[\Frac{2x_nr-|X|^2}{2Y_nr-|y|^2}\rigHt]^{\alphA/2}&\leQ&\frAc{Cr^{\ALpHa/2}}{|2Y_nr-|Y|^2|^{\aLPha/2}}\NOnUmbER\\&=&\frAc{Cr^{\alphA/2}}{|\sIgMa^2-2y_nr+Y_n^2|^{\alPHA/2}}\NOnumBer\\&=&\Frac{cr^{\alpHa/2}}{|(y_n-r)^2+\sigma^2-r^2|^{\alPha/2}}.\LabeL{B7}\enD{aligNed}$$

anD $$\fraC{1}{|x-Y|^{n+1}}\leq\Frac{C}{(1+|y|)^{N+1}}\leq\fRaC{C}{(1+|y'|)^{n+1}}=\frac{C}{(1+\sigma)^{n+1}}.\LabeL{b8}$$

For convEniEnCe oF eStimaTE, we ampLifY thE domain $d_{
&\leq&C_R \int_{\substack{ |y-x_r|>r\ \|y|\ leq  R,  \ , y_ n>0} }\left[\frac{2 x _nr- |x|^2}{2y_nr-|y|^2}\ri ght]^ {\ a lpha / 2} \frac {1}{|x- y |^ { n +1} }d y. \la be l {b 6}\en d{a ligned} $$

Here w e h av e used the f a ct  that the$\a lpha$-harmon icfuncti on  $u $  is b oun ded i n ther egion$$D_{R,r} =\ l eft\{y = (y',y_n ) \ le ft|\ right.|y-x_r|>r,~ | y| < R,~y_n>0\right \}.$$Th e  b o u nddep ends on $R $,  howe v er is i n de p e n den t  of $r$, sinc e $D_{R, r_ 1 } \ subset  D _{R ,  r_2}$  for$r _ 1>r _2$. For ea ch s uch fixed  opend omain $ D _{R,r}$ , thebou ndof t h e$\ alp ha $ -ha r mo nic fun ction $u $ca n bederi v e d from  th e in terio r smoothness( s ee,f orinsta nce [ @BKN ]and [ @FW])and t he  estimate up to  the  boundary  (se e [ @R S]).Set  $ y=( y', y_n), \ sigma=| y '|$ ,f o r  f ixed  $x$  and suf fi c i en tly larg e  $r$ ,  w eh ave

$$\ be gin {ali g n ed}
\ left [ \f rac{2x_n r-|x|^ 2 }{ 2y _nr-|y| ^2 }\righ t] ^{\ alp ha/2} & \leq &\frac {Cr^{\al pha/2 } }{|2y_nr-|y|^2 | ^{\alpha/2}}\ n on u m be r \\&= &\f rac{Cr^{\al pha/ 2 }}{| \sig m a^ 2-2 y _nr+y _n^2| ^{ \ al p ha/2}}\nonumber\\&= &\ frac{C r^{\a lpha/2}}{|(y_ n-r)^2+\si g m a ^2-r^2|^ {\al p ha / 2}}.\label{b7} \end{ aligned}$$ 
and $$\ frac{ 1}{|x-y| ^{n+1}}\l e q \frac{C} {(1 +|y |)^ {n+ 1 } }\ leq\frac{C}{( 1 + |y'| )^ {n+1}}= \fr ac{C}{( 1+\ sig ma) ^{n +1 }}.\label {b8}$$
Fo rco nv eni enceo f estima te , w eamp lifyt he dom ain $ D_{
&\leq&C_R \int_{\substack{|y-x_r|>r\\|y|\leq_R, \,_y_n>0}}\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}\frac{1}{|x-y|^{n+1}}dy.\label{b6}\end{aligned}$$

Here we have used_the fact_that_the $\alpha$-harmonic_function_$u$ is bounded_in the region_$$D_{R,r}=\left\{y=(y',y_n)\left|\right.|y-x_r|>r,~|y|<R,~y_n>0\right\}.$$ The bound depends_on $R$, however_is_independent of $r$, since $D_{R, r_1} \subset D_{R, r_2}$ for $r_1>r_2$. For each such_fixed_open domain_$D_{R,r}$,_the_bound of the $\alpha$-harmonic function_$u$ can be derived from_the interior_smoothness ( see, for instance [@BKN] and [@FW])_and_the estimate up_to the boundary ( see [@RS]).

Set  $y=(y',y_n), \sigma=|y'|$, for_fixed  $x$  and sufficiently large  $r$,_we have

$$\begin{aligned}
\left[\frac{2x_nr-|x|^2}{2y_nr-|y|^2}\right]^{\alpha/2}&\leq&\frac{Cr^{\alpha/2}}{|2y_nr-|y|^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|\sigma^2-2y_nr+y_n^2|^{\alpha/2}}\nonumber\\&=&\frac{Cr^{\alpha/2}}{|(y_n-r)^2+\sigma^2-r^2|^{\alpha/2}}.\label{b7}\end{aligned}$$

and $$\frac{1}{|x-y|^{n+1}}\leq\frac{C}{(1+|y|)^{n+1}}\leq\frac{C}{(1+|y'|)^{n+1}}=\frac{C}{(1+\sigma)^{n+1}}.\label{b8}$$

For_convenience_of_estimate, we amplify the_domain $D_{