非刚性常微分方程的数值解法通常会用四阶龙格库塔算法,其matlab函数对应ode45. 对于dy/dx = f(x,y),y(0)=y0. 其四阶龙格库塔公式如下: 对于通常计算,四阶已经够用,四阶以上函数f(x,y)计算工作量大大增加而精度提高较慢. 下面以龙格库塔法解洛伦兹方程为例: matlab代码如下: main.m: clear all; close all; clc; %系统龙格库塔法 [t,h] = ode45(@test_fun,[ ],[ ]); plot3(h(:,),h(:…