This subroutine will perform Gaussian elmination on the matrix that you pass to it. i.e., given A and b it can be used to find x, Ax = b