Dhage iteration method for approximating solutions of nonlinear differential equations with maxima,

7 years ago

Abstract

In this paper we study the initial value problem of first order nonlinear differential equations with maxima and discuss the existence and approximation of the solutions. The main result relies on the Dhage iteration method embodied in a recent hybrid fixed point theorem of Dhage (2014) in a partially ordered normed linear space. At the end, we give an example to illustrate the hypotheses and applicability of the abstract results of this paper.

Authors

B.C. Dhage
(Kasubai, Gurukul Colony, Ahmepur-413515, Dist. Latur, Maharashtra, India)

D. Otrocol
(Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy,
Technical University of Cluj-Napoca)

Keywords

Differential equations with maxima; Dhage iteration method; hybrid fixed point theorem; approximation of solutions

Cite this paper as:

B.C. Dhage, D. Otrocol, Dhage iteration method for approximating solutions of nonlinear differential equations with maxima, Fixed Point Theory, 19 (2018) no. 2, pp. 545-556.
DOI: 10.24193/fpt-ro.2018.2.43

[1] D.D. Bainov, S. Hristova, Differential Equations with Maxima, Chapman & Hall/CRC Pure and Applied Mathematics, 2011.

[2] B.C. Dhage, Quadratic perturbations of periodic boundary value problems of second order ordinary differential equations, Differ. Equ. & Appl., 2 (2010), 465-486.

[3] B.C. Dhage, Hybrid fixed point theory in partially ordered normed linear spaces and applications to fractional integral equations, Differ. Equ Appl., 5(2013), 155-184.

[4] B.C. Dhage, Partially condensing mappings in ordered normed linear spaces and applications to functional integral equations, Tamkang J. Math., 45(4)(2014), 397-426.

[5] B.C. Dhage, Nonlinear D-set-contraction mappings in partially ordered normed linear spaces and applications to functional hybrid integral equations, Malaya J. Mat., 3(1)(2015), 62-85.

[6] B.C. Dhage, Operator theoretic techniques in the theory of nonlinear hybrid differential equations, Nonlinear Anal. Forum, 20(2015), 15-31.

[7] B.C. Dhage, A new monotone iteration principle in the theory of nonlinear first order integrodifferential equations, Nonlinear Studies, 22(3)(2015), 397-417.

[8] B.C. Dhage, Some generalizations of a hybrid fixed point theorem in a partially ordered metric space and nonlinear functional integral equations, Differ. Equ Appl., 8(2016), 77-97.

[9] B.C. Dhage, S.B. Dhage, Approximating solutions of nonlinear first order ordinary differential equations, GJMS Special Issue for Recent Advances in Mathematical Sciences and Applications13, GJMS Vol., 2(2014), no. 2, 25-35.

[10] B.C. Dhage, S.B. Dhage, Approximating positive solutions of nonlinear first order ordinary quadratic differential equations, Cogent Mathematics, 2(2015), 1023671.

[11] B.C. Dhage, S.B. Dhage, S.K. Ntouyas, Approximating solutions of nonlinear hybrid differential equations, Appl. Math. Lett., 34(2014), 76-80.

[12] B.C. Dhage, S.B. Dhage, J.R. Graef, Dhage iteration method for initial value problems for nonlinear first order hybrid integrodifferential equations, J. Fixed Point Theory Appl., 17(2016), 309-325. 556 BAPURAO C. DHAGE AND DIANA OTROCOL

[13] S. Heikkil¨a and V. Lakshmikantham, Monotone Iterative Techniques for Discontinuous Nonlinear Differential Equations, Marcel Dekker Inc., New York 1994.

[14] J.J. Nieto, R. Rodriguez-Lopez, Contractive mappings theorems in partially ordered sets and applications to ordinary differential equations, Order, 22(2005), 223-239.

[15] D. Otrocol, Properties of the solutions of system of differential equations with maxima, via weakly Picard operator theory, Comm. in Applied Anal., 17(2013), no. 1, 99-107.

[16] D. Otrocol, Systems of functional differential equations with maxima, of mixed type, Electron. J. Qual. Theory Differ. Equ., 2014(2014), no. 5, 1-9.

[17] D. Otrocol, I.A. Rus, Functional-differential equations with ”maxima” via weakly Picard operators theory, Bull. Math. Soc. Sci. Math. Roumanie, 51(99)(2008), no. 3, 253-261.

[18] D. Otrocol, I.A. Rus, Functional-differential equations with maxima of mixed type argument, Fixed Point Theory, 9(2008), no. 1, 207-220.

[19] A. Petru¸sel, I.A. Rus, Fixed point theorems in ordered L-spaces, Proc. Amer. Math. Soc., 134(2006), 411-418.

[20] A.C.M. Ran, M.C. Reurings, A fixed point theorem in partially ordered sets and some applications to matrix equations, Proc. Amer. Math. Soc., 132(2004), 1435-1443.

PDF

https://doi.org/10.24193/fpt-ro.2018.2.43

About this paper

Journal

Fixed Point Theory

Publisher Name

House Book Science – Casa Cartii Stiinta, Cluj-Napoca, Romania

DOI

10.24193/fpt-ro.2018.2.43

Print ISSN

1583-5022

Online ISSN

2066-9208

MR

MR3821782

ZBL

1397.34108

Google Scholar

https://scholar.google.ro/scholar?oi=bibs&hl=en&cites=2571673060656038997

Fixed Point Theory, 19(2018), No.2, 545-556

http://www.math.ubbcluj.ro/^∼nodeacj/sfptcj.html

Dhage iteration method for approximating solutions of nonlinear differential equations with maxima

Bapurao C. Dhage^∗, Diana Otrocol^∗∗

