Return to Article Details On approximating the solutions of equations by the chord method and a method of Aitken-Steffensen type

On approximating the solutions of equations by the chord method and a method
of Aitken-Steffensen type

Adrian Diaconu§

January 14, 2013.

§Babeş-Bolyai University, Faculty of Mathematics and Computer Science, st. M. Kogălniceanu no. 1, 3400 Cluj-Napoca, Romania, e-mail: adiaconu@math.ubbcluj.ro.

In [ 13 ] we have studied the existence and the convergence of iterative methods that use generalized abstract divided differences (this notion being defined there). We have indicated a construction model for these differences as well. A special place has been given to the iterative method of the chord for which we have established a convergence theorem which in the same time ensures the existence of the solution of the considered equation. We have obtained the convergence order with the value 1+52. This value is inferior to 2, this last value representing the convergence order of the method of Newton-Kantorovich. This diminuation of the convergence order is the price to pay for the replacement of the Fréchet differential with the generalized abstract divided difference. In this paper we consider the issue of the improvement of the convergence order with respect to the method of Steffensen and Aitken-Steffensen or their generalizations through the method of the auxiliary sequences. This method will be presented in the paper together with the specification of the convergence order of the main sequence and the auxiliary sequences.

MSC. 65J15

Keywords. Convergence of the method of the Aitken-Steffensen type.

1 Introduction

In [ 13 ] we have studied certain approximation methods for the solutions of equations in linear normed spaces, methods that use the abstract divided differences or the generalized abstract divided differences. The main result concerned the study of the convergence of the iterative method of chord with the fixing of its convergence order.

Let us consider X,Y two linear normed spaces, denote by X:XR and Y:YR their norms respectively, and by θX and θY their null elements respectively. By (X,Y) we denote the set of the linear and continuous mappings defined from X to Y. The set (X,Y) is a linear normed space as well, if we define the norm :(X,Y)[0,+[, by U=suphX, hX=1U(h) for any U(X,Y). For the case of Y=R we denote by X the set (X,R), this set representing the space of real, linear and continuous functionals defined on the linear normed space X.

Let us consider now a set DX and a nonlinear mapping f:DY. Using this mapping we have the equation:

f(x)=θY.
1

We will study the approximation of its solutions.

In order to clarify the aforementioned notions we have the following definition:

Definition 1.1

Considering the nonlinear mapping f:DY together with the points x,yD, xy, any mapping Γf;x,y(X,Y) that verifies the equality:

Γf;x,y(xy)=f(x)f(y)
2

is called generalized abstract divided difference of the function f:DY at the points x,y.

In connection with the previous definition we have the following remark:

Remark 1.2

a) If we consider the theorem according to which in every linear normed space (X,X), for any aX{θX} there exists a linear and continuous functional uX such that u=1 and u(a)=aX. Therefore, for any x,yX with xy there exists the functional UxyX such that Uxy=1 and Uxy(xy)=xyX. At the same time there exists the functional UyxX such that Uyx=1 and Uyx(yx)=yxX as well. In the paper [ 13 ] there appears the mapping [x,y;f](X,Y), defined by the equality:

[x,y;f]h=Uxy(h)f(x)+Uyx(h)f(y)xyX
3

for any hX.
  This mapping verifies the equality (2) and it is called abstract divided difference of the nonlinear mapping f:DY at the points x,yD with xy. This mapping is a special case of generalized abstract divided difference.

  b) Let us suppose now that the space X is a space with a scalar product |:X×XR. Defining X:XR by xX=x|x, the space (X,X) is a linear normed space.

  For any x,yX with xy the functional UxyX from a) will be defined by:

Uxy(h)=h|xyxyX

for any hX. So, for the same elements x,yX with xy we have that the abstract divided difference [x,y;f](X,Y) is defined by:

[x,y;f]h=xy|h(f(x)f(y))xyX2

for any hX.

â–¡

The main result of the paper [ 13 ] concerns the convergence of the chord method for the approximation of a solution of the equation (1). This method consists in the consideration of an approximant sequence (xn)nND (here N=N{0} and this notation remains valid for the rest of this paper) and that for any nN verifies the equality:

Γf;xn1,xn(xn+1xn)+f(xn)=θY.
4

If we suppose the existence for a certain nN of the mapping Γf;xn1,xn1(Y,X) representing the inverse of the mapping Γf;xn1,xn(X,Y), for this number nN the equality (4) is equivalent to:

