求助各位大侠!

T=size(ret,1);
vsquare=var(ret);


for ilag = 1 : lag
    c = cov(ret(1:end - ilag), ret(ilag + 1:end));
    vsquare = vsquare + 2 * (1 - ilag / lag) * c(1, 2);
end

% for i=1:size(x,2)