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

摘要

正 四、程序流程图本算法程序实现过程为: 1、令tn+■=tn+h; 2、利用（5）式递推计算各阶差分,由（4）和（6）式同时得到Xn+■P和■; 3、将■代入微分方程（1）,得到一组隐式方程,用X(n■P作为初值预报,用Newton—Raphson迭代求得校正值Xn+1C即为tn+1处的数值解; 4、用（15）式计算截断误差EK,并与用户规定的最大容许误差Err比较,若超过最