Next Article in Journal
The Geometry of Dynamic Time-Dependent Best–Worst Choice Pairs
Next Article in Special Issue
Investigation of the Oscillatory Behavior of the Solutions of a Class of Third-Order Delay Differential Equations with Several Terms
Previous Article in Journal
A Framework for I*-Statistical Convergence of Fuzzy Numbers
Previous Article in Special Issue
On Some New Dynamic Hilbert-Type Inequalities across Time Scales
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Qualitative Outcomes on Monotone Iterative Technique and Quasilinearization Method on Time Scale

1
Department of Mathematics, Faculty of Sciences, Sakarya University, 54100 Sakarya, Turkey
2
Department of Mathematics, Faculty of Sciences, Gebze Technical University, 41400 Kocaeli, Turkey
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 640; https://doi.org/10.3390/axioms13090640
Submission received: 25 August 2024 / Revised: 13 September 2024 / Accepted: 14 September 2024 / Published: 19 September 2024
(This article belongs to the Special Issue Infinite Dynamical System and Differential Equations)

Abstract

:
In this paper, a nonlinear dynamic equation with an initial value problem (IVP) on a time scale is considered. First, applying comparison results with a coupled lower solution (LS) and an upper solution (US), we improved the quasilinearization method (QLM) for the IVP. Unlike other studies, we consider the LS and US pair of the seventh type instead of the natural type. It was determined that the solutions of the dynamic equation converge uniformly and monotonically to the unique solution of the IVP, and the convergence is quadratic. Moreover, we will use the delta derivative ( Δ γ ) instead of the classical derivative ( d γ ) in the proof because it studies a time scale. In the second part of the paper, we applied the monotone iterative technique (MIT) coupled with the LS and US, which is an effective method, proving a clear analytical representation for the solution of the equation when the relevant functions are monotonically non-decreasing and non-increasing. Then an example is given to illustrate the results obtained.

1. Introduction

A time scale is a system that combines discrete and continuous analysis. Stefan Hilger was the first to study this subject in 1990 [1]. Although it is a relatively new subject, it has become a field that has attracted more and more attention due to dynamic equations that allow differential equations and difference equations to be considered simultaneously. Some books, articles, and other studies on this subject are given with their references.
Kaymakçalan and Lawrence in [2] discussed the monotone iterative technique with the classical method for the IVP in a unified setting.
Akın et al. [3] used the QLM approximately on a time scale for the unique solution of the BVP that divides from above and below with the monotone convergent sequences of the LS and US.
Bhaskar and McRae [4] proved a fundamental theorem regarding the existence of coupled maximal and minimal solutions of dynamical systems.
West and Vatsala [5] examined some theorems for the type I pair of the LS and US. They obtained natural monotone sequences starting from the coupled LS and US of type I of the equation and observed that the results varied based on the iterative steps used to develop the arrays.
In their book, Lakshmikantham and Vatsala [6] presented a systematic development of a generalized QLM, illustrating the concepts and technical challenges encountered in the combined approach. They significantly increased the usefulness of the QLM, which has proven to be very effective in various research fields and applications.
Denton and Vatsala [7] presented comparison results of a nonlinear Riemann–Liouville fractional differential equation. The authors improved a monotone method for finite systems of q-order fractional equations using the coupled LS and US.
Ramírez and Vatsala [8] investigated a generalized MIT using the pair of the LS and US of Caputo fractional differential equations. They developed results that yield natural monotone sequences or convolute monotone sequences that converge uniformly and monotonically to minimal and maximal solutions.
Khavanin [9] extended the mixed monotony method to generate monotone sequences that converge to the unique solution of the IVP with a delay difference equation.
Daneev and Sizykh [10] proposed a new approach to real-time multifunctional automatic control systems depending on the combined use of dynamic technologies.
Jyoti and Singh [11] demonstrated a combined iterative approach based on the QLM and the Krasnoselskii–Mann approximation to evaluate solutions of nonlinear Dirichlet BVPs. In addition to this, the authors reduced the nonlinear problems to a set of linear equations using the QLM with Green’s function.
Wu and Shu [12] introduced an improved framework based on some important insights from geometry that is called geometric QLM. They proposed the fundamental theory of the QLM through the geometric features of convex areas to construct sequences.
Heydari et al. [13] presented a combination of quasilinearization to approximate the solution of nonlinear functional Volterra integral equations. They solved the linear integral equation, which is obtained from each iteration using the Legendre ordering method.
Izadi and Roul [14] offered a Vieta–Fibonacci matrix technique to find the solution of nonlinear and multiple singularity third-order Emden–Fowler equations. The authors applied the QLM to the underlying model problem to obtain an efficacious approach.
Verma and Urus [15] improved an MIT with lower and upper solutions for four-point nonlinear BVPs. They established an iterative design using quasilinearization and under convenient conditions, which proved that the resulting sequences converge uniformly to a solution in a given region.
Idiz et al. [16] investigated numerical solutions of fractional Lane–Emden-type equations that arise in astrophysics applications. They offered a numerical approach with the QLM using Legendre wavelets.
Yakar and Arslan [17] created new definitions for a causal extreme value problem with Riemann–Liouville fractional derivatives and examined the unique solution by combining the QLM.
Consider the following dynamic non-linear IVP:
u Δ = f γ , u + g γ , u + h γ , u u 0 = u 0 ,
in the particular case of f γ , u , g γ , u , h γ , u C r d J × R , R , where J = [ 0 , T ] .

2. Preliminaries

