function phi = multispike(m,n,inds,val)
dims = repmat([m],[1 n]);
phi = zeros(1,m^n);

phi(inds) = phi(inds) + val;

for v=n:-1:1
  xind = find(phi==v);
  for i=xind(:)'
    x = ind2coord(i,dims);
    Y = allrho1(x,m);
    Yind = coord2ind(Y,dims);
    badYi = find(v - phi(Yind) > 1);
    phi(Yind(badYi)) = 0*phi(Yind(badYi)) + v - 1;
  end
end
return