^∗Kasubai, Gurukul Colony, Ahmepur-413515, Dist. Latur, Maharashtra, India
E-mail: bcdhage@gmail.com
^∗∗“T. Popoviciu” Institute of Numerical Analysis, P.O.Box. 68-1, 400110 Cluj-Napoca, Romania E-mail: dotrocol@ictp.acad.ro

Abstract. In this paper we study the initial value problem of first order nonlinear differential equations with maxima and discuss the existence and approximation of the solutions. The main result relies on the Dhage iteration method embodied in a recent hybrid fixed point theorem of Dhage (2014) in a partially ordered normed linear space. At the end, we give an example to illustrate the hypotheses and applicability of the abstract results of this paper.

Key Words and Phrases: Differential equations with maxima, Dhage iteration method, hybrid fixed point theorem, approximation of solutions.

2010 Mathematics Subject Classification: 34A12, 34A45, 47H07, 47H10.

1. Introduction

The study of fixed point theorems for the contraction mappings in partially ordered metric spaces is initiated by Ran and Reurings [20] which are further continued by Nieto and Rodringuez-Lopez [13] and by Petruşel and Rus [19] and applied to boundary value problems of nonlinear first order ordinary differential equations and matrix equations for proving the existence results under certain monotonic conditions. Similarly, the study of hybrid fixed point theorems in a partially ordered metric spaces is initiated by Dhage [3, 4, 5] with applications to nonlinear differential and integral equation under weaker mixed conditions of nonlinearities. See Dhage [6, 7] and the references therein. In this paper we investigate the existence of approximate solutions of hybrid differential equations with maxima using the Dhage iteration method embodied in a hybrid fixed point theorem in a partially ordered spaces. We claim that the results of this paper are new to the theory of nonlinear differential equations with maxima.

Given a closed and bounded interval $J=[a,b]$ of the real line $\mathbb{R}$ for some $b>a\geq 0$ , we consider the following hybrid differential equation (in short HDE)