This part provides some basic concepts regarding the calculus on time scales and fundamental theorems needed in the next sections. We will introduce some useful tools for proofs.
Definition 1
([18]). Any nonempty closed subset of the real numbers is called a time scale. A set of real numbers, integers, and natural numbers, [ 4 , 5 ] [ 6 , 7 ] , are examples of a time scale. In contrast, sets of rational and irrational numbers are not a time scale. In addition to this, u Δ = d u d γ is the usual derivative if the time scale is taken as R.
Definition 2
([6]). A function ϰ is called LS of
u Δ = f γ , u u 0 = u 0
if
ϰ Δ f γ , ϰ , ϰ γ 0 u 0 ,
and similarly, ϱ is called US if
ϱ Δ f γ , ϱ , ϱ γ 0 u 0 .
Definition 3
([6]). Let ϰ and ϱ be rd-continuously differentiable functions such that ϰ ϱ on T. Then ϰ and ϱ are called a couple of seventh-type LS and US of (1) if
ϰ Δ f γ , ϱ + g γ , ϱ + h γ , ϱ ϰ γ 0 u 0
and
ϱ Δ f γ , ϰ + g γ , ϰ + h γ , ϰ ϱ γ 0 u 0
and called a couple of third-type LS and US if
ϰ Δ f γ , ϰ + g γ , ϰ + h γ , ϱ ϰ γ 0 u 0
ϱ Δ f γ , ϱ + g γ , ϱ + h γ , ϰ ϱ γ 0 u 0 .
Theorem 1
([6]). Let  ϰ , ϱ be the third-type LS and US of (1), respectively, and assume that
f γ , u 1 f γ , u 2 L u 1 u 2 ,
u 1 u 2 for some L > 0 whenever ϰ γ 0 ϱ γ 0 . Then this implies ϰ γ ϱ γ .
Theorem 2
([6]). Let  ϰ , ϱ be the LS and US of (1), respectively, such that ϰ γ ϱ γ . Then there exists a solution u γ of (1), satisfying ϰ γ u γ ϱ γ and ϰ γ 0 u γ 0 ϱ γ 0 .
Theorem 3
([19] Arzela–Ascoli theorem). Let f n γ be a sequence of functions defined on a compact set J, which is equicontinuous and uniformly bounded on J. Then there exists a subsequence f n , n = 1 , 2 , , which is uniformly convergent on J.

3. Results

3.1. The Quasilinearization Method for the Problem

