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