Step 1 corresponds to solve the linear system:
Step 2 is computed by the closed formula:
A continuation method over the parameter avoids to be “early” trapped by a local minima.
This is, we start with a quadratic function ( large) and we adjust the functional by systematically decreasing