function test_cy
[t,y]=ode45(@fn,[0 300],[1 0 0 0 0 6])

function dy=fn(t,y)
dy=zeros(6,1);
dy(1)=-0.05*y(1)*y(6)+0.11*y(2)*y(5);
dy(2)=0.05*y(1)*y(6)-0.11*y(2)*y(3)-0.215*y(2)*y(6)+1.228*y(3)*y(5);
dy(3)=0.215*y(2)*y(6)-1.228*y(3)*y(5)-0.242*y(3)*y(6)+0.007*y(4)*y(5);
dy(4)=0.242*y(3)*y(6)-0.007*y(4)*y(5);
dy(5)=0.5*y(1)*y(6)-0.11*y(2)*y(5)+0.215*y(2)*y(6)+1.228*y(3)*y(5)+0.242*y(3)*y(6)-0.007*y(4)*y(5);
dy(6)=-dy(5);

你的[0 1000]区间太大了好像
function test_cy
[t,y]=ode45(@fn,[0 100],[1 0 0 0 0 6])

function dy=fn(t,y)
dy=zeros(6,1);
dy(1)=-0.05*y(1)*y(6)+0.11*y(2)*y(5);
dy(2)=0.05*y(1)*y(6)-0.11*y(2)*y(3)-0.215*y(2)*y(6)+1.228*y(3)*y(5);
dy(3)=0.215*y(2)*y(6)-1.228*y(3)*y(5)-0.242*y(3)*y(6)+0.007*y(4)*y(5);
dy(4)=0.242*y(3)*y(6)-0.007*y(4)*y(5);
dy(5)=0.5*y(1)*y(6)-0.11*y(2)*y(5)+0.215*y(2)*y(6)+1.228*y(3)*y(5)+0.242*y(3)*y(6)-0.007*y(4)*y(5);
dy(6)=-dy(5);


有结果

我也是新手，交流交流

这个帖子帮到了我，必须提出来感谢！！！！！