In this article, we develop a new method to approximate numerically the fractional Laplacian of functions defined on \(\mathbb R\) , as well as some more general singular integrals. After mapping \(\mathbb R\) into a finite interval, we discretize the integral operator using a modified midpoint rule. The result of this procedure can be cast as a discrete convolution, which can be evaluated efficiently using the Fast-Fourier Transform (FFT). The method provides an efficient, second-order accurate approximation to the fractional Laplacian, without the need to truncate the domain. We first prove that the method gives a second-order approximation for the fractional Laplacian and other related singular integrals; then, we detail the implementation of the method using a fast convolution algorithm, and give numerical examples that support its efficacy and efficiency; finally, as an example of its applicability to an evolution problem, we employ the method for the discretization of the nonlocal part of the one-dimensional cubic fractional Schrödinger equation in the focusing case.
