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

    function dy = rigid(t,y)
dy = zeros(5,1);    % a column vector
%添加系数
dy(1) = y(1)+2*y(2)+y(4);
dy(2) = y(1)-y(2)*y(3)+y(4)-y(2);
dy(3) = y(1)*y(2)-y(3)+y(5)+y(3);
dy(4) = y(4)+y(2)-y(4);
dy(5) = y(5)-y(3)+y(5);