function [Q R] = Alg13(A) end