用于加性刚性常微分方程组的三阶非齐次方法
数值分析
2009-02-24 v1 经典分析与常微分方程
摘要
本文构造了一种求解加性分裂自治刚性常微分方程组的三阶方法。所构造的加性方法对隐式部分是L-稳定的,并允许使用雅可比矩阵的任意近似。与我们之前的论文相反,第四阶段是显式的。因此,所构造的方法由于第四阶段中间数值公式的L-稳定性而具有良好的稳定性,但每步的计算成本更低。基于局部误差和稳定性控制的自动步长选择得以执行。误差和稳定性控制的估计值是在没有显著额外计算成本的情况下获得的。数值实验表明了所实现积分算法的可靠性和效率。
引用
@article{arxiv.0902.3723,
title = {A non-homogeneous method of third order for additive stiff systems of ordinary differential equations},
author = {Evgeny Novikov and Anton Tuzov},
journal= {arXiv preprint arXiv:0902.3723},
year = {2009}
}
备注
12 pages