刚性常微分方程初值问题的一种数值解法及其程序实现(上)

摘要

本文提出了一种求解刚性常微分方程初值问题的数值解法。该法虽与Brayton等人提出的BDF方法及Bokhoven提出的PBD方法数学上完全等价,但把BDF及PBD从结构上统一了起来,并进一步减少了计算量。在该算法的程序设计中,采取了一系列自适应技术以防止不稳定的产生、保证计算精度及节省计算时间。最后,给出了程序框图及实验结果。