varband_path function

Solve main optimization problem along a path of lambda