trying to code in python using Gauss Elimination with partial pivot. I cannot get the partial pivot to work. Here is my code;

#YOU CAN USE EITHER OF THESE PIVOT FUNCTIONS NOT SURE WHICH IS BETTER OR CREATE A NEW ONE

def pivot1(): for i in range(n-1): p = np.argmax(np.abs(a[i:,i])) p = p + i a[[i, p], :] = a[[p, i], : ]

def pivot2 (): for i in range(n-1,0,-1): if(a[i-1,0] for j in range(0,n): c=a[i,j] a[i,j]=a[i-1,j] a[i-1,j]=c # Importing NumPy Library import numpy as np

n = 4

a = [[1, 2, 1, -1, 5], [3, 2, 4, 4, 16], [4, 4, 3, 4, 22], [2,0,1,5,15]]

#a = a.astype(float) a = np.array(a, dtype=np.float32)

print (a)

for i in range(n): pivot() for j in range(i+1, n): ratio = a[j][i]/a[i][i] for k in range(n+1): a[j][k] = a[j][k] - ratio * a[i][k] print(a)

x = np.zeros(n) x[n-1] = a[n-1][n]/a[n-1][n-1]

for i in range(n-2,-1,-1): x[i] = a[i][n] for j in range(i+1,n): x[i] = x[i] - a[i][j]*x[j] x[i] = x[i]/a[i][i]

print(' The answer for X: ') for i in range(n): print('X%d = %0.2f' %(i,x[i]), end = '\t')