On approximating the solutions of equations by the chord method and a method
of Aitken-Steffensen type\(^{\ast }\)
January 14, 2013.
\(^{\S }\)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 \(\tfrac {1+\sqrt{5}}{2}.\) 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 \(\left\Vert \cdot \right\Vert _{X}:X\rightarrow \mathbb {R}\) and \(\left\Vert \cdot \right\Vert _{Y}:Y\rightarrow \mathbb {R}\) their norms respectively, and by \(\theta _{X}\) and \(\theta _{Y}\) their null elements respectively. By \(\left( X,Y\right) ^{\ast }\) we denote the set of the linear and continuous mappings defined from \(X\) to \(Y.\) The set \(\left( X,Y\right) ^{\ast }\) is a linear normed space as well, if we define the norm \(\left\Vert \cdot \right\Vert :\left( X,Y\right) ^{\ast }\rightarrow \left[ 0,+\infty \right[ ,\) by \(\left\Vert U\right\Vert =\sup \nolimits _{h\in X,\ \left\Vert h\right\Vert _{X}=1}\left\Vert U\left( h\right) \right\Vert \) for any \(U\in \left( X,Y\right) ^{\ast }.\) For the case of \(Y=\mathbb {R}\) we denote by \(X^{\ast }\) the set \(\left( X,\mathbb {R}\right) ^{\ast },\) this set representing the space of real, linear and continuous functionals defined on the linear normed space \(X.\)
Let us consider now a set \(D\subseteq X\) and a nonlinear mapping \(f:D\rightarrow Y.\) Using this mapping we have the equation:
We will study the approximation of its solutions.
In order to clarify the aforementioned notions we have the following definition:
Considering the nonlinear mapping \(f:D\rightarrow Y\) together with the points \(x,y\in D,\) \(x\neq y,\) any mapping \(\Gamma _{f;x,y}\in \left( X,Y\right) ^{\ast }\) that verifies the equality:
is called generalized abstract divided difference of the function \(f:D\rightarrow Y\) at the points \(x,y.\)
In connection with the previous definition we have the following remark:
a) If we consider the theorem according to which in every linear normed space \(\left( X,\left\Vert \cdot \right\Vert _{X}\right) ,\) for any \(a\in X\diagdown \left\{ \theta _{X}\right\} \) there exists a linear and continuous functional \(u\in X^{\ast }\) such that \(\left\Vert u\right\Vert =1\) and \(u\left( a\right) =\left\Vert a\right\Vert _{X}.\) Therefore, for any \(x,y\in X\) with \(x\neq y\) there exists the functional \(U_{xy}\in X^{\ast }\) such that \(\left\Vert U_{xy}\right\Vert =1\) and \(U_{xy}\left( x-y\right) =\left\Vert x-y\right\Vert _{X}.\) At the same time there exists the functional \(U_{yx}\in X^{\ast }\) such that \(\left\Vert U_{yx}\right\Vert =1\) and \(U_{yx}\left( y-x\right) =\left\Vert y-x\right\Vert _{X}\) as well. In the paper [ 13 ] there appears the mapping \(\left[ x,y;f\right] \in \left( X,Y\right) ^{\ast },\) defined by the equality:
for any \(h\in X.\)
This mapping verifies the equality (2) and it is called abstract divided difference of the nonlinear mapping \(f:D\rightarrow Y\) at the points \(x,y\in D\) with \(x\neq y.\) 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 \(\left\langle \cdot |\cdot \right\rangle :X\times X\rightarrow \mathbb {R}.\) Defining \(\left\Vert \cdot \right\Vert _{X}:X\rightarrow \mathbb {R}\) by \(\left\Vert x\right\Vert _{X}=\sqrt{\left\langle x|x\right\rangle },\) the space \(\left( X,\left\Vert \cdot \right\Vert _{X}\right) \) is a linear normed space.
For any \(x,y\in X\) with \(x\neq y\) the functional \(U_{xy}\in X^{\ast }\) from a) will be defined by:
for any \(h\in X.\) So, for the same elements \(x,y\in X\) with \(x\neq y\) we have that the abstract divided difference \(\left[ x,y;f\right] \in \left( X,Y\right) ^{\ast }\) is defined by:
for any \(h\in X.\)
â–¡
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 \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) (here \(\mathbb {N}^{\ast }=\mathbb {N}\cup \left\{ 0\right\} \) and this notation remains valid for the rest of this paper) and that for any \(n\in \mathbb {N}\) verifies the equality:
If we suppose the existence for a certain \(n\in \mathbb {N}\) of the mapping \(\Gamma _{f;x_{n-1},x_{n}}^{-1}\in \left( Y,X\right) ^{\ast }\) representing the inverse of the mapping \(\Gamma _{f;x_{n-1},x_{n}}\in \left( X,Y\right) ^{\ast },\) for this number \(n\in \mathbb {N}\) the equality (4) is equivalent to:
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 \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) to a solution of the equation (1) the existence of which is also proved, is the following theorem:
We suppose that the following assumptions hold:
\(\left( X,\left\Vert \cdot \right\Vert _{X}\right) \) is a Banach space;
the nonlinear mapping \(f:D\rightarrow Y,\) for any \(x,y\in D\) with \(x\neq y\) admits a generalized abstract divided difference \(\Gamma _{f;x_{n-1},x_{n}}\in \left( X,Y\right) ^{\ast }\) and there exists a number \(L{\gt}0\) such that for any \(x,y,z\in D\) with \(x\neq y\) and \(y\neq z\) we have the following inequality:
\begin{equation} \left\Vert \Gamma _{f;x,y}-\Gamma _{f;y,z}\right\Vert \leq L\left\Vert x-z\right\Vert _{X}; \label{f6}\end{equation}6the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) is such that for any \(n\in \mathbb {N}\) we have \(x_{n-1}\neq x_{n}\) and the equality (4) is verified;
referring to the initial elements \(x_{0},x_{1}\in D\) of the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) we suppose the carrying out of the following conditions:
the mapping \(\Gamma _{f;x_{0},x_{1}}\in \left( X,Y\right) ^{\ast }\) is invertible and \(\Gamma _{f;x_{0},x_{1}}^{-1}\! \in \left( Y,X\right) ^{\ast };\)
there exist the numbers \(h_{0}\in \left] 0,1\right[ \) and \(B_{0}{\gt}0\) such that we have the following inequality:
\begin{equation} \big\Vert \Gamma _{f;x_{0},x_{1}}^{-1}\big\Vert \leq \tfrac {B_{0}}{1-h_{0}}; \label{f7}\end{equation}7if we note \(R_{0}=\left\Vert f\left( x_{0}\right) \right\Vert _{Y}\) and \(h_{1}=\tfrac {LB_{0}^{2}R_{0}}{\left( 1-h_{0}\right) ^{2}}\) there exists a number \(q\in \left[ \tfrac {1}{2},1\right[ \) such that we have the following inequality:
\begin{equation} d=\max \bigg\{ \tfrac {h_{0}}{\left( 1-q\right) ^{2}},\big[ \tfrac {h_{1}}{\left( 1-q\right) ^{2}}\big] ^{\frac{1}{\alpha }}\bigg\} <1, \label{f8}\end{equation}8where \(\alpha =\tfrac {1+\sqrt{5}}{2};\)
if \(\delta =\tfrac {\left( 1-q\right) ^{2}}{LB_{0}}\cdot \tfrac {d^{\alpha }}{1-d}\) and \(S\left( x_{0},\delta \right) =\left\{ \left. x\in X\right/ \left\Vert x-x_{0}\right\Vert _{X}\leq \delta \right\} \) the relation \(x_{1}\in S\left( x_{0},\delta \right) \subseteq D\) is true.
Then the following conclusions are true:
for any \(n\in \mathbb {N}\) we have that \(x_{n}\in S\left( x_{0},\delta \right) ,\) there exists the mapping \(\Gamma _{f;x_{n-1},x_{n}}^{-1}\in \left( Y,X\right) ^{\ast },\) and the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\) verifies the equality (5);
the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq X\) is convergent;
if \(\overline{x}=\lim \limits _{n\rightarrow \infty }x_{n},\) then \(\overline{x}\in S\left( x_{0},\delta \right) \) and \(f\left( \overline{x}\right) =\theta _{Y};\)
for any \(n\in \mathbb {N}^{\ast }\) the following estimate of the error of approximation is true:
\begin{equation} \left\Vert \overline{x}-x_{n}\right\Vert _{X}\leq \tfrac {\left( 1-q\right) ^{2}}{LB_{0}}\cdot \frac{d^{\alpha ^{n+1}}}{1-d^{\alpha ^{n}}}. \label{f9}\end{equation}9
The proof of this theorem is given in [ 13 ] .
Regarding this theorem we have the following remark:
From the conclusions of the theorem 2.1, corroborated with the fact that \(d{\lt}1,\) we deduce that there exists a number \(N\in \mathbb {N}\) such that for any number \(n\in \mathbb {N},\) \(n\geq N\) we have the following inequality:
This inequality indicates that the convergence order of the iterative method of the chord is \(\alpha =\tfrac {1+\sqrt{5}}{2}.\)â–¡
The inequality (10) is obvious from (9).
3 Certain remarks in connection with the previous result
We have the following statements:
Under the hypotheses of Theorem 2.1 the sequence of real and positive numbers \(\big( \big\Vert \Gamma _{f;x_{n-1},x_{n}}^{-1}\big\Vert \big) _{n\in \mathbb {N}}\) has an upper bound, for any \(n\in \mathbb {N}\) the following inequality taking place:
â–¡
Indeed, as \(\Vert \Gamma _{f;x_{n-1},x_{n}}^{-1}\Vert \leq B_{n},\) and using the recurrence relation of the sequence \(\left( B_{n}\right) _{n\in \mathbb {N}^{\ast }},\) we have \(\tfrac {B_{n}}{B_{n-1}}\leq \tfrac {1}{1-h_{n-1}}\) and so:
Evidently:
But:
therefore:
from where the inequality (11) is evident.
If in the hypotheses of Theorem 2.1 we choose \(q=\tfrac {1}{2},\) the values of the other constants from this theorem are:
and the inequality that expressed an upper bound of the error by which \(x_{n}\) approximates \(\overline{x}\) is:
For the upper bound of the sequence \(\big( \big\Vert \Gamma _{f;x_{n-1},x_{n}}^{-1}\big\Vert \big) _{n\in \mathbb {N}}\) for any \(n\in \mathbb {N},\) we have the following inequality:
â–¡
The statements are obvious:
For any \(n\in \mathbb {N}\) there exists the generalized abstract divided difference \(\Gamma _{f;x_{n},\overline{x}}\in \left( X,Y\right) ^{\ast },\) this mapping is invertible, so the mapping \(\Gamma _{f;x_{n},\overline{x}}^{-1}\in \left( Y,X\right) ^{\ast }\) exists and the sequence \(\big( \big\Vert \Gamma _{f;x_{n},\overline{x}}^{-1}\big\Vert \big) _{n\in \mathbb {N}}\) has an upper bound, more precisely, there exists a number \(p\in \mathbb {N}\) such that for any \(n\in \mathbb {N},\ n\geq p\) we have the inequality:
â–¡
It is obvious that \(\overline{x}\neq x_{n}\) for any \(n\in \mathbb {N}.\) Indeed, if there exists a number \(n_{0}\in \mathbb {N}\) such that \(\overline{x}=x_{n_{0}}\) then:
which is impossible.
As \(\overline{x}\neq x_{n}\) for any \(n\in \mathbb {N}\) we deduce the existence of the mapping \(\Gamma _{f;x_{n},\overline{x}}\in \left( X,Y\right) ^{\ast }\) representing the generalized abstract divided difference of the function \(f:D\rightarrow Y\) on the points \(x_{n}\) and \(\overline{x}.\)
In order to prove the invertibility of the mapping \(\Gamma _{f;x_{n},\overline{x}}\in \left( X,Y\right) ^{\ast }\) for any any \(n\in \mathbb {N}\) let us consider the following expression:
We have:
It is clear that:
As \(\lim \limits _{n\rightarrow \infty }\left\Vert \overline{x}-x_{n-1}\right\Vert _{X}=0,\) we deduce that there exists a number \(p\in \mathbb {N}\) such that for any \(n\in \mathbb {N},\ n\geq p,\) we have the inequality:
therefore \(\left\Vert W_{n}\right\Vert \leq \tfrac {1}{2}{\lt}1\) and so there exists the mapping \(\left( \mathbf{I}_{X}-W_{n}\right) ^{-1}\in \left( X,X\right) ^{\ast }\) and:
From the existence of the mapping \(\Gamma _{f;x_{n-1},x_{n}}^{-1}\in \left( Y,X\right) ^{\ast }\) and using the inequality \(\big\Vert \Gamma _{f;x_{n-1},x_{n}}^{-1}\big\Vert \leq B_{n}\) we deduce the existence of the mapping \(\Gamma _{f;x_{n},\overline{x}}^{-1}\in \left( Y,X\right) ^{\ast }\) by the equality:
again:
The sequence \(\left( \left\Vert \Gamma _{f;x_{n},\overline{x}}\right\Vert \right) _{n\in \mathbb {N}}\) also has an upper bound and for any \(n\in \mathbb {N}\) we have the inequality:
â–¡
Indeed, for any \(n\in \mathbb {N},\) we have:
There exists a number \(N\in \mathbb {N}\) such that for the same values \(n\in \mathbb {N}\) with \(n\geq N\) we have that:
and this inequality proves that the convergence order of this method is at least \(\alpha =\tfrac {1+\sqrt{5}}{2}.\) â–¡
Indeed, as \(\alpha ^{n+1}=\alpha ^{n}+\alpha ^{n-1}\) it is clear that \(\tfrac {d^{\alpha ^{n+1}}}{1-d^{\alpha ^{n}}}=d^{\alpha ^{n}}\cdot \tfrac {d^{\alpha ^{n-1}}}{1-d^{\alpha ^{n}}}\) and from the fact that \(\lim \limits _{n\rightarrow \infty }\tfrac {d^{\alpha ^{n-1}}}{1-d^{\alpha ^{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 \(\alpha =\tfrac {1+\sqrt{5}}{2}.\)
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:X\rightarrow X\) that verifies the inclusion \(Q\left( D\right) \subseteq D\) for a set \(D\subseteq X,\) generates the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D,\) starting from an arbitrary \(x_{0}\in D,\) by the verification for any \(n\in \mathbb {N}^{\ast }\) of the equality:
If for any \(n\in \mathbb {N}^{\ast }\) there exists the mapping \(\Gamma _{f;x_{n},Q\left( x_{n}\right) }^{-1}\in \left( Y,X\right) ^{\ast }\) the last equality is equivalent to:
A more general case is the one which use two mappings \(Q_{1},Q_{2}:X\rightarrow X\) which verify for \(i\in \left\{ 1,2\right\} \) the relations \(Q_{i}\left( D\right) \subseteq D.\) Starting from an arbitrary \(x_{0}\in D\) one build the sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) by the verification for any \(n\in \mathbb {N}^{\ast }\) of the equality:
equality which in the hypothesis of the existence of \(\Gamma _{f;Q_{1}\left( x_{n}\right) ,Q_{2}\left( x_{n}\right) }^{-1}\in \left( Y,X\right) ^{\ast }\) is equivalent to:
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:D\rightarrow Y\) together with \(Q_{1},Q_{2}:X\rightarrow X\) and one requests the verification of these conditions on every point of a set \(D\subseteq X.\)
We propose a more general frame and we will build an iterative process after as follows.
Let us consider a initial element \(x_{0}\in D.\) Besides the main sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\) \(\subseteq D\) we also use two auxiliary sequences \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D.\)
For these auxiliary sequences we request the existence of the numbers \(K_{1},\) \(K_{2},p,q{\gt}0\) such that for any \(n\in \mathbb {N}^{\ast }\) the following inequalities are verified:
Then, if for a number \(n\in \mathbb {N}^{\ast }\) we have available the elements \(y_{n},z_{n}\in D\) starting from \(x_{n}\in D,\) we will generate the new iterate \(x_{n+1}\in D\) by the following relation:
On account of the property of definition of the mapping \(\Gamma _{f;y_{n},z_{n}}\in \left( X,Y\right) ^{\ast }\) the equality (23) is equivalent with:
If for any \(n\in \mathbb {N}^{\ast }\) there exists the mapping \(\Gamma _{f;y_{n},z_{n}}^{-1}\in \left( Y,X\right) ^{\ast }\) we have:
In connection to the main sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\) and the auxiliary sequences \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) we have the following remarks. Here for a number \(k\in \mathbb {N}\) we denote by \(\left( X^{\left( k\right) },Y\right) ^{\ast }\) the set of mappings defined from \(X^{k}\) to \(Y\) that are \(k\)-linear and continuous. This set is a linear normed space as well with the norm \(\left\Vert \cdot \right\Vert :\left( X^{\left( k\right) },Y\right) ^{\ast }\rightarrow Y\) is defined for \(U\in \left( X^{\left( k\right) },Y\right) ^{\ast }\) by \(\left\Vert U\right\Vert =\sup \limits _{h_{i}\in X;\ \left\Vert h_{i}\right\Vert =1;\ i=\overline{1,k}}\left\Vert U\left( h_{1},...,h_{k}\right) \right\Vert _{Y}.\)
If the mapping \(f:D\rightarrow Y\) admits Fréchet differentials up to the order \(s-1,\) where \(s=\max \left\{ p,q\right\} \) and \(p,q\in \mathbb {N},\) the mapping \(f^{\left( s-1\right) }:D\rightarrow \left( X^{\left( s-1\right) },Y\right) ^{\ast }\) is a Lipschitz mapping, namely there exist the constants \(M^{\ast }{\gt}0\) such that for any \(x,y\in D\) the inequality:
is verified and for any \(n\in \mathbb {N}^{\ast }\) the inequalities:
then:
the mappings \(f^{\left( p-1\right) }\! :\! D\rightarrow \left( X^{\left( p-1\right) },Y\right) ^{\ast }\) and \(f^{\left( q-1\right) }:D\rightarrow \left( X^{\left( q-1\right) },Y\right) ^{\ast }\) are Lipschitz mappings, namely there exist the constants \(M_{1},M_{2}{\gt}0\) such that for any \(x,y\in D\) the following inequalities are true:
\begin{align} \Vert f^{\left( p-1\right) }\left( x\right) -f^{\left( p-1\right) }\left( y\right) \Vert & \leq M_{1}\Vert x-y\Vert _{X},\label{f49p}\\ \Vert f^{\left( q-1\right) }\left( x\right) -f^{\left( q-1\right) }\left( y\right) \Vert & \leq M_{2}\Vert x-y\Vert _{X};\nonumber \end{align}the sequences \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) verify the equalities (22) with:
\[ K_{1}=\alpha +\tfrac {M_{1}a^{p}}{p!},\ K_{2}=\beta +\tfrac {M_{2}b^{q}}{q!}. \]
â–¡
Indeed, from the hypothesis (26) we deduce that for any \(k\in \mathbb {N},\ k\leq s-1\) the mapping \(f^{\left( k\right) }:D\rightarrow \left( X^{\left( k\right) },Y\right) ^{\ast }\) is a Lipschitz mapping as well (we can use the well known theorem of Lagrange), particularly, the mappings \(f^{\left( p-1\right) }\) \(:D\rightarrow \left( X^{\left( p-1\right) },Y\right) ^{\ast }\) and \(f^{\left( q-1\right) }:D\rightarrow \left( X^{\left( q-1\right) },Y\right) ^{\ast }\) 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,y\in D\) the following inequalities are also verified:
(we can use Taylors’ formula for mappings).
On account of the fact that \(x_{n},y_{n}\) and \(z_{n}\in D,\) using the inequalities (27) and (29) we deduce that for any \(n\in \mathbb {N}^{\ast }\) the following inequalities are verified:
and similarly:
and these letter inequalities justify the statement of the present remark.
It is clear that if the first of the inequalities (22) is verified for any \(n\in \mathbb {N}\) with a certain \(K_{1}{\gt}0,\) this inequality is verified with any number \(K\geq \max \left\{ 1,K_{1}\right\} .\) The situation is identical regarding the second inequality from (22). In conclusion we can suppose that in these relations we have \(K_{1}=K_{2}=K\geq 1.\)
Identically, we can suppose that in the inequalities:
that are true for any \(n\in \mathbb {N},\) we can have \(b=a\geq 1.\)
In conclusion, for the main sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) together with the auxiliary sequences \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\) we can suppose that for any \(n\in \mathbb {N}^{\ast }\) we have \(y_{n}\neq z_{n}\) and there exist the numbers \(K,a\geq 1\) such that for any \(n\in \mathbb {N}^{\ast }\) the following inequalities are verified:
â–¡
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,q\geq 1,\) the following equation in \(x\) on the interval \(\left[ 0,+\infty \right[ :\)
We have the following remarks:
a) The equation (33) has an unique positive root and this root is \(\alpha \in \left] 0,1\right[ .\)
b) If \(\alpha \in \left] 0,1\right[ \) is the root of the equation (33) one verifies the following inequalities as well:
and these inequalities are equivalent to the following inequalities respectively:
â–¡
Indeed, let us consider the function \(\varphi :\left[ 0,+\infty \right[ \rightarrow \mathbb {R}\) defined by \(\varphi \left( x\right) =x^{p+q-1}+2x^{2}+2x-1.\) It is obvious that for any \(x\in \left[ 0,+\infty \right[ \) there exists the derivative \(\varphi ^{\prime }\left( x\right) \) at the point \(x\) and:
As it is obvious that for any \(x\in \left[ 0,+\infty \right[ \) we have that \(\varphi ^{\prime }\left( x\right) {\gt}0,\) therefore the function \(\varphi :\left[ 0,+\infty \right[ \rightarrow \mathbb {R}\) is a strictly increasing function, therefore an injective function, thus the equation \(\varphi \left( x\right) =0\) has at most one root.
As \(\varphi \left( 0\right) =-1,\) \(\varphi \left( 1\right) =4,\) this root exists indeed and it belongs to the interval \(\left] 0,1\right[ .\)
If \(\alpha \in \left] 0,1\right[ \) is the root of the equation (33) it is clear that \(2\alpha ^{2}+2\alpha -1=\)
\(=-\alpha ^{p+q-1}{\lt}0,\) whence \(\alpha ^{2}+2\alpha -1=-\alpha ^{p+q-1}-\alpha ^{2}{\lt}0\) and \(\alpha ^{2}+\alpha -1=\)
\(=-\alpha ^{p+q-1}-\alpha ^{2}-\alpha {\lt}0,\) therefore the relations (34) are true.
The inequality \(\alpha ^{2}+\alpha -1{\lt}0\) is equivalent to \(0{\lt}\alpha ^{2}{\lt}1-\alpha \) therefore to \(0{\lt}\tfrac {\alpha ^{2}}{1-\alpha }{\lt}1,\) the inequality \(\alpha ^{2}+2\alpha -1{\lt}0\) is equivalent to \(0{\lt}\alpha {\lt}1-\alpha -\alpha ^{2},\) therefore to \(0{\lt}\tfrac {\alpha }{1-\alpha -\alpha ^{2}}{\lt}1,\) again the inequality \(2\alpha ^{2}+2\alpha -1{\lt}0\) is equivalent to \(0{\lt}\alpha ^{2}{\lt}1-2\alpha -\alpha ^{2}\) therefore to \(0{\lt}\tfrac {\alpha ^{2}}{1-2\alpha -\alpha ^{2}}{\lt}1.\)
Let us consider now the numbers \(a,K,L,B_{0},R_{0}{\gt}0\) and the numbers \(p,q\geq 1\) and using these numbers we build the real number sequences \(\left( u_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( s_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( v_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( w_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( t_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( B_{n}\right) _{n\in \mathbb {N}^{\ast }}\) and \(\left( R_{n}\right) _{n\in \mathbb {N}^{\ast }}\) using the following recurrence relations:
It is obvious that this construction has a meaning if for any \(n\in \mathbb {N}^{\ast }\) we have that \(u_{n},\ v_{n},\ w_{n},\ t_{n}\in \mathbb {R\diagdown }\left\{ 1\right\} \) and \(B_{n},R_{n}{\gt}0.\)
It is clear that for any \(n\in \mathbb {N}^{\ast }\) we have:
as well.
Referring to the sequences that are defined by the relations (36) we have the following proposition:
If the following inequalities are verified:
(with the specification that for \(q=1\) the expression that has \(q-1\) in its denominator is \(+\infty ,\) and the same for the expression that has \(p-1\) in its denominator) and:
where \(\alpha \in \left] 0,1\right[ \) is the unique root of the equation (33), then for any \(n\in \mathbb {N}^{\ast }\) we have the following inequalities:
From the definition of \(d\) it is clear that:
therefore:
From \(B_{0}\leq \tfrac {1}{\sqrt{L}}K^{\frac{p-q+1}{2\left( q-1\right) }}\) we deduce that \(\tfrac {K^{\frac{p-q+1}{q-1}}}{LB_{0}^{2}}\geq 1,\) therefore as \(\tfrac {q-1}{p+q-1}\geq 0\) it is clear that \(\Big( \tfrac {K^{\frac{p-q+1}{q-1}}}{LB_{0}^{2}}\Big) ^{\frac{q-1}{p+q-1}}\geq 1{\gt}\alpha ,\) therefore:
therefore:
For \(q=1\) from the same definition of \(d\) we have:
namely \(u_{0}\leq \alpha \cdot \tfrac {\alpha }{K}\cdot d.\)
But \(\alpha {\lt}1\leq K,\) so \(\tfrac {\alpha }{K}{\lt}1,\) therefore \(u_{0}\leq \alpha d.\)
Identically if \(p{\gt}1\) we have that:
whence, in the same manner as in the case of \(u_{0},\) by inverting the roles of the numbers \(p\) and \(q,\) we deduce that \(s_{0}\leq \alpha d=\alpha d^{\left( p+q\right) ^{0}}{\lt}\alpha {\lt}1.\)
For \(p=1\) one can show the inequality \(s_{0}\leq \alpha d\) in the same manner as the inequality \(u_{0}\leq \alpha d\) for the case of \(q=1.\)
As \(v_{0}=\tfrac {a}{B_{0}}\cdot \tfrac {u_{0}s_{0}}{1-u_{0}}\) and using \(a\leq B_{0}\) and the inequalities concerning to \(u_{0}\) and \(s_{0},\) we have that:
the last inequality being the first from (35).
Also:
the last inequality being the second from (35), while:
the last inequality being the third from (35).
As \(\ u_{0},v_{0}\in \left] 0,1\right[ \) it is clear that \(v_{0},w_{0},t_{0}\in \left] 0,1\right[ \) as well.
Then we have:
therefore:
As:
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 \(n\in \mathbb {N}\) with \(n\leq k\) and we prove that they are also true for \(n=k+1.\)
For any \(i\in \mathbb {N}\) with \(i\leq k\) we have that:
From the equality that defined the sequence \(\left( B_{n}\right) _{n\in \mathbb {N}^{\ast }}\) we deduce that:
every fraction from those that multiply being greater than the unit.
Therefore \(B_{i}\geq B_{i-1}.\)
From this we actually have that \(B_{i}\geq B_{0},\) from where as \(2q-2\geq 0\) we have \(B_{i}^{2q-2}\geq B_{0}^{2q-2},\) namely \(\tfrac {1}{B_{i}^{2q-2}}\leq \tfrac {1}{B_{0}^{2q-2}}\) and so:
where:
If we multiply by \(C^{\frac{1}{p+q-1}}\) the two members of the inequality (43) we obtain that:
namely the inequality \(h_{i+1}\leq h_{i}^{p+q}\) if \(h_{i}=C^{\frac{1}{p+q-1}}u_{i}.\)
As \(i\in \left\{ 0,1,...,k\right\} \) we immediately deduce that:
Therefore we have that:
But \(\alpha \in \left] 0,1\right[ \) is the root of the equation (33), therefore:
therefore:
At the same time, if \(q{\gt}1\) we have that:
for the last inequality we have take into account (42).
For \(q=1\) we have:
So, from the inequality (44) we obtain that \(u_{k+1}\leq \alpha d^{\left( p+q\right) ^{k+1}}\) and as \(d{\lt}1,\) we obtain \(u_{k+1}\leq \alpha d^{\left( p+q\right) ^{k+1}}{\lt}\alpha {\lt}1.\)
One can obtain the inequality \(s_{k+1}\leq \alpha d^{\left( p+q\right) ^{k+1}}{\lt}\alpha {\lt}1\) in a similar manner with the one concerning \(u_{k+1}.\) It is only necessary to invert the roles of the numbers \(p\) and \(q.\)
We have:
Identically, as \(\tfrac {B_{i}}{B_{i-1}}\geq 1\) for any \(i\leq k\) we deduce that \(\tfrac {B_{k+1}}{B_{k}}\geq 1,\) therefore \(B_{k+1}\geq B_{k}\) and in fact \(B_{k+1}\geq B_{0},\) namely:
As \(a\leq B_{0}\) and \(d{\lt}1\) we have:
In the same manner we have:
and:
For the sequence \(\left( B_{n}\right) _{n\in \mathbb {N}^{\ast }}\) we have:
therefore:
Finally:
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 \(n\in \mathbb {N}^{\ast }.\) The proposition is proved.
6 The main result
We return to the issue of the convergence of the sequences \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }},\) \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\ \left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq D\subseteq X\) to the solution of the equation \(f\left( x\right) =\theta _{Y},\) where \(f:D\rightarrow Y.\)
We have the following fundamental result:
Suppose that the following assumptions hold:
The linear normed space \(\left( X,\left\Vert \cdot \right\Vert _{X}\right) \) is a Banach space;
The mapping \(f:D\rightarrow Y\) admits for any \(x,y\in D\) with \(x\neq y\) a generalized abstract divided difference \(\Gamma _{f;x,y}\in \left( X,Y\right) ^{\ast }\) and there exists a number \(L{\gt}0\) such that for any \(x,y,z\in D\) with \(x\neq y,\ y\neq z\) we have the inequality:
\[ \left\Vert \Gamma _{f;x,y}-\Gamma _{f;y,z}\right\Vert \leq L\left\Vert x-z\right\Vert _{X}; \]The main approximant sequence \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }}\) together with the secondary sequences \(\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }}\) and \(\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\) are such that for any \(n\in \mathbb {N}^{\ast }\) the following equality is fulfilled:
\begin{equation} \Gamma _{f;y_{n},z_{n}}\left( x_{n+1}-y_{n}\right) +f\left( y_{n}\right) =\theta _{Y} \label{f63p}\end{equation}45and the inequalities (32) with the constants \(a,K{\gt}0.\) We also have that \(f\left( y_{n}\right) ,f\left( z_{n}\right) \in Y\diagdown \left\{ \theta _{Y}\right\} ,\) \(y_{n}\neq z_{n}\) and we are in one of the following situations:
\(x_{n}\neq y_{n}\) and \(y_{n+1}\neq z_{n},\)
or
\(x_{n}\neq z_{n}\) and \(z_{n+1}\neq y_{n}.\)
The mapping \(\Gamma _{f;y_{0},z_{0}}\in \left( X,Y\right) ^{\ast }\) is invertible and \(\Gamma _{f;y_{0},z_{0}}^{-1}\in \left( Y,X\right) ^{\ast }.\)
Denoting:
\begin{align*} B_{0} & =\max \big\{ a,\big\Vert \Gamma _{f;y_{0},z_{0}}^{-1}\big\Vert \big\} , \\ R_{0}& =\left\Vert f\left( x_{0}\right) \right\Vert _{Y},\\ \overline{K} & =\max \big\{ K,\big( B_{0}\sqrt{L}\big) ^{\frac{2\left( q-1\right) }{p-q+1}},\big( B_{0}\sqrt{L}\big) ^{\frac{2\left( p-1\right) }{q-p+1}}\big\} , \\ d & =\tfrac {L\overline{K}B_{0}^{2}}{\alpha ^{2}}\cdot \max \Big\{ \tfrac {R_{0}^{p\left( p+q-1\right) }\overline{K}^{p-q+1}}{( LB_{0}^{2}) ^{q-1}},\tfrac {R_{0}^{q\left( p+q-1\right) }\overline{K}^{q-p+1}}{( LB_{0}^{2}) ^{p-1}}\Big\} , \\ \delta & =2aR_{0}+\tfrac {a\alpha ^{2}}{LB_{0}^{2}}\cdot \tfrac {d^{2}}{1-d^{2\left( p+q-1\right) }}+\tfrac {2\alpha }{L\overline{K}B_{0}}\cdot \tfrac {d}{1-d^{p+q-1}}, \end{align*}where \(\alpha \in \left] 0,1\right[ \) is the unique root of the equation (33), the conditions \(d{\lt}1\) and \(S\left( x_{0},\delta \right) =\left\{ \left. x\in X\right/ \left\Vert x-x_{0}\right\Vert _{X}\leq \delta \right\} \) are fulfilled.
Then the following conclusions are true:
for any \(n\in \mathbb {N}^{\ast }\) we have that \(x_{n},y_{n},z_{n}\in S\left( x_{0},\delta \right) ,\) there exists the mapping \(\Gamma _{f;y_{n},z_{n}}^{-1}\in \left( Y,X\right) ^{\ast }\) and:
\begin{equation} x_{n+1}=y_{n}-\Gamma _{f;y_{n},z_{n}}^{-1}f\left( y_{n}\right) =z_{n}-\Gamma _{f;y_{n},z_{n}}^{-1}f\left( z_{n}\right) ; \label{K}\end{equation}46the sequences \(\left( x_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( y_{n}\right) _{n\in \mathbb {N}^{\ast }},\left( z_{n}\right) _{n\in \mathbb {N}^{\ast }}\subseteq X\) are convergent to the limit \(\overline{x}\in S\left( x_{0},\delta \right) \) for that \(f\left( \overline{x}\right) =\theta _{Y};\)
for any \(n\in \mathbb {N}^{\ast }\) the following inequalities are fulfilled:
\begin{equation} \left\Vert x_{n+1}-x_{n}\right\Vert _{X}\leq \tfrac {a\alpha ^{2}}{LB_{0}^{2}}\cdot d^{2\left( p+q\right) ^{n-1}}+\tfrac {\alpha }{L\overline{K}B_{0}}\cdot d^{\left( p+q\right) ^{n}}; \label{K1}\end{equation}47\begin{align} \left\Vert x_{n}-\overline{x}\right\Vert _{X} \leq & \tfrac {a\alpha ^{2}}{LB_{0}^{2}}\cdot \tfrac {d^{2\left( p+q\right) ^{n-1}}}{1-d^{2\left( p+q\right) ^{n-1}\left( p+q-1\right) }}\label{K2}\\ & +\tfrac {\alpha }{L\overline{K}B_{0}}\cdot \tfrac {d^{\left( p+q\right) ^{n}}}{1-d^{\left( p+q\right) ^{n}\left( p+q-1\right) }};\nonumber \end{align}\begin{align} \max \big\{ \left\Vert y_{n}-\overline{x}\right\Vert _{X},\left\Vert z_{n}-\overline{x}\right\Vert _{X}\big\} \leq & \label{K3} \tfrac {a\alpha ^{2}}{LB_{0}^{2}}\cdot d^{2\left( p+q\right) ^{n-1}}\cdot \tfrac {2-d^{2\left( p+q\right) ^{n-1}\left( p+q-1\right) }}{1-d^{2\left( p+q\right) ^{n-1}\left( p+q-1\right) }}\nonumber \\ & +\tfrac {\alpha }{L\overline{K}B_{0}}\cdot \tfrac {d^{\left( p+q\right) ^{n}}}{1-d^{\left( p+q\right) ^{n}\left( p+q-1\right) }}.\medskip \nonumber \end{align}
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
- 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).