sigma = zeros(D,D,K);