The idea of the method is based on constructing a correction functional
by a general Lagrange multiplier and the multiplier is chosen in such a way that its correction solution is improved with
respect to the initial approximation or to the trial function.