function Bs = pivotGJ(B,i,j)
%function Bs = pivotGJ(B,i,j)
% Renvoie dans Bs le resultat du pivotage de la methode Gauss-Jordan
% appliquee a la matrice B avec le pivot B(i,j).

Bs=B; % pour dimensionner Bs et surtout tenir compte de l''anteriorite
      % en recuperant les colonnes de B deja modifiees.

[m,n] = size(B);
for ic=1:n
 Bs(:,ic) = B(:,ic) - ( B(i,ic)/B(i,j) ) * B(:,j);
end
Bs(i,:) = B(i,:)/B(i,j);
Bs(:,j) = zeros(m,1);
Bs(i,j) = 1;
