function P=transp(n,i,j)
%
% "Matriz de Transposición i-j"
%
% P = transp(n,i,j)
%
% Si A es una matriz de n renglones, PA es la matriz A con sus
% renglones i-esimo y j-esimo intercambiados.
% Si B es una matriz de n columnas, BP es la matriz B con sus
% columnas i-eisma y j-esima intercambiadas.
%
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% Revisión 1.0
% Sergio A. Pérez Ruiz
% Instituto Tecnológico Autónomo de México
% noviembre de 1996
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

if i<1 | j<1 | i>n | j>n
   disp('Indices incompatibles')
   return
end
P=eye(n);
row=P(i,:); P(i,:)=P(j,:); P(j,:)=row;
return