gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-22/rigid.m
function dy = rigid(t,y) dy = zeros(3,1); % 行向量 dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2);