function v=BidiaSVD(A)
[m,n]=size(A);
C=zeros(m+n,m+n);
C(1:n,n+1:m+n)=A';
C(n+1:m+n,1:n)=A;
v=eig(C);