Write a program in either Python or Julia to find the QR factorization of an m×n matrix. The program should take both m and n, as well as the entries of the matrix A, as inputs.
The output of the program should be an m×n matrix Q and an n×n matrix R such that A = QR. Q should be a matrix which has the property that Q'Q = I and R should be triangular.