d = diag(A) ;
a = abs(d) ;
f= sum(abs(A), 2) - a ;
i = find(f >= a) ;
A(i,i) = A(i,i) + diag(f(i)) ;
