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!

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&#039;t understand how to find the one correct root and reject all others. Bisection, Newton&#039;s and other methods seem to only work when there is one root in an interval. Would appreciate any help!

xgirlrogueim · Accepted Answer

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 &quot;close enough&quot; to the fixed point for Newton&#039;s method to work in practice.

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

Answer & Explanation

New Questions in Integral Calculus