
function W = reg_least_sq_path(X,y,Lambda,w_0)

m = size(X,2);
l = length(Lambda);
W = zeros(m,l);

w_t = w_0;
for t = [1:length(Lambda)]

  lambda_t = Lambda(t)
  w_t = reg_least_sq(X,y,lambda_t,w_t)
  W(:,t) = w_t;

end