xn+1=xnΓf;xn1,xn1f(xn)=xn1Γf;xn1,xn1f(xn1).
5

2 Theorem of convergence of the chord method

The main result with regard to the possibility that the equality (4) can be written under the form (5), together with the convergence of the sequence (xn)nND to a solution of the equation (1) the existence of which is also proved, is the following theorem:

Theorem 2.1

We suppose that the following assumptions hold:

  • (X,X) is a Banach space;

  • the nonlinear mapping f:DY, for any x,yD with xy admits a generalized abstract divided difference Γf;xn1,xn(X,Y) and there exists a number L>0 such that for any x,y,zD with xy and yz we have the following inequality:

    Γf;x,yΓf;y,zLxzX;
    6
  • the sequence (xn)nND is such that for any nN we have xn1xn and the equality (4) is verified;

  • referring to the initial elements x0,x1D of the sequence (xn)nND we suppose the carrying out of the following conditions:

    • the mapping Γf;x0,x1(X,Y) is invertible and Γf;x0,x11(Y,X);

    • there exist the numbers h0]0,1[ and B0>0 such that we have the following inequality:

      Γf;x0,x11B01h0;
      7
    • if we note R0=f(x0)Y and h1=LB02R0(1h0)2 there exists a number q[12,1[ such that we have the following inequality:

      d=max{h0(1q)2,[h1(1q)2]1α}<1,
      8

      where α=1+52;

    • if δ=(1q)2LB0dα1d and S(x0,δ)={xX/xx0Xδ} the relation x1S(x0,δ)D is true.

Then the following conclusions are true:

  • for any nN we have that xnS(x0,δ), there exists the mapping Γf;xn1,xn1(Y,X), and the sequence (xn)nN verifies the equality (5);

  • the sequence (xn)nNX is convergent;

  • if x=limnxn, then xS(x0,δ) and f(x)=θY;

  • for any nN the following estimate of the error of approximation is true:

    xxnX(1q)2LB0dαn+11dαn.
    9

The proof of this theorem is given in [ 13 ] .

Regarding this theorem we have the following remark:

Remark 2.2

From the conclusions of the theorem 2.1, corroborated with the fact that d<1, we deduce that there exists a number NN such that for any number nN, nN we have the following inequality:

xxnXdαnLB0.
10

This inequality indicates that the convergence order of the iterative method of the chord is α=1+52.â–¡

The inequality (10) is obvious from (9).

3 Certain remarks in connection with the previous result

We have the following statements:

Remark 3.1

Under the hypotheses of Theorem 2.1 the sequence of real and positive numbers (Γf;xn1,xn1)nN has an upper bound, for any nN the following inequality taking place:

Γf;xn1,xn1B0e(1q)d1dα1.
11

â–¡

Indeed, as Γf;xn1,xn1Bn, and using the recurrence relation of the sequence (Bn)nN, we have BnBn111hn1 and so:

BnB0(1h0)(1h1)...(1hn1).

Evidently:

1(1h0)(1h1)...(1hn1)[1nj=0n111hj]n==[1+1nj=0n1hj1hj]n[1+1n11qj=0n1hj]n.

But:

j=0n1hj(1q)2j=0n1dαj(1q)2d1dα1,

therefore:

1(1h0)(1h1)...(1hn1)[1+1n(1q)d1dα1]n<exp((1q)d1dα1),

from where the inequality (11) is evident.

Remark 3.2

If in the hypotheses of Theorem 2.1 we choose q=12, the values of the other constants from this theorem are:

d=max{4h0,(4h1)1α}<1, δ=dα4LB0(1d),
12

and the inequality that expressed an upper bound of the error by which xn approximates x is:

xxnXdαn+14LB0(1dαn).
13

For the upper bound of the sequence (Γf;xn1,xn1)nN for any nN, we have the following inequality:

Γf;xn1,xn1B0ed2(1dα1).
14

â–¡

The statements are obvious:

Remark 3.3

For any nN there exists the generalized abstract divided difference Γf;xn,x(X,Y), this mapping is invertible, so the mapping Γf;xn,x1(Y,X) exists and the sequence (Γf;xn,x1)nN has an upper bound, more precisely, there exists a number pN such that for any nN, np we have the inequality:

Γf;xn,x12B0e(1q)d1dα1.
15

â–¡

It is obvious that xxn for any nN. Indeed, if there exists a number n0N such that x=xn0 then:

xn0+1=xn0Γf;xn01,xn01f(xn0)=xn0,

which is impossible.

As xxn for any nN we deduce the existence of the mapping Γf;xn,x(X,Y) representing the generalized abstract divided difference of the function f:DY on the points xn and x.

In order to prove the invertibility of the mapping Γf;xn,x(X,Y) for any any nN let us consider the following expression:

Wn=Γf;xn1,xn1(Γf;xn1,xnΓf;xn,x)(X,X).

We have:

Γf;xn,x=Γf;xn1,xn(IXWn).

It is clear that:

WnBnLxxn1XB0Le(1q)d1dα1xxn1X.

As limnxxn1X=0, we deduce that there exists a number pN such that for any nN, np, we have the inequality:

B0Le(1q)d1dα1xxn1X12<1,

therefore Wn12<1 and so there exists the mapping (IXWn)1(X,X) and:

(IXWn)111Wn2.

From the existence of the mapping Γf;xn1,xn1(Y,X) and using the inequality Γf;xn1,xn1Bn we deduce the existence of the mapping Γf;xn,x1(Y,X) by the equality:

Γf;xn,x1=(IXWn)1Γf;xn1,xn1,

again:

Γf;xn,x1(IXWn)1Γf;xn1,xn12Bn2B0e(1q)d1dα1.

Remark 3.4

The sequence (Γf;xn,x)nN also has an upper bound and for any nN we have the inequality:

Γf;xn,xΓf;x0,x+(1q)2B0dα1d.
16

â–¡

Indeed, for any nN, we have:

Γf;xn,x=Γf;x0,x+k=1n(Γf;xk,xΓf;xk1,x)Γf;x0,x+k=1nΓf;xk,xΓf;xk1,xΓf;x0,x+Lk=1nxkxk1XΓf;x0,x+(1q)2B0k=1ndαkΓf;x0,x+(1q)2B0dα1d.

Remark 3.5

There exists a number NN such that for the same values nN with nN we have that:

xxnXdαnLB0,
17

and this inequality proves that the convergence order of this method is at least α=1+52. â–¡

Indeed, as αn+1=αn+αn1 it is clear that dαn+11dαn=dαndαn11dαn and from the fact that limndαn11dαn=0 the conclusion is obvious.

4 The Acceleration of the Convergence.
Iterative Methods of the Aitken-Steffensen type

The main conclusion of the introduction is the fact that the convergence order of the chord method is α=1+52.

At the same time it is well known that the convergence order of the Newton-Kantorovich method is 2, therefore greater. But the Newton-Kantorovich method uses the Fréchet differential instead of the divided difference.

One naturally thinks of improving the convergence order of the chord method, without renouncing at the divided difference in the favor of the Fréchet differential.

One way of doing this is using Steffensen’s method, that uses a mapping Q:XX that verifies the inclusion Q(D)D for a set DX, generates the sequence (xn)nND, starting from an arbitrary x0D, by the verification for any nN of the equality:

Γf;xn,Q(xn)(xn+1xn)+f(xn)=θY.
18

If for any nN there exists the mapping Γf;xn,Q(xn)1(Y,X) the last equality is equivalent to:

xn+1=xnΓf;xn,Q(xn)1f(xn).
19

A more general case is the one which use two mappings Q1,Q2:XX which verify for i{1,2} the relations Qi(D)D. Starting from an arbitrary x0D one build the sequence (xn)nND by the verification for any nN of the equality:

Γf;Q1(xn),Q2(xn)(xn+1Q1(xn))+f(Q1(xn))=θY,
20

equality which in the hypothesis of the existence of Γf;Q1(xn),Q2(xn)1(Y,X) is equivalent to:

xn+1=Q1(xn)Γf;Q1(xn),Q2(xn)1f(Q1(xn)).
21

A study of the convergence of this method, known as the iterative method of Aitken-Steffensen was made by Păvăloiu in [ 21 ] . The established result requests very strong conditions imposed on the mappings f:DY together with Q1,Q2:XX and one requests the verification of these conditions on every point of a set DX.

We propose a more general frame and we will build an iterative process after as follows.

Let us consider a initial element x0D. Besides the main sequence (xn)nN D we also use two auxiliary sequences (yn)nN,(zn)nND.

For these auxiliary sequences we request the existence of the numbers K1, K2,p,q>0 such that for any nN the following inequalities are verified:

f(yn)YK1f(xn)Yp,  f(zn)YK1f(xn)Yq.
22

Then, if for a number nN we have available the elements yn,znD starting from xnD, we will generate the new iterate xn+1D by the following relation:

Γf;yn,zn(xn+1yn)+f(yn)=θY.
23

On account of the property of definition of the mapping Γf;yn,zn(X,Y) the equality (23) is equivalent with:

Γf;yn,zn(xn+1zn)+f(zn)=θY.
24

If for any nN there exists the mapping Γf;yn,zn1(Y,X) we have:

xn+1=ynΓf;yn,zn1f(yn)=znΓf;yn,zn1f(zn).
25

In connection to the main sequence (xn)nN and the auxiliary sequences (yn)nN,(zn)nND we have the following remarks. Here for a number kN we denote by (X(k),Y) the set of mappings defined from Xk to Y that are k-linear and continuous. This set is a linear normed space as well with the norm :(X(k),Y)Y is defined for U(X(k),Y) by U=suphiX; hi=1; i=1,kU(h1,...,hk)Y.

Remark 4.1

If the mapping f:DY admits Fréchet differentials up to the order s1, where s=max{p,q} and p,qN, the mapping f(s1):D(X(s1),Y) is a Lipschitz mapping, namely there exist the constants M>0 such that for any x,yD the inequality:

f(s1)(x)f(s1)(y)MxyX,
26

is verified and for any nN the inequalities:

j=0p11j!f(j)(xn)(ynxn)jYαf(xn)Yp,j=0q11j!f(j)(xn)(znxn)jYβf(xn)Yq,ynxnXaf(xn)Y, znxnXbf(xn)Y,

then:

  • the mappings f(p1):D(X(p1),Y) and f(q1):D(X(q1),Y) are Lipschitz mappings, namely there exist the constants M1,M2>0 such that for any x,yD the following inequalities are true:

    f(p1)(x)f(p1)(y)M1xyX,f(q1)(x)f(q1)(y)M2xyX;
  • the sequences (xn)nN, (yn)nN,(zn)nND verify the equalities (22) with:

    K1=α+M1app!, K2=β+M2bqq!.

â–¡

Indeed, from the hypothesis (26) we deduce that for any kN, ks1 the mapping f(k):D(X(k),Y) is a Lipschitz mapping as well (we can use the well known theorem of Lagrange), particularly, the mappings f(p1) :D(X(p1),Y) and f(q1):D(X(q1),Y) are Lipschitz mappings and the rest of the conclusion j) is obvious especially the inequalities (28).

From these inequalities we deduce easily that for any x,yD the following inequalities are also verified:

f(y)j=0p11j!f(j)(x)(yx)jYM1p!yxXp,f(y)j=0q11j!f(j)(x)(yx)jYM2q!yxXq

(we can use Taylors’ formula for mappings).

On account of the fact that xn,yn and znD, using the inequalities (27) and (29) we deduce that for any nN the following inequalities are verified:

f(yn)Yf(yn)j=0p11j!f(j)(xn)(ynxn)jY+j=0p11j!f(j)(xn)(ynxn)jYM1p!ynxnXp+αf(xn)Yp(α+M1app!)f(xn)Yp,

and similarly:

f(zn)YM2q!znxnXq+βf(xn)Yq(β+M2bqq!)f(xn)Yq,

and these letter inequalities justify the statement of the present remark.

Remark 4.2

It is clear that if the first of the inequalities (22) is verified for any nN with a certain K1>0, this inequality is verified with any number Kmax{1,K1}. The situation is identical regarding the second inequality from (22). In conclusion we can suppose that in these relations we have K1=K2=K1.

Identically, we can suppose that in the inequalities:

ynxnXaf(xn)Y, znxnXbf(xn)Y,

that are true for any nN, we can have b=a1.
  

In conclusion, for the main sequence (xn)nND together with the auxiliary sequences (yn)nN,(zn)nND we can suppose that for any nN we have ynzn and there exist the numbers K,a1 such that for any nN the following inequalities are verified:

f(yn)YKf(xn)Yp,f(zn)YKf(xn)Yq,ynxnXaf(xn)Y,znxnXaf(xn)Y

â–¡

The statements are obvious.

5 The convergence of some auxiliary real number sequences

In connection with the enounced problem we consider, for the real numbers p,q1, the following equation in x on the interval [0,+[:

xp+q1+2x2+2x1=0.
33

We have the following remarks:

Remarks 5.1

a) The equation (33) has an unique positive root and this root is α]0,1[.
  b) If α]0,1[ is the root of the equation (33) one verifies the following inequalities as well:

α2+α1<0, α2+2α1<0, 2α2+2α1<0
34

and these inequalities are equivalent to the following inequalities respectively:

0<α21α<1, 0<α1αα2<1, 0<α212αα2<1.
35

â–¡

Indeed, let us consider the function φ:[0,+[R defined by φ(x)=xp+q1+2x2+2x1. It is obvious that for any x[0,+[ there exists the derivative φ(x) at the point x and:

φ(x)=(p+q1)xp+q2+4x+2.

As it is obvious that for any x[0,+[ we have that φ(x)>0, therefore the function φ:[0,+[R is a strictly increasing function, therefore an injective function, thus the equation φ(x)=0 has at most one root.

As φ(0)=1, φ(1)=4, this root exists indeed and it belongs to the interval ]0,1[.

If α]0,1[ is the root of the equation (33) it is clear that 2α2+2α1=
=αp+q1<0, whence α2+2α1=αp+q1α2<0 and α2+α1=
=αp+q1α2α<0, therefore the relations (34) are true.

The inequality α2+α1<0 is equivalent to 0<α2<1α therefore to 0<α21α<1, the inequality α2+2α1<0 is equivalent to 0<α<1αα2, therefore to 0<α1αα2<1, again the inequality 2α2+2α1<0 is equivalent to 0<α2<12αα2 therefore to 0<α212αα2<1.

Let us consider now the numbers a,K,L,B0,R0>0 and the numbers p,q1 and using these numbers we build the real number sequences (un)nN, (sn)nN, (vn)nN, (wn)nN, (tn)nN, (Bn)nN and (Rn)nN using the following recurrence relations:

un=LKBn2Rnp,sn=LKBn2Rnq,vn=aL2K2Bn3Rnp+q1un,wn=LKBn2Rnq(1un)(1vn),tn=aL2K2Bn3Rnp+q(1un)(1vn)(1wn),Bn+1=Bn(1un)(1vn)(1wn)(1tn),Rn+1=LK2Bn2Rnp+q.

It is obvious that this construction has a meaning if for any nN we have that un, vn, wn, tnR{1} and Bn,Rn>0.

It is clear that for any nN we have:

vn=aBnunsn1un,wn=sn(1un)(1vn),tn=vn(1vn)(1wn),Rn+1=unsnLBn2,

as well.

Referring to the sequences that are defined by the relations (36) we have the following proposition:

Proposition 5.2

If the following inequalities are verified:

aB01Lmin{Kpq+12(q1),Kqp+12(p1)}
38

(with the specification that for q=1 the expression that has q1 in its denominator is +, and the same for the expression that has p1 in its denominator) and:

d=LKB02α2max1p+q1{R0p(p+q1)Kpq+1(LB02)q1,R0q(p+q1)Kqp+1(LB02)p1}<1
39

where α]0,1[ is the unique root of the equation (33), then for any nN we have the following inequalities:

unαd(p+q)n<α<1,
snαd(p+q)n<α<1,vnα21αd2(p+q)n<α21α<1,wnα1αα2d(p+q)n<α1αα2<1,tnα212αα2d2(p+q)n<α212αα2<1,Bn+1Bn12α2α2,Rn+1α2LB02d2(p+q)n.

Proof â–¼
Let us consider first that q>1.

From the definition of d it is clear that:

LKB02R02α2[Kpq+1(LB02)q1]1p+q1d,

therefore:

u0αα(LB02Kpq+1q1)q1p+q1d.

From B01LKpq+12(q1) we deduce that Kpq+1q1LB021, therefore as q1p+q10 it is clear that (Kpq+1q1LB02)q1p+q11>α, therefore:

α(LB02Kpq+1q1)q1p+q1<1,
42

therefore:

u0αd=αd(p+q)0<α<1.

For q=1 from the same definition of d we have:

LKB02α2R0pKd,

namely u0ααKd.

But α<1K, so αK<1, therefore u0αd.

Identically if p>1 we have that:

LKB02R0qα2[Kqp+1(LB02)p1]1p+q1d,

whence, in the same manner as in the case of u0, by inverting the roles of the numbers p and q, we deduce that s0αd=αd(p+q)0<α<1.

For p=1 one can show the inequality s0αd in the same manner as the inequality u0αd for the case of q=1.

As v0=aB0u0s01u0 and using aB0 and the inequalities concerning to u0 and s0, we have that:

v0α21αd2=α21αd2(p+q)0<α21α<1,

the last inequality being the first from (35).

Also:

w0=s0(1u0)(1v0)αd(1α)(1α21α)==α1αα2d(p+q)0<α1αα2<1,

the last inequality being the second from (35), while:

t0=v0(1v0)(1w0)α21αd(1α21α)(1α1αα2)==α212αα2d(p+q)0<α212αα2<1,

the last inequality being the third from (35).

As  u0,v0]0,1[ it is clear that v0,w0,t0]0,1[ as well.

Then we have:

B1B0=1(1u0)(1v0)(1w0)(1t0)1(1α)(1α21α)(1α1αα2)(1α212αα2)=112α2α2,

therefore:

B1B012α2α2.

As:

R1=1LB02u0s0=α2LB02d2=α2LB02d2(p+q)0.

From the afore established relations we deduce that the properties (40)-(41) are true in the case of n=0.

We suppose that these properties are true for any nN with nk and we prove that they are also true for n=k+1.

For any iN with ik we have that:

ui+1=LKBi+12Ri+1p=LKBi2(12α2α2)2LpK2pBi2pRip(p+q)=Lp+1K2p+1Bi2p+2Rip(p+q)(12α2α2)2=Lp+1K2p+1(12α2α2)2(Bi2Rip)p+qBi2q2=Lp+1K2p+1(12α2α2)21Bi2q2(uiLK)p+q.

From the equality that defined the sequence (Bn)nN we deduce that:

BiBi1=11ui111vi111wi111si11,

every fraction from those that multiply being greater than the unit.

Therefore BiBi1.

From this we actually have that BiB0, from where as 2q20 we have Bi2q2B02q2, namely 1Bi2q21B02q2 and so:

ui+1L1qKpq+1(12α2α2)21B02(q1)uip+q=Cuip+q,
43

where:

C=Kpq+1(LB02)q11(12α2α2)2.

If we multiply by C1p+q1 the two members of the inequality (43) we obtain that:

C1p+q1ui+1(C1p+q1ui)p+q,

namely the inequality hi+1hip+q if hi=C1p+q1ui.

As i{0,1,...,k} we immediately deduce that:

h1h0p+q, h2h1p+qh0(p+q)2, ... ,hk+1hkp+q=h0(p+q)k+1.

Therefore we have that:

uk+1(1C)1p+q1(C1p+q1u0)(p+q)k+1.
44

But α]0,1[ is the root of the equation (33), therefore:

(12α2α2)2=α2(p+q1),

therefore:

C1p+q1u0=1α2[Kpq+1(LB0)q1]1p+q1LKB02R0p=LKB02α[R0p(p+q1)Kpq+1(LB0)q1]1p+q1d.

At the same time, if q>1 we have that:

(1C)1p+q1=(12α2α2)2p+q1[(LB02)q1Kpq+1]1p+q1=α2(LB02Kpq+1q1)q1p+q1<α;

for the last inequality we have take into account (42).

For q=1 we have:

(1C)1p+q1=(1C)1p=α2K<α.

So, from the inequality (44) we obtain that uk+1αd(p+q)k+1 and as d<1, we obtain uk+1αd(p+q)k+1<α<1.

One can obtain the inequality sk+1αd(p+q)k+1<α<1 in a similar manner with the one concerning uk+1. It is only necessary to invert the roles of the numbers p and q.

We have:

vk+1=aBk+1uk+1sk+11uk+1.

Identically, as BiBi11 for any ik we deduce that Bk+1Bk1, therefore Bk+1Bk and in fact Bk+1B0, namely:

vk+1aB0α2d2(p+q)k+11α.

As aB0 and d<1 we have:

vk+1α21αd2(p+q)k+1<α21α<1.

In the same manner we have:

wk+1=sk+1(1uk+1)(1vk+1)αd(p+q)k+1(1α)(1α21α)=α1αα2d(p+q)k+1<α1αα2<1

and:

tk+1=vk+1(1vk+1)(1wk+1)α21αd2(p+q)k+1(1α21α)(1α1αα2)=α212αα2d2(p+q)k+1<α212αα2<1.

For the sequence (Bn)nN we have:

Bk+2Bk+1=1(1uk+1)(1vk+1)(1wk+1)(1tk+1)1(1α)(1α21α)(1α1αα2)(1α212αα2)=112α2α2,

therefore:

Bk+2Bk+112α2α2.

Finally:

Rk+2=uk+1sk+1LBk+121LB0α2d2(p+q)k+1.

Therefore the inequalities (40) are also true for n=k+1.

On the basis of the principle of the mathematical induction these
inequalities are true for any nN. The proposition is proved.

6 The main result

We return to the issue of the convergence of the sequences (xn)nN, (yn)nN, (xn)nNDX to the solution of the equation f(x)=θY, where f:DY.

We have the following fundamental result:

Theorem 6.1

Suppose that the following assumptions hold:

  • The linear normed space (X,X) is a Banach space;

  • The mapping f:DY admits for any x,yD with xy a generalized abstract divided difference Γf;x,y(X,Y) and there exists a number L>0 such that for any x,y,zD with xy, yz we have the inequality:

    Γf;x,yΓf;y,zLxzX;
  • The main approximant sequence (xn)nN together with the secondary sequences (yn)nN and (zn)nN are such that for any nN the following equality is fulfilled:

    Γf;yn,zn(xn+1yn)+f(yn)=θY
    45

    and the inequalities (32) with the constants a,K>0. We also have that f(yn),f(zn)Y{θY}, ynzn and we are in one of the following situations:

    • xnyn and yn+1zn,

    or

    • xnzn and zn+1yn.

  • The mapping Γf;y0,z0(X,Y) is invertible and Γf;y0,z01(Y,X).

  • Denoting:

    B0=max{a,Γf;y0,z01},R0=f(x0)Y,K=max{K,(B0L)2(q1)pq+1,(B0L)2(p1)qp+1},d=LKB02α2max{R0p(p+q1)Kpq+1(LB02)q1,R0q(p+q1)Kqp+1(LB02)p1},δ=2aR0+aα2LB02d21d2(p+q1)+2αLKB0d1dp+q1,

    where α]0,1[ is the unique root of the equation (33), the conditions d<1 and S(x0,δ)={xX/xx0Xδ} are fulfilled.

Then the following conclusions are true:

  • for any nN we have that xn,yn,znS(x0,δ), there exists the mapping Γf;yn,zn1(Y,X) and:

    xn+1=ynΓf;yn,zn1f(yn)=znΓf;yn,zn1f(zn);
    46
  • the sequences (xn)nN,(yn)nN,(zn)nNX are convergent to the limit xS(x0,δ) for that f(x)=θY;

  • for any nN the following inequalities are fulfilled:

    xn+1xnXaα2LB02d2(p+q)n1+αLKB0d(p+q)n;
    47
    xnxXaα2LB02d2(p+q)n11d2(p+q)n1(p+q1)+αLKB0d(p+q)n1d(p+q)n(p+q1);
    max{ynxX,znxX}aα2LB02d2(p+q)n12d2(p+q)n1(p+q1)1d2(p+q)n1(p+q1)+αLKB0d(p+q)n1d(p+q)n(p+q1).\medskip

The proof of this result we will given in the following paper of this paperset.

Bibliography

1

I. K. Argyros, Polynomial Operator Equation in Abstract Spaces and Applications, CRC Press Boca Raton Boston London New York Washington D.C., 1998.

2

M. Balázs, On the chord method for solving nonlinear operational equations, Stud. Cerc. Mat., 20 (1968) no. 2, pp. 129–136, (in Romanian).

3

M. Balázs, M. and G. Goldner,G. Divided differences in Banach spaces and certain applications, Stud. Cerc. Mat., 21 (1969) no. 7, pp. 985–986 (in Romanian).

4

M. Balázs, On applying the chord method in solving nonlinear integral equations, Stud. Cerc. Mat., 23 (1971) no. 6, pp. 841–844, (in Romanian).

5

M. Balázs, and G. Goldner, On existence of divided differences in linear spaces, Rev. Anal. Numér. Théor. Approx., (1973) no.2, pp. 5–9. \includegraphics[scale=0.1]{ext-link.png}

6

J. G. P. Barnes, An algorithm for solving nonlinear equations based on the secant method, Cumput. J., (1965), pp. 55–72.

7

A. Diaconu, Interpolation dans les espaces abstaits. Méthodes itératives pour la résolution des équations opérationnelles obtenues par l’interpolation inverse (I), “Babeş-Bolyai" University, Faculty of Mathematics, Research Seminaries, Seminar of Functional Analysis and Numerical Methods, Preprint, (1981) no. 4, pp. 1–52, (in French).

8

A. Diaconu, Interpolation dans les espaces abstaits. Méthodes itératives pour la résolution des équations opérationnelles obtenues par l’interpolalation inverse (II), “Babeş-Bolyai" University, Faculty of Mathematics, Research Seminaries, Seminar of Functional Analysis and Numerical Methods, Preprint, (1984), no. 1, pp.  41–97, (in French).

9

A. Diaconu, Interpolation dans les espaces abstraits. Méthodes itératives pour la résolution des équations opérationnelles obtenues par l’interpolalation inverse (III), “Babeş-Bolyai" University, Faculty of Mathematics, Research Seminaries, Seminar of Functional Analysis and Numerical Methods, Preprint, (1985), no. 1, pp.  21–71, (in french).

10

A. Diaconu, Remarks on interpolation in certain linear spaces (I), Researches on Numerical Methods of Analysis and Optimization, (2000) 2-1(3), pp. 3–14.

11

A. Diaconu, Remarks on Interpolation in Certain Linear Spaces (II), Researches on Numerical Methods of Analysis and Optimization, 2-2(4), (2000), pp. 143–161.

12

A. Diaconu, On the interpolation in linear spaces, Mathematica Pannonica, 14/2 (2003), pp. 205–216.

13

A. Diaconu, On the chord method in linear normed spaces, Topics in Mathematics, Computer Science and Philosophy, A Festschrift for Wolfgang W. Breckner, pp.  73–94.

14

G. Goldner and M. Balázs, On the chord method and its modification for solving nonlinear operational equations, Stud. Cerc. Mat., 20 (1968) no.7, pp.  981–990. (in Romanian).

15

B. Janko and M. Balázs, On the chord method for solving nonlinear operational equations in normed spaces, Stud. Cerc. Mat., 19 (1967), no. 10, pp. 1433–1436, (in Romanian).

16

L. V. Kantorovich, Functional Analysis and Applied mathematics (in Russian), U. M. N. III, 5, (28)(1948), pp.  29-185.

17

L. V. Kantorovich, Functional Analysis (in Romanian), Editura Ştiinţifică şi Enciclopedică, Bucureşti, 1986.

18

V. L. Makarov and V. V. Hlobistov, Osnovı teorii polinomialnogo operatornogo interpolirovania, Institut Mathematiki H.A.H. Ukrain, Kiev, (1998) (in Russian).

19

I. Păvăloiu, Considerations regarding the iterative methods obtained by inverse interpolation, Stud. Cerc. Mat., 23 (1971), no. 10, pp.  1545–1549, (in Romanian).

20

I. Păvăloiu, Interpolation dans des espaces linéaires normés et application, Mathematica (Cluj), 12, (35), 1(1970), pp.  149–158.

21

I. Păvăloiu, Introduction to the theory of approximating the solutions of equations, Editura Dacia, Cluj-Napoca, 1976, (in Romanian).

22

P. M. Prenter, Lagrange and Hermite Interpolation in Banach Spaces, Journal of Approximation Theory, 4 (1971), pp. 419-432.

23

A. S. Sergeev, O metode hord, Sibirski Mat. Journal, 11, (2),(1961), pp.  282–289, (in Russian).

24

S. Ul’m, Ob interpolationih metodah resenia uravnenii v postransve Banaha, Izv. Acad. Nauk Estonskoi S.S.R. 1, (1963), pp.  24–30, (in Russian).

25

S. Ul’m, Ob odnom classe iterationnih metoda v postransve Ghilberta,Izv. Acad. Nauk Estonskoi S.S.R. 2,(1963) pp.  132–140, (in Russian).

26

S. Ul’m, Ob obobscennyh razdelennih raznostiah I., Izv. Acad. Nauk Estonskoi S.S.R. 16, 1,(1967), pp.  13–36, (in Russian).

27

S. Ul’m, Ob obobscennyh razdelennih raznostiah II., Izv. Acad. Nauk Estonskoi S.S.R. 16, 2,(1967), pp.  146–155, (in Russian).