function_y.m

关于病态系统的识别以及求解的方法。其中包括几种解法的求解精度和速度的比较。

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