function miscbound(n0)
global n
n=n0

%global optV
%optV = 1e6;

global maxr
maxr = 0;


while 1
  %x = rand(n,1); 
  x = rand(n,n);

  %y = fmincon(@sqrat,x,[],[],[],[],0*x,[]);
  %y = fminunc(@l2linf,x);
  y = fmincon(@l2linf,x,[],[],[],[],0*x,0*x+1);
  fprintf(' maxr = %1.9f \n',maxr);
end

return


function y = l2linf(x)
global n maxr

A = reshape(x,[n n]);
r = l2(A)/linf(A);

y = -r;

if r > maxr
  maxr = r;
  fprintf(' maxr = %1.9f \n',maxr);
end
return



function y = sqrat(x)
global optV optx

y = sum(x.^2)/sum(x)^2;

if y < optV
  optx = x;
  optV = y;
  fprintf(' optV = %1.9f \n',optV);
end
return