Metode Eliminasi Gauss-Jordan dengan Python 3.2

from numpy import*
#~~~~~~inisialisasi matrik augment~~~~~#
A = array([[1.,1.,0.,3.,4],
           [2.,1.,-1.,1.,1],
           [3.,-1.,-1.,2.,-3],
           [-1.,2.,3.,-1,4]])
print (‘Matriks A=’,A)
#===== METODE ELIMINASI GAUSS =========#
n=len(A)
#~~~~~~proses triangularisasi~~~~~~~~~~#
for k in range(0,n-1):
    for i in range(k+1,n):
        m=A[i][k]/A[k][k]
        for j in range(0,n+1):
            A[i][j]=A[i][j]-m*A[k][j]
#~~~~~~proses substitusi-mundur~~~~~~~~#
X = zeros((n,1))
X[n-1][0]=A[n-1][n]/A[n-1][n-1]
for j in range(n-2,-1,-1):
    S=0
    for i in range(j+1,n):
        S=S+A[j][i]*X[i][0]
        X[j][0]=(A[j][n]-S)/A[j][j]
#======================================#
print (‘X=’,X)
Exit mobile version