In the present paper, we consider an effective computational method to analyze a coupled dynamical system with Caputo–Fabrizio fractional derivative. The method is based on expanding the approximate solution into a symmetry Haar wavelet basis. The Haar wavelet coefficients are obtained by using the collocation points to solve an algebraic system of equations in mathematical physics. The error analysis of this method is characterized by a good convergence rate. Finally, some numerical examples are presented to prove the accuracy and effectiveness of this method.
