global A b c m n z nonbas bas myeps
 
%Define A, b, c, z here.
 
[m,n]=size(A);
nonbas = [1:n-m];
%vars = horzcat([1:n-m],zeros(1,n-n+m));
bas = [n-m+1:n];
myeps = .1^8;

