④解微分方程组的四阶龙格库塔法
  令h为计算步距,ih为离散的时间点。令
    
  则得到
    
  其中的四个n维向量由下式计算。
    
  解微分方程组的四阶龙格库塔法,是成熟的数值计算方法。方法的实质是对连续的微分方程进行离散化。