function InvA = InvLU_1(A)
% Inversa de A a partir de su fact. LU
%
 [m n]=size(A); InvA=zeros(m); b=zeros(m,1);
 [L U P] = lu(A);
 for i=1:m
    b(i)      = 1;
    InvA(:,i) = U\(L\(P*b));
    b(i)      = 0;
 end
end
