function [gama,L,R,B]=mincond_mu(A)
% gama = min cond(B), where B = L A R, L and R scaling diagonal matrices

[n,m]=size(A);
if n ~= m, error('Only square matriz!'), end

H=[zeros(n) inv(A);A zeros(n)];
[ub,D] = mu(H, ones(2*n,2), 'C');
L = diag(D(n+1:2*n));
R = inv(diag(D(1:n)));
gama = (ub(1)/ub(2))^2;
B = L*A*R;
