Next Article in Journal
Modified Proximal Algorithms for Finding Solutions of the Split Variational Inclusions
Previous Article in Journal
Properties of Spiral-Like Close-to-Convex Functions Associated with Conic Domains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Parametrization of Caputo-Type Fractional Differential Equations with Two-Point Nonlinear Boundary Conditions

Department of Mathematics, Faculty of Art and Science, Eastern Mediterranean University, Famagusta 99628, T. R. North Cyprus, 10 Mersin, Turkey
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(8), 707; https://doi.org/10.3390/math7080707
Submission received: 27 June 2019 / Revised: 2 August 2019 / Accepted: 5 August 2019 / Published: 6 August 2019
(This article belongs to the Section Mathematics and Computer Science)

Abstract

:
In this paper, we offer a new approach of investigation and approximation of solutions of Caputo-type fractional differential equations under nonlinear boundary conditions. By using an appropriate parametrization technique, the original problem with nonlinear boundary conditions is reduced to the equivalent parametrized boundary-value problem with linear restrictions. To study the transformed problem, we construct a numerical-analytic scheme which is successful in relation to different types of two-point and multipoint linear boundary and nonlinear boundary conditions. Moreover, we give sufficient conditions of the uniform convergence of the successive approximations. Also, it is indicated that these successive approximations uniformly converge to a parametrized limit function and state the relationship of this limit function and exact solution. Finally, an example is presented to illustrate the theory.

1. Introduction

In recent years, fractional differential equations subjected to different kind of boundary conditions have been studied such as periodic/anti-periodic, nonlocal, multipoint, and integral boundary conditions, see [1,2,3].
In this paper, we consider Caputo-type fractional differential equations with nonlinear boundary conditions. We apply the technique proposed in [4,5,6,7,8,9,10,11] for investigation and approximation of solutions of Caputo-type fractional differential equations with nonlinear boundary conditions. By using an appropriate parametrization technique, nonlinear boundary conditions are transformed to linear boundary conditions by using vector parameters. To study the transformed problem, we construct a numerical-analytic scheme which is successful in relation to different types of two-point and multipoint linear boundary and nonlinear boundary conditions [4,6,7,8,9,11,12].
According to the main idea of the numerical-analytic technique, certain types of successive approximations are constructed analytically. We give sufficient conditions for the uniform convergence of the successive approximations. Also, it is indicated that these successive approximations uniformly converge to a parametrized limit function and state the relationship of this limit function and exact solution. Finally, an example is discussed for illustration of the theory.

2. Background Material

In this section, some definitions of fractional calculus are presented which we use for the statement of the problem.
Definition 1.
Let f : 0 , R be a continuous function. Then
D 0 + q f ( t ) = 1 Γ ( n q ) d d t n 0 t ( t s ) n q 1 f ( s ) d s , n 1 < q < n , n N ,
is called the Riemann–Liouville fractional derivative of order q > 0 . Here Γ is the gamma function defined by Γ ( p ) = 0 e s s p 1 d s
Definition 2.
The Caputo derivative of order q for a function f : 0 , R can be written as
c D q f ( t ) = D 0 + q f ( t ) k = 0 n 1 t k k ! f ( k ) ( 0 ) , t > 0 , n 1 < q < n .
Remark 1.
Let f ( t ) C n 0 , . Then
c D q f ( t ) = 1 Γ ( n q ) 0 t f ( n ) ( s ) ( t s ) q + 1 n d s = I n q f ( n ) ( t ) , t > 0 , n 1 < q < n .

3. Statement of Fractional Differential Equation with Nonlinear Boundary Conditions and Identification of Parametrized Boundary-Value Problem

In this section, we state the Caputo-type fractional differential equation equipped with nonlinear boundary condition and we use the vector of parameters to reduce the nonlinear boundary conditions to the linear boundary condition.
Let us consider Caputo-type fractional differential equations with nonlinear boundary conditions
c D α x ( t ) = h ( t , x ( t ) ) , t 0 , T ,
A x ( 0 ) + B x ( T ) + g ( x ( 0 ) , x ( T ) ) = d , d R n ,
where c D α is the Caputo derivative of order α 0 , 1 , the functions h : 0 , T × D R , and g : D × D R n are continuous and the set D R n is closed and bounded domain. A and B are n × n matrices, det B 0 and d is a n dimensional vector.
By using appropriate parametrization technique [5], the given problem (1), (2) is reduced to certain parametrized two-point boundary conditions. To see that, we introduce the vectors of parameters
ω : = x ( 0 ) = ( ω 1 , ω 2 ω n ) T , ϕ : = x ( T ) = ϕ 1 , ϕ 2 , , ϕ n T ,
d ( ω , ϕ ) : = d g ( ω , ϕ ) .
and by using (4), the problem (1), (2) can be rewritten as follows:
c D α x ( t ) = h ( t , x ( t ) ) A x ( 0 ) + B x ( T ) = d ( ω , ϕ ) .

4. Conditions for Convergence of Successive Approximation

