我的 FEM 方案使用具有双线性形状函数的 4 节点四边形单元。我要解决的简单问题是。$\nabla ^2 f = 5$但随着节点数量的增加，解的图变得更平坦。我很困惑为什么会发生这种情况。下面的两个图显示了当我将网格的大小从 7x7 更改为大约。30x30，情节变得更加分散。边界条件也是边界上的$f = 0$

已修复，我的代码中有一些拼写错误