matlab -- exmple of Gram-Schmidt process

find the orthoggonal set of w1,w2,w3 , then save to v( , );


w1=[1 0 1]';
w2=[0 1 1]';
w3=[1 3 3]';

i=3;

w=zeros(i,i);
v=zeros(i,i);

w(:,1)=w1
w(:,2)=w2
w(:,3)=w3


v(:,1)=w1;

for j=2:i,
v(:,j)=w(:,j);
for k=1:j-1,
wdv=w(:,j)' * v(:,k)
v2=v(:,k)' * v(:,k)
v(:,j)=v(:,j)- vdw/v2 *v(:,k);
end
end