Some conditions are needed for studying of the successive approximation. In this study, the parametrized boundary-value problem (5) will be studied under the following conditions:
(A) 
The function h : 0 , T × R n R n satisfies the Lipschitz condition:
h ( t , u ) h ( t , v ) L u v
for all t 0 , T , u , v D , where L is a positive constant.
(B) 
Let
κ ( t ) = 2 t α Γ α + 1 1 t T α .
Then, κ ( t ) takes its maximum value at t = T 2 and
κ = T α 2 2 α 1 Γ α + 1 .
Define,
h = max t , x 0 , T × D h 1 2 + h 2 2
and a vector function δ : D × D R n is
δ ω , ϕ : = κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω ,
where I n is the n × n identity matrix and ω , ϕ D of the form (3). δ is the radius of a neighborhood C of the point ω D is defined as follows:
B ω , δ ω , ϕ : = x R n : x ω δ ω , ϕ for all ϕ D R n .
the set
D δ : = ω D : B ω , δ ω , ϕ D for all ϕ D
is nonempty.
(C) 
L κ < 1
where L is a positive constant and satisfies the inequality (6).
For studying of the solution of the parametrized boundary-value problem (5), we consider the sequence of functions x m which is defined by the iterative formula as follows:
x m ( t , ω , ϕ ) = ω + 1 Γ α 0 t t s α 1 h s , x m 1 s , ω , ϕ d s t T α 0 T T s α 1 h s , x m 1 s , ω , ϕ d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω
for t 0 , T , m = 1 , 2 , 3 where
x 0 ( t , ω , ϕ ) = x 01 , x 02 , . . , x 0 n T = ω D δ x m ( t , ω , ϕ ) = x m , 1 t , z , ϕ , x m , 2 t , z , ϕ , x m , n t , z , ϕ T
and ω , ϕ are considered to be parameters.
In addition, it is easy to see that the sequence of functions x m are satisfied linear parameterized boundary conditions (5) for all m 1 , ω D δ , ϕ D .
Now, we prove that the sequence of the functions (7) is uniformly convergent and show the relationship between this sequence of the functions and the limit function.
Theorem 1.
Assume that the parametrized boundary-value problem (5) satisfy the conditions A , B and C . Then for all fixed ϕ D and ω D δ , the following assertions are truė:
1. 
All functions of sequence (7) are continuous and satisfy the parametrized boundary conditions (5)
A x m ( 0 , ω , ϕ ) + B x m ( T , ω , ϕ ) = d ( ω , ϕ ) , m = 1 , 2 , 3
2. 
The sequence of functions (7) converges uniformly in t 0 , T as m to the limit function
x * ( t , ω , ϕ ) = lim m x m ( t , ω , ϕ ) .
3. 
The limit function x * satisfies the initial conditions
x * ( 0 , ω , ϕ ) = ω
and
A x * ( 0 , ω , ϕ ) + B x * ( T , ω , ϕ ) = d ( ω , ϕ )
.
4. 
The limit function (9) is the unique continuous solution of the integral equation
x ( t ) : = ω + 1 Γ α 0 t t s α 1 h s , x ( s ) d s t T α 0 T T s α 1 h s , x ( s ) d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω ,
or x ( t ) is the unique solution on the 0 , T of the Cauchy problem:
c D α x ( t ) = h ( t , x ( t ) ) + α Ω ω , ϕ , x ( 0 ) = ω
where
α Ω ω , ϕ = α T α 0 T T s α 1 h s , x * ( t , ω , ϕ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω .
5. 
Error estimation:
x * ( t , ω , ϕ ) x m ( t , ω , ϕ ) L κ m h κ + B 1 d ( ω , ϕ ) B 1 A + I n ω 1 1 L κ .
Proof. 
1.
Continuity of the sequence x m defined by (7) follows directly from the construction of sequence and by direct computation, it is easy to show that the sequence x m satisfies the parametrized boundary conditions (5).
2.
We prove that the sequence of functions is a Cauchy sequence in the Banach space C a , b , R n . At first, we need to show that x m ( t , ω , ϕ ) D for all ( t , ω , ϕ ) 0 , T × D δ × D , m N . We start from Equation (7). For m = 1 :
x 1 ( t , ω , ϕ ) = ω + 1 Γ α 0 t t s α 1 h s , x 0 s , ω , ϕ d s t T α 0 T T s α 1 h s , x 0 s , ω , ϕ d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω .
The Equation (13) can be written as follows:
x 1 ( t , ω , ϕ ) ω 1 Γ α 0 t t s α 1 t T α T s α 1 h s , ω d s + t T t T α T s α 1 h s , ω d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω : = I 1 + I 2 + I 3 .
We start from the estimation of I 1 :
I 1 1 Γ α 0 t t s α 1 t T α T s α 1 h d s = t T α T t α Γ α + 1 h ,
where the expression under the absolute value is nonnegative
1 t s 1 α t T α 1 t s 1 α t T α 1 T s 1 α .
Then, we estimate I 2 and I 3 :
I 2 1 Γ α t T t T α T s α 1 h s , ω d s = t T α T t α Γ α + 1 h
and
I 3 = t T α B 1 d ( ω , ϕ ) B 1 A + I n ω .
Substituting (15)–(17) into the relation (14) and we obtain the following result
x 1 ( t , ω , ϕ ) ω 2 t α Γ α + 1 1 t T α h + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω T α 2 2 α 1 Γ α + 1 h + B 1 d ( ω , ϕ ) B 1 A + I n ω = κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω = δ ω , ϕ .
Thus,
x 1 ( t , ω , ϕ ) D for ( t , ω , ϕ ) 0 , T × D δ × D .
By induction, it can be shown that all functions x m ( t , ω , ϕ ) defined by (7) also belong to the set D for all m = 1 , 2 , 3 , t 0 , T , ω D δ , ϕ D . To show that, we start with the difference between x m + 1 and x m :
x m + 1 ( t , ω , ϕ ) x m ( t , ω , ϕ ) = 1 Γ α ( 0 t t s α 1 h s , x m s , ω , ϕ ) h s , x m 1 s , ω , ϕ d s 0 T t T α T s α 1 × h s , x m s , ω , ϕ h s , x m 1 s , ω , ϕ d s
for m = 1 , 2 ,
Here, we denote the difference (19) by r m ( t , ω , ϕ ) as follows:
r m ( t , ω , ϕ ) : = x m ( t , ω , ϕ ) x m 1 ( t , ω , ϕ ) , for all m = 1 , 2 , 3 ,
We rewrite the inequality (18), by using (20) for m = 1 . Then, we obtain
r 1 ( t , ω , ϕ ) = x 1 ( t , ω , ϕ ) ω κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω .
Taking into account the Lipshitz condition A and the relation (21) for m = 2 into Equation (20), we get
r 2 ( t , ω , ϕ ) L Γ α 0 t t s α 1 t T α T s α 1 + t T t T α T s α 1 x 1 ( s , ω , ϕ ) ω d s = L Γ α 0 t t s α 1 t T α T s α 1
+ t T t T α T s α 1 d s r 1 ( s , ω , ϕ ) 2 L t α Γ α + 1 1 t T α κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω L κ 2 h + L κ B 1 d ( ω , ϕ ) B 1 A + I n ω .
Hence,
r 2 ( t , ω , ϕ ) L κ κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω .
Therefore, by using mathematical induction, we obtain the following inequality:
r m + 1 ( t , ω , ϕ ) L κ m κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω m = 0 , 1 , 2 . .
In view of (23) and by using triangular inequality we get
x m + j t , ω , ϕ x m t , ω , ϕ x m + j t , ω , ϕ x m + j 1 t , ω , ϕ + x m + j 1 t , ω , ϕ x m + j 2 t , ω , ϕ + + x m + 1 t , ω , ϕ x m t , ω , ϕ = r m + j t , ω , ϕ + r m + j 1 t , ω , ϕ + + r m + 1 t , ω , ϕ = i = 1 j r m + i t , ω , ϕ L κ m h κ + B 1 d ( ω , ϕ ) B 1 A + I n ω i = 1 L i 1 κ i 1 .
From the assumption C , it follows that
lim m 0 L κ m = 0 .
Hence, by (24), x m is Cauchy sequence and uniformly converges on 0 , T × D δ × D to a certain limit x * .
3.
Taking the limit in (8) as m , we see that x * satisfies the boundary conditions directly.
4.
By using contradiction, the uniqueness of the solution is shown. Assume that there are two limit functions such as x 1 * ( t , ω , ϕ ) and x 2 * ( t , ω , ϕ ) . Then, estimating the difference between x 1 * and x 2 *
x 1 * ( t , ω , ϕ ) x 2 * ( t , ω , ϕ ) L Γ α 0 t t s α 1 x 1 * ( s , ω , ϕ ) x 2 * ( s , ω , ϕ ) d s + 0 T T s α 1 x 1 * ( s , ω , ϕ , ψ ) x 2 * ( s , ω , ϕ ) d s L κ x 1 * x 2 *
Thus,
x 1 * x 2 * L κ x 1 * x 2 * .
It can be written
1 L κ x 1 * x 2 * 0 .
So, x 1 * x 2 * = 0 x 1 * x 2 * = 0 x 1 * = x 2 * .
5.
Passing to j in (24) we get
x 1 * ( t , ω , ϕ ) x 2 * ( t , ω , ϕ ) L κ m h κ + B 1 d ( ω , ϕ ) B 1 A + I n ω i = 1 L i 1 κ i 1 = L κ m h κ + B 1 d ( ω , ϕ ) B 1 A + I n ω 1 1 L κ .
 ☐
Remark 2.
If A = I n , B = I n , g x ( 0 ) , x ( T ) = 0 , d = 0 , boundary condition (2) becomes x ( 0 ) = x ( T ) . Please note that this problem was studied in [5].

5. Relationship between the Limit Function and the Solution of the Nonlinear Boundary-Value Problem

We consider the following equation
c D α x ( t ) = h ( t , x ) + ψ , t 0 , T
and
x 0 = ω ,
where ψ = c o l ψ 1 ψ n is the parameter of control.
Theorem 2.
Let ω D δ , ϕ D be arbitrarily defined vectors. Suppose that all conditions of Theorem 1 are satisfied. The solution x = t , ω , ϕ , ψ of the initial-value problem (25), (26) satisfies the boundary conditions (5) if and only if x = t , ω , ϕ , ψ coincides with the limit function x * = x * t , ω , ϕ , ψ of sequence (7). Moreover,
ψ = ψ ω , ϕ = α T α 0 T T s α 1 h s , x * ( t , ω , ϕ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω .
Proof. 
Sufficiency: The proof is similar to the proof of theorem in [9]. Necessity: We fix an arbitrary value ψ R n and assume that the problem
c D α x ( t ) = h ( t , x ) + ψ ¯ , t 0 , T
with initial condition x 0 = ω . The solution x ¯ = x ¯ ( t ) of the problem (28) satisfying the two-point boundary conditions (5):
A x ¯ ( 0 ) + B x ¯ ( T ) = d ( ω , ϕ ) .
Then, t 0 , T x ¯ is a solution of the integral equation
x ¯ ( t ) = ω + 1 Γ α 0 t t s α 1 h s , x ¯ ( s ) d s + t α ψ ¯ Γ α + 1 .
When t = T in (29), we get the following equation
x ( T ) = ω + 1 Γ α 0 T T s α 1 h s , x ( s ) d s + T α ψ ¯ Γ α + 1 .
Also,
x ¯ ( 0 ) = ω
and from the boundary conditions (5):
x ¯ ( T ) = B 1 d ( ω , ϕ ) A ω .
By using (30) and (31), we obtain
ψ ¯ = α T α 0 T T s α 1 h s , x ¯ ( s ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω .
Then, substituting (32) into the (29), we have
x ¯ ( t ) : = ω + 1 Γ α 0 t t s α 1 h s , x ¯ ( s ) d s t T α 0 T T s α 1 h s , x ¯ ( s ) d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω ,
Moreover, the limit function x * is a solution of the (25), (26) for ψ = ψ ω , ϕ of the form (27) and satisfies the boundary conditions (5).
x * ( t , ω , ϕ , ψ ) = ω + 1 Γ α 0 t t s α 1 h s , x * ( t , ω , ϕ , ψ ) d s + t α ψ Γ α + 1 .
Similarly, we start with the solution x * ( T , ω , ϕ , ψ ) of the integral equation:
x * ( T , ω , ϕ , ψ ) = ω + 1 Γ α 0 T T s α 1 h s , x * ( T , ω , ϕ , ψ ) d s + T α ψ Γ α + 1 .
Then, the limit function x * satisfies the following boundary conditions
A x * ( 0 , ω , ϕ , ψ ) + B x * ( T , ω , ϕ , ψ ) = d ( ω , ϕ )
with the initial condition
x * ( 0 , ω , ϕ , ψ ) = ω .
From (35), we obtain
x * ( T , ω , ϕ , ψ ) = B 1 d ( ω , ϕ ) A ω .
By using relations (34) and (36) we get
ψ ω , ϕ = α T α 0 T T s α 1 h s , x * ( s , ω , ϕ , ψ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω .
After substituting relation (37) into (33), we have
x * ( t , ω , ϕ , ψ ) : = ω + 1 Γ α 0 t t s α 1 h s , x * ( s , ω , ϕ , ψ ) d s t T α 0 T T s α 1 h s , x * ( s , ω , ϕ , ψ ) d s + t T α B 1 d ( ω , ϕ ) B 1 A + I n ω .
Taking the difference between x and x * , we get
x * ( t , ω , ϕ , ψ ) x ( t ) = 1 Γ α 0 t t s α 1 h s , x * ( s , ω , ϕ , ψ ) h s , x ( s ) d s t T α 0 T T s α 1 h s , x * s , ω , ϕ , ψ h s , x ( s ) d s .
Thus, we have the following inequalities between x * and x
x * ( t , ω , ϕ , ψ ) x ( t ) L Γ α 0 t t s α 1 x * ( s , ω , ϕ , ψ ) x ( s ) d s + 0 T T s α 1 x * ( s , ω , ϕ , ψ ) x ( s ) d s L κ x * x .
Thus,
x * x L κ x * x .
It can be written
1 L κ x * x 0 .
Therefore, we have
x * x = 0 x * x = 0 x * = x .
This means that the function x coincides with x * . Also, by using (32) and (37), we obtain ψ ω , ϕ = ψ ¯ . The theorem is proved. ☐
Theorem 3.
Assume that the conditions A , B and C are satisfied for the Caputo-type fractional differential Equation (1) with nonlinear boundary conditions (2). Then, x * · , ω * , ϕ * , ϕ * is a solution of the parametrized boundary-value problem (1), (5) if and only if ω * = ω 1 * , ω 2 * , . , ω n * and ϕ * = ϕ 1 * , ϕ 2 * , . , ϕ n * satisfy the system of determining algebraic or transcendental equations
Ω ω , ϕ = α T α 0 T T s α 1 h s , x * ( s , ω , ϕ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω = 0 ,
x * T , ω , ϕ = ϕ .
Proof. 
The result is obtained from the Theorem 2 and by observing that the differential Equation (11) coincides with (1) if and only if the couple ω * , ϕ * satisfies the equation
Ω ω * , ϕ * = 0 .
 ☐
The following assertion indicates the determining system of Equations (38), (39) shows all possible solution of the Caputo-type differential Equation (1) with nonlinear boundary conditions (2).
Remark 3.
Assume that all conditions of Theorem 1 are satisfied and there exist vectors ω D δ and ϕ D satisfying the system of determining Equations (38) and (39). Then the Caputo-type differential Equation (1) with nonlinear boundary conditions (2) have the solution x ( · ) such that
x ( 0 ) = ω , x ( T ) = ϕ .
Also, this solution has the following form
x ( t ) = x * t , ω , ϕ , t 0 , T ,
where x * is the limit function of sequence (7). Conversely, if the Caputo-type differential Equation (1) with nonlinear boundary conditions (2) has a solution x ( · ) , this solution necessarily has the form (40) and the system of determining Equations (38) and (39) is satisfied for
ω = x ( 0 ) , ϕ = x ( T ) .
Remark 4.
For some m 1 , a function Ω m : D × D R n is defined by the formula
Ω m ω , ϕ : = α T α 0 T T s α 1 h s , x m ( t , ω , ϕ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω ,
where ω and ϕ are given by (3). To study the solvability of the parametrized boundary-value problem (5), we consider the approximate determining system of algebraic equations of the form
Ω m ω , ϕ = α T α 0 T T s α 1 h s , x m ( t , ω , ϕ ) d s Γ α B 1 d ( ω , ϕ ) B 1 A + I n ω = 0 ,
x m T , ω , ϕ = ϕ ,
where x m is the vector function specified by the recurrence relation (7).

6. Example

Motivated by [5], we consider a system of Caputo-type fractional differential equations
c D α x 1 = x 2 = h 1 t , x 1 , x 2 c D α x 2 = 1 2 x 2 2 1 2 x 1 + t 8 x 2 + t 1 α 4 Γ 2 α + 2 t α + 1 + 1 16 Γ 2 + α = h 2 t , x 1 , x 2
with nonlinear boundary conditions
x 1 ( 0 ) + x 1 1 2 x 2 1 2 2 = 2 α + 1 + 1 2 α 8 Γ α + 2 1 64 , x 2 ( 0 ) + x 1 1 2 x 2 1 2 = 2 α + 1 2 α 8 Γ α + 2 1 8
The pair of the functions
x 1 * = 2 t α + 1 + 1 8 Γ α + 2 , x 2 * = t 4
are the exact solution of the Caputo-type fractional differential Equation (43) with nonlinear boundary conditions (44). Then, the nonlinear boundary conditions can be shown by the form of matrix vectors as follows:
A x 0 + B x 1 2 + g x ( 0 ) , x 1 2 = d ,
where
A = 1 0 0 1 , B = 1 0 1 1 , d = 2 α + 1 + 1 2 α 8 Γ α + 2 1 64 2 α + 1 2 α 8 Γ α + 2 1 8 , g x ( 0 ) , x 1 2 = [ x 2 1 2 ] 2 0 .
Also, det ( B ) = 1 0 .
Then, new parameters are introduced as follows:
x ( 0 ) = ω : = ω 1 ω 2 , x 1 2 = ϕ : = ϕ 1 ϕ 2
By using the parameters in (46), the nonlinear boundary condition (44) can be written in the following form:
A x ( 0 ) + B x 1 2 = d g ( ω , ϕ ) .
Thus,
d ( ω , ϕ ) = d g ( ω , ϕ ) = 2 α + 1 + 1 2 α 8 Γ α + 2 1 64 + ϕ 2 2 2 α + 1 2 α 8 Γ α + 2 1 8 .
By using (47), the nonlinear boundary conditions (44) are transformed to the linear conditions as follows:
A x ( 0 ) + B x 1 2 = d ( ω , ϕ ) .
The conditions of convergence of successive approximations ( A ) , ( B ) , and ( C ) are checked. At first, the domain D is defined as follows:
D = x 1 , x 2 : x 1 1 , x 2 3 4 , t 0 , 0.5 .
Then, the first condition A which is related with Lipschitz condition is satisfied as follows:
L = max 0 , 1 , 1 / 2 , 7 / 8
Thus,
L = 1 .
Then,
κ = 0.2143
and
h 1.6207
are obtained for α = 0.9 . The vector δ ω , ϕ is stated as follows:
δ ω , ϕ : = κ h + B 1 d ( ω , ϕ ) B 1 A + I n ω . 0.3473 + 0.0565 + 2 ϕ 2 4 6 ϕ 2 2 ( 0.111833 + ω 1 ) 0.9866 ω 1 + 5 ω 1 2
Therefore, the condition ( C ) is satisfied.
Thus, it is verified that all needed conditions are fulfilled. Hence, we can proceed with the procedure of the numerical-analytic scheme described above. Therefore, the sequence of approximate solutions are constructed. For the Caputo-type boundary-value problem (43), (48) the successive approximations (7) have the following form:
x m , 1 t , ω , ϕ : = ω 1 + 1 Γ α 0 t t s α 1 h 1 s , x m 1 , 1 s , ω , ϕ , x m 1 , 1 s , ω , ϕ d s t T α 0 T T s α 1 h 1 s , x m 1 , 1 s , ω , ϕ , x m 1 , 1 s , ω , ϕ d s + t T α 2 α + 1 + 1 2 α 8 Γ α + 2 1 64 + ϕ 2 2 2 ω 1 ,
x m , 2 t , ω , ϕ : = ω 2 + 1 Γ α 0 t t s α 1 h 2 s , x m 1 , 1 s , ω , ϕ , x m 1 , 1 s , ω , ϕ d s t T α 0 T T s α 1 h 2 s , x m 1 , 1 s , ω , ϕ , x m 1 , 1 s , ω , ϕ d s + t T α 1 8 Γ ( α + 2 ) + 7 64 + ϕ 2 2 ω 1 , where α = 0.9 .
Then, by using the program Mathematica, we get following results for α = 0.9 .
Iteration 1: We start from the approximate system of algebraic Equations (41) and (42) for m = 1 . Then, the approximate system has the following solutions:
ω 1 = ω 11 = 0.0656973365195 ,
ω 2 = ω 12 = 0.00219529679272 ,
ϕ 1 = ϕ 11 = 0.179133148137 ,
ϕ 2 = ϕ 12 = 0.239437851344 .
Substituting (50)–(53) into the equations of x 1 , 1 and x 1 , 2 . Then, we obtain x 1 , 1 t and x 1 , 2 t . Figure 1 shows the graphic of x 1 , 1 t and x 1 t . On the other hand, Figure 2 indicates the graphic of x 1 , 2 t and x 2 t .
Also, for the first iteration, the maximum deviations of the exact solution are
max t 0 , 1 x 1 * ( t ) x 11 ( t ) 0.02893
max t 0 , 1 x 2 * ( t ) x 12 ( t ) 0.01547
Similarly, we use (41) and (42) to find the unknown parameters for each iteration. Also, for each iteration the solutions of approximate systems are so close with (50)–(53). Therefore, for the next iterations, components of exact and approximate solutions are shown by figures and with their maximum errors.
Iteration 50: The graphs of the first and second components of the exact and approximate (in the fifth iteration) solutions are shown in Figure 3 and Figure 4 respectively.
The following inequalities are related to the maximum deviation of the exact solution with its fifth approximations.
max t 0 , 1 x 1 * ( t ) x 501 ( t ) 0.02096
max t 0 , 1 x 2 * ( t ) x 502 ( t ) 0.01744
Iteration 100: The graphs of the first and second components of the exact and approximate (in the hundredth iteration) solutions are shown in Figure 5 and Figure 6 respectively.
For the hundredth approximation the maximum deviations of the exact solution are
max t 0 , 1 x 1 * ( t ) x 1001 ( t ) 0.01311
max t 0 , 1 x 2 * ( t ) x 1002 ( t ) 0.01471
Iteration 150: The graphs of the first and second components of the exact and approximate (in the one hundred and fifth iteration) solutions are shown on Figure 7 and Figure 8 respectively.
The following inequalities are related with the maximum deviations of the exact solution with its one hundred and fifth approximations.
max t 0 , 1 x 1 * ( t ) x 1501 ( t ) 0.008333
max t 0 , 1 x 2 * ( t ) x 1502 ( t ) 0.01077
Iteration 200: The graphs of the first and second components of the exact and approximate (in the two hundredth iteration) solutions are shown in Figure 9 and Figure 10, respectively.
For the two hundredth iteration, the maximum deviations of the exact solution are
max t 0 , 1 x 1 * ( t ) x 2001 ( t ) 0.00487
max t 0 , 1 x 2 * ( t ) x 2002 ( t ) 0.006097
Iteration 250: The graphs of the first and second components of the exact and approximate (in the two hundredth and fifth iteration) solutions are shown in Figure 11 and Figure 12, respectively.
The following inequalities are related with the maximum deviations of the exact solution with its two hundredth and fifth approximations
max t 0 , 1 x 1 * ( t ) x 2501 ( t ) 0.0004704
max t 0 , 1 x 2 * ( t ) x 2502 ( t ) 0.0006233
Iteration 300: The graphs of the first and second components of the exact and approximate (in the three hundredth iteration) solutions are shown in Figure 13 and Figure 14, respectively.
The following inequalities are related with the maximum deviations of the exact solution with its three hundredth approximations.
max t 0 , 1 x 1 * ( t ) x 3001 ( t ) 0.00007809
max t 0 , 1 x 2 * ( t ) x 3002 ( t ) 0.00006241
Iteration 364: The graphs of the first and second components of the exact and approximate (in the three hundred and sixty-fourth iteration) solutions are shown in Figure 15 and Figure 16, respectively.
The following inequalities are related with the maximum deviations of the exact solution with its three hundred and sixty-fourth approximations.
max t 0 , 1 x 1 * ( t ) x 1001 ( t ) 1.209 × 10 6
max t 0 , 1 x 2 * ( t ) x 1002 ( t ) 5.813 × 10 6
The results of the 1st component of the exact solution and its 1st approximation are compared for some t values in Table 1 with errors.
The results of the second component of the exact solution and its 1st approximation are compared for some t values in Table 2 with errors.
The results of the 1st component of the exact solution and its last approximation are compared for some t values in Table 3 with errors.
Note that, for example, 1.563 × 10−8 means to multiply 1.563 by 0.00000001. The results of the second component of the exact solution and its last approximation are compared for some t values in Table 4 with errors.

7. Conclusions

In this article, we studied Caputo-type fractional differential equations with parametrized boundary conditions. To study of the solution of the Caputo-type fractional differential equation, successive approximations are considered. It is shown that these successive approximations are uniformly convergent and the relationship between the limit function and the exact solution of the boundary-value problem is stated.

Author Contributions

The authors have made the same contribution. All authors read and approved the final manuscript.

Funding

The research was completed under the research project BAPC-04-19-02 of the Eastern Mediterranean University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmad, B.; Ntouyas, S.K.; Tariboon, J. Nonlocal fractional-order boundary value problems with generalized Riemann-Liouville integral boundary conditions. J. Comp. Anal. Appl. 2017, 23, 1281–1296. [Google Scholar]
  2. Ahmad, B. Existence of solutions for fractional differential equations of order q∈(2, 3] with anti-periodic boundary conditions. J. Appl. Math. Comput. 2010, 34, 385–391. [Google Scholar] [CrossRef]
  3. Mahmudov, N.; Emin, S. Fractional-order boundary value problems with Katugampola fractional integral conditions. Adv. Differ. Equ. 2018, 2018, 81. [Google Scholar] [CrossRef] [Green Version]
  4. Marynets, K. On construction of the approximate solution of the special type integral boundary-value problem. Electron. J. Qual. Theory Differ. Equ. 2016, 6, 1–14. [Google Scholar] [CrossRef]
  5. Ronto, M.I.; Marynets, K.V. On the parametrization of boundary-value problems with two-point nonlinear boundary conditions. Nonlinear Oscil. 2012, 14, 379–413. [Google Scholar] [CrossRef]
  6. Ronto, M.; Marynets, K.V.; Varha, J.V. Further results on the investigation of solutions of integral boundary value problems. Tatra Mt. Math. Publ. 2015, 63, 247–267. [Google Scholar] [CrossRef] [Green Version]
  7. Ronto, A.; Ronto, M. Periodic successive approximations and interval halving. Miskolc Math. Notes 2012, 13, 459–482. [Google Scholar] [CrossRef]
  8. Ronto, M.; Samoilenko, A.M. Numerical-Analytic Methods in the Theory of Boundary-Value Problems; World Scientific: Singapore, 2000. [Google Scholar]
  9. Fečkan, M.; Marynets, K. Approximation approach to periodic BVP for fractional differential systems. Eur. Phys. J. 2017, 226, 3681–3692. [Google Scholar] [CrossRef]
  10. Fečkan, M.; Marynets, K. Approximation approach to periodic BVP for mixed fractional differential systems. J. Comput. Appl. Math. 2018, 339, 208–217. [Google Scholar] [CrossRef]
  11. Ronto, M.; Marynets, K. Parametrization for Nonlinear problems with integral boundary conditions. Electron. J. Qual. Theory Differ. Equ. 2012, 99, 1–23. [Google Scholar] [CrossRef]
  12. Batelli, F.; Fečkan, M. Hand Book of Differential Equations: Ordinary Differential Equations; Elsevier: Amsterdam, The Netherlands, 2008. [Google Scholar]
Figure 1. The first component of the exact solution and its first approximation.
Figure 1. The first component of the exact solution and its first approximation.
Mathematics 07 00707 g001
Figure 2. The second component of the exact solution and its first approximation.
Figure 2. The second component of the exact solution and its first approximation.
Mathematics 07 00707 g002
Figure 3. The first component of the exact solution and its fifth approximation.
Figure 3. The first component of the exact solution and its fifth approximation.
Mathematics 07 00707 g003
Figure 4. The second component of the exact solution and its fifth approximation.
Figure 4. The second component of the exact solution and its fifth approximation.
Mathematics 07 00707 g004
Figure 5. The first component of the exact solution and its hundredth approximation.
Figure 5. The first component of the exact solution and its hundredth approximation.
Mathematics 07 00707 g005
Figure 6. The second component of the exact solution and its hundredth approximation.
Figure 6. The second component of the exact solution and its hundredth approximation.
Mathematics 07 00707 g006
Figure 7. The first component of the exact solution and its one hundred and fifth approximation.
Figure 7. The first component of the exact solution and its one hundred and fifth approximation.
Mathematics 07 00707 g007
Figure 8. The second component of the exact solution and its one hundred and fifth approximation.
Figure 8. The second component of the exact solution and its one hundred and fifth approximation.
Mathematics 07 00707 g008
Figure 9. The first component of the exact solution and its two hundredth approximation.
Figure 9. The first component of the exact solution and its two hundredth approximation.
Mathematics 07 00707 g009
Figure 10. The second component of the exact solution and its two hundredth approximation.
Figure 10. The second component of the exact solution and its two hundredth approximation.
Mathematics 07 00707 g010
Figure 11. The first component of the exact solution and its two hundredth and fifth approximation.
Figure 11. The first component of the exact solution and its two hundredth and fifth approximation.
Mathematics 07 00707 g011
Figure 12. The second component of the exact solution and its two hundredth and fifth approximation.
Figure 12. The second component of the exact solution and its two hundredth and fifth approximation.
Mathematics 07 00707 g012
Figure 13. The first component of the exact solution and its three hundredth approximation.
Figure 13. The first component of the exact solution and its three hundredth approximation.
Mathematics 07 00707 g013
Figure 14. The second component of the exact solution and its three hundredth approximation.
Figure 14. The second component of the exact solution and its three hundredth approximation.
Mathematics 07 00707 g014
Figure 15. The first component of the exact solution and its three hundred and sixty-fourth approximation.
Figure 15. The first component of the exact solution and its three hundred and sixty-fourth approximation.
Mathematics 07 00707 g015
Figure 16. The second component of the exact solution and its three hundred and sixty-fourth approximation.
Figure 16. The second component of the exact solution and its three hundred and sixty-fourth approximation.
Mathematics 07 00707 g016
Table 1. Comparing approximated and first component of exact Solution with error.
Table 1. Comparing approximated and first component of exact Solution with error.
tExact SolutionApproximated SolutionErrorRelative Errors
0.10.063750.077040.013290.2085
0.20.06750.088380.020880.3093
0.30.073750.099730.025980.3523
0.40.08250.11110.028570.3467
0.50.093750.12240.028670.3056
0.60.10750.13380.026260.2447
0.70.12380.14510.021350.1721
0.80.14250.15640.013950.0975
0.90.16380.16780.004040.0244
10.18750.17910.0083670.0448
Table 2. Comparing Approximated and second component of exact solution with error.
Table 2. Comparing Approximated and second component of exact solution with error.
tExact SolutionApproximated SolutionErrorRelative Errors
0.10.0250.019920.0050820.2032
0.20.050.042150.0078470.1570
0.30.0750.064640.010360.1381
0.40.10.087490.012510.1251
0.50.1250.11080.014160.1136
0.60.150.13480.015180.1013
0.70.1750.15950.015460.0886
0.80.20.18510.014870.0745
0.90.2250.21170.013280.0591
10.250.23940.010560.0424
Table 3. Comparing Approximated and first component of exact solutions with error.
Table 3. Comparing Approximated and first component of exact solutions with error.
tExact SolutionApproximated SolutionErrorRelative Errors
0.10.063750.063751.563 × 10−82.418 × 10−7
0.20.06750.06756.25 × 10−89.2593 × 10−7
0.30.073750.073751.406 × 10−71.9064 × 10−6
0.40.08250.08250.000000253.0303 × 10−6
0.50.093750.093753.906 × 10−74.1664 × 10−6
0.60.10750.10755.625 × 10−75.2326 × 10−6
0.70.12380.12387.656 × 10−76.1842 × 10−6
0.80.14250.14250.0000017.0175 × 10−6
0.90.16380.16380.0000012667.7289 × 10−6
10.18750.18750.0000015638.3360 × 10−6
Table 4. Comparing Approximated and second component of exact Solutions with error.
Table 4. Comparing Approximated and second component of exact Solutions with error.
tExact SolutionApproximated SolutionErrorRelative Errors
0.10.0250.0250.0000006252.5000 × 10−5
0.20.050.050.000001252.5000 × 10−5
0.30.0750.0750.0000018752.5000 × 10−5
0.40.10.10.00000252.5000 × 10−5
0.50.1250.1250.0000031252.5000 × 10−5
0.60.150.150.000003752.5000 × 10−5
0.70.1750.1750.0000043752.5000 × 10−5
0.80.20.20.0000052.5000 × 10−5
0.90.2250.2250.0000056252.5000 × 10−5
10.250.250.000006252.5000 × 10−5

Share and Cite

MDPI and ACS Style

Mahmudov, N.I.; Emin, S.; Bawanah, S. On the Parametrization of Caputo-Type Fractional Differential Equations with Two-Point Nonlinear Boundary Conditions. Mathematics 2019, 7, 707. https://doi.org/10.3390/math7080707

AMA Style

Mahmudov NI, Emin S, Bawanah S. On the Parametrization of Caputo-Type Fractional Differential Equations with Two-Point Nonlinear Boundary Conditions. Mathematics. 2019; 7(8):707. https://doi.org/10.3390/math7080707

Chicago/Turabian Style

Mahmudov, Nazım I., Sedef Emin, and Sameer Bawanah. 2019. "On the Parametrization of Caputo-Type Fractional Differential Equations with Two-Point Nonlinear Boundary Conditions" Mathematics 7, no. 8: 707. https://doi.org/10.3390/math7080707

APA Style

Mahmudov, N. I., Emin, S., & Bawanah, S. (2019). On the Parametrization of Caputo-Type Fractional Differential Equations with Two-Point Nonlinear Boundary Conditions. Mathematics, 7(8), 707. https://doi.org/10.3390/math7080707

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