\left\{\begin{array}[l]{l}x^{\prime}(t)=f(t,x(t))+g\left(t,\underset{a\leq\xi\leq t}{\max}\,x(\xi)\right),\\ x(a)=\alpha_{0}\in\mathbb{R},\end{array}\right.

(1.1)

for all $t\in J=[a,b]$ and $f,g~:~J\times\mathbb{R}\rightarrow\mathbb{R}$ are continuous functions.

By a solution of equation (1.1) we mean a differentiable function $x\in C(J,\mathbb{R})$ that satisfies equation (1.1), where $C(J,\mathbb{R})$ is the space of continuous real-valued functions defined on $J$ .

Differential equations with maxima are often met in the applications, for instance in the theory of automatic control. Numerous results on existence and uniqueness, asymptotic stability as well as numerical solutions have been obtained. To name a few, we refer the reader to [1, 14, 15, 16, 17] and the references therein. The HDE (1.1) is a linear perturbation of first type of nonlinear differential equations. The details of different types of perturbation appears in Dhage [2]. The special cases of the HDE (1.1) in the form

\left\{\begin{array}[]{l}x^{\prime}(t)=f(t,x(t)),\ t\in J,\\[5.69054pt] x(a)=\alpha_{0},\end{array}\right.

(1.2)

and

\left\{\begin{array}[]{l}x^{\prime}(t)=g\left(t,\underset{a\leq\xi\leq t}{\max}\,x(\xi)\right),\ t\in J,\\[8.53581pt] x(a)=\alpha_{0},\end{array}\right.

(1.3)

have already been discussed in the literature for different aspects of the solutions using usual Picard iteration method. See Bainov and Hristova [1] and the references therein for the details. In this paper we discuss the HDE (1.1) for existence and approximation of solutions via a new approach based upon the Dhage iteration method. In consequence, we obtain the existence and approximation results for HDEs (1.2) and (1.3) as special cases which are also new to the literature.

In the following section we give some preliminaries and the key tool that will be used for proving the main result of this paper.

2. Preliminaries

Throughout this paper, unless otherwise mentioned, let $(E,\preceq,\|\cdot\|)$ denote a partially ordered normed linear space. Two elements $\,x\,$ and $\,y\,$ in $E$ are said to be comparable if either the relation $x\preceq y$ or $y\preceq x$ holds. A non-empty subset $C$ of $E$ is called a chain or totally ordered if all the elements of $C$ are comparable. It is known that $E$ is regular if $\{x_{n}\}$ is a nondecreasing (resp. nonincreasing) sequence in $E$ such that $x_{n}\to x^{*}$ as $n\to\infty$ , then $x_{n}\preceq x^{*}$ (resp. $x_{n}\succeq x^{*}$ ) for all $n\in\mathbb{N}$ . The conditions guaranteeing the regularity of $E$ may be found in Heikkilä and Lakshmikantham [18] and the references therein.

We need the following definitions (see Dhage [3, 4, 5] and the references therein) in what follows.

Definition 2.1.

A mapping $\mathcal{T}:E\to E$ is called isotone or nondecreasing if it preserves the order relation $\preceq$ , that is, if $x\preceq y$ implies $\mathcal{T}x\preceq\mathcal{T}y$ for all $x,y\in E$ . Similarly, $\mathcal{T}$ is called nonincreasing if $x\preceq y$ implies $\mathcal{T}x\succeq\mathcal{T}y$ for all $x,y\in E$ . Finally, $\mathcal{T}$ is called monotonic or simply monotone if it is either nondecreasing or nonincreasing on $E$ .

Definition 2.2.

A mapping $\mathcal{T}:E\to E$ is called partially continuous at a point $a\in E$ if for $\epsilon>0$ there exists a $\delta>0$ such that $\|\mathcal{T}x-\mathcal{T}a\|<\epsilon$ whenever $x$ is comparable to $a$ and $\|x-a\|<\delta$ . $\mathcal{T}$ called partially continuous on $E$ if it is partially continuous at every point of it. It is clear that if $\mathcal{T}$ is partially continuous on $E$ , then it is continuous on every chain $C$ contained in $E$ .

Definition 2.3.

A non-empty subset $S$ of the partially ordered Banach space $E$ is called partially bounded if every chain $C$ in $S$ is bounded. An operator $\mathcal{T}$ on a partially normed linear space $E$ into itself is called partially bounded if $\mathcal{T}(E)$ is a partially bounded subset of $E$ . $\mathcal{T}$ is called uniformly partially bounded if all chains $C$ in $\mathcal{T}(E)$ are bounded by a unique constant.

Definition 2.4.

A non-empty subset $S$ of the partially ordered Banach space $E$ is called partially compact if every chain $C$ in $S$ is a relatively compact subset of $E$ . A mapping $\mathcal{T}:E\to E$ is called partially compact if $\mathcal{T}(E)$ is a partially relatively compact subset of $E$ . $\mathcal{T}$ is called uniformly partially compact if $\mathcal{T}$ is a uniformly partially bounded and partially compact operator on $E$ . $\mathcal{T}$ is called partially totally bounded if for any bounded subset $S$ of $E$ , $\mathcal{T}(S)$ is a partially relatively compact subset of $E$ . If $\mathcal{T}$ is partially continuous and partially totally bounded, then it is called partially completely continuous on $E$ .

Remark 2.1.

Suppose that $\mathcal{T}$ is a nondecreasing operator on $E$ into itself. Then $\mathcal{T}$ is a partially bounded or partially compact if $\,\mathcal{T}(C)\,$ is a bounded or relatively compact subset of $E$ for each chain $C$ in $E$ .

Definition 2.5.

The order relation $\preceq$ and the metric $\,d\,$ on a non-empty set $E$ are said to be compatible if $\{x_{n}\}$ is a monotone sequence, that is, monotone nondecreasing or monotone nonincreasing sequence in $E$ and if a subsequence $\{x_{n_{k}}\}$ of $\{x_{n}\}$ converges to $x^{*}$ implies that the original sequence $\{x_{n}\}$ converges to $x^{*}$ . Similarly, given a partially ordered normed linear space $(E,\preceq,\|\cdot\|)$ , the order relation $\preceq$ and the norm $\|\cdot\|$ are said to be compatible if $\preceq$ and the metric $d$ defined through the norm $\|\cdot\|$ are compatible. A subset $S$ of $E$ is called Janhavi if the order relation $\preceq$ and the metric $d$ or the norm $\|\cdot\|$ are compatible in it. In particular, if $S=E$ , then $E$ is called a Janhavi metric or Janhavi Banach space.

Clearly, the set $\mathbb{R}$ of real numbers with usual order relation $\leq$ and the norm defined by the absolute value function $|\cdot|$ has this property. Similarly, the finite dimensional Euclidean space $\mathbb{R}^{n}$ with usual componentwise order relation and the standard norm possesses the compatibility property and so is a Janhavi Banach space.

Definition 2.6.

An upper semi-continuous and monotone nondecreasing function $\psi:\mathbb{R}_{+}\to\mathbb{R}_{+}$ is called a $\mathcal{D}$ -function provided $\psi(0)=0$ . An operator $\mathcal{T}:E\to E$ is called partially nonlinear $\mathcal{D}$ -contraction if there exists a $\mathcal{D}$ -function $\psi$ such that

\|\mathcal{T}x-\mathcal{T}y\|\leq\psi\big(\|x-y\|\big)

(2.1)

for all comparable elements $x,y\in E$ , where $0<\psi(r)<r$ for $r>0$ . In particular, if $\psi(r)=k\,r$ , $k>0$ , $\mathcal{T}$ is called a partial Lipschitz operator with a Lischitz constant $k$ and moreover, if $0<k<1$ , $\mathcal{T}$ is called a partial linear contraction on $E$ with a contraction constant $k$ .

The Dhage iteration method embodied in the following applicable hybrid fixed point theorem of Dhage [4] in a partially ordered normed linear space is used as a key tool for our work contained in this paper. The details of a Dhage iteration method is given in Dhage [6, 7, 8], Dhage et al. [11, 12] and the references therein.

Theorem 2.1 (Dhage [4]).

Let $\big(E,\preceq,\|\cdot\|\big)$ be a regular partially ordered complete normed linear space such that every compact chain $C$ of $E$ is Janhavi. Let $\mathcal{A},\mathcal{B}:E\to E$ be two nondecreasing operators such that

(a)

$\mathcal{A}$ is partially bounded and partially nonlinear $\mathcal{D}$ -contraction,
(b)

$\mathcal{B}$ is partially continuous and partially compact, and
(c)

there exists an element $x_{0}\in E$ such that $x_{0}\preceq\mathcal{A}x_{0}+\mathcal{B}x_{0}\,$ or $\,x_{0}\succeq\mathcal{A}x_{0}+\mathcal{B}x_{0}$ .

Then the operator equation $\mathcal{A}x+\mathcal{B}x=x$ has a solution $x^{*}$ in $E$ and the sequence $\{x_{n}\}$ of successive iterations defined by $x_{n+1}=\mathcal{A}x_{n}+\mathcal{B}x_{n}$ , n=0,1,…, converges monotonically to $x^{*}$ .

Remark 2.2.

The condition that every compact chain of $E$ is Janhavi holds if every partially compact subset of $E$ possesses the compatibility property with respect to the order relation $\preceq$ and the norm $\|\cdot\|$ in it.

Remark 2.3.

We remark that hypothesis (a) of Theorem 2.1 implies that the operator $\mathcal{A}$ is partially continuous and consequently both the operators $\mathcal{A}$ and $\mathcal{B}$ in the theorem are partially continuous on $E$ . The regularity of $E$ in above Theorem 2.1 may be replaced with a stronger continuity condition of the operators $\mathcal{A}$ and $\mathcal{B}$ on $E$ which is a result proved in Dhage [3, 4].

3. Main Result

In this section, we prove an existence and approximation result for the HDE (1.1) on a closed and bounded interval $J=[a,b]$ under mixed partial Lipschitz and partial compactness type conditions on the nonlinearities involved in it. We place the HDE (1.1) in the function space $C(J,\mathbb{R})$ of continuous real-valued functions defined on $J$ . We define a norm $\|\,\cdot\,\|$ and the order relation $\,\leq\,$ in $C(J,\mathbb{R})$ by

\|x\|=\displaystyle\sup_{t\in J}|x(t)|

(3.1)

and

x\leq y\iff x(t)\leq y(t)\quad\mbox{for all }\,\,t\in J.

(3.2)

Clearly, $C(J,\mathbb{R})$ is a Banach space with respect to above supremum norm and also partially ordered w.r.t. the above partially order relation $\leq$ . It is known that the partially ordered Banach space $C(J,\mathbb{R})$ is regular and lattice so that every pair of elements of $E$ has a lower and an upper bound in it. The following useful lemma concerning the Janhavi subsets of $C(J,\mathbb{R})$ follows immediately from the Arzelá-Ascoli theorem for compactness.

Lemma 3.1.

Let $\big(C(J,\mathbb{R}),\leq,\|\cdot\|\big)$ be a partially ordered Banach space with the norm $\|\cdot\|$ and the order relation $\leq$ defined by (3.1) and (3.2) respectively. Then every partially compact subset of $C(J,\mathbb{R})$ is Janhavi.

Proof.

The proof of the lemma is well-known and appears in the papers of Dhage [6, 7], Dhage and Dhage [9, 10], Dhage et al. [12] and so we omit the details. ∎

We need the following definition in what follows.

Definition 3.1.

A differentiable function $u\in C(J,\mathbb{R})$ is said to be a lower solution of the equation (1.1) if it satisfies

\left\{\begin{array}[c]{l}u^{\prime}(t)\leq f(t,u(t))+g\left(t,\underset{a\leq\xi\leq t}{\max}\,u(\xi)\right),\\ u(a)\leq\alpha_{0},\end{array}\right.

(*)

for all $t\in J$ . Similarly, a differentiable function $v\in C(J,\mathbb{R})$ is called an upper solution of the HDE (1.1) if the above inequality is satisfied with reverse sign.

We consider the following set of assumptions in what follows:

(H₁)

There exist constants $\lambda>0,\mu>0$ with $\lambda\geq\mu$ such that

$0\leq[f(t,x)+\lambda x]-[f(t,y)+\lambda y]\leq\mu(x-y)$

for all $t\in J$ and $x,y\in\mathbb{R}$ , $x\geq y$ .
(H₂)

There exists a constant $M>0$ such that $\left|g(t,x)\right|\leq M,$ for all $t\in J,\ x\in\mathbb{R};$
(H₃)

$g(t,x)$ is nondecreasing in $x$ for each $t\in J$ .
(H₄)

HDE (1.1) has a lower solution $u\in C(J,\mathbb{R}).$

Now we consider the following HDE

\left\{\begin{array}[c]{l}x^{\prime}(t)+\lambda x(t)=\widetilde{f}(t,x(t))+g\left(t,\underset{a\leq\xi\leq t}{\max}\,u(\xi)\right),\\ x(a)=\alpha_{0},\end{array}\right.

(3.3)

for all $t\in J=[a,b],\$ where $\widetilde{f},g~:~J\times\mathbb{R}\rightarrow\mathbb{R}$ and $\widetilde{f}(t,x)=f(t,x)+\lambda x,\ \lambda>0$ .

Remark 3.1.

A differentiable function $\ u\in C(J,\mathbb{R})$ is a solution of the equation (3.3) if and only if it is a solution of the equation (1.1) defined on $J.$

We also consider the following condition in what follows.

(H₅)

There exists a constant $K>0$ such that $\left|\widetilde{f}(t,x)\right|\leq K,$ for all $t\in J\$ and $x\in\mathbb{R};$

Lemma 3.2.

Suppose that the hypotheses $(H_{2}),(H_{3})$ and $(H_{5})$ hold. Then a function $x\in C(J,\mathbb{R})$ is a solution of the HDE (3.3) if and only if it is a solution of the nonlinear integral equation

\begin{split}x(t)&=\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,x(s))\,ds\\ &\quad\quad\quad+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds,\end{split}

(3.4)

for all $t\in J$ .

Theorem 3.1.

Suppose that hypotheses $(H_{1})-(H_{5})$ hold. Then the HDE (1.1) has a solution $x^{\ast}$ defined on $J$ and the sequence $\{x_{n}\}\,$ of successive approximations defined by

\begin{gathered}x_{0}=u,\\ x_{n+1}(t)=\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,x_{n}(s))\,ds\\ \quad\quad\quad+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds,\end{gathered}

(3.5)

for all $t\in J$ , converges monotonically to $x^{\ast}$ .

Proof.

Set $E=C(J,\mathbb{R})$ . Then, in view of Lemma 3.1, every compact chain $C$ in $E$ possesses the compatibility property with respect to the norm $\|\cdot\|$ and the order relation $\leq$ so that every compact chain $C$ is Janhavi in $E$ .

Define two operators $\mathcal{A}$ and $\mathcal{B}$ on $E$ by

\mathcal{A}x(t)=\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,x(s))\,ds,\ t\in J,

(3.6)

and

\mathcal{B}x(t)=e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds,\ t\in J.

(3.7)

From the continuity of the integral, it follows that $\mathcal{A}$ and $\mathcal{B}$ define the operators $\mathcal{A},\mathcal{B}:E\rightarrow E$ . Applying Lemma 3.2, the HDE (1.1) is equivalent to the operator equation

\mathcal{A}x(t)+\mathcal{B}x(t)=x(t),\ t\in J.

Now, we show that the operators $\mathcal{A}$ and $\mathcal{B}$ satisfy all the conditions of Theorem 2.1 in a series of following steps.

Step I: $\mathcal{A}$ and $\mathcal{B}$ are nondecreasing on $E.$

Let $x,y\in E$ be such that $x\geq y.$ Then by hypothesis $(H_{1})$ , we get

	$\displaystyle\mathcal{A}x(t)$	$\displaystyle=\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,x(s))\,ds$
		$\displaystyle\geq\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,y(s))\,ds$
		$\displaystyle=\mathcal{A}y(t),$

