LAUSR

This paper deal with the numerical method, based on the operational matrices of the Haar wavelet orthonormal functions approach to approximate solutions to a class of coupled systems of time-fractional order partial differential equations (FPDEs.). By introducing the fractional derivative of the Caputo sense, to avoid the tedious calculations and to promote the study of wavelets to beginners, we use the integration property of this method with the aid of the aforesaid orthogonal matrices which convert the coupled system under some consideration into an easily algebraic system of Lyapunov or Sylvester equation type. The advantage of the present method, including the simple computation, computer-oriented, which requires less space to store, time-efficient, and it can be applied for solving integer (fractional) order partial differential equations. Some specific and illustrating examples have been given; figures are used to show the efficiency, as well as the accuracy of the, achieved approximated results. All numerical calculations in this paper have been carried out with MATLAB.