Metode Gauss Eliminasi menggunakan Python

#Gauss Eliminasi
import numpy as np
def gaussEliminasi(a,b):
n = len(b)
# Eliminasi fase
for k in range (0,n-1):
                for i in range(k+1,n):
                        if a[i,k]!=0.0:
                                lam = a[i,k]/a[k,k]
                                a[i,k+1:n] = a[i,k+1:n] - lam*a[k,k+1:n]
                                b[i] = b[i]-lam*b[k]
#back substitusion
for k in range (n-1,-1,-1):
b[k] = (b[k] - np.dot(a[k,k+1:n],b[k+1:n]))/a[k,k]
return b

a = np.array([[30.0 -20.0, 0.01],\
     [-20.0, 30.0, -10.0],\
     [0.0, -10.0, -10.0]])
b = np.array ([19.6, 29.4, 24.5])

aOrig = a.copy()
#bOrig = b.copy()
x = gaussEliminasi(a,b)
#det = np.proud(np.diagonal(a))

print('x =\n',x)
#print('\ndet =',det)
print('\ncheck: A*x =\n',np.dot(aOrig,x))
input('\nPress return to exit')


Hasil eksekusi