for all $t\in J$ .

Next, we show that the operator $\mathcal{B}$ is also nondecreasing on $E$ . Let $x,y\in E$ be such that $x\geq y.$ Then $x(t)\geq y(t)$ for all $t\in J$ . Since $y$ is continuous on $[a,t]$ , there exists a $\xi^{*}\in[a,t]$ such that $y(\xi^{*})=\displaystyle\max_{a\leq\xi\leq t}y(\xi)$ . By definition of $\leq$ , one has $x(\xi^{*})\geq y(\xi^{*})$ . Consequently, we obtain

\displaystyle\max_{a\leq\xi\leq t}x(\xi)\geq x(\xi^{*})\geq y(\xi^{*})=\displaystyle\max_{a\leq\xi\leq t}y(\xi).

Now, using hypothesis (H₃), it can be shown that the operator $\mathcal{B}$ is also nondecreasing on $E.$

Step II: $\mathcal{A}$ is partially bounded and partially contraction on $E.$

Let $x\in E$ be arbitrary. Then by $(H_{5})$ we have

	$\displaystyle\left\|\mathcal{A}x(t)\right\|$	$\displaystyle\leq\left\|\alpha_{0}e^{-\lambda t}\right\|+e^{-\lambda t}\int_{a}^{t}e^{-\lambda s}\left\|\widetilde{f}(s,x(s))\right\|\,ds$
		$\displaystyle\leq\left\|\alpha_{0}\right\|+K\int_{a}^{b}e^{\lambda s}\,ds$
		$\displaystyle\leq\left\|\alpha_{0}\right\|+e^{\lambda a}K(b-a),$

