Solve the problem that involves iteration

Click For Summary
SUMMARY

The discussion focuses on solving an iterative problem involving the function defined by the equation ##\dfrac{e^x}{4x^2-1}= -2##. Participants analyze the convergence of the iteration method, particularly around the asymptote at ##x=0.5##, and the failure of the sign-change method between ##x=0## and ##x=1## due to discontinuity. The second derivative test confirms that the fixed point ##α## cannot be found as it leads to negative logarithmic values. The conversation emphasizes the importance of understanding the limitations of numerical methods in the context of discontinuous functions.

PREREQUISITES
  • Understanding of iterative methods in numerical analysis
  • Familiarity with calculus concepts, including derivatives and asymptotes
  • Knowledge of the sign-change rule for root-finding
  • Experience with logarithmic functions and their properties
NEXT STEPS
  • Study the limitations of the sign-change rule in numerical methods
  • Explore Newton-Raphson method for faster convergence in root-finding
  • Investigate the behavior of functions with discontinuities and their impact on numerical solutions
  • Learn about fixed-point iteration and stability analysis in iterative methods
USEFUL FOR

Mathematicians, engineering students, and anyone involved in numerical analysis or solving equations using iterative methods will benefit from this discussion.

chwala
Gold Member
Messages
2,828
Reaction score
420
Homework Statement
see attached.
Relevant Equations
iterative techniques
1686668177209.png


part (a)

Asymptote at ##x=0.5##

part (b)

##\dfrac{e^x}{4x^2-1}= -2##

##e^x=2-8x^2##

##2e^x= 4-16x^2##

##16x^2=4-2e^x##

##x^2= \dfrac{4-2e^x}{16}##

##x=\dfrac{\sqrt{4-2e^x}}{4}##

part (c)

...

##x_{2}=0.2851##

##x_{3}=0.2894##

##x_{4}=0.2881##

##x_{5}=0.2885##

##x_{6}=0.2884##

##x_{7}=0.2884##

##α=0.2884##

For part(d)

not sure here, but i checked directly with

##F(x_n) = \ln (2-8x^2_n)## ...and noted that after a few iterations we were ending up with natural log of negative numbers thus ##α## cannot be found.

on a different approach, using second derivative,

##F^{'}(x) = \left[\dfrac{1}{2-8x^2} × -16x\right]##

##F^{'}(x) = \dfrac{8x}{4x^2-1}##

##F^{'}(0.3) = -3.75##

##F^{''}(x)=- \left[\dfrac{8+32x^2}{(4x^2-1)^2}\right]##

This will always be negative irrespective of ##x## value thus ##α## cannot be found.
 
Last edited:
Physics news on Phys.org
chwala said:
Homework Statement: see attached.
Relevant Equations: iterative techniques

View attachment 327810

part (a)

Asymptote at ##x=0.5##
This doesn't explain why the sign change method doesn't work for x = 0 and x = 1.
chwala said:
part (b)

##\dfrac{e^x}{4x^2-1}= -2##

##e^x=2-8x^2##

##2e^x= 4-16x^2##

##16x^2=4-2e^x##

##x^2= \dfrac{4-2e^x}{16}##

##x=\dfrac{\sqrt{4-2e^x}}{4}##
This looks to be correct, but with quite a few unnecessary steps. All you need to do is to get the ##x^2## term on one side and the term with ##e^x## on the other.
Also, there is a positive and a negative solution for x.
chwala said:
<snip>
Part c looks fine.
For part(d)

not sure here, but i checked directly with

##F(x_n) = \ln (2-8x^2_n)## ...and noted that after a few iterations we were ending up with ##\ln ## of negative numbers thus ##α## cannot be found.
...
##F^{'}(x) = \left[\dfrac{1}{2-8x^2} × -16x\right]##

##F^{'}(x) = \dfrac{8x}{4x^2-1}##

##F^{'}(0.3) = -3.75##

##F^{''}(x)=- \dfrac{8+32x^2}{(4x^2-1)^2}## will always be negative irrespective of ##x## value thus ##α## cannot be found.
 
Last edited:
For part d, I think there might be a typo in the problem statement.
"By considering F'(0.3)..." seems irrelevant to me, but if ##x_0 = 0.3## is ##x_1 = F(x_0)## in the domain of F?
 
Last edited:
chwala said:
on a different approach, using second derivative,
<snip>

This will always be negative irrespective of x value thus α cannot be found.
I don't see how this answers the question asked in part d. In my previous post I mentioned that the problem author might have made a typo, which in my view sent you off on a wild goose chase.
 
@Mark44 I am informed question is just okay (correct)... one needs to check the modulus of the derivative and establish the fact that there is no convergence to a particular value...unless that's not correct boss.
 
chwala said:
@Mark44 I am informed question is just okay (correct)
Who informed you of this? I don't see how F'(0.3) is relevant to the problem at all.
chwala said:
one needs to check the modulus of the derivative and establish the fact that there is no convergence to a particular value...
Based on what is given in the problem, the above seems to me to be just handwaving.
 
The iteration in (d) is the (local) inverse of the iteration in (c), so if one moves you towards the root, the other must move you away from it.

\alpha is an unstable fixed point of the iteration x_{n+1} = F(x_n) if |F&#039;(\alpha)| &gt; 1. Here we are asked to look at F&#039;(0.3), which is close to, but not at, the fixed point. However by continuity of F&#039; we might be able to conclude that |F&#039;(\alpha)| &gt; 1.
 
  • Informative
