What will the following code snippet output for the LU factorization of matrix A?

import numpy as np

A = np.array([[2, 1], [4, 3]])

n = A.shape[0]

X = np.eye(n)

Y = A.copy()

for i in range(n):

    for j in range(i + 1, n):

        X[j, i] = Y[j, i] / Y[i, i]

        Y[j] = Y[j] - X[j, i] * Y[i]

print(X)

print(Y)