Abstract
The main focus of this paper was to find the approximate solution of a class of second-order multi-pantograph delay differential equations with singularity. We used the shifted version of Vieta-Lucas polynomials with some symmetries as the main base to develop a collocation approach for solving the aforementioned differential equations. Moreover, an error bound of the present approach by using the maximum norm was computed and an error estimation technique based on the residual function is presented. Finally, the validity and applicability of the presented collocation scheme are shown via four numerical test examples.