La solución fundamental de la ecuación de difusión fraccional en espacio-tiempo multidimensional se estudia aplicando el principio de subordinación, que establece una relación con la función gaussiana clásica. Se obtienen representaciones integrales en términos de funciones de Mittag-Leffler para la solución fundamental y el núcleo de subordinación. Las representaciones integrales obtenidas se utilizan para la evaluación numérica de la solución fundamental para diferentes valores de los parámetros.