Bilateral approximations of the roots of scalar equations by Lagrange-Aitken-Steffensen method of order three

Abstract

We study the monotone convergence of two general methods of Aitken-Steffenssen type. These methods are obtained from the Lagrange inverse interpolation polynomial of degree two, having controlled nodes. The obtained results provide information on controlling the errors at each iteration step.

Author

Ion Păvăloiu
(Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy)

Keywords

Aitken-Steffenssen methods; Lagrange inverse interpolation

PDF

PDF-LaTeX file (on the journal website).

Cite this paper as:

I. Păvăloiu, Bilateral approximations of the roots of scalar equations by Lagrange-Aitken-Steffensen method of order three, Rev. Anal. Numér. Théor. Approx., 35 (2006) no. 2, pp. 173-182. https://doi.org/10.33993/jnaat352-843

About this paper

Print ISSN

1222-9024

Online ISSN

2457-8126

References

[1] Balazs, M., A bilateral approximating method for finding the real roots of real equations, Rev. Anal. Numer. Theor. Approx., 21 (2), pp. 111-117, 1992.

[2] Casulli, V., Trigiante, D.,  The convergence order for iterative multipoint procedures, Calcolo, 13 (1), pp. 25-44, 1997.

[3] Costabile, F., Gualtieri, I.M., Luceri, R., A new iterative method for the computation of the solution of nonlinear equations, Numer. Algorithms, 28, pp. 87-100, 2001.

[4] Frontini, M., Hermite interpolation and a new iterative method for the computation of the roots of non-linear equations, Calcolo, 40, pp. 109-119, 2003.

[5] Grau, M., An improvement to the computing of nonlinear equation solutions, Numer. Algorithms., 34, pp. 1-12, 2003.

[6] Ostrowski, A., Solution of Equations in Euclidian and Banach Spaces, Academic Press, New York and London, 1973.

[7] Păvăloiu I., Optimal efficiency index for iterative methods of interpolatory type, Computer Science Journal of Moldova, 1 (5), pp. 20-43, 1997.

[8] Păvăloiu I.,  Approximation of the roots of equations by Aitken-Steffensen-type monotonic sequences, Calcolo, 32 (1-2), pp. 69-82, 1995.

[9] Păvăloiu I., Optimal problems concerning interpolation methods of solution of equations, Publications de L’Institut Mathematique, 52 (66), pp. 113-126, 1992.

[10] Păvăloiu I., Optimal effiency index of a class of Hermite iterative methods, with two steps, Rev. Anal. Numer. Theor. Approx., 29 (1), pp. 83-89, 2000.

[11] Păvăloiu I., Local convergence of general Steffensen type methods, Rev. Anal. Numer. Theor. Approx., 33 (1), pp. 79-86, 2004.

[12] Păvăloiu I.,  and Pop, N., Interpolation and Applications, Risoprint, Cluj-Napoca, 2005 (in Romanian).

[13] Păvăloiu I., On a Steffensen-Hermite-type Method for approximating the solution of nonlinear equations, Rev. Anal. Numer. Theor. Approx., 25 1, pp. 87-94, 2006.

[14] Păvăloiu I.,  Bilateral approximation of solutions of equations by order-three Steffensen type methods, Studia Univ. “Babeș-Bolyai”, Mathematica, Vol. LI, no. 3, pp. 105-114, 2006.

[15] Traub, J.F., Iterative Methods for Solutions of Equations, Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1964.

[16] Turowicz, B. A., Sur les derivees d’ordre superieur d’une function inverse, Ann. Polon. Math., 8, pp. 265-269, 1960.

Paper (preprint) in HTML form

BILATERAL APPROXIMATIONS OF THE ROOTS OF SCALAR EQUATIONS BY LAGRANGE-AITKEN-STEFFENSEN METHOD OF ORDER THREE∗

BILATERAL APPROXIMATIONS OF THE ROOTS OF SCALAR EQUATIONS BY LAGRANGE-AITKEN-STEFFENSEN METHOD
OF ORDER THREE

Ion Păvăloiu
(Date: January 11, 2006.)
Abstract.

