求解stiff常微分方程组只计算函数值的一类L-稳定显式单步法

摘要

一 引 言 讫今为止,解stiff常微分方程组的初值问题 y’=f(x,y),y(x_0)=η,x∈[a,b],η,y,f∈R^m,x_0=α,的绝大多数可行的数值方法都需要计算函效f(x,y)的Jacobi矩阵αf/αy及对与αf/αy有关的某矩阵作LU分解.当m很大时,其计算量和存贮量都是惊人的.而显式的Runge-Kutta方法,虽然只计算函数f的直,但其绝对稳定域是有界的,故不宜于解stiff方程.