把微分代数方程组转化为微分方程组时需要加Mass质量矩阵，options=odeset('mass',M).
那么在微分代数方程组加入时滞后，再转化为时滞微分代数方程组，使用Mass质量矩阵，options还是使用odeset()函数调用么，我用了ddeset()函数调用，即
options=ddeset('mass',M),运行时是说mass无法识别