#def LU(A):
n=len(A)
L,U=np.zeros([n,n]),np.zeros([n,n])

#for i in range (1):
i=0
Lt=A[:,i]/A[i,i]
Ut= A[i]
A=A-np.outer(Lt,Ut)
L[:,i]=Lt
U[i]=Ut
print (f"Matrix L is : \n {L}")
print (f"Matrix U is : \n{U}")
print (f"Matrix A is : \n{A}")