function K1 = getK1(K,m,n)
dims = repmat([m],[1 n]);
dims1 = repmat([m],[1 n-1]);

K1 = zeros(m^(n-1),1);

for xind2n = 1:m^(n-1)
  x2n = ind2coord(xind2n,dims1);
  for x1 = 1:m
    xind = coord2ind([x1 x2n],dims);
    K1(xind2n) = K1(xind2n) + K(xind);
  end
end
return