% -x    -2L <= 0
% -x -y -2L <= 0
%    -y -2L <= 0
%        -L <= 0
A = [ -1, -1,  0,  0;
       0, -1, -1,  0;
      -2, -2, -2, -1 ];

% -x    -2L <= 0
%  x +y -2L <= 0
%    -y -2L <= 0
%        -L <= 0
%A = [ -1,  1,  0,  0;
%       0,  1, -1,  0;
%      -2, -2, -2, -1 ];

% -x        <= 0
%  x +y +2L <= 0
%    -y     <= 0
%        -L <= 0
%A = [ -1,  1,  0,  0;
%       0,  1, -1,  0;
%       0,  2,  0, -1 ];

[U,E,V] = PCC_convert(A);
U

for k=1:size(U,2)
    if abs(U(3,k)) >= eps
        U(:,k) = U(:,k)/abs(U(3,k));
    end
end
U

