On the Lagged Diffusivity Method for the solution of nonlinear finite difference systems