I have a differntial equation du/dx=f(x,u), where f(x,u)=d psi/dx+p(u−psi(x))^4, where psi is a given function of x independent of u.



Answered question


I have a differntial equation d u / d x = f ( x , u ), where
f ( x , u ) = d ψ d x + p ( u ψ ( x ) ) 4 ,
where ψ is a given function of x independent of u.
The backward Euler Method written for step i is
u i + 1 = u i + h f ( x i + 1 , u i + 1 ) .
So, at every step we have a 4th degree polynomial to solve for u i + 1 , which can have more than one root. I am told to employ a root-finding algorithm. But I don't understand how to find the one correct root and reject all others. Bisection, Newton's and other methods seem to only work when there is one root in an interval. Would appreciate any help!

Answer & Explanation



Beginner2022-09-07Added 13 answers

If h is small and f is Lipschitz you will have that u i + 1 is close to u i . Hence, your first guess for the Newton algorithm ( u i itself) will be always "close enough" to the fixed point for Newton's method to work in practice.

Do you have a similar question?

Recalculate according to your conditions!

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?