for all $t\in J$ . Taking the supremum over $t,$ we obtain $\left\|\mathcal{A}x(t)\right\|\leq\left|\alpha_{0}\right|+e^{\lambda a}K(b-a),$ so $\mathcal{A}$ is a bounded operator on $E$ . This implies that $\mathcal{A}$ is partially bounded on $E.$

Let $x,y\in E$ be such that $x\geq y.$ Then by $(H_{1})$ we have

	$\displaystyle\left\|\mathcal{A}x(t)-\mathcal{A}y(t)\right\|$	$\displaystyle\leq\left\|e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\left[\widetilde{f}(s,x(s))-\widetilde{f}(s,y(s))\right]\,ds\right\|$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\mu\left\|x(s)-y(s)\right\|\,ds$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\lambda\left\|x(s)-y(s)\right\|\,ds$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}\frac{d}{ds}e^{\lambda s}\left\\|x-y\right\\|\,ds$
		$\displaystyle\leq(1-e^{-\lambda a})\left\\|x-y\right\\|,$

for all $t\in J$ . Taking the supremum over $t$ , we obtain $\left\|\mathcal{A}x-\mathcal{A}y\right\|\leq L\left\|x-y\right\|,$ for all $x,y\in E$ with $x\geq y.$ Hence $\mathcal{A}$ is a partially contraction on $E$ and which also implies that $\mathcal{A}$ is partially continuous on $E.$

Step III: $\mathcal{B}$ is partially continuous on $E.$

Let $\{x_{n}\}_{n\in\mathbb{N}}$ be a sequence in a chain $C$ such that $x_{n}\rightarrow x,$ for all $n\in\mathbb{N}$ . Then

	$\displaystyle\underset{n\rightarrow\infty}{\lim}\mathcal{B}x_{n}(t)$	$\displaystyle=\underset{n\rightarrow\infty}{\lim}e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds$
		$\displaystyle=e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\left[\underset{n\rightarrow\infty}{\lim}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\right]\,ds$
		$\displaystyle=e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds$
		$\displaystyle=\mathcal{B}x(t),$

for all $t\in J$ . This shows that $\mathcal{B}x_{n}$ converges to $\mathcal{B}x$ pointwise on $J$ .

Now we show that $\{\mathcal{B}x_{n}\}_{n\in\mathbb{N}}$ is an equicontinuous sequence of functions in $E.\$ Let $t_{1},t_{2}\in J$ with $t_{1}<t_{2}$ . We have

	$\displaystyle\left\|\mathcal{B}x_{n}(t_{2})-\mathcal{B}x_{n}(t_{1})\right\|$	$\displaystyle=\left\|e^{-\lambda t_{2}}\int_{a}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right.$
		$\displaystyle\quad\quad\quad\left.-e^{-\lambda t_{1}}\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\leq\left\|\left(e^{-\lambda t_{2}}-e^{-\lambda t_{1}}\right)\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\quad+\left\|e^{-\lambda t_{2}}\int_{t_{1}}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\rightarrow 0\quad\text{ as }\quad t_{2}\rightarrow t_{1},$

uniformly for all $n\in\mathbb{N}.$ This shows that the convergence $\mathcal{B}x_{n}\rightarrow\mathcal{B}x$ is uniform and hence $\mathcal{B}$ is partially continuous on $E$ .

Step IV: $\mathcal{B}$ is partially compact operator on $E.$

