For solving sets of linear equations, Gauss-Jordan elimination produces both the solution of the equations for one or more right-hand side vectors $b$, and also the matrix inverse $A^{−1}$.

The subroutine used here is adapted from Numerical Recipies. You can find the subroutine and relevant files in the code folder.

Example

$$x + y + 2 = 5$$

$$2x + 3y + 5z = 8$$

$$4x + 5z = 2$$

This can be written in agumented matrix form as follows:

The two matrices can be then use as input in to the gauss-jordan subroutine to find the solution for $x$, $y$ and $z$.

• lieq.f - Main Routine to input the matrices
• gauss-jordan.f - Gauss-Jordan Elimination Subroutine