function [Mp,Vbp,Vnp] = pivotGJTR(M,Vb,Vn,i,j)
%function [Mp,Vbp,Vnp] = pivotGJTR(M,Vb,Vn,i,j)
% Renvoie dans [Mp,Vbp,Vnp] le resultat du pivotage de la methode du simplexe
% avec le pivot M(i,j), sous sa forme tableau reduit.

%
% echange des variables associees au pivot
Vbp = Vb;
Vnp = Vn;
Vbp(i-1) = Vn(j-1);
Vnp(j-1) = Vb(i-1);
%
% pivotage
Mp = M;
%
il = [1:i-1, i+1:size(M,1)];
for ic=1:size(M,2)
 Mp(il,ic) = M(il,ic) - ( M(i,ic)/M(i,j) ) * M(il,j);
end
Mp(i,:)  =   M(i,:)/M(i,j);
Mp(il,j) = -M(il,j)/M(i,j);
Mp(i,j)  = 1/M(i,j);
