function s=SVDFranIt(A);
[n,m]=size(A);
if n>m
    B=BiDiag(A);
    Bhat=B(1:m,1:m);
else
    A=A';
    B=BiDiag(A);
    Bhat=B(1:n,1:n);
end
s=svdFrancis(Bhat);
end