mod_grams.m

various matlab programs to slove various problems

function [Q,R] = mod_grams(A)
[m,n]=size(A);
Q=zeros(m,n);   Q(1:m,1) = A(1:m,1);   R=zeros(n);   R(1,1)=1;
for k = 1:n
 R(k,k) = norm (A(1:m,k));    Q(1:m,k) = A(1:m,k)/R(k,k);
 for j=k+1:n
    R (k,j) = Q (1:m,k)' * A(1:m,j);
    A (1:m,j) = A (1:m,j) - Q(1:m,k)*R(k,j);
 end
end