function phi = randphi_hlp(m,n)
phi = zeros(m^n,1);
combs = generate_combs(m^n,n);
nc = size(combs,1);
rp = randperm(nc);
rc = rp(1);
com = combs(rc,:);
phi = cumsum(com');
return