How can one compute the *Neumann spectral fractional Laplacian* of power function, $(-\Delta)^s x^{\alpha}$, with $\alpha >0$, in an interval $(0,1)$. I'm only aware of the formula in the whole space.