我试图用mathematica去求解如下的反应扩散方程，因为是新手，有很多不明白的地方，具体模型如下：
其中区域一维Omega=（0，pi），d，r, gamma, alpha都是常数，可以取值为：
我的程序代码是这样的：
solution = NDSolve[{D[m[x, t], t] == D[m[x, t], x, x] + m[x, t] (0.5 a[x, t] - 1/(1 + m[x, t])),
0.5 D[a[x, t], t] == D[a[x, t], x, x] + 0.1 (1 - a[x, t]) - m[x, t] a[x, t], m[x, 0] == 0.1250 + 0.1 Cos[2 x], a[x, 0] == 0.444 - 0.1 Cos[2 x], Derivative[1, 0][m][0, t] == 0, Derivative[1, 0][m][Pi, t] == 0, Derivative[1, 0][a][0, t] == 0 Derivative[1, 0][a][Pi, t] == 0}, {m, a}, {x, 0, Pi}, {t, 0, 800}]
然后就出现了这样的提示：
之后不知道如何处理了，关键画图不知道如何画。
同时我想问下mathematica是否能处理变量带有时滞的情况。
求吧里大牛不吝赐教！多谢！

1. True的问题看这里：https://tieba.baidu.com/p/4850533592
2. 某个地方掉了个逗号，加上即可顺利求解。
3. 如果我记得不错，NDSolve还求不了含时滞的
偏微分方程。