我想解决下面的问题

$$\begin{equation} \begin{aligned} \eta u-\Delta u &=f, &\text{in $\Omega$}\\ \end{aligned} \end{equation}$$

在哪里$\Omega=(0,1)\times (0,1)$, 我的边界是这样的$\frac{\partial u}{\partial n} = 0$如果$y=0$要么$1$， 和$u=0$如果$x=0$要么$1$. 使用幽灵节点可以处理诺伊曼边界的低阶离散化，但我面临的问题是我想在诺伊曼边界中包含四个角点而不是在狄利克雷边界中，这给我带来了问题，因为我不能在拐角处使用内部方程。

我希望我的问题很清楚。