mathematica数值求解非线性偏微分方程组

想数值解解这样一个非线性偏微分方程，mma里显示错误，有人能告诉我应该怎么改吗🙏🙏🙏