function [A, ipivtr, ipivtc, iflag] = lu_pc( A );
% 
% LU_CP  computes the LU--decomposition with complete
% pivoting of a matrix A 
% 
% Usage
%       [A, ipivt, iflag] = lu_pp( A )
% 
% input:
%        A:     the n by n matrix A
% 
% output:
%        A:      the LU-decomposition of  A
% 
%        ipivtr: information on row exchanges
%                ipivtr(k) = i means that in step k 
%                (entries k:n of) rows k and i were interchanged
%
%        ipivtc: information on column exchanges
%                ipivtc(k) = i means that in step k columns
%                k and i were interchanged
% 
%        iflag: error flag
%               iflag = 0  Row reductions could be performed,
%                          A is upper triangular
%               iflag = 1  dimension of A or of b is not correct
% 
%
% Your name
% Jan 10, 2008

%%%%% insert code


%end of lu_cp.m