rigid.m

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

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);