时滞微分方程的求解。许多动力系统随时间的演化不仅依赖于系统当前的状态，而且依赖于系统过去

某一时刻或若千个时刻的状态，这样的系统被称为时滞动力系统。时滞非线性动力系统有着比用常微分方程所描述的动力系统更加丰富的动力学行为，例如，一阶的自治时滞非线性系统就可能出现混沌运动。时滞微分方程的一般形式为

式中：T≥0为时滞常数。在Matlab中提供了命令dde23来直接求解时滞微分方程。其调用格式为801=dde23(ddefun，lags，history，tspan，options)，

其中，ddfun为描述时滞微分方程的函数;lags为时滞常数向量;history为描述t≤to时的状态变量值的函数;tspan为求解的时间区间;options为求解器的参数设置。该函数的返回值sol是结构体数据，其中sol.x成员变量为时间向量l，sol.y成员变量为各个时刻的状态向量构成的矩阵，其每一个行对应着一个状态变量的取值。求解如下时滞微分方程组：

已知，在i≤0时，x(t)=5，x2(t)=0，x(1)=1，试求该方程组在[0，40]上的数值解。