Function:     INVMAT

Purpose:      INVMAT is a routine for inverting a matrix. The algorithm used
              is the Gauss-Jordan algorithm described in Stoer, Numerische
              Matematik, 1 Teil.

Category:     MATH

File:         invmat.shl

Author:       K.G. Begeman

Use:          INTEGER INVMAT( MATRIX ,    Input/Output    REAL ARRAY
                              MATDIM ,       Input        INTEGER
                              DECDIM )       Input        INTEGER


              INVMAT       Returns 0 on success, 1 when matrix cannot be
                           inverted.
              MATRIX       On input contains matrix to be inverted,
                           on output the inverted matrix.
              MATDIM       Actual size of matrix.
                           Note: maximum value for MATDIM is 512.
              DECDIM       Size of MATRIX as declared in main program if
                           MATRIX was declared as a two dimensional array.
                           If MATRIX was declared one dimensional, DECDIM
                           should be equal to MATDIM.

Updates:      Jul  14, 1991: KGB, Document created.