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:X→Xis 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+1→Xwhich 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 eachx∈X.

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

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

Orx0,x1,…,xk∈Xare given elements.

Let us designate byD⊂Xa 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+1ai⁢r⁢(ini,ini),

Orai∈ℝ,ai≥0,i=1,2,…,k+1,for each (in1,in2,…,ink+1), (in1,in2,…,ink+1)∈Dk+1.

Consider the equation:

(5) a1+a2⁢t+…+ak+1⁢tk=tk+1.

We know that in the case whereai≥0  for eachi=1,2,…,k+1 and where there exists at least one numberi1,1≤i1≤k+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=0∞provided 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={x∈X:r⁢(x,x0)≤a1−t1}⊆D,Ort1is the positive root of equation ( 5 ) and

    a≥max⁡{r⁢(x0,x1),1t1⁢r⁢(x1,x2),…,1t1k−1⁢r⁢(xk−1,xk)};

then we have the following properties:

  • j)

    xn∈Sfor eachn=0,1,…;

  • jj)

    the sequel(xn)n=0∞is convergent and if we denote byx∗the limit of the sequence(xn)n=0∞, SOx∗∈SAndx∗is the unique solution of equation ( 1 ) of the sphereS;

  • jjj)

    we have the following inequalities:

    r⁢(x∗,xn)≤a⁢t1n1−t1,for each ⁢n=0,1,…
Demonstration.

From conclusion ii) it follows thatr⁢(xi,xi+1)≤a⁢t1ifor eachi=0,1,…,k−1. Fori=1,2,…,kwe have the following inequalities:

(8) r⁢(xi,x0) ≤r⁢(x0,x1)+r⁢(x1,x2)+…+r⁢(xi−1,xi)
≤a+a⁢t1+…+a⁢t1i−1
<a1−t1

from which we deduce thatxi∈S,Fori=1,2,…,k.

Subsequently we assume thatxi∈SAndr⁢(xi−1,xi)≤a⁢ti−1,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=1−k+1ai⁢r⁢(xn+i−1,xn+i)≤a⁢(a1⁢t1n+…+ak+1⁢t1k+n)=a⁢t1n+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+2∈S. 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+s−1,xn+s)
≤a⁢(t1n+t1n+1+…+t1n+s−1)≤a⁢t1n1−t1

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

Let us designate byx∗=limn→∞xn. So by passing to limit in the inequalities ( 9 ) with s→∞, on a:

(10) r⁢(xn,x∗)≤a⁢t1n1−t1,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+1ai⁢r⁢(xn+i−1,x∗)
≤a⁢t1n+k+11−t1+∑i=1k+1ai⁢a⁢t1n+i−11−t1
=a1−t1⁢[t1n+k+1+t1n⁢∑i=1k+1ai⁢t1i−1]
=2⁢a⁢t1n+k+11−t1.

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 elementx∗⊂Sis the solution to equation ( 1 ). ∎

Let us now show thatx∗is the unique solution of equation ( 1 ) which belongs to the sphereS. Suppose instead that equation ( 1 ) has at least two solutionsx∗Andand∗ 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+1ai⋅r⁢(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+1→Xwhich 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=0∞provided 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∗,…,xk∗such that the following conditions are met:

(13) r⁢(x0∗,x1∗) ≤a+2⁢e1−b
r⁢(x1∗,x2∗) ≤a⁢t1+2⁢e1−b
…⁢…⁢…⁢…⁢…⁢…..
r⁢(xk−1∗,xk∗) ≤a⁢t1k−1+2⁢e1−b.

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

r⁢(xk∗,xk+1∗)≤a⁢t1k+2⁢e1−b.

Let us also assume that the set:

(14) S∗={x∈X:r⁢(x,x0)≤a⁢(2−t1)1−t1+e1−b}

is contained inD.

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

(15) r⁢(x0,x0∗) ≤a+e1−b
r⁢(x1,x1∗) ≤a⁢t1+e1−b
…⁢…⁢…⁢…⁢…⁢…⁢…
r⁢(xk,xk∗) ≤a⁢t1k+e1−b

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

(16) r⁢(xn+k+1,xn+k+1∗)≤a⁢t1n+k+1+e1−b
(17) r⁢(xn+k+1∗,xn+k+2∗)≤a⁢t1n+k+1+2⁢e1−b

And

(18) r⁢(x∗,xn+k+1∗)≤a⁢(2−t1)1−t1⁢t1n+k+1+e1−b.

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+a1⁢r⁢(x0,x0∗)+a2⁢r⁢(x1,x1∗)+…⁢ak+1⁢r⁢(xk,xk∗)
≤e+e1−b⁢(a1+…+ak+1)+a⁢(a1+a2⁢t1+…+ak+1⁢t1k)
=e1−b+a⁢t1k+1

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

(20) r⁢(xk+1∗,x0) ≤r⁢(xk+1∗,xk+1)+r⁢(xk+1,x0)
≤a1−t1+e1−b+a⁢t1k+1≤a+a1−t1+e1−b
=a⁢(2−t1)1−t1+e1−b.

Generally, if we assume that

xn∗,xn+1∗,…,xn+k∗∈S∗

and that we have the following inequalities:

r⁢(xn+i,xn+1∗)≤a⁢t1n+1+e1−b

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+a1⁢r⁢(xn,xn∗)+…+ak+1⁢r⁢(xn+k,xn+k∗)
≤e+a⁢t1n⁢(a1+a2⁢t1+…+ak+1⁢t1k)+e⁢b1−b
e+a⁢t1n+k+1+e⁢b1−b=a⁢t1n+k+1+e1−b

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+1∗∈S∗.

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∗))
≤2⁢e+2⁢e⁢b1−b+a⁢(a1+a2⁢t1+…+ak+1⁢t1k)
=2⁢e1−b+a⁢t1k+1,

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

Now assuming that we have the following inequalities:

r⁢(xn+i∗,xn+i+1∗)≤2⁢e1−b+a⁢t1n+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∗)
≤a⁢t1n+k+11−t1+a⁢t1n+k+1+e1−b
=a⁢(1−t1)1−t1⁢t1n+k+1+e1−b

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 wholeD⊂X,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=a⁢x+b=G⁢(x),x∈ℝand the approximate equation x∗=a∗⁢x∗+b∗=G∗⁢(x∗).

Because|G⁢(x)−G∗⁢(x)|=|a−a∗|⋅|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