We study the monotone convergence of two general methods of Aitken-Steffenssen type. These methods are obtained from the Lagrange inverse interpolation polynomial of degree two, having controlled nodes. The obtained results provide information on controlling the errors at each iteration step.

Key words and phrases:
Aitken-Steffenssen methods, Lagrange inverse interpolation.
1991 Mathematics Subject Classification:
65H05.
This work has been supported by MEdC-ANCS under grant 2CEEX-06-11-96.
“Tiberiu Popoviciu” Institute of Numerical Analysis, P.O. Box. 68-1, Cluj-Napoca,Romania, e-mail: pavaloiu@ictp.acad.ro.

1. Introduction

It is well known that the Steffensen, Aitken, and Aitken-Steffensen methods are obtained from the inverse Lagrange interpolation polynomial of degree one, with controlled nodes [9][13], [7]. Consider the equation

(1.1) f(x)=0

where

f:[a,b]R,a,bR,a<b.

We also consider the following three equations, each of them equivalent with equation (1.1):

(1.2) x=g(x),g:[a,b][a,b]

and

(1.3) x = g1(x),g1:[a,b][a,b],
x = g2(x),g2:[a,b][a,b].

The Steffensen method in given by relations

(1.4) xn+1=xnf(xn)[xn,g(xn);f],x0[a,b],n=0,1,2,,

and, analogously, the Aitken method is of the following form:

(1.5) xn+1=g1(xn)f(g1(xn))[g1(xn),g2(xn);f],x0[a,b],n=0,1,2,.

Finally, the Aitken-Steffensen method is given by the relations:

(1.6) xn+1=g1(xn)f(g1(xn))[g1(xn),g2(g1(xn));f],x0[a,b],n=0,1,2,.

The order of convergence of all three methods (1.4)–(1.6) is at least two, and this order depends on the functions g and g1, g2 respectively. Essentially, the methods (1.4)–(1.6) are obtained from the method of chord where the interpolation nodes depend on the functions g and respectively g1,g2. In papers [2], [9], [10], [13] some conditions had to be considered in order that all the three methods (1.4)–(1.6) generate two sequences (un)n0 and (vn)n0 with the properties:

α) The sequence (un)n0 is increasing and the sequence (vn)n1 is decreasing:

β)limun=limvn=x¯, where x¯ is the root of equation (1.1), x¯ [a,b].

Practically, such sequences are very interesting, because by inequalities

max{x¯un,vnx¯}vnun,n=0,1,2,,

the errors of approximation at every step of iteration may be controlled.

Let a1,a2, a3[a,b] be three nodes of interpolation, and let b1, b2, b3 the values of the function f, i.e.

b1=f(a1),b2=f(a2),b3=f(a3).

Suppose that the function f:[a,b]F, is bijective where F=f([a,b]).

Then there exists f1:F[a,b] and the following equality holds [13], [15]:

f1(y) = a1+[b1,b2;f1](yb1)
+[b1,b2,b3;f1](yb1)(yb2)
+[y,b1,b2,b3;f1](yb1)(yb2)(yb3),

for every yF.

If x¯[a,b] is the root of equation (1.1), then x¯=f1(0), and by (1) one obtains the following representation of x¯:

x¯ = a1[b1,b2;f1]b1+[b1,b2,b3;f1]b1b2
[0,b1,b2,b3;f1]b1b2b3.

By relations (see [13])

[b1,b2;f1]=1[a1,a2;f]

and

[b1,b2,b3;f1]=[a1,a2,a3;f][a1,a2;f][a1,a3;f][a2,a3;f],

using (1), one obtains the following approximation for x¯:

(1.9) a4=a1f(a1)[a1,a2;f][a1,a2,a3;f]f(a1)f(a2)[a1,a2;f][a1,a3;f][a2,a3;f]

and the error:

(1.10) x¯a4=[0,b1,b2,b3;f1]f(a1)f(a2)f(a3).

If fC3([a,b]) and f(x)0, x[a,b], then f1C3(F) and the following equality holds (see [13], [17]):

(1.11) [f1(y)]′′′=3[f′′(x)]2f(x)f′′′(x)[f(x)]5,

where y=f(x).

Using the mean value formula for divided differences (see [13]) it follows that there exists ηF such that

