On approximating the solutions of equations in metric spaces

Abstract

Let \(\left( X,\rho \right)\) be a complete metric space, \(f:X\rightarrow X\) a nonlinear mapping. In order to solve the equation \(x=f\left( x\right) \) we consider a multistep method \[x_{n+k+1}=G(x_{n},x_{n+1},…,x_{n+k}), \quad n=1,2,… \] generated by a mapping \(G:X^{k+1}\rightarrow X\), whose diagonal restriction coincides with \(f\): \(G(x,…,x)=f(x)\). Under Lipschitz assumption on \(G\) we determine the algebraic equation whose unique positive solution leads to the convergence order of the iterations. We also study the case when the operator \(G\) replaced by an approximation of it.

Authors

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

Title

Original title (in French)

Sur l’approximation des racines des equations dans une espace métrique

English translation of the title

On approximating the solutions of equations in metric spaces

Keywords

multistep iterative methods; convergence; successive approximations

PDF

Cite this paper as:

I. Păvăloiu, Sur l’approximation des racines des equations dans une espace métrique, Seminar on functional analysis and numerical methods, Preprint no. 1 (1989), pp. 95-104 (in French).

About this paper

Journal

Seminar on functional analysis and numerical methods,
Preprint

Publisher Name

“Babes-Bolyai” University,
Faculty of Mathematics,
Research Seminars

DOI

Not available yet.

References

[1] I. Pavaloiu, I., Serb, Sur des methodes iteratives optimales, Research Seminars, Seminar on Functional Analysis and Numerical Methods, Preprint Nr.1 (1983), 175–182.

[2] I.A. Rus, An iterative method for the solution of the equation x = f (x, x, . . . , x), Anal. Num´er. Theor. Approx., 10 (1981), 95–100.

[3] Weinischke, J.H., Uber eine klasse von Iterationverfahren, Numerische Mathematik 6 (1964), 395–404.

Paper (preprint) in HTML form

"Babeş-Bolyai" University

Faculty of Mathematics and Physics

Research Seminars

Seminar on Functional Analysis and Numerical Methods

Preprint Nr.1, 1989, pp.95-104

On approximating the solutions of equations in metric spaces

Ion Pavaloiu

In this note we study the evaluation of errors that arise during the numerical resolution of equations in metric spaces using certain multi-step iteration methods.

Consider a metric space(X,r)complete and the following equation:

(1) x=f(x),

Orf:XXis any operator.

Let us designate byXk+1=X×X××X,that is to say the Cartesian product of the setXwith himselfk+1times.

For the resolution of equation ( 1 ) we consider the applicationG:Xk+1Xwhich we assume is restricted to the diagonal of spaceXk+1coincides with the operatorf, that's to say:

(2) G(x,x,,x)=f(x),

for eachxX.

Consider the following(xn)n=0provided by the following iteration process:

(3) xn+k+1=G(xn,xn+1,,xn+k),n=0,1,2,

Orx0,x1,,xkXare given elements.

Let us designate byDXa bounded set in this space. We assume that on the setDk+1the operatorGverifies the Lipschitz condition, that is:

(4) r(G(in1,in2,,ink+1),G(in1,in2,,ink+1))i=1k+1air(ini,ini),

Orai,ai0,i=1,2,,k+1,for each (in1,in2,,ink+1), (in1,in2,,ink+1)Dk+1.

Consider the equation:

(5) a1+a2t++ak+1tk=tk+1.

We know that in the case whereai0  for eachi=1,2,,k+1 and where there exists at least one numberi1,1i1k+1for which ai1>0, this equation has only one roott1real and positive.

If in addition:

(6) b=i=1k+1ai<1

then this root verifies the inequality

(7) 0<t1<1.

Regarding the convergence of the sequence(xn)n=0provided by method ( 3 ) we have the following theorem:

Theorem 1 .

If the applicationGand the initial elements x0,x1,,xkmeet the following conditions:

  • i)

    L’application Gfulfills condition ( 4 ) where ai,i=1,2,,k+1, satisfy the relation ( 6 );

  • ii)

    The wholeS={xX:r(x,x0)a1t1}D,Ort1is the positive root of equation ( 5 ) and

    amax{r(x0,x1),1t1r(x1,x2),,1t1k1r(xk1,xk)};

