用matlab可以解如下方程组吗？
 

初次遇到这样的方程，困扰很久，不知道如何求解，希望各位高手赐教，谢谢！：）


da/dt=da^2/dx^2+f(a,u,j) (1) 其中f(a,u,j)=-a+exp(a)-exp(-j+a)+12.5*u

dj/dt=-j-10*(j+exp(a)-1) (2)

du/dt=-u-2*(u-a) (3)

边界条件：

da(0,t)/dx=da(L,t)/dx=0 变量x的取值范围[0,L]

我现在只会用pdepe函数解方程(1), 但加了 (2)， (3)两个约束方程后，就不知道如何编程了。