(1.12) [0,b1,b2,b3;f1]=(f1)′′′(η)6.

Because f is bijective, it follows that there exists ξ[a,b] such that η=f(ξ), and by (1.11) and (1.12) one obtains:

(1.13) [0,b1,b2,b3;f1]=3[f′′(ξ)]2f(ξ)f′′′(ξ)6[f(ξ)]5.

By (1.9), if one considers particular nodes a1,a2,a3 it is possible to obtain different methods of Steffensen type, of Aitken type or of Aitken-Steffensen type.

Let xn[a,b] be an approximation of the root x¯ of equation (1.1).

If one considers a1=xn, a2=g(xn), a3=g(g(xn)), then it follows the following method of Steffensen type [9], [10], [13]:

xn+1 = xnf(xn)[xn,g(xn);f]
[xn,g(xn),g(g(xn));f]f(xn)f(g(xn))[xn,g(xn);f][xn,g(g(xn));f][g(xn),g(g(xn));f],

x0[a,b],n=0,1,2,.

If a1=xn, a2=g1(xn), a3=g2(g1(xn)), then one obtains the following method of Aitken-Steffensen type:

xn+1 = xnf(xn)(xn,g1(xn);f]
[xn,g1(xn),g2(g1(xn));f]f(xn)f(g1(xn))[xn,g1(xn);f][xn,g2(g1(xn));f][g1(xn),g2(g1(xn));f],

n=0,1,, x0[a,b], and finally, for a1=xn,a2=g1(xn), a3=g2(xn) one obtains the following method of Aitken type:

xn+1 = xnf(xn)[xn,g1(xn);f]
[xn,g1(xn),g2(xn);f]f(xn)f(g1(xn))[xn,g1(xn);f][xn,g2(xn);f][g1(xn),g2(xn);f],
n = 0,1,,x0[a,b].

Using the symmetry of Lagrange polynomial with respect to nodes, by permutations of a1, a2, a3 in methods (1)–(1), one obtains the same results for xn+1.

In [15] conditions are given in order that method (1) generates sequences approximating the root of equation (1.1) bilaterally. In this paper we study methods (1) and (1) and we obtain same conditions in order that the sequences generated by there methods are bilateral approximations of the root x¯ of equations (1.1).

2. The convergence of Aitken-Steffensen method

In the following we study the method (1) and we search the conditions on the sequences (xn)n0, (g1(xn))n0 and (g2(g1(xn))n0 generated from this method in order that they are monotonic sequences, bilaterally approximating the root x¯ of equation (1.1).

We need the following hypothesis:

  1. a)

    g1 is increasing on [a,b];

  2. b)

    g2 is continuous and decreasing on [a,b];

  3. c)

    equation (1.1) has a solution x¯[a,b] and g1(x¯)=g2(x¯)=x¯,

  4. d)

    function f is in C3([a,b]), and for every x[a,b] the following relation is fulfilled:

    (2.1) 3[f′′(x)]2f(x)f′′′(x)<0;
  5. e)

    function g1 satisfies inequality

    |g1(x)g1(x¯)|L|xx¯|,x[a,b],

    where 0<L<1.

Concerning the convergence of sequence (xn)n0 generated by (1), the following Theorem holds:

Theorem 2.1.

Let x0[a,b] and f,g1,g2 verify the following conditions:

  1. i)1

    f is increasing on [a,b];

  2. ii)1

    f is convex on [a,b];

  3. iii)1

    functions g1,g2 and f verify hypotheses a)e);

  4. iv)1

    x0>x¯ and g2(g1(x0))a.

Then the sequences (xn)n0, (g1(xn))n0 and (g2(g1(xn))n0 generated by (1) verify the properties:

  1. j)1

    sequences (xn)n0 and (g1(xn))n0 are decreasing and bounded from below by x¯ ;

  2. jj)1

    sequence (g2(g1(xn))n0 is increasing and bounded from above by x¯;

  3. jjj)1

    at every iteration step the following inequalities hold:

    (2.2) xnx¯xng2(g1(xn)),n=0,1,;
  4. jv)1

    limxn=limg1(xn)=limg2(g1(xn)=x¯.

Proof.

By hypotheses a) and x0>x¯ it follows g1(x0)>g1(x¯), i.e g1(x0)>x¯. Hypothesis e) implies g1(x0)<x0. From 𝐢𝐯1) it follows g2(g1(x0)a and from b) and g1(x0)>x¯ one obtains g2(g1(x1))<x¯. Now using 𝐢1) it follows:

