function zungrq import Complex; input Complex RQ[:, :] "Square or rectangular matrix"; input Complex tau[:] "elementary reflectors"; output Complex Q[size(RQ, 1), size(RQ, 2)] "matrix Q"; end zungrq;