function F = fillsupsuma(m,n)
global myopt
myopt = optimset('Diagnostics','off');
myopt = optimset('Display','off');

F = zeros(m^n,1);
for xi=1:m^n
  x = i2c(xi,m,n);
  x = x/m;
  %x = (x - m/2)/m;
  a = quadprog(eye(n),0*x',[],[],x,1,[],[],[],myopt);
  a = a/sqrt(sum(a.^2));
  F(xi) = x*a;
end

return


function i=c2i(c,m,n)
dims = repmat(m,[1 n]);
i = coord2ind(c,dims);
return

function c=i2c(i,m,n)
dims = repmat(m,[1 n]);
c = ind2coord(i,dims);
return