How calculate programs/packages like Matlab, Python/scipy, ...the inverse jacobian for multidimensional Newton's method? x_(n+1)=xn−(J(x_n)^(−1)∗f(x_n))

Vorbeckenuc

Vorbeckenuc

Answered question

2022-07-18

How calculate programs/packages like Matlab, Python/scipy, ...the inverse jacobian for multidimensional Newton's method?
x n + 1 = x n ( J ( x n ) 1 f ( x n )

Answer & Explanation

Vartavk

Vartavk

Beginner2022-07-19Added 11 answers

Most numerical packages give you the option of either computing the Jacobian yourself and passing it to the solver, or of numerically approximating it with a finite difference scheme. I imagine that in general while performing Newton's method or other methods expressed in terms of an inverse Jacobian, these packages do not actually compute the inverse for reasons of stability. Instead, they solve the linear system J ( x n ) x n + 1 = J ( x n ) x n f ( x n ) for x n + 1 at each time step.
The Jacobian is not always invertible; in order to use Newton's method the Jacobian must be invertible, though. You can see this in one dimension, with a function R R with a critical point which is not an extremum.

Do you have a similar question?

Recalculate according to your conditions!

New Questions in Analysis

Ask your question.
Get an expert answer.

Let our experts help you. Answer in as fast as 15 minutes.

Didn't find what you were looking for?