%script pivgrow
for n=1:10
    m=2^n
    A=-ones(m);
    A=eye(m)+tril(A,-1);
    A(:,m)=1;
    [L,U,piv]=GEpiv(A);
    norm(L*U-A)
end
