×

求解刚性和非刚性ODE的单个代码。 (英语) Zbl 0601.65055号

在这篇有趣的文章中,作者描述了一个对具有初始条件(*)\(y'=f(x,y),y(a)=c\),\(y\ In{mathbb{R}}^s)的一阶常微分方程组的刚性和非刚性问题都有效的代码。策略如下:集成是从可变阶Adams PECE方法开始的。如果在这部分积分过程中,满足某些测试,则假设系统是刚性的,然后引入隐式后向微分多步公式(BDF)。这项工作的背景是,许多实际刚性问题都包含一个非刚性子系统。如果可以分离(*)的非刚性子系统,则可以保留亚当斯公式的其余部分。作者提出了系统分区的自动测试(*)。如果问题(*)很难解决,并且作者的测试可以正确识别瞬态区域,那么在初始范围内使用Adams代码将进一步提高效率。对于非刚性问题,该算法将像通常的变阶Adams代码一样执行。
审核人:H.阿德

MSC公司:

65升05 常微分方程初值问题的数值方法
65L20英寸 常微分方程数值方法的稳定性和收敛性
34A34飞机 非线性常微分方程和系统
PDF格式BibTeX公司 XML格式引用
全文: 内政部