我正在尝试使用教科书：Nodal Discontinuous Galerkin Methods Algorithms, Analysis, and Applications，让一些代码适用于一维泊松方程。

我使用以下公式（对于同质情况）：

为了解释不均匀的情况，教科书提出以下建议：

其中由齐次情况给出。当我尝试在上面等式的右侧实现添加的项时，我得到以下响应，这似乎几乎是正确的；但是，第一项和最后一项与它们指定的 bc 有很大的不同：$\mathcal A$

我使用以下代码实现了额外的 rhs 术语：

e = zeros((basisdegree+1),numelements);
enp = e; enp(end) = 1; e1 = e; e1(1) = 1;
extraterm = Dr*(enp*b - e1*a);

wheree只是一个与解决方案大小相同的向量u，其中 1 位于第一个和最后一个元素中。Dr的齐次方程相同。$\mathcal{M}^{-1} \mathcal{S}$