“…where F 2 is the Gegenbauer wavelet coefficients vector. Now by substituting Equations (26), (32)-(34), (41) and (42) into Equations (24) and (32), (33), (35), and (41)-(43) into Equation (25), respectively, then using operational matrices of integration, we get the residuals functions R 1 (x, t) and R 2 (x, t) for this system as follows:…”