www.gusucode.com > 数学建模竞赛：露天矿的开采问题，运用lingo与matlab编程，内附代码 > 数学建模竞赛：露天矿的开采问题，运用lingo与matlab编程，内附代码/work4/aaa.m

    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（时间，值（行,列）