We will present a convenient, practical method of generating examples of matrices for which the Gaussian elimination process can be done in integer arithmetic, resulting in an echelon form with integer entries. The method can be refined so that the inverse has integer entries. The method is based on the L-U factorization, and we will also give an elementary justification for that factorization