<p>This study extends the Haar wavelet collocation method (HWCM) to solve the systems of fractional pantograph delay differential equations (FPDDEs) with initial value problems (IVP). From a numerical solution perspective, the pantograph delay operational matrix of Haar wavelets was derived. A comprehensive and clearly articulated methodology was presented, accompanied by flowcharts and pseudocode, to systematically outline the implementation process of HWCM. Additionally, the existence and uniqueness of the solutions were examined and the convergence of the proposed algorithm was analyzed. The effectiveness of the algorithm was assessed through four numerical experiments in which various error values and experimental convergence orders of the obtained solutions were computed and compared with the results from other recent methods. The presented data demonstrate the accuracy and efficiency of the proposed methodology. The HWCM shows extensive applicability, is suitable for both linear and nonlinear cases, and maintains manageable computational complexity with increasing number of equations or pantograph-type delay functions.</p>