function gmin = condmin(Gf)
% Minimize cond.no. with both input and output scaling
% WARNING: May be inaccurate because of inaccurate mu-calculation
%
[dum1,n,m,dum2]=minfo(Gf);
zero = 0*eye(n);
H = abv( sbs(zero, minv(Gf)), sbs(Gf, zero));
blk0=[1 1]; blk1=blk0;
for k=2:n
  blk1 = [blk1; blk0];
end
blk2=blk1; blk=[blk1; blk2];
mu1=mu(H,blk,'UC9'); gmin1=sel(mu1,':',1);
gmin=mmult(gmin1,gmin1);
