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.