function alpha_best = line_search_one(xs,solvec,num_iters,ts,ys);
	alpha = 1;
	factor = 0.95;
	resvec = build_residual(xs,ts,ys);
	
	alphas = zeros(num_iters,1);
	residual_norms = zeros(num_iters,1);
	alphas(1) = alpha;
	for i=1:num_iters
		if i>1
			alpha = factor*alpha;
		end
		alphas(i) = alpha;
	
    	xnew = xs - alpha*solvec;
		resvec = build_residual(xnew,ts,ys);
		residual_norms(i) = norm(resvec);
	end

	[min_val, min_ind] = min(residual_norms);
	alpha_best = alphas(min_ind);	
end
