带有时滞项的微积分方程组，大神可以看一下吗？

含时滞的微分方程怎么求特征方程

在数学领域中,时滞微分方程, 或延时微分方程(dde) 是一类微分方程, 其中未知函数的在确定时刻的导数由先前时刻函数所决定.

m-文件rigid.m第二行dy=zeros(3,

说y1(0)=0,y2(0)=1,y3(0)=1自变量0?

解微程

y1'=y2y3

y2'=-y1y3

y3'=-0.5y1y2

y1(0)=0,y2(0)=1,y3(0)=1

1、建立m-文件rigid.m：

function dy=rigid(t,y)

dy=zeros(3,1);

dy(1)=y(2)*y(3);

dy(2)=-y(1)*y(3);

dy(3)=-0.51*y(1)*y(2);

2、取t0=0,tf=12,输入命令做图：

[T,Y]=ode45('rigid',[0 12],[0 1 1]);

plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')

m-文件rigid.m第二行dy=zeros(3,

三阶，因为导数最高是第三次求导，不用公式，就是阶数的定义，（一个方程如果最高次导是(y')^3，那它也是一阶方程，不是三阶，因为阶数是看最高几次求导的），只要看式子里最高是几次导就可以判断了，大佬，我也是学生来的，加个私信以后有问题可以互相讨论