dde卷积积分

Question

我正在尝试使用dde求解该微分方程：

我对积分项有疑问。在下面显示的代码中，tau是延迟，但是我不能仅指定一个常量值，因为它也在从0到t的整数中。有人知道如何处理吗？谢谢！

function sol = exer_3

   sol = dde23(@exer3f,tau,[0; 0],[0, 10]);
   figure
   plot(sol.x,sol.y)


function v = exer3f(t,y,Z)

    k = 125; m = 5; F = 1; w = 8;
    c=@(t)exp(-t^2);
    ylag = Z(:,1);
    v = zeros(2,1);
    v(1)=y(2);
    v(2) = -(k*y(1) - F*cos(w*t) + integral(@(tau)c(tau).*ylag(1), 0, t,'ArrayValued',true))./m;

我还将这个问题发布在MathWorks主页上。