clear;
clc;
A=[1 0 0 0 
   1 1 0 0 
   1 1 1 0
   1 1 1 1];
V(:,1)=[1,0,0,0]';
H=zeros(5,4);
for i=1:4
    w=A*V(:,i);
    for j=1:i
        H(j,i)=V(:,j)'*A*V(:,i);
        w=w-H(j,i)*V(:,j);
    end
    H(i+1,i)=norm(w);
    V(:,i+1)=w/H(i+1,i);
end
V=V(:,1:4)
H=H(1:4,:)