f(x0)>0,f(g1(x0))>0

and

f(g2(g1(x0))<0,

and by considering 𝐢1) and 𝐢𝐢1) for n=0 in (1) it follows that x1<x0. Using (1.13), by hypotheses d) and (1.10) for a4=x1 and b1=f(x0),b2=f(g1(x0)),b3=f(g2(g1(x0)) it follows x¯x1<0, i.e x1>x¯. By hypotheses a) and x1<x0 it follows g1(x1)<g1(x0) and then, by b) one obtains g2(g1(x1))g2(g1(x0)), and g2(g1(x1))<x¯. Let xm, m be an element of sequence (xn)n0 generated by (1) and suppose that xm>x¯. Then one obtains the relations:

(2.3) g2(g1(xm))<g2(g1(xm+1))<x¯<g1(xm+1)<xm+1<g1(xm)<xm.

Inequality xm+1<g1(xm) follows from equality:

xm+1 = xmf(xm)[xm,g1(xm);f]
[xm,g1(xm),g2(g1(xm));f]f(xm)f(g1(xm)[xm,g1(xm);f][xm,g2(g1(xm));f][g1(xm),g2(g1(xm));f]
= g1(xm)f(g1(xm)[xm,g1(xm);f]
[xm,g1(xm),g2(g1(xm));f]f(g1(xm))f(xm)[xm,g1(xm);f][xm,g2(g1(xm));f][g1(xm),g2(g1(xm));f]

and from the hypothesis of the theorem.

From relations (2.3) it follows (2.2) and conclusions 𝐣1) and 𝐣𝐣1). Conclusion 𝐣𝐯1) is obvious. The theorem is proved. ∎

Remark 2.2.

If function f is concave and decreasing on [a,b], then function h:[a,b] defined by h(x)=f(x) in convex and increasing. The relation (2.1) is also verified for h. Consequently the sequence generated by (1) for function h, verifies all the conditions of Theorem 2.1, and then the conclusions of this theorem hold.

An analogous proof with that from Theorem 2.1 is valid for the following:

Theorem 2.3.

Let x0[a,b] and f,g1,g2 verify the following conditions:

  1. i)2

    function f is increasing on [a,b];

  2. ii)2

    function f is concave on [a,b];

  3. iii)2

    functions g1,g2 and f verify the hypotheses a)e);

  4. iv)2

    x0<x¯ and g2(g1(x0)b.

Then sequences (xn)n0, (g1(xn))n0 and (g2(g1(xn0))) generated by (1) have the following properties:

  1. j)2

    sequences (xn)n0 and (g1(xn))n0 are increasing and bounded from above by x¯;

  2. jj)2

    sequence g2(g1(xn)) is decreasing and bounded from below by x¯;

  3. jjj)2

    the following relations hold:

    (2.4) x¯xng2(g1(xn))xn,n=0,1,2,,
  4. jv)2

    limxn=limg1(xn)=limg2(g1(xn))=x¯.

Remark 2.4.

If function f is decreasing and convex, then h1:[a,b] defined by h1(x)=f(x) is increasing and concave, h1verifies all hypotheses of Theorem 2.3 and consequently all the conclusions 𝐣2)𝐣𝐯2) hold.

3. Convergence of Aitken type method

In order to study the convergence of the sequences generated by (1) we must suppose that functions f,g1 and g2 verify hypotheses a)-e) of section 2.

The following theorem holds:

Theorem 3.1.

Let x0[a,b] and the functions f,g1 and g2 verify the conditions:

  1. i)3

    f is increasing an [a,b];

  2. ii)3

    f is convex an [a,b];

  3. iii)3

    f, g1 and g2 verifies the hypotheses a)e);

  4. iv)3

    x0>x¯ and g2(x0)>a.

