MATLAB求解拉普拉斯代码-ZetaTrap2D:二维曲线上拉普拉斯、亥姆霍兹和斯托克斯势的局部校正梯形规则

MATLAB求解拉普拉斯代码二维边界积分算子的高阶局部校正梯形规则 这是论文随附的“zeta 校正正交”的 MATLAB 代码： B. Wu 和 PG Martinsson，Zeta 校正：构建奇异积分算子的校正梯形规则的新方法。 （2020 年，arxiv） 它包含修正的梯形规则，这些规则对于平滑闭合曲线上的拉普拉斯、亥姆霍兹和斯托克斯层势是稳定且高阶准确的。 作者：吴博伟，2020/7，更新于2021/1 还包含从 Alex Barnett 的包修改的支持功能 另见表面积分的相关代码。 依赖的注意事项： 最多 42 个订单的 Zeta 更正使用预先计算的权重。 （它们对应于kapur_rokhlin_sep_log.m函数的输入k <= 21 ）对于更高阶的更正，需要来自 Symbolic Math Toolbox 的vpa函数。 （此存储库中的测试文件不需要 Symbolic Toolbox。） 例子 以下是使用高阶 zeta 正交求解斯托克斯和亥姆霍兹外边界值问题的测试文件的一些示例输出。 收敛性根据边界上使用的点数显示。 主要测试文件说明： test_sto2d_bie.m