我需要在上找到的数值解，其中。$-\Delta U = f$$\Omega = [0,1]^2$$U|_{\partial \Omega} = 0$

我找到了一个方法：POISSONFDin NMLibforOctave with signature: [U,X,Y]=POISSONFD(A,C,B,D,NX,NY,FUN,BOUND)，

但我不明白如何指定绑定。就像是：

BOUND = inline('x * y * (y - 1) * (x - 1)', 'x', 'y').

这将在边界上提供良好的价值（我希望）。我可以用吗？还是会变坏？或者我应该使用其他包？