Also, can you explain what each line of the code does?

Write a user-defined MATLAB function that solves a system of n linear equations, [a] [x] = [b], with the Gauss-Jordan method. the program should include pivoting in which the pivot row is switched with the row that has a pivot element with the largest absolute numerical value. For the function name and arguments use x = Gauss Jordan (a, b), where a is the matrix of coefficients, b is the right-hand-side column of constants, and x is the solution