For most orthogonal systems and their corresponding Fourier series, the study of the almost everywhere convergence for functions in Lp requires very complicated research, harder than in the case of the mean convergence. For instance, for trigonometric series, the almost everywhere convergence for functions in L2 is the celebrated Carleson theorem, proved in 1966 (and extended to Lp by Hunt in 1967). In this paper, we take the system jn (x) = sqrt(2 ( + 2 n + 1)) J + 2 n + 1 (x) x- - 1, n = 0, 1, 2, ... (with J being the Bessel function of the first kind and of the order ), which is orthonormal in L2 ((0, ), x2 + 1 d x), and whose Fourier series are the so-called Fourier-Neumann series. We study the almost everywhere convergence of Fourier-Neumann series for functions in Lp ((0, ), x2 + 1 d x) and we show that, surprisingly, the proof is relatively simple (inasmuch as the mean convergence has already been established). © 2009 Elsevier B.V. All rights reserved.