rigid.m

来自「病态系统的配套的拉普拉斯反变换：从状态方程到状态方程的解」· M 代码 · 共 5 行

function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) =-21*y(1)+19*y(2)-20*y(3);
dy(2) = 19*y(1)-21*y(2)+20*y(3);
dy(3) =40*y(1)-40*y(2)-40*y(3);