function [u,v] = lrsvecs(A,sigma)
		
[m,n] = size(A);

[v,s] = eigs(A'*A - sigma^2*eye(n), 1, 'sm');
v = v/norm(v);

u = A*v;
u = u/norm(u);
				
return;