**7. Conclusions**

In this paper, we proposed a hybrid method based on finite difference and Haar wavelets approximations. The scheme is applied for the numerical solution of (1 + 1)- and (1 + 2)-dimensional time fraction partial differential equations. The accuracy and applicability of the scheme is validated through some test problems. The tabulated data and graphical solution show that the scheme works very well for time fractional problems.

**Author Contributions:** Conceptualization, A.G. and S.H.; Methodology, A.G.; Software, A.G.; Validation, S.H., M.H. and M.A.J.; Formal Analysis, A.G.; Investigation, M.A.J.; Resources, P.K.; Writing–Original Draft Preparation, A.G.; Writing–Review and Editing, M.H.; Visualization, P.K.; Supervision, S.H. and P.K.; Project Administration, P.K.; Funding Acquisition, P.K.

**Funding:** The project was supported by the Center of Excellence in Theoretical and Computational Science (TaCS-CoE), Faculty of Science, King Mongkut's University of Technology Thonburi (KMUTT).

**Acknowledgments:** The authors are thankful to anonymous reviewers for their fruitful suggestion which improved the quality of the manuscript. Also we are thankful for the financial support of the Center of Excellence in Theoretical and Computational Science(TaCS-CoE), Faculty of Science, King Mongkut's University of Technology Thonburi (KMUTT).

**Conflicts of Interest:** The authors declare no conflict of interest.
