HAHA VODOX

import numpy as np

ifile = open("in1.txt","r")

#file theke input

print("The number of the unknown variables is: ")
num = int(ifile.readline())
print(num)
A=[]
#print("Matrix A is: ")
for i in range(num):
    A_read = ifile.readline().rstrip().split(" ")
    A.append([float(x) for x in A_read ])


B=[]
#print("Matrix B is: ")
for i in range(num):
    B_read = ifile.readline().rstrip().split(" ")
    B.append([float(x) for x in B_read ])

def LUDecomposition(matrix,size):
    lower = [[0 for i in range(size)] for i in range(size)]
    upper = [[0 for i in range(size)] for i in range(size)]

    for i in range(size):

        for j in range(i,size):
            sumU = 0
            for k in range(i):
                sumU += (lower[i][k] * upper[k][j])
            upper[i][j] = matrix[i][j] - sumU
        for j in range(i,size):
            if (i == j):
                lower[i][i] = 1
            else:
                sumL = 0
                for k in range(i):
                    sumL += (lower[j][k] * upper[k][i])
                lower[j][i] = (matrix[j][i] - sumL)/upper[i][i]





    print(upper)
    print(lower)

LUDecomposition(A,num)