Then sequences (xn)n0, (g1(xn))n0 and (g2(xn))n0 generated by (1) have the properties:

  1. j)3

    sequences (xn)n0 and (g1(xn))n0 are decreasing and bounded from below by x¯ ;

  2. jj)3

    sequence (g2(xn))n0 is increasing and bounded from above by x¯;

  3. jjj)3

    the following relations hold,

    xnx¯xng2(xn),n=0,1,;
  4. jv)3

    limxn=limg1(xn)=limg2(xn)=x¯.

Proof.

Let be xm[a,b], xm>x¯ and g2(xm)>a, where m. By a) it follows that g1(xm)>x¯ and by b) g2(xm)<x¯. Using e) one obtains g1(xm)<xm. By hypothesis 𝐢3) and 𝐢𝐢3), and by the above relations and (1), for n=m it follows xm+1<xm.

For a1=xm, a2=g1(xm) and a3=g2(xm) in (1.10), by (1.13) and hypothesis d) it follows that xm+1>x¯. Observe that xm+1 in (1) may by expressed in the following way:

xm+1 = g1(xm)f(g1(xm))[g1(xm),xm;f]
[xm,g1(xm),g2(xm);f]f(g1(xm))f(xm)[xm,g1(xm);f][xm,g2(xm);f][g1(xm),g2(xm);f]

and then xm+1<g1(xm). By relation xm+1<xm it follows that g2(xm+1)>g2(xm). Consequently it follows that for every m, the following relations hold:

g2(xm)<g2(xm+1)<x¯<g1(xm+1)<xm+1<g1(xm)<xm.

By these relations conclusions 𝐣𝐣𝐣3) and 𝐣𝐯3) also follow. ∎

Remark 3.2.

If function f is decreasing and concave, then function h(x)=f(x), x[a,b] verifies all the hypothesis of Theorem 3.1 and consequently the sequence (xn)n1 generated by (1) verifies all the conclusions in Theorem 3.1.

Analogously, the following result can be proved

Theorem 3.3.

Let x0[a,b] and functions f,g1,g2 have the following properties:

  1. i4)

    f is increasing on [a,b];

  2. ii)4

    f is concave on [a,b];

  3. iii)4

    f,g1, and g2 verify hypothesis a)-e);

  4. iv)4

    x0<x¯ and g2(x0)b.

Then sequences (xn)n0, (g1(xn))n0 and (g2(xn))n0 generated by (1) verify the properties:

  1. j)4

    (xn)n0 and (g1(xn))n0 are increasing and bounded from above by x¯;

  2. jj)4

    sequence (g2(xn))n0 is decreasing and bounded from below by x¯;

  3. jjj)4

    the following relations hold:

    x¯xng2(xn)xn,n=0,1,2,,
  4. jv)4

    limxn=limg1(xn)=limg2(xn)=x¯.

Remark 3.4.

If function f is decreasing and convex, then function h(x) =f(x), x[a,b] is increasing and concave. It follows that sequence (xn)n0 generated by (1) verifies the hypotheses and all the conclusions of Theorem 3.3.

4. The determination of the auxiliary functions

In the following, for every situation concerning the monotonicity and convexity of function f, functions g1 and g2 can be determined such that conditions a), b), c) and e) should be verified. In the following we thoroughly present the case in which function f is increasing and convex.

Supposing that f(x)>0 for every x[a,b], one considers the functions:

(4.1) g1(x)=xf(x)f(b),

and

(4.2) g2(x)=xf(x)f(a).

Then

g1(x)=1f(x)f(b)0

for every x[a,b] and, consequently g1 is increasing. Analogously

g2(x)=1f(x)f(a)0

for every x[a,b] and then g2 is a decreasing function. It follows that g1 and g2 verify hypotheses a) and b). The hypothesis e) is also verified because g1(x)<1, for every x[a,b]. In iv1) the inequality g2(g1(x0)a holds. Because g1(x0)<x0 and x0>x¯, the above inequality is verified if g2(x0)a. This follows, because g2(g1(x0))>g2(x0) (g2 is a decreasing functions). This means that the following relation must hold:

x0f(x0)f(a)a

i.e.

x0a+f(x0))f(a).

