function phi1 = concatmuphi(phi,mu,m,n)

maxs = 0;
dims = repmat([m],[1 n]);

phi1 = zeros(1,m^n);

for xind = 1:m^n
  x = ind2coord(xind,dims);
  xind1 = coord2ind(x(2:end),dims(2:end));
  phi1(xind) = mu(x(1)) + phi(xind1);
end
return