In the current study, a numerical scheme based on the Haar wavelet is proposed to solve a coupled system of fractional-order integral-differential equations. The proposed method is to derive the operational matrix of fractional-order integration, and that is used to transform the main problem to a system of algebraic equations. Additionally, the convergence analysis theorem of this system is rigorously established and the numerical results show that the proposed method is practicable and effective for solving such kinds of problem. (C) 2019 InternationalAssociation forMathematics andComputers in Simulation (IMACS). Published by Elsevier B.V. All rights reserved.