Let $C$ be an arbitrary chain in $E$ . We show that $\mathcal{B}(C)$ is uniformly bounded and equicontinuous set in $E$ . First we show that $\mathcal{B}(C)$ is uniformly bounded. Let $y\in\mathcal{B}(C)$ be any element. Then there is an element $x\in C$ such that $y=\mathcal{B}x.$ By hypothesis (H₂)

	$\displaystyle\left\|y(t)\right\|$	$\displaystyle=\left\|\mathcal{B}x(t)\right\|$
		$\displaystyle=\left\|e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds\right\|$
		$\displaystyle\leq\int_{a}^{t}e^{\lambda s}\left\|g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,\right\|\,ds$
		$\displaystyle\leq\int_{a}^{b}e^{\lambda b}M\,ds$
		$\displaystyle\leq e^{\lambda b}M(b-a)=r,$

for all $t\in J$ . Taking the supremum over $t$ we obtain $\left\|y\right\|\leq\left\|\mathcal{B}x\right\|\leq r,$ for all $y\in\mathcal{B}(C)$ . Hence $\mathcal{B}(C)$ is uniformly bounded subset of $E$ . Next we show that $\mathcal{B}(C)$ is an equicontinuous set in $E$ . Let $t_{1},t_{2}\in J$ , with $t_{1}<t_{2}$ . Then, for any $y\in\mathcal{B}(C),$ one has

	$\displaystyle\big\|y(t_{2})$	$\displaystyle-y(t_{1})\big\|=\left\|\mathcal{B}x(t_{2})-\mathcal{B}x(t_{1})\right\|$
		$\displaystyle=\left\|e^{-\lambda t_{2}}\int_{a}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds-e^{-\lambda t_{1}}\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds\right\|$

	$\displaystyle\leq\left\|\left(e^{-\lambda t_{2}}-e^{-\lambda t_{1}}\right)\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds\right\|$
	$\displaystyle\quad\quad+\left\|e^{-\lambda t_{2}}\int_{t_{1}}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds\right\|$
	$\displaystyle\rightarrow 0\quad\text{ as }\quad t_{1}\rightarrow t_{2}$

uniformly for all $y\in\mathcal{B}(C).$ This shows that $\mathcal{B}(C)$ is an equicontinuous subset of $E$ . So $\mathcal{B}(C)$ is a uniformly bounded and equicontinuous set of functions in $E$ and hence it is compact in view of Arzelá-Ascoli theorem. Consequently $\mathcal{B}:E\rightarrow E$ is a partially compact operator of $E$ into itself.

Step V: $u$ satisfies the inequality $u\leq\mathcal{A}u+\mathcal{B}u$ .

By hypothesis $(H_{4})$ the equation (1.1) has a lower solution $u$ defined on $J$ . Then we have

\left\{\begin{array}[c]{l}u^{\prime}(t)\leq f(t,u(t))+g\left(t,\underset{a\leq\xi\leq t}{\max}\,u(\xi)\right),\ t\in J,\\ u(a)\leq\alpha_{0}.\end{array}\right.

(3.8)

Adding $\lambda u(t)$ on both sides of the first inequality in (3.8), we obtain

u^{\prime}(t)+\lambda u(t)\leq f(t,u(t))+\lambda u(t)+g\left(t,\underset{a\leq\xi\leq t}{\max}\,u(\xi)\right),\ t\in J.

Again, multiplying the above inequality by $e^{\lambda t},$

\left(e^{\lambda t}u(t)\right)^{\prime}\leq e^{\lambda t}\widetilde{f}(t,u(t))+e^{\lambda t}g\left(t,\underset{a\leq\xi\leq t}{\max}\,u(\xi)\right).

(3.9)

A direct integration of (3.9) from $a$ to $t$ yields

\begin{split}u(t)&\leq\alpha_{0}e^{-\lambda t}+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\widetilde{f}(s,u(s))ds\\ &\quad\quad+e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,u(\xi)\right)ds,\end{split}

(3.10)

for $t\in J$ . From definitions of the operators $\mathcal{A}$ and $\mathcal{B}$ it follows that $u(t)\leq\mathcal{A}u(t)+\mathcal{B}u(t)$ , for all $t\in J$ . Hence $u\leq\mathcal{A}u+\mathcal{B}u$ . Thus $\mathcal{A}$ and $\mathcal{B}$ satisfy all the conditions of Theorem 2.1 and we apply it to conclude that the operator equation $\mathcal{A}x+\mathcal{B}x=x$ has a solution. Consequently the integral equation and the equation (1.1) has a solution $x^{\ast}$ defined on $J$ . Furthermore, the sequence $\{x_{n}\}_{n=0}^{\infty}$ of successive approximations defined by (3.5) converges monotonically to $x^{\ast}$ . This completes the proof. ∎

Remark 3.2.

The conclusion of Theorem 3.1 also remains true if we replace the hypothesis (H₄) with the following one.

(H ${}^{\prime}_{4}$ )

The HDE (1.1) has an upper solution $v\in C(J,\mathbb{R})$ .

Remark 3.3.

We note that if the HDE (1.1) has a lower solution $u$ as well as an upper solution $v$ such that $u\leq v$ , then under the given conditions of Theorem 3.1 it has corresponding solutions $x_{*}$ and $x^{*}$ and these solutions satisfy $x_{*}\leq x^{*}$ . Hence they are the minimal and maximal solutions of the HDE (1.1) in the vector segment $[u,v]$ of the Banach space $E=C(J,\mathbb{R})$ , where the vector segment $[u,v]$ is a set of elements in $C(J,\mathbb{R})$ defined by

[u,v]=\{x\in C(J,\mathbb{R})\mid u\leq x\leq v\}.

