sigma(:,:,k) = sigma(:,:,k)./Nk(k);