The main purpose of this work is to provide a spectral method based on the operational matrices of the Legendre polynomials for solving neutral multi-pantograph equations. We analyze the convergence properties of the proposed method. It is shown that for the multi-pantograph equations of neutral type, the spectral method yields the exponential order of convergence. Some examples are given to demonstrate the high precision, fast computation and good performance of the new scheme.
