Para una ecuación de difusión fraccional en el espacio con un coeficiente de superdifusión no lineal y con la presencia de un efecto de retraso, se construye un método numérico de rejilla. Se utilizan procedimientos de interpolación y extrapolación para tener en cuenta el retraso funcional. En cada paso de tiempo, el algoritmo se reduce a resolver un sistema lineal con una matriz principal que tiene dominancia diagonal. Se demuestra la convergencia del método en la norma máxima. Se presentan los resultados de experimentos numéricos con retrasos constantes y variables.