Recently, the nonlinearity continuation method has been used to numerically solve boundary value problems for steady-state Richards equation. The can be considered as a predictor-corrector procedure with simplest form which applied date having trivial, zeroth-order predictor. In this article, effect of more sophisticated predictor technique is examined. Numerical experiments are performed finit...