rigid.m

来自「简单的解常微分方程的matlab程序！压缩包」· M 代码 · 共 6 行

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