aaa.m

来自「数学建模竞赛：露天矿的开采问题」· M 代码 · 共 3 行

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5 1e-5 1e-5]);%精度
[T,Y] = ode45(@rigid,[0:.001:1],[0 1 1 1 1],options);% 时间步长  初始值
 plot(T,Y(:,1:3))%plot（时间，值（行,列）