This is because the order relation $\leq$ defined by (3.2) is equivalent to the order relation defined by the order cone $\mathcal{K}=\{x\in C(J,\mathbb{R})\mid x\geq\theta\}$ which is a closed set in $C(J,\mathbb{R})$ .

In the following we illustrate our hypotheses and the main abstract result for the validity of conclusion.

Example 3.1.

We consider the following HDE

\left\{\begin{array}[]{l}x^{\prime}(t)=\arctan x(t)-x(t)+\tanh\left(\underset{0\leq\xi\leq t}{\max}\,x(\xi)\right),\ t\in J=[0,1],\\ x(0)=1.\end{array}\right.

(3.11)

Here $f(t,x)=\arctan x(t)-x(t)$ and $g(t,x)=\tanh x$ . The functions $f$ and $g$ are continuous on $J\times\mathbb{R}.$ Next, we have

0\leq\arctan x(t)-\arctan y(t)\leq\frac{1}{\xi^{2}+1}(x-y),

for all $x,y\in\mathbb{R},x>\xi>y$ . Therefore $\lambda=1>\frac{1}{\xi^{2}+1}=\mu.$ Hence the function $f$ satisfies the hypothesis $(H_{1}).$ Moreover, the function $\widetilde{f}(t,x)=\arctan x(t)$ is bounded on $J\times\mathbb{R}$ with bound $K=\pi/2$ , so that the hypothesis $(H_{5})$ is satisfied. The function $g$ is bounded on $J\times\mathbb{R}$ by $M=1,$ so $(H_{2})$ holds. The function $g(t,x)$ is increasing in $x$ for each $t\in J$ , so the hypothesis $(H_{3})$ is satisfied. The HDE (3.11) has a lower solution $u(t)=-2t+1$ , $t\in[0,1]$ . Thus all hypothesis of Theorem 3.1 are satisfied and hence the HDE (3.11) has a solution $x^{\ast}$ defined on $J$ and the sequence $\{x_{n}\}_{n=0}^{\infty}$ defined by

	$\displaystyle x_{0}$	$\displaystyle=-2t+1,$
	$\displaystyle x_{n+1}(t)$	$\displaystyle=e^{-t}+e^{-t}\int_{0}^{t}e^{s}\arctan x_{n}(s)\,ds$
		$\displaystyle\quad+e^{-t}\int_{0}^{t}e^{s}\tanh\left(\underset{0\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds$

for each $t\in J$ , converges monotonically to $x^{\ast}.$

Remark 3.4.

Finally while concluding, we mention that the study of this paper may be extended with appropriate modifications to the nonlinear hybrid differential equation with maxima,

\left\{\begin{array}[l]{l}x^{\prime}(t)=f\left(t,x(t),\max\limits_{a\leq\xi\leq t}\,x(\xi)\right)+g\left(t,x(t),\max\limits_{a\leq\xi\leq t}\,x(\xi)\right),\\[11.38109pt] x(a)=\alpha_{0}\in\mathbb{R},\end{array}\right.

(3.12)

for all $t\in J=[a,b]$ , where $f,g~:~J\times\mathbb{R}\times\mathbb{R}\to\mathbb{R}$ are continuous functions. When $g\equiv 0$ , the differential equation (3.12) reduces to the nonlinear differential equations with maxima,

\left\{\begin{array}[l]{l}x^{\prime}(t)=f\left(t,x(t),\max\limits_{a\leq\xi\leq t}\,x(\xi)\right),\,\,\,t\in J,\\[11.38109pt] x(a)=\alpha_{0}\in\mathbb{R},\end{array}\right.

(3.13)

which is studied in Otrocol and Rus [16] for existence and uniqueness theorem via Picard iterations under strong Lipschitz condition. Therefore, the obtained results for differential equation (3.12) with maxima via Dhage iteration method will include the existence and approximation results for the differential equation with maxima (3.13) under weak partial Lipschitz condition.

Acknowledgment

The authors are thankful to Prof. I.A. Rus (Romania) and Dr. Shyam B. Dhage (India) for giving some suggestions for the improvement of this paper.

References

[1] D.D. Bainov, S. Hristova, Differential equations with maxima, Chapman & Hall/CRC Pure and Applied Mathematics, 2011.
[2] B.C. Dhage, Quadratic perturbations of periodic boundary value problems of second order ordinary differential equations, Differ. Equ. & Appl., 2 (2010), 465-486.
[3] B.C. Dhage, Hybrid fixed point theory in partially ordered normed linear spaces and applications to fractional integral equations, Differ. Equ Appl., 5(2013), 155-184.
[4] B.C. Dhage, Partially condensing mappings in ordered normed linear spaces and applications to functional integral equations, Tamkang J. Math., 45 (4)(2014), 397-426.
[5] B.C. Dhage, Nonlinear $\mathcal{D}$ -set-contraction mappings in partially ordered normed linear spaces and applications to functional hybrid integral equations, Malaya J. Mat., 3(1)(2015), 62-85.
[6] B.C. Dhage, Operator theoretic techniques in the theory of nonlinear hybrid differential equations, Nonlinear Anal. Forum, 20(2015), 15-31.
[7] B.C. Dhage, A new monotone iteration principle in the theory of nonlinear first order integro-differential equations, Nonlinear Studies, 22(3)(2015), 397-417.
[8] B.C. Dhage, Some generalizations of a hybrid fixed point theorem in a partially ordered metric space and nonlinear functional integral equations, Differ. Equ Appl., 8(2016), 77-97.
[9] B.C. Dhage, S.B. Dhage, Approximating solutions of nonlinear first order ordinary differential equations, GJMS Special issue for Recent Advances in Mathematical Sciences and Applications-13, GJMS Vol., 2(2014), No. 2, 25-35.
[10] B.C. Dhage, S.B. Dhage, Approximating positive solutions of nonlinear first order ordinary quadratic differential equations, Cogent Mathematics, (2015), 2: 1023671.
[11] B.C. Dhage, S.B. Dhage, S.K. Ntouyas, Approximating solutions of nonlinear hybrid differential equations, Appl. Math. Lett., 34(2014), 76-80.
[12] B.C. Dhage, S.B. Dhage, J.R. Graef, Dhage iteration method for initial value problems for nonlinear first order hybrid integrodifferential equations, J. Fixed Point Theory Appl., 17(2016): DOI 10.1007/s11784-015-0279-3.
[13] J.J. Nieto, R. Rodriguez-Lopez, Contractive mappings theorems in partially ordered sets and applications to ordinary differential equations, Order, 22(2005), 223-239.
[14] D. Otrocol, Properties of the solutions of system of differential equations with maxima, via weakly Picard operator theory, Communications in Applied Analysis, 17(2013), No. 1, 99-107.
[15] D. Otrocol, Systems of functional differential equations with maxima, of mixed type, Electron. J. Qual. Theory Differ. Equ., 2014 (2014), No. 5, 1-9.
[16] D. Otrocol, I.A. Rus, Functional-differential equations with “maxima” via weakly Picard operators theory, Bull. Math. Soc. Sci. Math. Roumanie, 51(99)(2008), No. 3, 253-261.
[17] D. Otrocol, I.A. Rus, Functional-differential equations with maxima of mixed type argument, Fixed Point Theory, 9 (2008), No. 1, pp. 207-220.
[18] S. Heikkilä and V. Lakshmikantham, Monotone Iterative Techniques for Discontinuous Nonlinear Differential Equations, Marcel Dekker inc., New York 1994.
[19] A. Petruşel, I.A. Rus, Fixed point theorems in ordered L-spaces, Proc. Amer. Math. Soc., 134(2006), 411-418.
[20] A.C.M. Ran, M.C. Reurings, A fixed point theorem in partially ordered sets and some applications to matrix equations, Proc. Amer. Math. Soc., 132 (2004), 1435-1443.

Received: ; Accepted:

2018

Approximation with arbitrary order by certain linear positive operators

March 13, 2023

Theoretical models of hematopoietic cell dynamics related to bone marrow transplantation

May 4, 2022

Analysis of a planar differential system arising from hematology

April 26, 2022

	$\displaystyle\left\|\mathcal{A}x(t)\right\|$	$\displaystyle\leq\left\|\alpha_{0}e^{-\lambda t}\right\|+e^{-\lambda t}\int_{a}^{t}e^{-\lambda s}\left\|\widetilde{f}(s,x(s))\right\|\,ds$
		$\displaystyle\leq\left\|\alpha_{0}\right\|+K\int_{a}^{b}e^{\lambda s}\,ds$
		$\displaystyle\leq\left\|\alpha_{0}\right\|+e^{\lambda a}K(b-a),$

	$\displaystyle\left\|\mathcal{A}x(t)-\mathcal{A}y(t)\right\|$	$\displaystyle\leq\left\|e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\left[\widetilde{f}(s,x(s))-\widetilde{f}(s,y(s))\right]\,ds\right\|$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\mu\left\|x(s)-y(s)\right\|\,ds$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}e^{\lambda s}\lambda\left\|x(s)-y(s)\right\|\,ds$
		$\displaystyle\leq e^{-\lambda t}\int_{a}^{t}\frac{d}{ds}e^{\lambda s}\left\\|x-y\right\\|\,ds$
		$\displaystyle\leq(1-e^{-\lambda a})\left\\|x-y\right\\|,$

	$\displaystyle\left\|\mathcal{B}x_{n}(t_{2})-\mathcal{B}x_{n}(t_{1})\right\|$	$\displaystyle=\left\|e^{-\lambda t_{2}}\int_{a}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right.$
		$\displaystyle\quad\quad\quad\left.-e^{-\lambda t_{1}}\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\leq\left\|\left(e^{-\lambda t_{2}}-e^{-\lambda t_{1}}\right)\int_{a}^{t_{1}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\quad+\left\|e^{-\lambda t_{2}}\int_{t_{1}}^{t_{2}}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x_{n}(\xi)\right)\,ds\right\|$
		$\displaystyle\rightarrow 0\quad\text{ as }\quad t_{2}\rightarrow t_{1},$

	$\displaystyle\left\|y(t)\right\|$	$\displaystyle=\left\|\mathcal{B}x(t)\right\|$
		$\displaystyle=\left\|e^{-\lambda t}\int_{a}^{t}e^{\lambda s}g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,ds\right\|$
		$\displaystyle\leq\int_{a}^{t}e^{\lambda s}\left\|g\left(s,\underset{a\leq\xi\leq s}{\max}\,x(\xi)\right)\,\right\|\,ds$
		$\displaystyle\leq\int_{a}^{b}e^{\lambda b}M\,ds$
		$\displaystyle\leq e^{\lambda b}M(b-a)=r,$

Dhage iteration method for approximating solutions of nonlinear differential equations with maxima,

Abstract

Authors

Keywords

Cite this paper as:

References

PDF

About this paper

Journal

Publisher Name

DOI

Print ISSN

Online ISSN

MR

ZBL

Google Scholar

Dhage iteration method for approximating solutions of nonlinear differential equations with maxima

1. Introduction

2. Preliminaries

Definition 2.1.

Definition 2.2.

Definition 2.3.

Definition 2.4.

Remark 2.1.

Definition 2.5.

Definition 2.6.

Theorem 2.1 (Dhage [4]).

Remark 2.2.

Remark 2.3.

3. Main Result

Lemma 3.1.

Proof.

Definition 3.1.

Remark 3.1.

Lemma 3.2.

Theorem 3.1.

Proof.

Remark 3.2.

Remark 3.3.

Example 3.1.

Remark 3.4.

Acknowledgment

References

Related Posts