Let
Ω 0 = γ , u : ϰ 0 γ u ϱ 0 γ γ T
where ϰ 0 γ , ϱ 0 γ C r d J , R . Under some convenient conditions by using comparison theorems with the QLM, the monotone sequences that converge to the solution of (1) were obtained and showed the convergence rate.
Theorem 4.
Suppose that the following hypotheses hold:
(A1) Let ϰ 0 γ , ϱ 0 γ C r d J , R be coupled with the seventh-type LS and US of (1) such that ϰ 0 γ ϱ 0 γ on J.
(A2) Assume that f u u γ , u 0 , g u u γ , u + ϕ u u γ , u 0 whenever ϕ u u exists and ϕ u u γ , u > 0 and h u u γ , u 0 on Ω 0 , where f γ , u , g γ , u , h γ , u C r d 2 J × R , R and
(A3) f u ( γ , u ) 0 , g u ( γ , u ) 0 and h u ( γ , u ) 0 .
Then there exist monotone sequences ϰ n γ and ϱ n γ , which converge to the unique solution of (1) uniformly and monotonically, and the convergence is quadratic.
Proof of Theorem 4.
In view of the condition (A2), the Lipschitz conditions below are satisfied. For ϰ 0 γ u 2 u 1 ϱ 0 γ ,
L u 1 u 2 f γ , u 1 f γ , u 2 L u 1 u 2 L u 1 u 2 g γ , u 1 g γ , u 2 L u 1 u 2 L u 1 u 2 h γ , u 1 h γ , u 2 L u 1 u 2 , L > 0
Since f u u γ , u 0 , g u u γ , u + ϕ u u γ , u 0 and h u u γ , u 0 , the following inequalities, which are very important for the proof, can be written, respectively. For u v ,
f γ , u f γ , v + f u γ , v u v ,
g γ , u g γ , v + g u γ , v + ϕ u γ , v u v ϕ γ , u ϕ γ , v
and
h γ , u h γ , v + h u γ , u u v .
Consider auxiliary linear IVP as follows:
u Δ = w γ , ϰ 0 , ϱ 0 ; v = f γ , ϱ 0 + f u γ , ϱ 0 v ϱ 0 + g γ , ϱ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 v ϱ 0 ϕ γ , v ϕ γ , ϱ 0 + h γ , ϱ 0 + h u γ , ϱ 0 v ϱ 0 , u 0 = u 0 ,
v Δ = W γ , ϰ 0 , ϱ 0 ; u = f γ , ϰ 0 + f u γ , ϰ 0 u ϰ 0 + g γ , ϰ 0 + g u γ , ϰ 0 + ϕ u γ , ϰ 0 u ϰ 0 ϕ γ , u ϕ γ , ϰ 0 + h γ , ϰ 0 + h u γ , ϰ 0 u ϰ 0 , v 0 = u 0
Due to the condition (A1) in the hypothesis, it is clear that ϰ 0 Δ f γ , ϱ 0 + g γ , ϱ 0 + h γ , ϱ 0 = w γ , ϰ 0 , ϱ 0 ; ϱ 0 . Namely, ϰ 0 is an LS for (7). Similarly, it is known to be ϱ 0 Δ f γ , ϰ 0 + g γ , ϰ 0 + h γ , ϰ 0 , and using the inequalities given in (4)–(6), we can write
ϱ 0 Δ f γ , ϱ 0 + f u γ , ϱ 0 ϰ 0 ϱ 0 + g γ , ϱ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϰ 0 ϱ 0 ϕ γ , ϰ 0 ϕ γ , ϱ 0 + h γ , ϱ 0 + h u γ , ϰ 0 ϰ 0 ϱ 0 f γ , ϱ 0 + f u γ , ϱ 0 ϰ 0 ϱ 0 + g γ , ϱ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϰ 0 ϱ 0 ϕ γ , ϰ 0 ϕ γ , ϱ 0 + h γ , ϱ 0 + h u γ , ϱ 0 ϰ 0 ϱ 0 = w γ , ϰ 0 , ϱ 0 ; ϰ 0
Hence, ϱ 0 is a US for (7). Then there exists a unique solution of (7) called ϰ 1 γ such that ϰ 0 γ ϰ 1 γ ϱ 0 γ with the help of Theorem 2. Similarly, to guarantee the solution ϱ 1 γ , we will show that ϰ 0 and ϱ 0 are the LS and US of (8), respectively. Since ϰ 0 and ϱ 0 are the seventh-type LS and US, we have ϰ 0 Δ f γ , ϱ 0 + g γ , ϱ 0 + h γ , ϱ 0 . When the inequalities obtained from (4)–(6) is used, the following expression is obtained:
ϰ 0 Δ f γ , ϰ 0 + f u γ , ϰ 0 ϱ 0 ϰ 0 + g γ , ϰ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϱ 0 ϰ 0 ϕ γ , ϱ 0 ϕ γ , ϰ 0 + h γ , ϰ 0 + h u γ , ϰ 0 ϱ 0 ϰ 0 = f γ , ϰ 0 + f u γ , ϱ 0 ϱ 0 ϰ 0 + g γ , ϰ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϱ 0 ϰ 0 ϕ γ , ϱ 0 ϕ γ , ϰ 0 + h γ , ϰ 0 + h u γ , ϰ 0 ϱ 0 ϰ 0 = W γ , ϰ 0 , ϱ 0 ; ϱ 0
and
ϱ 0 Δ f γ , ϰ 0 + g γ , ϰ 0 + h γ , ϰ 0 = W γ , ϰ 0 , ϱ 0 ; ϰ 0 .
Thus, ϰ 0 and ϱ 0 are the LS and US of (8); then there exists a unique solution of (8) called ϱ 1 such that ϰ 0 γ ϱ 1 γ ϱ 0 γ . Now we have to show that ϰ 1 ϱ 1 . By using (7),
ϰ 1 Δ = w γ , ϰ 0 , ϱ 0 ; ϱ 1 = f γ , ϱ 0 + f u γ , ϱ 0 ϱ 1 ϱ 0 + g γ , ϱ 0 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϱ 1 ϱ 0 ϕ γ , ϱ 1 ϕ γ , ϱ 0 + h γ , ϱ 0 + h u γ , ϱ 0 ϱ 1 ϱ 0
can be written. If the inequalities obtained from (4)–(6) are substituted into the equation above, we obtain
ϰ 1 Δ f γ , ϱ 1 + f u γ , ϱ 0 ϱ 0 ϱ 1 + f u γ , ϱ 0 ϱ 1 ϱ 0 + g γ , ϱ 1 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϱ 0 ϱ 1 ϕ γ , ϱ 0 ϕ γ , ϱ 1 + g u γ , ϱ 0 + ϕ u γ , ϱ 0 ϱ 1 ϱ 0 ϕ γ , ϱ 1 ϕ γ , ϱ 0 + h γ , ϱ 1 + h u γ , ϱ 1 ϱ 0 ϱ 1 + h u γ , ϱ 0 ϱ 1 ϱ 0 f γ , ϱ 1 + g γ , ϱ 1 + h γ , ϱ 1 .
Similarly, from (8),
ϱ 1 Δ = W γ , ϰ 0 , ϱ 0 ; ϰ 1 = f γ , ϰ 0 + f u γ , ϰ 0 ϰ 1 ϰ 0 + g γ , ϰ 0 + g u γ , ϰ 0 + ϕ u γ , ϰ 0 ϰ 1 ϰ 0 ϕ γ , ϰ 1 ϕ γ , ϰ 0 + h γ , ϰ 0 + h u γ , ϰ 0 ϰ 1 ϰ 0 .
By using (4)–(6),
ϱ 1 Δ f γ , ϰ 1 + f u γ , ϰ 1 ϰ 0 ϰ 1 + f u γ , ϰ 0 ϰ 1 ϰ 0 + g γ , ϰ 1 + g u γ , ϰ 1 + ϕ u γ , ϰ 1 ϰ 0 ϰ 1 ϕ γ , ϰ 0 ϕ γ , ϰ 1 + g u γ , ϰ 0 + ϕ u γ , ϰ 0 ϰ 1 ϰ 0 ϕ γ , ϰ 1 ϕ γ , ϰ 0 + h γ , ϰ 1 + h u γ , ϰ 0 ϰ 0 ϰ 1 + h u γ , ϰ 0 ϰ 1 ϰ 0 = f γ , ϰ 1 + g γ , ϰ 1 + h γ , ϰ 1 .
Since ϰ 1 , ϱ 1 are the seventh-type LS and US with ϰ 1 ( 0 ) ϱ 1 ( 0 ) , then by Theorem 1, we obtain ϰ 1 ( γ ) ϱ 1 ( γ ) . Consequently, the inequality below can be written as
ϰ 0 γ ϰ 1 γ ϱ 1 γ ϱ 0 γ .
Contemplate the following initial value problems to obtain the next level:
u Δ = w γ , ϰ 1 , ϱ 1 ; v = f γ , ϱ 1 + f u γ , ϱ 1 v ϱ 1 + g γ , ϱ 1 + g u γ , ϱ 1 + ϕ u γ , ϱ 1 v ϱ 1 ϕ γ , v ϕ γ , ϱ 1 + h γ , ϱ 1 + h u γ , ϱ 1 v ϱ 1 , u 0 = u 0 ,
v Δ = W γ , ϰ 1 , ϱ 1 ; u = f γ , ϰ 1 + f u γ , ϰ 1 u ϰ 1 + g γ , ϰ 1 + g u γ , ϰ 1 + ϕ u γ , ϰ 1 u ϰ 1 ϕ γ , u ϕ γ , ϰ 1 + h γ , ϰ 1 + h u γ , ϰ 1 u ϰ 1 , v 0 = u 0 .
By (7), it is clear that
ϰ 1 Δ f γ , ϱ 1 + g γ , ϱ 1 + h γ , ϱ 1 = w γ , ϰ 1 , ϱ 1 ; ϱ 1 ,
and with (4)–(6), we have
ϱ 1 Δ f γ , ϰ 1 + g γ , ϰ 1 + h γ , ϰ 1 f γ , ϱ 1 + g u γ , ϱ 1 ϰ 1 ϱ 1 + g γ , ϱ 1 + g u γ , ϱ 1 + ϕ u γ , ϱ 1 ϰ 1 ϱ 1 ϕ γ , ϰ 1 ϕ γ , ϱ 1 + h γ , ϱ 1 + h u γ , ϱ 1 ϰ 1 ϱ 1 = w γ , ϰ 1 , ϱ 1 ; ϰ 1 .
Therefore, ϰ 1 and ϱ 1 are the LS and US that are the seventh type of (9). In a similar way, it can be shown that ϰ 1 and ϱ 1 are the coupled seventh-type lower and upper solutions of (10). Based on Theorem 2, there exist the unique solutions ϰ 2 and ϱ 2 of (9) and (10), respectively, such that ϰ 1 γ ϰ 2 γ ϱ 1 γ and ϰ 1 γ ϱ 2 γ ϱ 1 γ . Therefore, it can be obtained as follows:
ϰ 0 γ ϰ 1 γ ϰ 2 γ ϱ 2 γ ϱ 1 γ ϱ 0 γ
If this process continues like this, we can write
ϰ 0 γ ϰ 1 γ ϰ 2 γ · · · ϰ n γ ϱ n γ · · · ϱ 2 γ ϱ 1 γ ϱ 0 γ .
Each element of ϰ n γ and ϱ n γ sequences is respectively a solution of the following system:
ϰ n + 1 Δ = w γ , ϰ n , ϱ n ; ϱ n + 1 , ϰ n + 1 0 = u 0
ϱ n + 1 Δ = W γ , ϰ n , ϱ n ; ϰ n + 1 , ϱ n + 1 0 = u 0 .
Because the solutions ϰ 0 γ and ϱ 0 γ are bounded, the sequences ϰ n γ and ϱ n γ are also uniformly bounded on T. Moreover, it can be easily shown that ϰ 0 γ and ϱ 0 γ are equicontinuous. According to Theorem 3, the sequences converge to the unique solution of (1) uniformly. Now, it will be shown that this convergence is quadratic. For this purpose, define
p n γ = u γ ϰ n γ 0 p n 0 = 0 ,
q n γ = ϱ n γ u γ 0 q n 0 = 0 .
If we take the delta derivative of both sides of (11) defined above and take g γ , u + ϕ γ , u = G γ , u then we have
p n Δ = [ f γ , u f γ , ϱ n 1 ] [ f u γ , ϱ n 1 ] ϱ n ϱ n 1 + [ G γ , u G γ , ϱ n 1 ] G u γ , ϱ n 1 ϱ n ϱ n 1 + ϕ γ , ϱ n ϕ γ , u g u γ , ϱ n 1 + ϕ u γ , ϱ n 1 ϱ n ϱ n 1 + [ h γ , u h γ , ϱ n 1 ] h u γ , ϰ n 1 ϱ n ϱ n 1 .
If the definitions of ϰ n and ϱ n are taken into account and applying the mean value theorem then the above expression can be written as
p n Δ = f u γ , a ϱ n 1 u [ f u γ , ϱ n 1 ] ϱ n ϱ n 1 G u γ , b ϱ n 1 u + ϕ u γ , c ϱ n u G u γ , ϱ n 1 ϱ n ϱ n 1 h u γ , d ϱ n 1 u h u γ , ϰ n 1 ϱ n ϱ n 1 f 1 u ( γ , u ) q n 1 [ f 1 u γ , ϱ n 1 ] q n q n 1 G u γ , u q n 1 + ϕ u γ , u q n G u γ , ϱ n 1 q n q n 1 f 3 u γ , ϱ n 1 q n 1 f 3 u γ , ϰ n 1 q n q n 1
where u < a < ϱ n 1 , u < b < ϱ n 1 , u < c < ϱ n ,   u < d < ϱ n 1 . Therefore,
p n Δ q n 1 [ f u γ , ϱ n 1 f u ( γ , u ) ] + q n 1 [ G u γ , ϱ n 1 G u γ , u ] + q n 1 [ h u γ , ϰ n 1 h u γ , ϱ n 1 ] + q n [ f u γ , ϱ n 1 + ϕ u γ , u g u γ , u ϕ u γ , u h u γ , ϰ n 1 q n 1 2 f u u γ , δ + q n 1 2 G u u γ , ζ ( q n 1 ) h u u ( γ , η ) ϱ n 1 ϰ n 1 + q n [ f u γ , ϱ n 1 g u γ , u h u γ , ϰ n 1 ] = q n 1 2 f u u γ , δ + q n 1 2 G u u γ , ζ ( q n 1 ) h u u ( γ , η ) q n 1 + p n 1 + q n [ f u γ , ϱ n 1 g u γ , u h u γ , ϰ n 1 ] = q n 1 2 [ f u u γ , δ + G u u γ , ζ h u u ( γ , η ) ] + q n [ f u γ , ϱ n 1 g u γ , u h u γ , ϰ n 1 ] ( p n 1 ) ( q n 1 ) h u u ( γ , η ) ,
where u < δ , ζ < ϱ n 1 , ϰ n 1 < η < ϱ n . Since continuous functions are bounded in a closed set, we can write
p n Δ q n 1 2 A + q n B ( p n 1 ) ( q n 1 ) C ,
where A, B, and C are positive constants and
f u u γ , δ + G u u γ , ζ h u u ( γ , η ) A , f u γ , ϱ n 1 + g u γ , u + h u γ , ϰ n 1 B , h u u ( γ , η ) C .
When Cauchy inequality is applied to the term ( p n 1 ) ( q n 1 ) C , then
p n Δ q n 1 2 A + q n B p n 1 2 2 C q n 1 2 2 C = p n 1 2 2 C + q n 1 2 [ A C 2 ] + q n B
is obtained. Similarly, taking the delta derivative of both sides of (12),
q n Δ = ϱ n Δ u Δ = [ f γ , ϰ n 1 f γ , u ] + f u γ , ϱ n 1 ϰ n ϰ n 1 + [ g γ , ϰ n 1 + ϕ γ , ϰ n 1 ] [ g γ , u + ϕ γ , u ] + [ ϕ γ , u ϕ γ , ϰ n ] + g u γ , ϰ n 1 + ϕ u γ , ϰ n 1 ϰ n ϰ n 1 + [ h γ , ϰ n 1 h γ , u ] + h u γ , ϰ n 1 ϰ n ϰ n 1
can be written. Similarly, by applying the above procedures, we obtain
q n Δ p n 1 f u u γ , δ ϱ n 1 ϰ n 1 + p n [ f u γ , ϱ n 1 G u γ , ϰ n 1 + ϕ u γ , u h u γ , ϰ n 1 ] p n 1 h u u ( γ , η ) u ϰ n 1 = p n 1 2 [ f u u γ , δ h u u ( γ , η ) ] + ( p n 1 ) ( q n 1 ) f u u γ , δ + p n [ f u γ , ϱ n 1 G u γ , ϰ n 1 + ϕ u γ , u h u γ , ϰ n 1 ] .
Hence, one can see
q n Δ p n 1 2 K + ( p n 1 ) ( q n 1 ) L + p n M ,
where K, L,and M are positive constants and
f u u γ , δ h u u ( γ , η ) K , f u u γ , δ L , f u γ , ϱ n 1 + G u γ , ϰ n 1 + ϕ u γ , u + h u γ , ϰ n 1 M .
If Cauchy inequality is implemented to the term ( p n 1 ) ( q n 1 ) L , then
q n Δ p n 1 2 K + p n 1 2 2 L + q n 1 2 2 L + p n M = p n 1 2 [ K + L 2 ] + q n 1 2 2 L + p n M .
When (13) and (14) are considered together, it can be written as
r n Δ P r n + R r n 1 2 ,
which is linear in r n , where r n = p n q n , P = 0 B M 0 , R = C 2 A + C 2 K + L 2 L 2 . Now, if Gronwall’s inequality is applied and integrated from 0 to γ , then it can be obtained that
r n e P γ r n 0 + 0 γ R r n 1 2 e P s Δ s = 0 γ R r n 1 2 e P γ s Δ s .
Consequently, we have
max J r n max J R r n 1 2 0 γ e P γ s Δ s max J r n 1 2 R e P γ 1 P max J r n 1 2 R e P γ P .
That is,
max J u ϰ n R e P γ P max J u ϰ n 2 , max J u ϱ n R e P γ P max J ϱ n u 2 .
This result shows that the sequences ϰ n γ and ϱ n γ converge to the unique solution u γ of (1) quadratically, and this completes the proof. □

3.2. The Monotone Iterative Technique for the Problem

Let us consider the following IVP:
u Δ = f γ , u + g γ , u h γ , u u 0 = u 0 .
It has been determined that the sequences obtained with the monotone iteration technique converge to the extremal solutions of the problem (15) uniformly and monotonically.
Theorem 5.
(A1) Let ϰ 0 γ , ϱ 0 γ C r d T , R be coupled with the third-type LS and US of (15) such that ϰ 0 γ ϱ 0 γ on J.
(A2) Suppose that the functions f γ , u , g γ , u , h γ , u are non-decreasing with respect to u. Then there exist the monotone sequences ϰ n γ and ϱ n γ converging uniformly and monotonically to the extremal solutions ρ (minimal solution) and r (maximal solution), respectively, which are solutions of (15).
Proof of Theorem 5.
We note that since ϰ 0 γ , ϱ 0 γ are the couple of the third-type LS and US, then
ϰ 0 Δ f γ , ϰ 0 + g γ , ϰ 0 h γ , ϱ 0 , ϰ 0 0 u 0 , ϱ 0 Δ f γ , ϱ 0 + g γ , ϱ 0 h γ , ϰ 0 , ϱ 0 0 u 0 .
Let us consider the following systems:
ϰ m + 1 Δ = f γ , ϰ m + g γ , ϰ m h γ , ϱ m , ϰ m + 1 0 = u 0 , ϱ m + 1 Δ = f γ , ϱ m + g γ , ϱ m h γ , ϰ m , ϱ m + 1 0 = u 0 .
Now, it will be shown that ϰ 0 ϰ 1 . For this purpose, define p = ϰ 0 ϰ 1 . It is trivial that p 0 = ϰ 0 0 ϰ 1 0 0 . If we take the delta derivative of both sides, then we have
p Δ = ϰ 0 Δ ϰ 1 Δ f γ , ϰ 0 + g γ , ϰ 0 h γ , ϱ 0 f γ , ϰ 0 g γ , ϰ 0 + h γ , ϱ 0 = 0 .
Due to p Δ 0 with p 0 0 , then based on the comparison theorem, p γ 0 . Hence, ϰ 0 ϰ 1 . Similarly, we show ϱ 1 ϱ 0 . Let q = ϱ 1 ϱ 0 , q 0 = ϱ 1 0 ϱ 0 0 0 . Taking the delta derivative, we obtain
q Δ = ϱ 1 Δ ϱ 0 Δ f γ , ϱ 0 + g γ , ϱ 0 h γ , ϰ 0 f γ , ϱ 0 + g γ , ϱ 0 h γ , ϰ 0 = 0 .
Therefore, we obtain q Δ 0 , p 0 0 , and by the comparison theorem, q 0 = ϱ 1 0 ϱ 0 0 0 ϱ 1 ϱ 0 . Now, it will be shown that ϰ 1 ϱ 1 . For this, let p = ϰ 1 ϱ 1 ; then p 0 = ϰ 1 0 ϱ 1 0 = 0 . Taking the delta derivative of both sides and using the properties of the functions f , g , h , we can write the following:
p Δ = ϰ 1 Δ ϱ 1 Δ = f γ , ϰ 0 + g γ , ϰ 0 h γ , ϱ 0 f γ , ϱ 0 g γ , ϱ 0 + h γ , ϰ 0 0 .
Since p Δ 0 , p 0 0 , again by the comparison theorem, p γ 0 . Therefore, ϰ 1 ϱ 1 . Now, let us show ϰ 1 u ϱ 1 . Let q = u ϱ 1 and p = ϰ 1 u with q 0 = u 0 ϱ 1 0 = 0 and p 0 = ϰ 1 0 u 0 = 0 . Since the functions are non-decreasing, we obtain
q Δ = u Δ ϱ 1 Δ = f γ , u + g γ , u h γ , u f γ , ϱ 0 g γ , ϱ 0 + h γ , ϰ 0 0 .
and
p Δ = ϰ 1 Δ u Δ = f γ , ϰ 0 + g γ , ϰ 0 h γ , ϱ 0 f γ , u g γ , u + h γ , u 0 .
Hence,
ϰ 0 ϰ 1 u ϱ 1 ϱ 0 .
We will use mathematical induction to generalize this. For k = 1 , it is clear that ϰ k 1 ϰ k u ϱ k ϱ k 1 is true. Assume that, for k > 1 , ϰ k 1 ϰ k u ϱ k ϱ k 1 is satisfied. We will show that ϰ k ϰ k + 1 u ϱ k + 1 ϱ k is also satisfied. For this, we define q = ϱ k + 1 ϱ k and p = ϰ k ϰ k + 1 . We take the delta derivative on both sides, respectively. Then we have
q Δ = ϱ k + 1 Δ ϱ k Δ f γ , ϱ k + g γ , ϱ k h γ , ϰ k f γ , ϱ k 1 + g γ , ϱ k 1 h γ , ϰ k 1 0
and
p Δ = ϰ k Δ ϰ k + 1 Δ f γ , ϰ k 1 + g γ , ϰ k 1 h γ , ϱ k 1 f γ , ϰ k g γ , ϰ k + h γ , ϱ k 0 .
Therefore, we can say that ϱ k + 1 ϱ k and ϰ k ϰ k + 1 . Following similar steps, it can be shown that u ϱ k + 1 and ϰ k + 1 u . Therefore, the following conclusion is reached:
ϰ 0 ϰ 1 · · · ϰ m ϰ m + 1 u ϱ m + 1 ϱ m · · · ϱ 1 ϱ 0 .
To show that these sequences converge uniformly, we must show that they are uniformly bounded and equicontinuous. Since ϰ 0 and ϱ 0 are bounded, then there exists L 1 , L 2 > 0 such that, for every γ N , ϰ 0 L 1 and ϱ 0 L 2 . Since, for every m > 0 , ϰ 0 ϱ m ϱ 0 , then it can be written that 0 ϱ m ( γ ) ϰ 0 ( γ ) ϱ 0 ( γ ) ϰ 0 ( γ ) . Therefore, ϱ m ( γ ) M , M > 0 and ϱ m ( γ ) are uniformly bounded. Similarly, it can be shown that ϰ m ( γ ) is uniformly bounded. Now, we will show that ϱ m ( γ ) is equicontinuous. For this, we will find any δ > 0 when γ 1 γ 2 < δ , ϰ m ( γ 1 ) ϰ m ( γ 2 ) < ε is satisfied. For 0 γ 1 γ 2 , m > 0
ϱ m ( γ 1 ) ϱ m ( γ 2 ) = u 0 + 0 γ 1 f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s u 0 0 γ 2 f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s = γ 2 γ 1 f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s γ 2 γ 1 f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s .
Since ϱ m ( γ ) and ϰ m ( γ ) are uniformly bounded and the functions f ( γ , u ) ,   g ( γ , u ) , h ( γ , u ) are bounded, then there exists Γ ˜ such that
ϰ m ( γ 1 ) ϰ m ( γ 2 ) γ 2 γ 1 f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s Γ ˜ γ 1 γ 2 .
If we choose δ = ε Γ ˜ , then ϱ m ( γ ) is equicontinuous. Similarly, it can be shown that ϰ m ( γ ) is equicontinuous. By Theorem 3, the sequences ϰ m k ( γ ) and ϱ m k ( γ ) converge to the extremal solutions ρ γ and r γ of (15) uniformly and monotonically. Finally, it will be shown that ρ γ and r γ are coupled with third-type minimal and maximal solutions of (15). From previous findings, we have
lim m ϰ m ( γ ) = ρ ( γ ) , lim m ϱ m ( γ ) = r ( γ ) .
If we integrate the system (16) from 0 to γ , then it can be obtained that
ϰ m + 1 ( γ ) = u 0 + 0 γ f ( s , ϰ m ( s ) ) + g ( s , ϰ m ( s ) ) h ( s , ϱ m ( s ) ) Δ s , ϱ m + 1 ( γ ) = u 0 + 0 γ f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s
and take the limit while m ,
lim m ϰ m + 1 ( γ ) = lim m u 0 + 0 γ f ( s , ϰ m ( s ) ) + g ( s , ϰ m ( s ) ) h ( s , ϱ m ( s ) ) Δ s ρ = u 0 + 0 γ f ( s , ρ ( s ) ) + g ( s , ρ ( s ) ) h ( s , r ( s ) ) Δ s .
After taking the delta derivative of both sides, then we obtain
ρ Δ = f γ , ρ + g γ , ρ h γ , r .
Similarly, it can be written that
lim m ϱ m + 1 ( γ ) = lim m u 0 + 0 γ f ( s , ϱ m 1 ( s ) ) + g ( s , ϱ m 1 ( s ) ) h ( s , ϰ m 1 ( s ) ) Δ s r = u 0 + 0 γ 1 f ( s , r ( s ) ) + g ( s , r ( s ) ) h ( s , ρ ( s ) ) Δ s
and
r Δ = f γ , r + g γ , r h γ , ρ .
Consequently, it can be obtained that
ρ Δ = f γ , ρ + g γ , ρ h γ , r , ρ 0 = u 0 , r Δ = f γ , r + g γ , r h γ , ρ , r 0 = u 0
and
ϰ 0 ϰ 1 · · · ρ u r · · · ϱ 1 ϱ 0 .
Therefore, ρ γ and r γ are coupled with third-type minimal and maximal solutions of (15). This completes the proof.

3.3. An Example

For the particular case of the system (15), consider the following dynamic initial value problem:
u Δ = u 2 4 + u 2 2 u 2 3 ,
where f γ , u = u 2 4 , g γ , u = u 2 2 , h γ , u = u 2 3 , γ 0 , 1 . Let ϰ 0 γ = 1 , ϱ 0 γ = 1 2 . Therefore, the following inequalities are satisfied:
ϰ 0 Δ = 0 f γ , ϰ 0 + g γ , ϰ 0 h γ , ϱ 0 = 11 12 ,
and
ϱ 0 Δ = 0 f γ , ϱ 0 + g γ , ϱ 0 h γ , ϰ 0 = 17 96 .
Hence, ϰ 0 γ , ϱ 0 γ are coupled with the third-type LS and US of (15). If the functions ϰ 1 γ , ϱ 1 γ are selected as ϰ 1 γ = 11 12 γ 1 and ϱ 1 γ = 17 96 γ + 1 2 , then
ϰ 0 γ ϰ 1 γ ϱ 1 γ ϱ 0 γ , γ 0 , 1
is satisfied. If we continue in this way, ϰ 2 γ , ϰ 3 γ , · · · , ϰ n γ and ϱ 2 γ , ϱ 3 γ , · · · , ϱ n γ can be obtained. Therefore, all the conditions of Theorem 5 are satisfied.

4. Conclusions

In this work, we studied a unique solution by combining techniques from the generalized QLM, using comparison results with coupled seventh-type lower and upper solutions, and improved the QLM to given nonlinear differential equations on a time scale. Under convenient conditions, it was observed that the monotone sequences converge to the unique solution of the original problem uniformly and monotonically. Furthermore, we observed that this convergence is quadratic. Additionally, we applied a monotone iterative technique with coupled third-type lower and upper solutions when the functions involved were monotonically non-decreasing. The sequences obtained converge to the extremal solutions of the problem uniformly and monotonically. Finally, we calculated a numerical solution for a dynamic IVP based on the MIT.

Author Contributions

All authors (Ş.Ç., Y.Y. and C.Y.) contributed to this article. All authors have read and agreed to the published version of the manuscript.

Funding

This publication has been supported by Sakarya University and TUBITAK (the Scientific and Technological Research Council of Türkiye).

Data Availability Statement

Data are contained within this published article.

Acknowledgments

The first author would like to thank TUBITAK (the Scientific and Technological Research Council of Türkiye) for the financial support during his PhD studies.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
QLMquasilinearization method
LSlower solution
USupper solution
IVPinitial value problem

References

  1. Hilger, S. Analysis on measure chains—A unified approach to continuous and discrete calculus. Results Math. 1990, 18, 18–56. [Google Scholar] [CrossRef]
  2. Kaymakçalan, B.; Lawrence, B. Coupled solutions and monotone iterative techniques for some nonlinear IVPs on time scales. Nonlinear Anal. Real World Appl. 2003, 4, 245–259. [Google Scholar] [CrossRef]
  3. Akın-Bohner, E.; Atici, F.M. Quasilinearization approach for two point nonlinear boundary value problems on time scales. Rocky Mt. J. Math. 2005, 35, 142–149. [Google Scholar] [CrossRef]
  4. Bhaskar, T.G.; McRae, F.A. Monotone iterative techniques for nonlinear problems involving the difference of two monotone functions. Appl. Math. Comput. 2002, 133, 187–192. [Google Scholar] [CrossRef]
  5. West, I.H.; Vatsala, A.S. Generalized monotone iterative methodfor initial value problems. Appl. Math. Lett. 2004, 17, 1231–1237. [Google Scholar] [CrossRef]
  6. Lakshmikantham, V.; Vatsala, A.S. Generalized Quasilinearization for Nonlinear Problems; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013; 440p. [Google Scholar]
  7. Denton, Z.; Vatsala, A.S. Monotone iterative technique for finite systems of nonlinear riemann- Liouville fractional differential equations. Opusc. Math. 2011, 31, 327–339. [Google Scholar] [CrossRef]
  8. Ramírez, J.D.; Vatsala, A.S. Generalized monotone iterative technique for Caputo fractional differential equation with periodic boundary condition via initial value problem. Int. J. Differ. Equ. 2012, 2012, 842813. [Google Scholar] [CrossRef]
  9. Khavanin, M. The Method of Mixed Monotony and First Order Delay Differential Equations. Pro Math. 2012, 26, 87–96. [Google Scholar]
  10. Daneev, A.V.; Sizykh, V.N. Approximate-Optimal Synthesis of Operational Control Systems for Dynamic Objects on the Basis of Quasilinearization and Sufficient Optimality Conditions. J. Comput. Syst. Sci. Int. 2022, 61, 918–934. [Google Scholar] [CrossRef]
  11. Jyoti; Singh, M. An iterative technique for a class of Dirichlet nonlinear BVPs: Troesch’s problem. Comput. Appl. Math. 2023, 42, 163. [Google Scholar] [CrossRef]
  12. Wu, K.; Shu, C.W. Geometric quasilinearization framework for analysis and design of bound-preserving schemes. SIAM Rev. 2023, 65, 1031–1073. [Google Scholar] [CrossRef]
  13. Zare, F.; Heydari, M.; Loghmani, G.B. Quasilinearization-based Legendre collocation method for solving a class of functional Volterra integral equations. Asian-Eur. J. Math. 2023, 16, 2350078. [Google Scholar] [CrossRef]
  14. Izadi, M.; Roul, P. A new approach based on shifted Vieta- Fibonacci-quasilinearization technique and its convergence analysis for nonlinear third-order Emden–Fowler equation with multi-singularity. Commun. Nonlinear Sci. Numer. Simul. 2023, 117, 106912. [Google Scholar] [CrossRef]
  15. Verma, A.; Urus, N. Well ordered monotone iterative technique for nonlinear second order four point dirichlet bvps. Math. Model. Anal. 2022, 27, 59–77. [Google Scholar] [CrossRef]
  16. İdiz, F.; Tanoğlu, G.; Aghazadeh, N. A numerical method based on Legendre wavelet and quasilinearization technique for fractional Lane-Emden type equations. Numer. Algorithms 2024, 95, 181–206. [Google Scholar] [CrossRef]
  17. Yakar, C.; Arslan, M. Quasilinearization method for causal terminal value problems involving Riemann-Liouville fractional derivatives. Electron. J. Differ. Equ. 2019, 2019, 1–11. [Google Scholar]
  18. Bohner, M.; Peterson, A. Dynamic Equations on Time Scales: An Introduction with Applications; Birkhäuser: Boston, MA, USA, 2001. [Google Scholar]
  19. Green, J.W.; Valentine, F.A. On the Arzelà-Ascoli Theorem. Math. Mag. 1961, 34, 199–202. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Çetin, Ş.; Yılmaz, Y.; Yakar, C. Qualitative Outcomes on Monotone Iterative Technique and Quasilinearization Method on Time Scale. Axioms 2024, 13, 640. https://doi.org/10.3390/axioms13090640

AMA Style

Çetin Ş, Yılmaz Y, Yakar C. Qualitative Outcomes on Monotone Iterative Technique and Quasilinearization Method on Time Scale. Axioms. 2024; 13(9):640. https://doi.org/10.3390/axioms13090640

Chicago/Turabian Style

Çetin, Şahap, Yalçın Yılmaz, and Coşkun Yakar. 2024. "Qualitative Outcomes on Monotone Iterative Technique and Quasilinearization Method on Time Scale" Axioms 13, no. 9: 640. https://doi.org/10.3390/axioms13090640

APA Style

Çetin, Ş., Yılmaz, Y., & Yakar, C. (2024). Qualitative Outcomes on Monotone Iterative Technique and Quasilinearization Method on Time Scale. Axioms, 13(9), 640. https://doi.org/10.3390/axioms13090640

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop