function phi = fillmask(y,m)
n = length(y);
phi = zeros(1,m^n);
dims = repmat([m],[1 n]);
for xind = 1:m^n
  x = ind2coord(xind,dims);
  phi(xind) = rho(x,y);
end
return