Gaussian Elimination Using Python

The goal of Gaussian elimination is to convert a matrix to its row echelon form, $ref$. We discussed the steps of the algorithm in detail in the previous lesson. In this lesson, we’ll translate those steps into a Python program.

Implementation

The code below is an implementation of the Gaussian elimination algorithm in Python. The Elimination function takes a two-dimensional np.array as input. It performs the elimination while printing each step and returns the $ref$ of the input array as output.