function [y] = GS(A)
y(:,1)=A(:,1)/norm(A(:,1));
Num=0;
[Ahang,Alie]=size(A);
for j=2:Alie
sum(:,1)=zeros(1,Ahang);
for i=1:j-1
h(i,j)=A(:,j)'*y(:,i);
sum(:,1) =sum(:,1)+h(i,j)*y(:,i);
Num=Num+1;
end
y(:,j)=A(:,j)-sum(:,1);
y(:,j)=y(:,j)/norm(y(:,j));
end
disp('Num=');
disp(Num)
