Gauss.py

import sys
print("->>[̲̅J̲̅.̲̅L̲̅υ̲̅c̲̅α̲̅s̲̅™̲̅]<<-\n")
matriz=[]
termos_ind=[]
n=int(input("Digite a ordem do sistema: "))
resultado=[0]*n
for i in range(n):
    vet=[]
    for j in range(n):
        vet.append(float(input("digite o coeficiente da icógnita "+str(i+1)+str(j+1)+": ")))
    matriz.append(vet)
    
for i in range(n):
    termos_ind.append(float(input("digite o "+str(i+1)+"º valor do vetor dos termos independentes: ")))
for i in range(n):
    matriz[i].append(termos_ind[i])
res=int(input("Escolha um tipo de substituição: \n1-progressiva(Triangular inferior) 2-retroativa(Triangular superior): "))

if res!=1 and res!=2:
    while(res!=1 and res!=2):
        res=int(input("Escolha um tipo de substituição: \n1-progressiva(Triangular inferior) 2-retroativa(Triangular superior): "))
        
if res==2:
    for i in range(n-1):#n-1 porque ta comparando de duas em duas listas/ zerando a linha abaixo da diag principal
        ciclo=n
        pivo=matriz[i][i]
        while(pivo==0):
            matriz.append(matriz[i])
            del(matriz[i])#mudança de linha joga a linha do pivo para ultima linha
            pivo=matriz[i][i]
            ciclo=ciclo-1
            if pivo!=0: break
            if ciclo==0: #coluna toda zerada
                for xa in matriz:
                    del(xa[i]) #deleta coluna zerada
                print("Sistema incompatível")
                sys.exit(0)
                break       
        for j in range(i+1,n):
            coeficiente=matriz[j][i]/pivo
            for k in range(n+1):
                matriz[j][k]=matriz[j][k]-coeficiente*matriz[i][k]
        for i in matriz:
            for j in i:
                print(j,end=" | ")#matriz com duas casas decimais arredondadas
            print()
        print("\n\n")
    termos_ind2=[]
    for i in range(n):
        termos_ind2.append(matriz[i][n])
        del(matriz[i][n])
    resultado[n-1]=termos_ind2[n-1]/matriz[n-1][n-1]
    for i in range(n-1,-1,-1):
        soma=0
        for j in range(1+i,n):
            soma+=matriz[i][j]*resultado[j]
        if matriz[i][i]==0:
            print("Sistema incompatível")
            sys.exit(0)
        else:
            resultado[i]=(termos_ind2[i]-soma)/matriz[i][i]
    print("\n")
    print("Vetor solução:")
    for i in resultado:
        print("|"+str("%f"%i)+"|")
        
else:
    for i in range(n-1,0,-1): #zerando a linha acima da diag principal
        pivo=matriz[i][i]
        if pivo==0:
            print("Pivo não pode ser zero!")
            sys.exit(0)
        else:
            for j in range(i-1,-1,-1):
                coeficiente=matriz[j][i]/pivo
                for k in range(n+1):
                    matriz[j][k]=matriz[j][k]-coeficiente*matriz[i][k]
                for t in matriz:
                    for u in t:
                        print("%f"%u,end=" | ")
                    print("\n")
                print("\n\n")
    print("\n")
    print("Vetor solução:")
    for i in range(n):
        print("|"+str("%f"%(matriz[i][n]/matriz[i][i]))+"|")