fixed point function (nonlinear equation)

71 Views Asked by At

here's the following problem, I'm trying to find a real root by fixed-point iteration method but I can't find a properly $g(x)$ that meets the condition that $|g'(x_0)|<1$. Well, my nonlinear equation is $f(x)=1.08^x-\frac{125}{81}=0$ so I add +x in both sides, $g(x)=1.08^x-\frac{125}{81}+x=x$ but when I do $g'(x) =d/dx(1.08^x - 125/81 + x) = 0.076961*1.08^x + 1$ and I try to solve inequality to know a initial value $x_0$ I got this after set this equation $g'(x)=0.076961*1.08^x + 1 < 1$: $Im(x) = (π (2 n + 1))/(3 log(3) - 2 log(5))$,$n \in Z $ . PS: My real root of $x \approx 5.637457293$.

Thanks in advance.

1

There are 1 best solutions below

0
On BEST ANSWER

Calling

$$ f(x) = a^x-b $$

we have

$$ f(x)= f(x_k) + f'(x_k)(x-x_k) + O(|x-x_k|^2) $$

so with $x$ near a root we have

$$ f(x)\approx 0 = f(x_k)+f'(x_k)(x-x_k) $$

thus

$$ x = x_{k+1} = x_k-\frac{f(x_k)}{f'(x_k)} $$

now calling $g(x) = x - \frac{f(x)}{f'(x)}$ we have

$$ x_{k+1} = g(x_k)\\ x_k = g(x_{k-1}) $$

and

$$ x_{k+1}-x_k = g(x_k)-g(x_{k-1}) $$

but

$$ g(x_k) = g(x_{k-1})+g'(x_{k-1})(x_k-x_{k-1})+O(|x_k-x_{k-1}|^2) $$

or

$$ x_{k+1}-x_k = g'(x_{k-1})(x_k-x_{k-1})+O(|x_k-x_{k-1}|^2) $$

now assuming small $O(|x_k-x_{k-1}|^2)$ we can assure convergence if $|g'(x)| < 1$ at the root vicinity.

In our case we have

$$ g(x) = x-\frac{a^x-b}{a^x\ln a} $$

and for convergence (sufficient) we need $|g'(x)| = |1-b a^{-x}|< 1$