dA(t)/dt=a*b^1/2*exp[(1-i*c)*F(t-T)/2-(1-i*d)*G(t-T)/2]A(t-T)-a*A(t);
dB(t)/dt=-e*B(t)+f*D(t)*[1-B(t)]-g*B(t)-exp(-G(t))*[exp(F(t))-1]*abs(A(t))^2;
dC(t)/dt=-h*C(t)+j*E(t)*[1-C(t)]-k*C(t)-s*[1-exp(-G(t)]*abs(A(t))^2;
dD(t)/dt=D0-m*D(t)-2*f*D(t)*[1-B(t)]+2*g*B(t);
dE(t)/dt=E0-n*E(t)-2*j*E(t)*[1-C(t)]+2*k*C(t);
其中
F(t)=2*u*[2*B(t)-1]
G(t)=2*v*[2*C(t)-1]
其中变量都是时间t的函数，这里牵扯到时间延迟F(t-T),G(t-T),A(t-T), 其中T为恒定量，abs(A(t))^2为模的平方，其他a,b...都为常数，求这些函数随时间的变化关系曲线，能否用的DDE23 来求？我写了下面代码，运行起来也没发现错误，但没有关于T的周期性，感觉不对，那位高手给指点一下，十分感谢！！！
clear
sol=dde23(@DDE_SL,[5,5,5],[0;0;0;60;5;3.2;6],[0 20]); 
figure
plot(sol.x,sol.y(1,:));
figure
plot(sol.x,sol.y(2,:));
figure
plot(sol.x,sol.y(3,:));

function dydt=DDE_SL(t,y,Z)    
a=10;
b=0.2;
c=2;
d=2;
e=0.01;
h=0.01;
f=2;
j=2;
g=0.1;
k=0.1;
s=10;
m=0.01;
n=1;
D0=5;
E0=6;
u=1.6;
v=3;
dydt=[a*b^(1/2)*exp((1-i*c)*Z(1)/2-(1-i*d)*Z(2)/2)*Z(3)-a*y(1);
    -e*y(2)+g*y(4)*(1-y(2))-j*y(2)-exp(-y(7))*(exp(y(6))-1)*(abs(y(1)))^2;
    -f*y(3)+h*y(5)*(1-y(3))-k*y(3)-s*(1-exp(-y(7)))*(abs(y(1)))^2;
    D0-m*y(4)-2*f*y(4)*(1-y(2))+2*g*y(2);
    E0-n*y(5)-2*j*y(5)*(1-y(3))+2*k*y(3);
    4*u*y(2);
    4*v*y(3);];