then we have the following properties:

  • j)

    xnSfor eachn=0,1,;

  • jj)

    the sequel(xn)n=0is convergent and if we denote byxthe limit of the sequence(xn)n=0, SOxSAndxis the unique solution of equation ( 1 ) of the sphereS;

  • jjj)

    we have the following inequalities:

    r(x,xn)at1n1t1,for each n=0,1,
Demonstration.

From conclusion ii) it follows thatr(xi,xi+1)at1ifor eachi=0,1,,k1. Fori=1,2,,kwe have the following inequalities:

(8) r(xi,x0) r(x0,x1)+r(x1,x2)++r(xi1,xi)
a+at1++at1i1
<a1t1

from which we deduce thatxiS,Fori=1,2,,k.

Subsequently we assume thatxiSAndr(xi1,xi)ati1,for eachi=1,2,,n+k+1. From this inequality we deduce that:

r(xn+k+2,xn+k+1)r(G(xn,xn+1,,xn+k),G(xn+1,,xn+k+1))
i=1k+1air(xn+i1,xn+i)a(a1t1n++ak+1t1k+n)=at1n+k+1,

where we used the fact thatt1verifies equation ( 5 ).

From the above inequalities, and in a manner similar to that used in the demonstration of relation ( 8 ), we deduce thatxn+k+2S. From the principle of mathematical induction it therefore follows that all the elements of the sequence(xn)n=0 belong to the sphereS.

On the completeness of spaceXand inequalities

(9) r(xn+s,xn) r(xn,xn+1)++r(xn+s1,xn+s)
a(t1n+t1n+1++t1n+s1)at1n1t1

for eachs=1,2,,it follows that the following(xn)n=0is convergent.

Let us designate byx=limnxn. So by passing to limit in the inequalities ( 9 ) with s, on a:

(10) r(xn,x)at1n1t1,n=0,1,

Let us now demonstrate thatG(x,x,,x)=x.

In fact we have:

r(x,G(x,x,,x))
r(xn+k+1,x)+r(xn+k+1,G(x,x,,x))
r(xn+k+1,x)+r(G(xn,xn+1,,xn+k),G(x,x,,x))
r(xn+k+1,x)+i=1k+1air(xn+i1,x)
at1n+k+11t1+i=1k+1aiat1n+i11t1
=a1t1[t1n+k+1+t1ni=1k+1ait1i1]
=2at1n+k+11t1.

This inequality is verified for each n, from which it follows that we have the equalityx=G(x,x,,x), but this relation and relation ( 2 ) imply that the elementxSis the solution to equation ( 1 ). ∎

Let us now show thatxis the unique solution of equation ( 1 ) which belongs to the sphereS. Suppose instead that equation ( 1 ) has at least two solutionsxAndand inS. So fromx=f(x)Andand=f(and), and from ( 2 ) it follows that:

r(x,and) =r(G(x,x,,x),G(and,and,,and))
i=1k+1air(x,and)

from which, if we take into account ( 6 ) we obtain

0<r(x,and)<r(x,and).

This last inequality being impossible, it follows that the assumption is false, therefore that equation ( 1 ) has a unique solution.

Let us then consider the applicationG:Xk+1Xwhich fills with the applicationGthe following condition:

(11) r(G(in1,in2,,ink+1),G(in1,in2,,ink+1))e,e>0

for each(in1,in2,,ink+1)Dk+1.

To solve equation ( 1 ) we will replace the sequence(xn)n=0provided by the relation ( 3 ), subsequently(xn)n=0 provided by the following iterative process:

(12) xn+k+1=G(xn,xn+1,,xn+k),n=0,1,,

where we assume that the initial elements have been chosenx0,x1,,xksuch that the following conditions are met:

(13) r(x0,x1) a+2e1b
r(x1,x2) at1+2e1b
..
r(xk1,xk) at1k1+2e1b.

Regarding the relationshipxk+1=G(x0,x1,,xk)we assume that we have the following condition

r(xk,xk+1)at1k+2e1b.

Let us also assume that the set:

(14) S={xX:r(x,x0)a(2t1)1t1+e1b}

is contained inD.

Furthermore we assume that the initial elementsx0,x1,,xkof method ( 12 ) fill with the elementsx0,x1,,xkof ( 3 ) the following conditions:

