We provide a novel approach to the numerical solution of the family of nonlocal elliptic equations $(-\Delta)^su=f$ in $\Omega$, subject to some homogeneous boundary conditions $\mathcal{B}(u)=0$ on $\partial \Omega$, where ...