Likes   Reactions: chwala
More rigourously, if x_n = 0.3 + \xi_n then \xi_n satisfies <br /> \xi_{n+1} + 0.3 = F(0.3 + \xi_n) = F(0.3) + \xi_nF&#039;(0.3) + \dots which to leading order is of the form <br /> \xi_{n+1} = \alpha + \beta\xi_n which has solution <br /> \xi_n = \left(\xi_0 + \frac{\alpha}{\beta - 1}\right)\beta^n - \frac{\alpha}{\beta - 1}. We can see that if |\beta| &gt; 1 then \xi_n \to \infty and if |\beta| &lt; 1 then \xi_n \to -\alpha/(\beta - 1).
 
  • Informative
Likes   Reactions: Mark44 and chwala
Mark44 said:
This doesn't explain why the sign change method doesn't work for x = 0 and x = 1.
This looks to be correct, but with quite a few unnecessary steps. All you need to do is to get the ##x^2## term on one side and the term with ##e^x## on the other.
Also, there is a positive and a negative solution for x.
For part (a),
...sign change between ##f(0)## and ##f(1)## will imply that a root exists between ##x=0## and ##x=1##. In our case,

##f(0)=1##

and

##f(1)=2.91##

there is no sign change thus no root exists between the two points.
 
Last edited:
  • #10
chwala said:
For part (a),
...sign change between ##f(0)## and ##f(1)## will imply that a root exists between ##x=0## and ##x=1##. In our case,

##f(0)=1## and ##f(1)=2.91##

there is no sign change thus no root exists between the two points.
Look at the graph you posted as part of the problem statement! Clearly there is a root between x = 0 and x = 1. What you were supposed to do was to explain why the sign-change rule isn't working here, and not conclude that there is no root, which is completely at odds with the graph.
 
  • Like
Likes   Reactions: chwala
  • #11
The sign method will fail with a false negative if there are an even number of simple roots in the interval, or an odd number of simple roots and a vertical asymptote at which the function changes sign, or a double root. It will fail with a false positive if the interval contains an odd number of vertical asymptotes at which the function changes sign.
 
  • Like
Likes   Reactions: chwala
  • #12
pasmith said:
The sign method will fail with a false negative if there are an even number of simple roots in the interval, or an odd number of simple roots and a vertical asymptote at which the function changes sign, or a double root. It will fail with a false positive if the interval contains an odd number of vertical asymptotes at which the function changes sign.
@pasmith
thanks, ...is this the correct way to reason on this question? The term false negative, false positive...are relatively new to me. Thanks mate.
 
  • #13
chwala said:
For part (a),
...sign change between ##f(0)## and ##f(1)## will imply that a root exists between ##x=0## and ##x=1##. In our case,

##f(0)=1##

and

##f(1)=2.91##

there is no sign change thus no root exists between the two points.

Ok, would it be correct if i handle part (a) as follows;

...there is no sign change at ##f(0)## and ##f(1)## but a root does exist between ##x=0## and ##x=1## implying that the sign-rule change does not work for this problem.
 
  • #14
chwala said:
Ok, would it be correct if i handle part (a) as follows;

...there is no sign change at ##f(0)## and ##f(1)## but a root does exist between ##x=0## and ##x=1## implying that the sign-rule change does not work for this problem.
No, not correct. Again, they are asking why the sign-change rule doesn't work, not just to state that it doesn't work based on looking at the graph.

The textbook you're working from should state exactly what the sign-change rule is, as well as what limitations it has. Does it require that the function in question be a polynomial (in which case it is continuous everywhere) or does it merely require that the function be continuous?
 
  • #15
Mark44 said:
No, not correct. Again, they are asking why the sign-change rule doesn't work, not just to state that it doesn't work based on looking at the graph.

The textbook you're working from should state exactly what the sign-change rule is, as well as what limitations it has. Does it require that the function in question be a polynomial (in which case it is continuous everywhere) or does it merely require that the function be continuous?
The sign-change rule does not work because the function is not continous between ##x=0## and ##x=1##.

The rule is appropriate for continous functions in small intervals.
 
  • #16
chwala said:
The sign-change rule does not work because the function is not continous between ##x=0## and ##x=1##.

The rule is appropriate for continous functions in small intervals.
Now you're on the right track. Can you state the sign-change rule and any restrictions that are placed on its use?
 
  • #18
What restriction is given in the page at the link you included? Does the function in this thread satisfy this restriction?
 
  • #19
Mark44 said:
What restriction is given in the page at the link you included? Does the function in this thread satisfy this restriction?
Our function is discontinous...there is an asymptote and that is why the sign-change rule fails. I hope i am getting your question correctly.

The other scenario on where the sign-change rule fails is when the interval is too large giving a possibility of many roots occurring; in particular,
even number of roots may mean that roots are missed entirely or odd number of roots may mean that not all roots are identified.
 
  • #20
That's what I was looking for.
 
  • Like
Likes   Reactions: chwala
  • #21
For part (c) my solution was not correct, i should have taken more significant figures...

part (c)

...

##x_{2}=0.28507##

##x_{3}=0.28943##

##x_{4}=0.28817##

##x_{5}=0.28853##

##x_{6}=0.28843##

##x_{7}=0.28846##

##x_{8}=0.28845##
##α=0.2885##
 

Similar threads

  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
10
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
Replies
4
Views
2K
  • · Replies 6 ·
Replies
6
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
1
Views
2K
Replies
4
Views
2K