Я пытаюсь решить систему нелинейных числовых DE на mathematica 12, но я получаю сообщение: 
Lx = Ly = 10^-7;
NDSolve[{
mz[x,y](D[my[x,y],{x,2}]+D[my[x,y],{y,2}])-
my[x,y](D[mz[x,y],{x,2}]+D[mz[x,y],{y,2}])==0,
mz[x,y](D[mx[x,y],{x,2}]+D[mx[x,y],{y,2}])-
mx[x,y](D[mz[x,y],{x,2}]+D[mz[x,y],{y,2}])==0,
my[x,y](D[mx[x,y],{x,2}]+D[mx[x,y],{y,2}])-
mx[x,y](D[my[x,y],{x,2}]+D[my[x,y],{y,2}])==0,
(D[mx[x,y],x]/.x->0)==0,
(D[mx[x,y],x]/.x->Lx)==0,
(D[mx[x,y],y]/.y->0)==0,
(D[mx[x,y],y]/.y->Ly)==0,
(D[my[x,y],x]/.x->0)==0,
(D[my[x,y],x]/.x->Lx)==0,
(D[my[x,y],y]/.y->0)==0,
(D[my[x,y],y]/.y->Ly)==0,
(D[mz[x,y],x]/.x->0)==0,
(D[mz[x,y],x]/.x->Lx)==0,
(D[mz[x,y],y]/.y->0)==0,
(D[mz[x,y],y]/.y->Ly)==0,
{mx[x, y],my[x, y],mz[x, y]}, {x, 0, Lx}, {y, 0, Ly}]
Я не понимаю, почему это так? Потому что мои граничные условия линейны.