我想在 Matlab 中实现 P2（甚至 P3...）有限元。我需要能够控制刚性和质量矩阵的构造（$\int \nabla \varphi_i\cdot \nabla \varphi_j$和$\int \varphi_i \varphi_j$) （因为我想在之后解决一些精度计算问题，并且可能使用 Advanpix 超越机器精度......）

如果已经有一些很好的 Matlab 代码可以做到这一点，请指点我。如果没有，是否有任何好的参考资料可以让我学会正确地做到这一点？

如果可能的话，我想避免手动进行所有计算，因为我确信其他人之前已经完成了它们:) ...（找到正确的基函数，明确计算相邻基函数乘积的所有积分等。 .)