rigid.m

自己编的四节龙格琨塔程序

%example
function dy = rigid(t,y,u)
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)*u;