The Help page of
FindRoot says: "by default,
FindRoot uses Newton's method (Newton-Raphson) to solve a nonlinear system". But I still have the following puzzles:
How will the Jacobian be calculated in each iteration: symbolically or using finite difference approximation only?
Is there any damping method of the step being used when the Jacobian at any step is invertible?