(15) r(x0,x0) a+e1b
r(x1,x1) at1+e1b
r(xk,xk) at1k+e1b

Using these assumptions and further using the fact thatG verifies the hypotheses of Theorem 1 , we will demonstrate that for eachnwe have the following inequalities:

(16) r(xn+k+1,xn+k+1)at1n+k+1+e1b
(17) r(xn+k+1,xn+k+2)at1n+k+1+2e1b

And

(18) r(x,xn+k+1)a(2t1)1t1t1n+k+1+e1b.

Indeed, by using relations ( 15 ) we obtain the following relations:

(19) r(xk+1,xk+1)=
=r(G(x0,x1,,xk),G(x0,x1,,xk))
r(G(x0,x1,,xk),G(x0,x1,,xk))+
+r(G(x0,x1,,xk),G(x0,,xk))
e+a1r(x0,x0)+a2r(x1,x1)+ak+1r(xk,xk)
e+e1b(a1++ak+1)+a(a1+a2t1++ak+1t1k)
=e1b+at1k+1

that is to say that inequality ( 16 ) is verified, forn=0.Let us now show thatxk+1S. Using inequalities ( 19 ) and ( 8 ) we obtain:

(20) r(xk+1,x0) r(xk+1,xk+1)+r(xk+1,x0)
a1t1+e1b+at1k+1a+a1t1+e1b
=a(2t1)1t1+e1b.

Generally, if we assume that

xn,xn+1,,xn+kS

and that we have the following inequalities:

r(xn+i,xn+1)at1n+1+e1b

for eachi=1,2,,k,then in a manner similar to that used to demonstrate inequality ( 19 ) we obtain:

(21) r(xn+k+1,xn+k+1) e+a1r(xn,xn)++ak+1r(xn+k,xn+k)
e+at1n(a1+a2t1++ak+1t1k)+eb1b
e+at1n+k+1+eb1b=at1n+k+1+e1b

from which it follows that we have the inequalities ( 16 ) for eachn. From relations ( 21 ) and ( 8 ), in a manner similar to that used in the deduction of relation ( 20 ), it follows thatxn+k+1S.

Let us now demonstrate that the inequalities ( 17 ) are verified. If we consider the relations ( 13 ) we have the following inequalities:

(22) r(xk+1,xk+2) r(G(x0,x1,,xk),G(x1,,xk+1))
r(G(x0,x1,,xk),G(x0,x1,xk))
+r(G(x0,x1,,xk),G(x1,x2,,xk+1))
+r(G(x1,x2,,xk+1),G(x1,x2,,xk+1))
2e+2eb1b+a(a1+a2t1++ak+1t1k)
=2e1b+at1k+1,

so we obtained the inequality ( 17 ) forn=0.

Now assuming that we have the following inequalities:

r(xn+i,xn+i+1)2e1b+at1n+i,

Fori=0,1,,k,in a manner similar to that employed in the deduction of inequality ( 22 ) we obtain inequalities ( 17 ).

For inequalities ( 18 ) we have

r(x,xn+k+1) r(x,xn+k+1)+r(xn+k+1,xn+k+1)
at1n+k+11t1+at1n+k+1+e1b
=a(1t1)1t1t1n+k+1+e1b

that is, inequalities ( 18 ).

Inequalities ( 18 ) provide an assessment of the distance between the exact solution of equation ( 1 ) and its approximation, obtained using the iteration method ( 12 ).

Noticed .

If we notice that the wholeDX,on which condition ( 4 ) is fulfilled, is itself a metric space, it follows that we can consider theorem 1 as a special case of the main result contained in the work [ 3 ] . From our point of view the importance of this theorem consists in the fact that the setDbeing limited it allows easy choice of applicationG, application that fills with the applicationGthe condition ( 11 ). Let us imagine for example the trivial case of the equation x=ax+b=G(x),xand the approximate equation x=ax+b=G(x).

Because|G(x)G(x)|=|aa||x|, it follows that the difference|G(x)G(x)|is not bounded, so even in the case0<a<1condition ( 11 ) cannot be fulfilled in the whole space.

Bibliography


Institute of Computing

Post Office 1

C.P. 68

3400 Cluj-Napoca

Romania



This paper is in final form and no version of it is or will be submitted for publication elsewhere.

1989

Related Posts