On the chord method

Abstract

Let \(X_{1},X_{2}\) be two Banach spaces, \(f:X_{1}\rightarrow X_{2}\) a nonlinear mapping and consider the chord method for solving the equation \(f\left(x\right) =0\): \[x_{n+1}=x_n-[x_{n-1},x_{n};f]^{-1}f(x_n), \quad n=1,…\] Under some simple conditions on the divided differences of order one of \(f\), of the form \[\|[y, u; f] − [x, y; f]\| ≤ l_1 \|x − u\| ^p + l_2 \|x − y\|^p + l_3 \|y − u\|^p\] we show that the chord method converge to the solution. We obtain error estimations and determine the convergence order.

Authors

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

Keywords

chord method; equations in Banach spaces; error estimation; convergence order

PDF

Cite this paper as:

I. Păvăloiu, On the chord method, Bul. Ştiinţ. Univ. Baia Mare, Seria B. Fasc. Mat.-Fiz., 7 (1991), pp. 61-66.

About this paper

Journal

Bul. Ştiinţ. Univ. Baia Mare, Seria B. Fasc. Mat.-Fiz., 7 (1991)

Publisher Name

Bul. Ştiinţ. Univ. Baia Mare

Print ISSN

Not available yet.

Online ISSN

Not available yet.

References

[1] Argyros, I. K., The secant method and fixed points on nonlinear operators, Mh. Math.,  106 (1988), 85–94.

[2] Balazs, M. ¸si Goldner, G., Observatii asupra diferentelor divizate si asupra metodei coardei. Revista de analiza numerica si teoria aproximatiei vol. 3 (1974) fasc. 1, 19–30.

[3] Pavaloiu, I., Remarks on the secant method for the solution of nonlinear operatorial equations, Research Seminars, Seminar on Mathematical Analysis, Preprint no. 7, (1991), pp. 127 132.

[4] Pavaloiu, I., Introducere in teoria aproximarii solutiilor ecuatiilor. Ed. Dacia, ClujNapoca, 1976.

[5] Schmidt, I. W., Eine Ubertagungen der Regula Falsi auf Gleichungen, in ”Banachraumen” I ZAMM, 48, 1–8 (1963).

[6] Schmidt, I. W., Eine Ubertagungen der Regula Falsi auf Gleichungen, in ”Banachraumen” II 97–110 (1963).

Paper (preprint) in HTML form

On the chord method

Bul.St.Univ.Baia Mare

Seria B, Matematică-Informatică, vol.VII (1991) Nr.1-2, 61-66


On the chord method

I. Păvăloiu (Cluj-Napoca)

In the paper [1], I.K. Argyros considers as divided difference of the mapping f:X1X2, where X1andX2 are Banach spaces, a linear mapping [x,y;f](X1,X2) which fulfils the following conditions:

  • (a)

    [x,y;f](yx)=f(y)f(x), for every x,yD, where DX1 is a subset of X1;

  • (b)

    there exist the real constants l10, l20, l30 and p(0,1] such that for every x,y,uD the following inequality holds:

    [y,u;f][x,y;f]l1xup+l2xyp+l3yup.

In [1] the hypothesis that the equation:

(1) f(x)=0

admits a simple solution x in adopted, and conditions for the convergence of the sequence (xn)n0 generated by the chord method:

(2) xn+1=xn[xn1,xn;f]1f(xn),n=1,2,,x0,x1D0

are given.

In a recent paper [2] there is shown that, with the hypotheses considered in [1], the convergence speed of the sequence generated by (2) and the error estimation are featured by the inequality:

(3) xxn+1αd0t1n+1,

where α is a precised constant, 0<d0<1 and t1 is the positive root of the equation t2tp.

We shall admit further down that the divided difference operator fulfils the conditions (a) and (b), and search for supplementary conditions in order to make equation (1) admit a solution x into a precised domain D0 and the sequence (xn)n0 generated by (2) converge to this solution.

Observe firstly that the identity:

(4) xn[xn1,xn;f]1f(xn)=xn1[xn1,xn;f]1f(xn1)

holds for every n=1,2, with the hypothesis that the linear mapping [xn1,xn;f] admits an inverse mapping.

The following identity

(5) f(xn+1)= f(xn)+[xn1,xn;f](xn+1xn)
+([xn,xn+1;f][xn1,xn;f])(xn+1xn),n=1,2,

