function normval = estimate_spectral_norm(A,maxiters)
	m = size(A,1);
	n = size(A,2);
	
	v = randn(n,1);
	v = v/norm(v);

	for k=1:maxiters
		w = A'*(A*v);
		normval = norm(w)/norm(v);
		v = w/norm(w);
	end
	normval = sqrt(normval);
end
