考虑化学反应动力学模型，设三种化学物质的浓度随时间变化的函数为y1(t)，y2(t)，y(t)，则浓度由下列方程给出
其中k1和k2是两个反应的速度常数，假定初始浓度为y1(0)=y2(0)=y3(0)=1．取k1=1,分别用k2=10,100,1000进行试验．对每个k2，分别用四阶R-K方法，四阶阿当姆斯预测-校正法及梯形法求解，针对不同步长，比较各种方法的精度和稳定性．从t=0开始计算到近似稳定状态或可以明显看出解不稳定或方法无效为止．