Resumen
In this paper, we present a robust algorithm to solve numerically a family of two-dimensional fractional integro differential equations. The Haar wavelet method is upgraded to include in its construction the Laplace transform step. This modification has proven to reduce the accumulative errors that will be obtained in case of using the regular Haar wavelet technique. Different examples are discussed to serve two goals, the methodology and the accuracy of our new approach.