mu(k,:) = mu(k,:) + res(k,i).*X(k,:);