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:X1→X2, where X1andX2 are Banach spaces, a linear mapping [x,y;f]∈ℒ⁢(X1,X2) which fulfils the following conditions:

  • (a)

    [x,y;f]⁢(y−x)=f⁢(y)−f⁢(x), for every x,y∈D, where D⊆X1 is a subset of X1;

  • (b)

    there exist the real constants l1≥0, l2≥0, l3≥0 and p∈(0,1] such that for every x,y,u∈D the following inequality holds:

    ‖[y,u;f]−[x,y;f]‖≤l1⁢‖x−u‖p+l2⁢‖x−y‖p+l3⁢‖y−u‖p.

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)n≥0 generated by the chord method:

(2) xn+1=xn−[xn−1,xn;f]−1⁢f⁢(xn),n=1,2,…,x0,x1∈D0

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) ‖x∗−xn+1‖≤α⁢d0t1n+1,

where α is a precised constant, 0<d0<1 and t1 is the positive root of the equation t2−t−p.

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)n≥0 generated by (2) converge to this solution.

Observe firstly that the identity:

(4) xn−[xn−1,xn;f]−1⁢f⁢(xn)=xn−1−[xn−1,xn;f]−1⁢f⁢(xn−1)

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

The following identity

(5) f⁢(xn+1)= f⁢(xn)+[xn−1,xn;f]⁢(xn+1−xn)
+([xn,xn+1;f]−[xn−1,xn;f])⁢(xn+1−xn),n=1,2,…

Let B>0,α>0,0<d0<1, and x0,x∈X1. Consider the sphere

(6) U={x∈X1:‖x−x0‖≤B⁢α⁢d01−d0t1−1}

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

(7) t2−t−p=0

The following theorem holds:

Theorem 1.

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

  • (1)

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

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

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

    ‖x∗−xn‖≤B⁢α⁢d0t1n1−d0t1n⁢(t1−1).
Proof.

From (2) for n=1 we deduce:

‖x2−x1‖≤B⁢‖f⁢(x1)‖≤B⁢α⁢d0t1

from which, taking also into account iii. it follows

‖x2−x0‖ ≤‖x2−x1‖+‖x1−x0‖
≤B⁢α⁢d0t1+B⁢α⁢d0
≤B⁢α⁢d0⁢(1+d0t1−1)
<B⁢α⁢d01−d0t1−1

from which it results that x2∈U.

Using the fact that x2∈U, the identities (4) and (5), and the inequality a), we obtain

‖f⁢(x2)‖ ≤‖x2−x1‖⁢(l1⁢‖x2−x0‖p+l2⁢‖x1−x0‖p+l3⁢‖x2−x1‖p)
≤B⁢‖f⁢(x1)‖⁢(l1⁢Bp⁢‖f⁢(x0)‖p+l2⁢‖x1−x0‖p+l3⁢Bp⁢‖f⁢(x1)‖p)
≤B⁢α⁢d0t1⁢(l1⁢Bp⁢αp⁢d0p+l2⁢Bp⁢αp⁢d0p+l3⁢Bp⁢αp⁢d0t1)
≤Bp+1⁢αp+1⁢d0t1+p⁢(l1+l2+l3⁢d0p⁢(t1−1))
=Bp+1⁢αp+1⁢(l1+l2+l3⁢d0p⁢(t1−1))⁢d0t12−p≤α⁢d0t12

since αp⁢Bp+1⁢(l1+l2+l3⁢d0p⁢(t1−1))≤αp⁢Bp+1⁢(l1+l2+l3)<1.

From the above inequality follows therefore:

‖f⁢(x2)‖≤α⁢d0t12

Suppose by induction that:

  • (a’)

    xi∈U,i=0,1,…,k;

  • (b’)

    ‖f⁢(xi)‖≤α⁢d0t1i,i=1,2,…,k.

Then, for xk+1we have:

‖xk+1−x0‖≤‖xk+1−xk‖+‖xk−xk−1‖+…+‖x1−x0‖
≤B⁢‖f⁢(xk)‖+B⁢‖f⁢(xk−1)‖+…+B⁢α⁢d0
≤B⁢α⁢d0t1k+B⁢α⁢d0t1k−1+…+B⁢α⁢d0
=B⁢α⁢d0⁢(1+d0t1−1+d0t12−1+…+d0t1k−1)
≤B⁢α⁢d0⁢(1+d0t1−1+d02⁢(t1−1)+…+d0k⁢(t1−1))≤B⁢α⁢d01−d0t1−1

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

‖f⁢(xk+1)‖ ≤Bp+1⁢αp+1⁢(l1+l2+l3⁢dop⁢t1k−1⁢(t1−1))⁢d0t1k−1⁢(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)n≥0 is fundamental.

Indeed, for every n,s∈ℕ we have:

‖xn+s−xn‖ ≤∑k=nn+s−1‖xk+1−xk‖≤∑k=nn+s−1B⁢‖f⁢(xk)‖≤B⁢α⁢∑k=nn+s−1d0t1k
=B⁢α⁢d0t1n⁢∑k=nn+s−1d0t1k−t1n=B⁢α⁢d0t1n⁢∑k=nn−s−1d0t1n⁢(t1k−n−1)
≤B⁢α⁢d0t1n⁢∑k=nn+s−1d0t1n⁢(k−n)⁢(t1−1)=B⁢α⁢d0t1n⁢∑k=nn+s−1(d0t1n⁢(t1−1))k−n
≤B⁢α⁢d0t1n1−d0t1n⁢(t1−1).

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

‖xn+s−xn‖≤B⁢α⁢d0t1n1−d0t1n⁢(t1−1)

follows the inequality:

‖x∗−xn‖≤B⁢α⁢d0t1n1−d0t1n⁢(t1−1).

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