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

    function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);