This present investigation is contemplated to provide Legendre spectral collocation method for solving multi-Pantograph delay boundary value problems (BVPs). In this regard, an equivalent integral form of such BVPs has been considered. The proposed method is based on Legendre–Gauss collocation nodes and Legendre–Gauss quadrature rule. Convergence analysis associated to the presented scheme has been provided to show its applicability theoretically. Some numerical examples are given to demonstrate the efficiency, accuracy, and versatility of our method. Numerical results confirm the theoretical predictions and are superior with respect to several recent numerical methods including Hermite collocation approach, Laguerre collocation technique and the reproducing kernel method.