En este documento, presentamos un algoritmo robusto para resolver numéricamente una familia de ecuaciones diferenciales integro-fraccionarias bidimensionales. El método de la onda Haar se mejora para incluir en su construcción el paso de la transformada de Laplace. Esta modificación ha demostrado reducir los errores acumulativos que se obtendrían en caso de usar la técnica regular de la onda Haar. Se discuten diferentes ejemplos para cumplir dos objetivos, la metodología y la precisión de nuestro nuevo enfoque.