Let B>0,α>0,0<d0<1, and x0,xX1. Consider the sphere

(6) U={xX1:xx0Bαd01d0t11}

where t1=1+1+4p2 that is, the positive root, the equation:

(7) t2tp=0

The following theorem holds:

Theorem 1.

If the divided difference [x,y;f] fulfils the conditions (a) and (b) for every x,yU and the following hypotheses:

  • (1)

    the mapping [x,y;f] admits a bounded inverse mapping for every x,yU, namely there exists a constant B>0 such that [x,y;f]1B

  • (ii)
    α=1B(1+p)/p(l1+l2+l3)1/p;
  • (iii)
    x1x0Bαd0,f(x0)αd0,f(x1)αd0t1

    are also fulfilled, then equation (1) has at least one solution xU and the sequence (xn)n0 generated by (2) converges to x, the convergence speed and the error estimation being featured by the inequality:

    xxnBαd0t1n1d0t1n(t11).
Proof.

From (2) for n=1 we deduce:

x2x1Bf(x1)Bαd0t1

from which, taking also into account iii. it follows

x2x0 x2x1+x1x0
Bαd0t1+Bαd0
Bαd0(1+d0t11)
<Bαd01d0t11

from which it results that x2U.

Using the fact that x2U, the identities (4) and (5), and the inequality a), we obtain

f(x2) x2x1(l1x2x0p+l2x1x0p+l3x2x1p)
Bf(x1)(l1Bpf(x0)p+l2x1x0p+l3Bpf(x1)p)
Bαd0t1(l1Bpαpd0p+l2Bpαpd0p+l3Bpαpd0t1)
Bp+1αp+1d0t1+p(l1+l2+l3d0p(t11))
=Bp+1αp+1(l1+l2+l3d0p(t11))d0t12pαd0t12

since αpBp+1(l1+l2+l3d0p(t11))αpBp+1(l1+l2+l3)<1.

From the above inequality follows therefore:

f(x2)αd0t12

Suppose by induction that:

  • (a’)

    xiU,i=0,1,,k;

  • (b’)

    f(xi)αd0t1i,i=1,2,,k.

Then, for xk+1we have:

xk+1x0xk+1xk+xkxk1++x1x0
Bf(xk)+Bf(xk1)++Bαd0
Bαd0t1k+Bαd0t1k1++Bαd0
=Bαd0(1+d0t11+d0t121++d0t1k1)
Bαd0(1+d0t11+d02(t11)++d0k(t11))Bαd01d0t11

from which follows that xk+1U. Proceeding now for xk+1, as in the case of x2, we obtain:

f(xk+1) Bp+1αp+1(l1+l2+l3dopt1k1(t11))d0t1k1(t1+p)
Bp+1αp+1(l1+l2+l3)d0t1k+1αd0t1k+1

It results therefore that the relations (a’) and (b’) hold for every i.

Now we shall show that the sequence (xn)n0 is fundamental.

Indeed, for every n,s we have:

xn+sxn k=nn+s1xk+1xkk=nn+s1Bf(xk)Bαk=nn+s1d0t1k
=Bαd0t1nk=nn+s1d0t1kt1n=Bαd0t1nk=nns1d0t1n(t1kn1)
Bαd0t1nk=nn+s1d0t1n(kn)(t11)=Bαd0t1nk=nn+s1(d0t1n(t11))kn
Bαd0t1n1d0t1n(t11).

By the last inequality and the fact that 0<d0<1 and t1>1 follows that the sequence (xn)n2 is fundamental. For s, from the inequality:

xn+sxnBαd0t1n1d0t1n(t11)

follows the inequality:

xxnBαd0t1n1d0t1n(t11).

In [1] Argyros showed that if the divided difference [x,y;f] fulfils the conditions (a) and (b) then f is Fréchet differentiable and [x,x;f]=f(x). From this fact follows that the mapping f is continuous on B: hence at limit for n in the inequality:

f(xn)αd0t1n,

one obtains:

f(x)0

from which results f(x)=0. With this the theorem is entirely proved ∎

Remark 2.

In [5], [6] Schmidt imposes in the divided difference conditions similar to the conditions (a) and (b) given by Argyros in [1], but for p=1. The same conditions are reproduced in [2], too.

References

Institutul de Calcul

Academia Română

Cluj-Napoca

1991

Related Posts