Nonsingular Matrices
Our algorithm for calculating matrix inverses called for augmentation by
identity matrices. We can accomplish this with the
syntax A.augment(identity_matrix(A.nrows()), subdivide=True).
Of course, sage is capable of calculating inverses directly. The appropriate
syntax is A.inverse().