Because f(x0)>0 and f(a)>0, for the veridicity of last inequality it is sufficient that

a+f(x0)f(a)x¯,

where x¯ is the root of equation (1.1). This last inequality may be realized if x0 is sufficiently close to x¯.

Consequently, the hypotheses of theorems 2.1 and 3.1 are realized for g1 and g2 considered above. The other cases may be similarly analyzed.

5. Order of convergence

In the following we prove that every method (1) and (1) have the order of convergence three. The order of convergence for method (1) was treated in [15], and this order is at least three. Assume the following hypotheses:

α) function g2 verifies relation

|g2(x)g2(x¯)|p|xx¯|,

for every x[a,b], where p>0, p;

β) m|f(x)|M, for every x[a,b], where m>0 and M>0 are real numbers.

γ)|3[f(x)]2f(x)f′′′(x)|q, for every x[a,b], where q>0, q .

In hypotheses α),β),γ) and e), by (1.10) and (1.13), for sequence (xn)n0 generated by (1), one obtains:

|x¯xn+1|qM3L2p6m5|x¯xn|3,n=0,1,2,

and this means that the order of convergence for method (1) is at least three.

With the same hypotheses, for tree sequence (xn)n0 generated by (1) it follows:

|x¯xn+1|qM3Lp6m5|x¯xn|3,n=0,1,2,

i.e. the order of convergence of (1) is also at least three.

References

  • [1]
  • [2] Balázs, M., A bilateral approximating method for finding the real roots of real equations, Rev. Anal. Numér. Théor. Approx., 21 (2), pp. 111–117, 1992.
  • [3] Casulli, V., Trigiante, D., The convergence order for iterative multipoint procedures, Calcolo, 13 (1), pp. 25–44, 1997.
  • [4] Costabile, F., Gualtieri, I. M., Luceri, R., A new iterative method for the computation of the solution of nonlinear equations, Numer. Algorithms, 28, pp. 87–100, 2001.
  • [5] Frontini, M., Hermite interpolation and a new iterative method for the computation of the roots of non-linear equations, Calcolo, 40, pp. 109–119, 2003.
  • [6] Grau, M., An improvement to the computing of nonlinear equation solutions, Numer. Algorithms., 34, pp. 1–12, 2003.
  • [7] Ostrowski, A., Solution of Equations in Euclidian and Banach Spaces, Academic Press, New York and London, 1973.
  • [8] Păvăloiu, I., Optimal efficiency index for iterative methods of interpolatory type, Computer Science Journal of Moldova, 1 (5), pp. 20–43, 1997.
  • [9] Păvăloiu, I., Approximation of the roots of equations by Aitken-Steffensen-type monotonic sequences, Calcolo, 32 (1–2), pp. 69–82, 1995.
  • [10] Păvăloiu, I., Optimal problems concerning interpolation methods of solution of equations, Publications de L’Institut Mathématique, 52 (66), pp. 113–126, 1992.
  • [11] Păvăloiu, I., Optimal effiency index of a class of Hermite iterative methods, with two steps, Rev. Anal. Numér. Théor. Approx., 29 (1), pp. 83–89, 2000.
  • [12] Păvăloiu, I., Local convergence of general Steffensen type methods, Rev. Anal. Numér. Théor. Approx., 33 (1), pp. 79–86, 2004.
  • [13] Păvăloiu, I. and Pop, N., Interpolation and Applications, Risoprint, Cluj-Napoca, 2005 (in Romanian).
  • [14] Păvăloiu, I., On a Steffensen-Hermite-type Method for approximating the solution of nonlinear equations, Rev. Anal. Numér. Théor. Approx., 25 1, pp. 87–94, 2006.
  • [15] Păvăloiu, I., Bilateral approximation of solutions of equations by order-three Steffensen type methods, Studia Univ. “Babeş-Bolyai”, Mathematica, Vol. LI, no. 3, pp. 105–114, 2006.
  • [16] Traub, J. F., Iterative Methods for Solutions of Equations, Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1964.
  • [17] Turowicz, B. A., Sur les derivées d’ordre supérieur d’une function inverse, Ann. Polon. Math., 8, pp. 265–269, 1960.
  • [18]
2006

Related Posts