Implementing an arc-length method for a robust approach in solving systems of nonlinear equations