using python without using any powerful functions besides import math and copy.

Write a function defined as: def GaussElim(Aaug):

Purpose of GaussElim(Aaug): use the Gauss Elimination with Back Substitution method to solve the set of linear equations expressed in matrix form as Ax =b. Both A and b are contained in the function argument Aaug.

Aaug: an augmented matrix containing [A | b ] having N rows and N+1 columns, where N is the number of equations in the set.

return value: x an array of size N, containing the solution to Ax = b

Notes: Do NOT use partial pivoting in this program. (I want everyone to get the exact same answer, and I dont want you working quite that hard).

Write and call a main() function that uses your GaussElim function to calculate and print the solution to the same sets of linear equations given.

MyA = [[4, -1, -1, 3], [-2, -3, 1, 9], [-1, 1, 7, -6]] using initial guesses of all zeros. Perform 22 iterations.

anotherA = [[4, 3, 1, -1, 2], [2, -5, 0, -2, -3], [-3, 3, -6, 1, 5], [0, 1, 4, 8, -2]] using initial guesses of all zeros. Perform 3 iterations.