function zgeqrf import Complex; input Complex A[:, :] "Square or rectangular matrix"; output Complex Q[size(A, 1), size(A, 2)] "Square upper triangular matrix"; output Complex R[size(A, 2), size(A, 2)] "Square upper triangular matrix"; end zgeqrf;