function [b, iflag] = lu_cp_sl( A, b, ipivtr, ipivtc );
% 
% LU_CP_SL  computes the solution of a linear system
% A x = b using the LU-decomposition with complete
% pivoting computed by LU_CP.
% 
% Usage:
%        [b, iflag] = lu_cp_sl( A, b, ipivtr, ipivtc )
% 
% input:
%        A:      the LU-decomposition of  A  computed by lu_cp
%
%        ipivtr: row permutation information for the LU-decomposition of  A  
%                computed by lu_cp.m
% 
%        ipivtc: column permutation  information for the LU-decomposition 
%                of  A  computed by lu_cp.m
% 
%        b:     the right hand side b. b can have more than one column. 
%               In this case a system with multiple right hand sides is solved.
%  
% output:
%        b:     the solution of the linear system if iflag = 0
%
%        iflag: error flag
%               iflag = 0  A is nonsingular, solution x is computed.
%               iflag = 1  dimension of A or of b is not correct
%               iflag > 1  zero diagonal element of U
%                          detected in row iflag+1
%  
%
% Your name
% Jan 10, 2008

%%%%% insert code



% end of lu_cp_sl