Simple rational approximation (SRA) is a subset of interpolating methods using rational functions. Especially, SRA interpolates a given function with a specific rational function whose poles and zeros are simple, which means that there is no multiplicity in poles and zeros. Sometimes, it only implies simple poles. The main application of SRA lies in finding the zeros of secular functions. A divide-and-conquer algorithm to find the eigenvalues and eigenvectors for various kinds of matrices is well known in numerical analysis. In a strict sense, SRA implies a specific interpolation using simple rational functions as a part of the divide-and-conquer algorithm. Since such secular functions consist of a series of rational functions with simple poles, SRA is the best candidate to interpolate the zeros of the secular function. Moreover, based on previous researches, a simple zero that lies between two adjacent poles can be considerably well interpolated by using a two-dominant-pole rational function as an approximating function.
One-point third-order iterative method: Halley's formula The origin of the interpolation with rational functions can be found in the previous work done by Edmond Halley. Halley's formula is known as one-point third-order iterative method to solve f ( x ) = 0 {\displaystyle \,f(x)=0} by means of approximating a rational function defined by
h ( z ) = a z + b + c . {\displaystyle h(z)={\frac {a}{z+b}}+c.}
We can determine a, b, and c so that
h ( i ) ( x ) = f ( i ) ( x ) , i = 0 , 1 , 2. {\displaystyle h^{(i)}(x)=f^{(i)}(x),\qquad i=0,1,2.}
Then solving h ( z ) = 0 {\displaystyle \,h(z)=0} yields the iteration
x n + 1 = x n − f ( x n ) f ′ ( x n ) ( 1 1 − f ( x n ) f ″ ( x n ) 2 ( f ′ ( x n ) ) 2 ) . {\displaystyle x_{n+1}=x_{n}-{\frac {f(x_{n})}{f'(x_{n})}}\left({\frac {1}{1-{\frac {f(x_{n})f''(x_{n})}{2(f'(x_{n}))^{2}}}}}\right).}
This is referred to as Halley's formula. This geometrical interpretation h ( z ) {\displaystyle h(z)} was derived by Gander (1978), where the equivalent iteration also was derived by applying Newton's method to
g ( x ) = f ( x ) f ′ ( x ) = 0. {\displaystyle g(x)={\frac {f(x)}{\sqrt {f'(x)}}}=0.}
We call this algebraic interpretation g ( x ) {\displaystyle g(x)} of Halley's formula.
One-point second-order iterative method: Simple rational approximation Similarly, we can derive a variation of Halley's formula based on a one-point second-order iterative method to solve f ( x ) = α ( ≠ 0 ) {\displaystyle \,f(x)=\alpha (\neq 0)} using simple rational approximation by
h ( z ) = a z + b . {\displaystyle h(z)={\frac {a}{z+b}}.}
Then we need to evaluate
h ( i ) ( x ) = f ( i ) ( x ) , i = 0 , 1. {\displaystyle h^{(i)}(x)=f^{(i)}(x),\qquad i=0,1.}
… excerpt ends here. Continue reading the full article.
