Next Article in Journal
Hopf-like Bifurcation Analysis of a Fractional-Order Tumor-Lymphatic Model Involving Two Time Delays
Previous Article in Journal
On Quasi-Subordination for Bi-Univalency Involving Generalized Distribution Series
Previous Article in Special Issue
Positive Periodic Solutions for a First-Order Nonlinear Neutral Differential Equation with Impulses on Time Scales
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions

1
Applied Technology College, Soochow University, Suzhou 215325, China
2
Department of Mathematics, School of Electronics & Information Engineering, Taizhou University, Taizhou 318000, China
*
Authors to whom correspondence should be addressed.
Symmetry 2024, 16(6), 774; https://doi.org/10.3390/sym16060774
Submission received: 16 May 2024 / Revised: 14 June 2024 / Accepted: 17 June 2024 / Published: 20 June 2024

Abstract

:
The Hadamard fractional derivative and integral are important parts of fractional calculus which have been widely used in engineering, biology, neural networks, control theory, and so on. In addition, the periodic boundary conditions are an important class of symmetric two-point boundary conditions for differential equations and have wide applications. Therefore, this article considers a class of nonlinear Hadamard fractional coupling ( p 1 , p 2 ) -Laplacian systems with periodic boundary value conditions. Based on nonlinear analysis methods and the contraction mapping principle, we obtain some new and easily verifiable sufficient criteria for the existence and uniqueness of solutions to this system. Moreover, we further discuss the generalized Ulam–Hyers (GUH) stability of this problem by using some inequality techniques. Finally, three examples and simulations explain the correctness and availability of our main results.

1. Introduction

In this article, we consider the following nonlinear Hadamard fractional coupling ( p 1 , p 2 ) -Laplacian system with periodic boundary value conditions
  H D 1 + α 1 [ Φ p 1 (   H D 1 + β 1 u 1 ( t ) ) ] = f 1 ( t , u 1 ( t ) , u 2 ( t ) ) , t ( 1 , e ] ,   H D 1 + α 2 [ Φ p 2 (   H D 1 + β 2 u 2 ( t ) ) ] = f 2 ( t , u 1 ( t ) , u 2 ( t ) ) , t ( 1 , e ] , u 1 ( 1 ) = u 1 ( e ) ,   H D 1 + β 1 u 1 ( 1 ) =   H D 1 + β 1 u 1 ( e ) , u 2 ( 1 ) = u 2 ( e ) ,   H D 1 + β 2 u 2 ( 1 ) =   H D 1 + β 2 u 2 ( e ) .
where 1 < α i 2 , 1 < β i 2 , p i > 1 ,   H D 1 + * represents the Hadamard fractional derivative of * -order and Φ p i ( x ) = | x | p i 2 x are p i -Lapacian operators with inverses Φ p i 1 = Φ q i , provided that 1 p i + 1 q i = 1 , f i C ( [ 1 , e ] × R 2 , R ) , i = 1 , 2 .
Hadamard [1] first proposed Hadamard fractional calculus in 1892. Hadamard fractional calculus has some obvious differences compared to Riemann–Liouville fractional calculus. The most direct manifestation is that the kernel function k H ( t , s ) = ( log t s ) α 1 of Hadamard fractional calculus is different from the kernel function k R ( t , s ) = ( t s ) α 1 of Riemann–Liouville fractional calculus. Moreover, μ > 0 , k H ( μ t , μ s ) = k H ( t , s ) and k R ( μ t , μ s ) = μ α 1 k R ( t , s ) k R ( t , s ) , which is also different. For more details on Hadamard fractional calculus, we refer the reader to [2,3,4,5,6] and the references therein. In recent years, the study of Hadamard fractional differential equations has attracted the attention of many scholars, mainly focusing on the existence, stability and approximation of solutions (see [7,8,9,10,11,12,13,14,15,16,17,18,19]). For example, Huang et al. [9] applied a nonlinear alternative of Leray–Schauder to study the existence of solutions to a nonlinear coupled Hadamard fractional system.
To describe the turbulence problem in porous media, Leibenson [20] proposed the following p-Laplacian differential equation model in 1983
Φ p ( x ( t ) ) = f ( t , x ( t ) ) , t ( 0 , 1 ] ,
where the p-Laplacian operator is defined by Φ p ( z ) = | z | p 2 z ( p > 1 ) ; its inverse is Φ p 1 = Φ q when 1 p + 1 q = 1 . Because it can describe the basic mechanical structure of turbulence problems, many scholars have begun to focus on the dynamics of nonlinear fractional differential equations with the p-Laplacian. In recent years, many excellent results (see [21,22,23,24,25,26]) have been obtained in the study of nonlinear fractional differential equations with the p-Laplacian. For example, Zhao [21] studied a nonlinear Hadamard fractional differential equation with the p-Laplacian. He defined two different distances in a metric space to discuss the solvability, approximation and stability of this system. Compared to this article, our system (1) is a coupled system of equations, and we study the solvability and stability of the system in a Banach space. In terms of norm definition and estimation, coupled systems are more complex and more difficult to solve than the single equation. In [26], Li et al. studied the existence of at least triple positive solutions of fractional-delay differential equations with a p-Laplacian operator using the Avery–Peterson theorem.
Hyers and Ulam [27,28] raised the Ulam–Hyers (UH) stability in the 1940s. On this basis, many UH-type stability concepts have been proposed successively, such as the generalized UH stability, Ulam–Hyers–Rassias stability and generalized Ulam–Hyers–Rassias stability. Recently, the study of UH-type stability has still been very popular, especially for fractional-order differential systems. There are many papers dealing with UH-type stability of nonlinear fractional differential systems (see [21,22,25,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45], among others). However, the UH-type stability of Hadamard fractional differential coupling systems is rarely studied because the study of the former is much more difficult than that of a single differential equation. To the best of our knowledge, there are no papers focusing on nonlinear Hadamard fractional coupling Laplacian systems with periodic boundary value conditions, which are an interesting and challenging problem.
The purpose of this article is to investigate the solvability, GUH stability and simulation of solutions for system (1). Our main contributions are as follows: (i) Since there are few papers dealing with nonlinear Hadamard fractional coupling Laplacian systems with periodic boundary value conditions, we first consider system (1) to fill this gap. (ii) We obtain some novel sufficient conditions for the existence, uniqueness and GUH stability of solutions to system (1). (iii) Based on the integral equation and differential equation, we have obtained the numerical solution and simulation of system (1) using the appropriate ODE toolbox in MATLAB.
The structure of this paper is as follows. Section 2 mainly introduces some necessary information and lemmas about Hadamard fractional calculus. In Section 3, some sufficient conditions on the existence and uniqueness of a solution are obtained by the contraction mapping principle. We prove the GUH stability of Hadamard fractional coupling Laplacian system (1) by applying nonlinear analysis methods and inequality techniques in Section 4. Section 5 provides the numerical solutions and simulations for three examples, illustrating the correctness and validity of our main results. Finally, we present a simple conclusion in Section 6.

2. Preliminaries

In this section, we introduce some important concepts and lemmas of Hadamard fractional calculus.
Definition 1
([3]). For a > 0 , the left-sided Hadamard fractional integral of order α > 0 for a function f : [ a , ) R is defined by
  H J a + α f ( t ) = 1 Γ ( α ) a t log t s α 1 f ( s ) d s s ,
provided the integral exists, where Γ ( α ) = 0 t α 1 e t d t and log ( · ) = log e ( · ) .
Definition 2
([3]). For a , α > 0 and f C n [ a , ) , the α-order left-sided Hadamard fractional derivative is defined by
  H D a + α f ( t ) = 1 Γ ( n α ) t d d t n a t log t s n α 1 f ( s ) d s s ,
where n 1 < α n , n = [ α ] + 1 and [ · ] is the Gaussian function.
Lemma 1
([3]). Let a , b , α > 0 and f C n ( a , b ) L 1 ( a , b ) , then
  H J a + α ( H D a + α f ( t ) ) = f ( t ) + i = 1 n d i log t a α i ,
where d 1 , d 2 , , d n are some real constants, and n = [ α ] + 1 .
Lemma 2
([21]). Let p > 1 . The p-laplacian operator Φ p ( z ) = | z | p 2 z admits the properties as follows:
( i )
If z 0 , then Φ p ( z ) = z p 1 , and Φ p ( z ) is increasing with respect to z;
( ii )
For all z , w R , Φ p ( z w ) = Φ p ( z ) Φ p ( w ) ;
( iii )
If 1 p + 1 q = 1 , then Φ q [ Φ p ( z ) ] = Φ p [ Φ q ( z ) ] = z , for all z R ;
( iv )
For all z , w 0 , z w Φ q ( z ) = Φ q ( w ) ;
( v )
0 z Φ q 1 ( w ) 0 Φ q ( z ) w ;
( vi )
| Φ q ( z ) Φ q ( w ) | ( q 1 ) M ¯ q 2 | z w | , q 2 , 0 z , w M ¯ ;
( vii )
| Φ q ( z ) Φ q ( w ) | ( q 1 ) M ̲ q 2 | z w | , 1 < q < 2 , z , w M ̲ 0 .
Lemma 3
([46]). Let X be a Banach space and E X be closed. If T : E E is contracted, then T admits a unique fixed point x * E .
Lemma 4.
Assume that 1 < α i 2 , 1 < β i 2 , and p i > 1 are some constants, f i C ( [ 1 , e ] × R 2 , R ) , i = 1 , 2 . Then, BVP (1) is equivalent to the following integral equation
u 1 ( t ) = B 1 u ( e ) ( log t ) β 1 1 + B 1 u ( t ) , u 2 ( t ) = B 2 u ( e ) ( log t ) β 2 1 + B 2 u ( t ) .
where B i u ( t ) =   H J 1 + β i ( Φ q i [ A i u ( e ) ( log t ) α i 1 + A i u ( t ) ] ) , A i u ( t ) =   H J 1 + α i f i ( t , u 1 ( t ) , u 2 ( t ) ) , i = 1 , 2 .
We can obtain Equation (2) using the method in reference [21], so we omitted the proof of Lemma 4.

3. Existence and Uniqueness of Solutions

We take X = C ( [ 1 , e ] , R ) × C ( [ 1 , e ] , R ) . For all u = ( u 1 , u 2 ) X , define the norm u = ( u 1 , u 2 ) = max { sup t [ 1 , e ] | u 1 ( t ) | , sup t [ 1 , e ] | u 2 ( t ) | } , then ( X , · ) is a Banach space. Subsequently, we will inquire into the existence and stability of (1.1) on ( X , · ) . In addition, we need the following underlying assumptions in the whole paper.
( H 1 )
1 < α i 2 , 1 < β i 2 and 1 < p i 2 are some constants, f i C ( [ 1 , e ] × R 2 , R ) , i = 1 , 2 ;
( H 2 )
There are two constants M i > 0 ( i = 1 , 2 ) , such that for all t [ 1 , e ] , ( u 1 , u 2 ) R 2 , 0 f i ( t , u 1 ( t ) , u 2 ( t ) ) M i ;
( H 3 )
There exist some continuous functions l i j ( t ) 0 ( i , j = 1 , 2 ) , such that for all t [ 1 , e ] , u 1 , u 2 , u ¯ 1 , u ¯ 2 R , | f i ( t , u 1 , u 2 ) f i ( t , u ¯ 1 , u ¯ 2 ) | l i 1 ( t ) | u 1 u ¯ 1 | + l i 2 ( t ) | u 2 u ¯ 2 | ;
( H 4 )
max { κ 1 , κ 2 } < 1 , where κ i = 4 ( q i 1 ) Γ ( α i + β i + 1 ) M i Γ ( α i + 1 ) q i 2 j = 1 2 l i j e , l i j e = max 1 t e { l i j ( t ) } , i = 1 , 2 .
Theorem 1.
If ( H 1 ) ( H 4 ) are true, then system (1) has a unique solution ( u 1 * ( t ) , u 2 * ( t ) ) X .
Proof. 
For all u = ( u 1 ( t ) , u 2 ( t ) ) X , based on Lemma 4, we define the vector operator L : X X as
L ( u 1 ( t ) , u 2 ( t ) ) = ( L 1 ( u 1 ( t ) , u 2 ( t ) ) , L 2 ( u 1 ( t ) , u 2 ( t ) ) ) ,
where
L 1 ( u 1 ( t ) , u 2 ( t ) ) = B 1 u ( e ) ( log t ) β 1 1 + B 1 u ( t ) , L 2 ( u 1 ( t ) , u 2 ( t ) ) = B 2 u ( e ) ( log t ) β 2 1 + B 2 u ( t ) .
By ( H 2 ) , we have
0   H J 1 + α i f i ( t , u 1 ( t ) , u 2 ( t ) ) M i Γ ( α i + 1 ) , t [ 1 , e ] , i = 1 , 2 .
According to ( H 2 ) and A i u ( t ) =   H J 1 + α i [ f i ( t , u 1 ( t ) , u 2 ( t ) ) ] ( i = 1 , 2 ) , for all t [ 1 , e ] , we know that
[ A i u ( t ) ( log t ) α i 1 ] = 1 Γ ( α i ) 1 t 1 log s log t α i 1 f i ( s , u 1 ( s ) , u 2 ( s ) ) d s s = ( α i 1 ) Γ ( α i ) 1 t 1 t 1 log s log t α i 2 log s ( log t ) 2 f i ( s , u 1 ( s ) , u 2 ( s ) ) d s s > 0 ,
which implies that A i u ( t ) ( log t ) α i 1 ( i = 1 , 2 ) is monotonically increasing with respect to t [ 1 , e ] . So, we have
0 < A i u ( e ) ( log t ) α i 1 A i u ( t ) < M i Γ ( α i + 1 ) , i = 1 , 2 .
For all u = ( u 1 , u 2 ) , u ¯ = ( u ¯ 1 , u ¯ 2 ) X , t [ 1 , e ] , we derive from (vi) in Lemma 2, ( H 3 ) and (5) that
  H J 1 + β i | A i u ( t ) A i u ¯ ( t ) | =   H J 1 + β i |   H J 1 + α i [ f i ( t , u 1 ( t ) , u 2 ( t ) ) ]   H J 1 + α i [ f i ( t , u ¯ 1 ( t ) , u ¯ 2 ( t ) ) ] |   H J 1 + α i + β i | f i ( t , u 1 ( t ) , u 2 ( t ) ) f i ( t , u ¯ 1 ( t ) , u ¯ 2 ( t ) ) |   H J 1 + α i + β i [ l i 1 ( t ) | u 1 ( t ) u ¯ 1 ( t ) | + l i 2 ( t ) | u 2 ( t ) u ¯ 2 ( t ) | ] 1 Γ ( α i + β i + 1 ) j = 1 2 l i j e · u u ¯ , i = 1 , 2 ,
and
| B i u ( t ) B i u ¯ ( t ) | = |   H J 1 + β i Φ q i [ A i u ( e ) ( log t ) α i 1 + A i u ( t ) ] Φ q i [ A i u ¯ ( e ) ( log t ) α i 1 + A i u ¯ ( t ) ] |   H J 1 + β i | Φ q i [ A i u ( e ) ( log t ) α i 1 A i u ( t ) ] Φ q i [ A i u ¯ ( e ) ( log t ) α i 1 A i u ¯ ( t ) ] | ( q i 1 ) M i Γ ( α i + 1 ) q i 2   H J 1 + β i | ( A i u ( e ) A i u ¯ ( e ) ) ( log t ) α 1 1 ( A i u ( t ) A i u ¯ ( t ) ) | ( q i 1 ) M i Γ ( α i + 1 ) q i 2   H J 1 + β i | A i u ( e ) A i u ¯ ( e ) | +   H J 1 + β i | A i u ( t ) A i u ¯ ( t ) | 2 ( q i 1 ) Γ ( α i + β i + 1 ) M i Γ ( α i + 1 ) q i 2 j = 1 2 l i j e u u ¯ , i = 1 , 2 .
It follows from (4) and (9) that
| L 1 ( u 1 ( t ) , u 2 ( t ) ) L 1 ( u ¯ 1 ( t ) , u ¯ 2 ( t ) ) | = | ( B 1 u ( e ) B 1 u ¯ ( e ) ) ( log t ) β 1 1 ( B 1 u ( t ) B 1 u ¯ ( t ) ) | 4 ( q 1 1 ) Γ ( α 1 + β 1 + 1 ) M 1 Γ ( α 1 + 1 ) q 1 2 j = 1 2 l 1 j e u u ¯ = κ 1 u u ¯ ,
and
| L 2 ( u 1 ( t ) , u 2 ( t ) ) L 2 ( u ¯ 1 ( t ) , u ¯ 2 ( t ) ) | = | ( B 2 u ( e ) B 2 u ¯ ( e ) ) ( log t ) β 2 1 ( B 2 u ( t ) B 2 u ¯ ( t ) ) | 4 ( q 2 1 ) Γ ( α 2 + β 2 + 1 ) M 2 Γ ( α 2 + 1 ) q 2 2 j = 1 2 l 2 j e u u ¯ = κ 2 u u ¯ .
By (10) and (11), we obtain
L ( u 1 ( t ) , u 2 ( t ) ) L ( u ¯ 1 ( t ) , u ¯ 2 ( t ) ) max { κ 1 , κ 2 } · u u ¯ .
Due to 0 < max { κ 1 , κ 2 } < 1 , (12) means that the operator L : X X is contractive. Hence, it follows from Lemma 3 that L has a unique fixed point ( u 1 * ( t ) , u 2 * ( t ) ) X , which is the unique solution of system (1). The proof is complete. □

4. GUH Stability

In this section, we mainly discuss the GUH stability of problem (1) using direct analysis methods for u = ( u 1 , u 2 ) X and ϵ > 0 . Consider the following inequality
  H D 1 + α 1 [ Φ p 1 ( H D 1 + β 1 u 1 ( t ) ) ] f 1 ( t , u 1 ( t ) , u 2 ( t ) ) ϵ , t ( 1 , e ] ,   H D 1 + α 2 [ Φ p 2 ( H D 1 + β 2 u 2 ( t ) ) ] f 2 ( t , u 1 ( t ) , u 2 ( t ) ) ϵ , t ( 1 , e ] , u 1 ( 1 ) = u 1 ( e ) ,   H D 1 + β 1 u 1 ( 1 ) =   H D 1 + β 1 u 1 ( e ) , u 2 ( 1 ) = u 2 ( e ) ,   H D 1 + β 2 u 2 ( 1 ) =   H D 1 + β 2 u 2 ( e ) .
Definition 3.
Suppose that for all ϵ > 0 and u = ( u 1 , u 2 ) X satisfying (13), there is a unique u ¯ = ( u ¯ 1 , u ¯ 2 ) X satisfying (1) such that
u ( t ) u ¯ ( t ) ω ϵ ,
where ω > 0 is a constant. Then, problem (1) is said to be Ulam–Hyers stable.
Definition 4.
Suppose that for all ϵ > 0 and u = ( u 1 , u 2 ) X satisfying (13), there is a unique u ¯ = ( u ¯ 1 , u ¯ 2 ) X satisfying (1) such that
u ( t ) u ¯ ( t ) ϖ ( ϵ ) ,
where ϖ C ( R , R + ) with ϖ ( 0 ) = 0 . Then, problem (1) is said to be GUH stable.
Remark 1.
u = ( u 1 , u 2 ) X is a solution of inequality (13) iff there exist φ = ( φ 1 , φ 2 ) X such that
(1)
| φ 1 ( t ) | ϵ and | φ 2 ( t ) | ϵ , 1 < t e ;
(2)
  H D 1 + α 1 [ Φ p 1 (   H D 1 + β 1 u 1 ( t ) ) ] = f 1 ( t , u 1 ( t ) , u 2 ( t ) ) + φ 1 ( t ) , 1 < t e ;
(3)
  H D 1 + α 2 [ Φ p 2 (   H D 1 + β 2 u 2 ( t ) ) ] = f 2 ( t , u 1 ( t ) , u 2 ( t ) ) + φ 2 ( t ) , 1 < t e ;
(4)
u 1 ( 1 ) = u 1 ( e ) ,   H D 1 + β 1 u 1 ( 1 ) =   H D 1 + β 1 u 1 ( e ) , u 2 ( 1 ) = u 2 ( e ) ,   H D 1 + β 2 u 2 ( 1 ) =   H D 1 + β 2 u 2 ( e ) .
Theorem 2.
If ( H 1 ) ( H 4 ) are satisfied, then problem (13) is GUH stable.
Proof. 
On the basis of Lemma 4 and Remark 1, the inequality (13) is solved by
u 1 ( t ) = B 1 u φ ( e ) ( log t ) β 1 1 + B 1 u φ ( t ) , u 2 ( t ) = B 2 u φ ( e ) ( log t ) β 2 1 + B 2 u φ ( t ) .
where B i u φ ( t ) =   H J 1 + β i ( Φ q i [ A i u φ ( e ) ( log t ) α i 1 + A i u φ ( t ) ] ) , A i u φ ( t ) =   H J 1 + α i [ f i ( t , u 1 ( t ) , u 2 ( t ) ) + φ i ( t ) ] , i = 1 , 2 .
According to Theorem 1 and Lemma 4, the unique solution u * = ( u 1 * , u 2 * ) X of system (1) is
u 1 * ( t ) = B 1 u * ( e ) ( log t ) β 1 1 + B 1 u * ( t ) , u 2 * ( t ) = B 2 u * ( e ) ( log t ) β 2 1 + B 2 u * ( t ) .
where B i u * ( t ) =   H J 1 + β i ( Φ q i [ A i u * ( e ) ( log t ) α i 1 + A i u * ( t ) ] ) , A i u * ( t ) =   H J 1 + α i [ f i ( t , u 1 * ( t ) , u 2 * ( t ) ) ] , i = 1 , 2 .
Similar to (8) and (9), we have
  H J 1 + β i | A i u φ i ( t ) A i u * ( t ) | =   H J 1 + β i |   H J 1 + α i [ f i ( t , u 1 ( t ) , u 2 ( t ) ) + φ i ( t ) ]   H J 1 + α i [ f i ( t , u 1 * ( t ) , u 2 * ( t ) ) ] |   H J 1 + α i + β i [ | f i ( t , u 1 ( t ) , u 2 ( t ) ) f i ( t , u 1 * ( t ) , u 2 * ( t ) ) | + | φ i ( t ) | ]   H J 1 + α i + β i [ l i 1 ( t ) | u 1 ( t ) u 1 * ( t ) | + l i 2 ( t ) | u 2 ( t ) u 2 * ( t ) | + | φ i ( t ) | ] 1 Γ ( α i + β i + 1 ) j = 1 2 l i j e · u u ¯ + ϵ , i = 1 , 2 ,
and
| B i u φ i ( t ) B i u * ( t ) | = |   H J 1 + β i Φ q i [ A i u φ i ( e ) ( log t ) α i 1 + A i u φ i ( t ) ] Φ q i [ A i u * ( e ) ( log t ) α i 1 + A i u * ( t ) ] |   H J 1 + β i | Φ q i [ A i u φ i ( e ) ( log t ) α i 1 A i u φ i ( t ) ] Φ q i [ A i u * ( e ) ( log t ) α i 1 A i u * ( t ) ] | ( q i 1 ) M i + ϵ Γ ( α i + 1 ) q i 2   H J 1 + β i | ( A i u φ i ( e ) A i u * ( e ) ) ( log t ) α i 1 ( A i u φ i ( t ) A i u * ( t ) ) | ( q i 1 ) M i + ϵ Γ ( α i + 1 ) q i 2   H J 1 + β i | A i u φ i ( e ) A i u * ( e ) | +   H J 1 + β i | A i u φ i ( t ) A i u * ( t ) | 2 ( q i 1 ) Γ ( α i + β i + 1 ) M i + ϵ Γ ( α i + 1 ) q i 2 j = 1 2 l i j e · u u ¯ + ϵ , i = 1 , 2 .
By the same manner as (10) and (11), we apply (14)–(17) to obtain
| u i ( t ) u i * ( t ) | = | ( B i u φ i ( e ) B i u * ( e ) ) ( log t ) β i 1 ( B i u φ i ( t ) B i u * ( t ) ) | | B i u φ i ( e ) B i u * ( e ) | + | B i u φ i ( t ) B i u * ( t ) | 4 ( q i 1 ) Γ ( α i + β i + 1 ) M i + ϵ Γ ( α i + 1 ) q i 2 j = 1 2 l i j e · u u ¯ + ϵ = κ i ( ϵ ) u u * + λ i ( ϵ ) ϵ , i = 1 , 2 ,
where κ i ( ϵ ) = 4 ( q i 1 ) Γ ( α i + β i + 1 ) M i + ϵ Γ ( α i + 1 ) q i 2 j = 1 2 l i j e , λ i ( ϵ ) = 4 ( q i 1 ) Γ ( α i + β i + 1 ) M i + ϵ Γ ( α i + 1 ) q i 2 , i = 1 , 2 .
We know from (18) that
u u * max { κ 1 ( ϵ ) , κ 2 ( ϵ ) } u u * + max { λ 1 ( ϵ ) , λ 2 ( ϵ ) } ϵ .
For all ϵ > 0 ( ϵ small enough), condition ( H 4 ) ensures that 0 < max { κ 1 ( ϵ ) , κ 2 ( ϵ ) } < 1 . Then, it follows from (19) that
u u * max { λ 1 ( ϵ ) , λ 2 ( ϵ ) } ϵ 1 max { κ 1 ( ϵ ) , κ 2 ( ϵ ) } = ϖ ( ϵ ) .
Obviously, when ϵ > 0 , ϖ ( ϵ ) > 0 and ϖ ( 0 ) = 0 . Therefore, we know from (20) and Definition 4 that problem (1) is GUH stable. The proof is complete. □

5. Three Examples and Simulations

This section provides three examples to verify the validity and correctness of our main results.
Example 1.
In (1), we take α 1 = 1.2 , α 2 = 1.4 , β 1 = 1.8 , β 2 = 1.6 , p 1 = 3 2 , p 2 = 7 4 , f 1 ( t , u 1 , u 2 ) = 2 + cos ( 2 t ) 20 [ 5 π 4 + arctan ( u 1 ( t ) ) + arctan ( u 2 ( t ) ) ] , f 2 ( t , u 1 , u 2 ) = 4 + sin ( 2 t ) 30 [ 7 π 4 + arctan ( u 1 ( t ) ) + 2 arctan ( u 2 ( t ) ) ] . By a simple calculation, we have
π 80 < f 1 ( t , u 1 , u 2 ) < 27 π 80 , π 80 < f 2 ( t , u 1 , u 2 ) < 13 π 24 ,
| f 1 ( t , u 1 , u 2 ) f 1 ( t , u ¯ 1 , u ¯ 2 ) | < 2 + cos ( 2 t ) 20 [ | u 1 u ¯ 1 | + | u 2 u ¯ 2 | ] , | f 2 ( t , u 1 , u 2 ) f 2 ( t , u ¯ 1 , u ¯ 2 ) | < 4 + sin ( 2 t ) 30 | u 1 u ¯ 1 | + 4 + sin ( 2 t ) 15 | u 2 u ¯ 2 | .
Thus, the conditions ( H 1 ) ( H 3 ) are true. Consequently, M 1 = 27 π 80 , M 2 = 13 π 24 , l 11 e = l 12 e 3 20 , l 21 e 1 6 , l 22 e 1 3 . When p 1 = 3 2 , p 2 = 7 4 , then q 1 = 3 > 2 , q 2 = 7 3 > 2 , and
κ 1 = 4 ( q 1 1 ) Γ ( α 1 + β 1 + 1 ) M 1 Γ ( α 1 + 1 ) q 1 2 j = 1 2 l 1 j e 0.3849 < 1 , κ 2 = 4 ( q 2 1 ) Γ ( α 2 + β 2 + 1 ) M 2 Γ ( α 2 + 1 ) q 2 2 j = 1 2 l 2 j e 0.4936 < 1 .
Thus, ( H 4 ) holds. From Theorems 1 and 2, we claim that system (1) in Example 1 has a unique solution, which is GUH stable.
Remark 2.
In Example 1, α i , β i ( i = 1 , 2 ) are all rational numbers. α 1 = 1.2 is close to 1, β 1 = 1.8 is close to 2, and α 2 = 1.4 and β 2 = 1.6 are close to 1.5. To further verify the correctness of our results and the sensitivity of the numerical simulation to the parameters, we set α i and β i close to 2 in Example 2.
Example 2.
In (1), assume α 1 = 1.9 , α 2 = 1.95 , β 1 = 1.95 , β 2 = 1.9 , p 1 , p 2 , f 1 ( t , u 1 , u 2 ) and f 2 ( t , u 1 , u 2 ) are the same as in Example 1. Then, conditions ( H 1 ) ( H 3 ) hold. We simply perform some calculations to obtain
κ 1 = 4 ( q 1 1 ) Γ ( α 1 + β 1 + 1 ) M 1 Γ ( α 1 + 1 ) q 1 2 j = 1 2 l 1 j e 0.0725 < 1 , κ 2 = 4 ( q 2 1 ) Γ ( α 2 + β 2 + 1 ) M 2 Γ ( α 2 + 1 ) q 2 2 j = 1 2 l 2 j e 0.1337 < 1 .
Thus, ( H 4 ) also holds. It follows from Theorems 1 and 2 that system (1) in Example 2 has a unique solution, which is GUH stable.
Remark 3.
In Examples 1 and 2, p 1 and p 2 are all rational numbers. To further verify the correctness of our results and the sensitivity of the numerical simulation to the parameters, we chose some irrational numbers p 1 and p 2 satisfying 1 < p 1 , p 2 < 2 in Example 3.
Example 3.
In (1), assume p 1 = 3 , p 2 = 2.5 , α 1 , α 2 , β 1 , β 2 , f 1 ( t , u 1 , u 2 ) and f 2 ( t , u 1 , u 2 ) are the same as in Example 1. Then, conditions ( H 1 ) ( H 3 ) hold. When p 1 = 3 and p 2 = 2.5 , then q 1 = 2.366 > 2 and q 2 = 2.7208 > 2 . A simple computation gives that
κ 1 = 4 ( q 1 1 ) Γ ( α 1 + β 1 + 1 ) M 1 Γ ( α 1 + 1 ) q 1 2 j = 1 2 l 1 j e 0.2694 < 1 , κ 2 = 4 ( q 2 1 ) Γ ( α 2 + β 2 + 1 ) M 2 Γ ( α 2 + 1 ) q 2 2 j = 1 2 l 2 j e 0.7197 < 1 .
Thus, ( H 4 ) also holds. We conclude from Theorems 1 and 2 that system (1) in Example 3 has a unique solution, which is GUH stable.
To complete the numerical simulations of Examples 1–3, we present a concise algorithm below. Let ( v 1 ( t ) , v 2 ( t ) ) = (   H D 1 + β 1 u 1 ( t ) ,   H D 1 + β 2 u 2 ( t ) ) ; then, system (1) can be rewritten as
u 1 ( t ) = 1 Γ ( β 1 ) 1 e ( log e s ) β 1 1 v 1 ( s ) d s s · ( log t ) β 1 1 1 t ( log t s ) β 1 1 v 1 ( s ) d s s , u 2 ( t ) = 1 Γ ( β 2 ) 1 e ( log e s ) β 2 1 v 2 ( s ) d s s · ( log t ) β 2 1 1 t ( log t s ) β 2 1 v 2 ( s ) d s s , v 1 ( t ) = [ 1 Γ ( α 1 ) 1 t ( log t s ) α 1 1 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 Γ ( α 1 ) 1 e ( log e s ) α 1 1 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 1 1 ] q 1 1 , v 2 ( t ) = [ 1 Γ ( α 2 ) 1 t ( log t s ) α 2 1 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 Γ ( α 2 ) 1 e ( log e s ) α 2 1 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 2 1 ] q 2 1 .
Taking the derivative on both sides of system (21), we get
d u 1 ( t ) d t = β 1 1 Γ ( β 1 ) t 1 e ( log e s ) β 1 1 v 1 ( s ) d s s · ( log t ) β 1 2 1 t ( log t s ) β 1 2 v 1 ( s ) d s s , d u 2 ( t ) d t = β 2 1 Γ ( β 2 ) t 1 e ( log e s ) β 2 1 v 2 ( s ) d s s · ( log t ) β 2 2 1 t ( log t s ) β 2 2 v 2 ( s ) d s s , d v 1 ( t ) d t = ( q 1 1 ) [ 1 Γ ( α 1 ) 1 t ( log t s ) α 1 1 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 Γ ( α 1 ) 1 e ( log e s ) α 1 1 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 1 1 ] q 1 2 × α 1 1 Γ ( α 1 ) t [ 1 t ( log t s ) α 1 2 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 e ( log e s ) α 1 1 f 1 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 1 2 ] , d v 2 ( t ) d t = ( q 2 1 ) [ 1 Γ ( α 2 ) 1 t ( log t s ) α 2 1 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 Γ ( α 2 ) 1 e ( log e s ) α 2 1 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 2 1 ] q 2 2 × α 2 1 Γ ( α 2 ) t [ 1 t ( log t s ) α 2 2 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s 1 e ( log e s ) α 2 1 f 2 ( s , u 1 ( s ) , u 2 ( s ) ) d s s · ( log t ) α 2 2 ] .
We can use the appropriate ODE toolbox in MATLAB to complete numerical simulations of (22). Based on the above analysis, we employ the ODE113 toolbox in MATLAB R2017b on the three examples to present numerical simulations of their solutions and GUH stability. The simulation diagrams are shown in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6.

6. Conclusions

The Hadamard fractional derivative and integral are important parts of fractional calculus which have been widely used in engineering, biology, neural networks, control theory, and so on. Compared with a single Hadamard fractional system, Hadamard fractional coupling systems have a more complex structure and a wide range of applications. As far as we know, there are no works dealing with nonlinear Hadamard fractional coupling Laplacian systems with periodic boundary value conditions. So, we investigated system (1) to fill this gap. In this article, we consider a class of nonlinear Hadamard fractional coupling ( p 1 , p 2 ) -Laplacian systems with periodic boundary value conditions. Based on nonlinear analysis methods and contraction mapping principle, we obtain some new and easily verifiable sufficient criteria of the existence, uniqueness and GUH stability of solutions of system (1). Examples 1–3 and their simulations demonstrate the correctness and availability of our main results. Meanwhile, Figure 2, Figure 4 and Figure 6 also indicate that the solution of (1) is sensitive and dependent on parameters p i , α i and β i , i = 1 , 2 . In addition, inspired by recently published papers [9,47,48], we will investigate the Lyapunov stability of fractional differential equations, the coincidence theory of fractional differential equations, and fractional differential equations involving fractional derivative impulses in the future.

Author Contributions

X.L.: conceptualization, methodology, investigation, writing—original draft, and writing—review and editing; K.Z.: conceptualization, methodology, investigation and formal analysis; H.X.: conceptualization, methodology, investigation and project administration. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the talent introduction and research initiation fund of the Applied Technology College of Soochow University and the research start-up funds for high-level talents of Taizhou University.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data were used to support this study.

Acknowledgments

The authors sincerely thank the editors and reviewers for their help and useful suggestions to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hadamard, J. Essai sur l’étude des fonctions données par leur développment de Taylor. J. Math. Pures Appl. 1892, 8, 101–186. [Google Scholar]
  2. Ahmad, B.; Alsaedi, A.; Ntouyas, S.; Tariboon, J.J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  3. Kilbas, A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  4. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus: Models and Numerical Methods; World Scientific: Singapore, 2012. [Google Scholar]
  5. Miller, K.; Ross, B. An introduction to the Fractional Calculus and Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  6. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014. [Google Scholar]
  7. Aljoudi, S.; Ahmad, B.; Nieto, J.; Alsaedi, A. A coupled system of Hadamard type sequential fractional differential equations with coupled strip conditions. Chaos Soliton. Fract. 2016, 91, 39–46. [Google Scholar] [CrossRef]
  8. Benchohra, M.; Bouriah, S.; Graef, J. Boundary value problems for nonlinear implicit Caputo-Hadamard-type fractional differential equations with impulses. Mediterr. J. Math. 2017, 14, 206. [Google Scholar] [CrossRef]
  9. Huang, H.; Zhao, K.; Liu, X. On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Math. 2022, 7, 19221–19236. [Google Scholar] [CrossRef]
  10. Mouy, M.; Boulares, H.; Alshammari, S.; Alshammari, M.; Laskri, Y.; Mohammed, W.W. On averaging principle for Caputo-Hadamard fractional stochastic differential pantograph equation. Fractal Fract. 2023, 7, 31. [Google Scholar] [CrossRef]
  11. Ortigueira, M.; Bohannan, G. Fractional scale calculus: Hadamard vs. Liouville. Fractal Fract. 2023, 7, 296. [Google Scholar] [CrossRef]
  12. Zhao, K. Existence and UH-stability of integral boundary problem for a class of nonlinear higher-order Hadamard fractional Langevin equation via Mittag–Leffler functions. Filomat 2023, 37, 1053–1063. [Google Scholar] [CrossRef]
  13. Nyamoradi, N.; Ahmad, B. Hadamard fractional differential equations on an unbounded domain with integro-initial conditions. Qual. Theor. Dyn. Syst. 2024, 23, 183. [Google Scholar] [CrossRef]
  14. Alruwaily, Y.; Venkatachalam, K.; El-hady, E. On some impulsive fractional integro-differential equation with anti-periodic conditions. Fractal Fract. 2024, 8, 219. [Google Scholar] [CrossRef]
  15. Hammad, H.; Aydi, H.; Kattan, D. Integro-differential equations implicated with Caputo-Hadamard derivatives under nonlocal boundary constraints. Phys. Scr. 2024, 99, 025207. [Google Scholar] [CrossRef]
  16. Kassim, M. Convergence to logarithmic-type functions of solutions of fractional systems with Caputo-Hadamard and Hadamard fractional derivatives. Fract. Calc. Appl. Anal. 2024, 27, 281–318. [Google Scholar] [CrossRef]
  17. Zhang, K.; Xu, J. Solvability for a system of Hadamard-type hybrid fractional differential inclusions. Demonstr. Math. 2023, 56, 20220226. [Google Scholar] [CrossRef]
  18. Ciftci, C.; Deren, F. Analysis of p-Laplacian Hadamard fractional boundary value problems with the derivative term involved in the nonlinear term. Math. Method. Appl. Sci. 2023, 46, 8945–8955. [Google Scholar] [CrossRef]
  19. Rafeeq, A.; Thabet, S.; Mohammed, M.; Kedim, I.; Vivas-Cortez, M. On Caputo-Hadamard fractional pantograph problem of two different orders with Dirichlet boundary conditions. Alex. Eng. J. 2024, 86, 386–398. [Google Scholar] [CrossRef]
  20. Leibenson, L. General problem of the movement of a compressible uid in a porous medium. Izv. Akad. Nauk Kirg. SSR Ser. Biol. Nauk 1983, 9, 7–10. [Google Scholar]
  21. Zhao, K. Solvability, Approximation and Stability of Periodic Boundary Value Problem for a Nonlinear Hadamard Fractional Differential Equation with p-Laplacian. Axioms 2023, 12, 733. [Google Scholar] [CrossRef]
  22. Zhao, K. Study on the stability and its simulation algorithm of a nonlinear impulsive ABC-fractional coupled system with a Laplacian operator via F-contractive mapping. Adv. Contin. Discret. Models 2024, 2024, 5. [Google Scholar] [CrossRef]
  23. Rao, S.; Ahmadini, A. Multiple positive solutions for system of mixed Hadamard fractional boundary value problems with (p1,p2)-Laplacian operator. AIMS Math. 2023, 8, 14767–14791. [Google Scholar] [CrossRef]
  24. Alkhazzan, A.; Al-Sadi, W.; Wattanakejorn, V.; Khan, H.; Sitthiwirattham, T.; Etemad, S.; Rezapour, S. A new study on the existence and stability to a system of coupled higher-order nonlinear BVP of hybrid FDEs under the p-Laplacian operator. AIMS Math. 2022, 7, 14187–14207. [Google Scholar] [CrossRef]
  25. Khan, H.; Chen, W.; Sun, H. Analysis of positive solution and Hyers-Ulam stability for a class of singular fractional differential equations with p-Laplacian in Banach space. Math. Methods Appl. Sci. 2018, 41, 3430–3440. [Google Scholar] [CrossRef]
  26. Li, S.; Zhang, Z.; Jiang, W. Multiple positive solutions for four-point boundary value problem of fractional delay differential equations with p-Laplacian operator. Appl. Numer. Math. 2021, 165, 348–356. [Google Scholar] [CrossRef]
  27. Ulam, S. A Collection of Mathematical Problems. Interscience Tracts in Pure and Applied Mathmatics; Interscience: New York, NY, USA, 1906. [Google Scholar]
  28. Hyers, D. On the stability of the linear functional equation. Proc. Natl. Acad. Sci. USA 1941, 27, 2222–2240. [Google Scholar] [CrossRef] [PubMed]
  29. Ali, Z.; Zada, A.; Shah, K. On Ulam’s stability for a coupled systems of nonlinear implicit fractional differential equations. Bull. Malays. Math. Sci. Soc. 2019, 42, 2681–2699. [Google Scholar] [CrossRef]
  30. Ahmad, M.; Zada, A.; Ghaderi, M.; George, R.; Rezapour, S. On the existence and stability of a neutral stochastic fractional differential system. Fractal Fract. 2022, 6, 203. [Google Scholar] [CrossRef]
  31. Ahmadova, A.; Mahmudov, N. Ulam-Hyers stability of Caputo type fractional stochastic neutral differential equations. Stat. Probabil. Lett. 2021, 168, 108949. [Google Scholar] [CrossRef]
  32. Li, M.; Wang, J. Existence results and Ulam type stability for conformable fractional oscillating system with pure delay. Chaos Soliton. Fract. 2022, 161, 112317. [Google Scholar] [CrossRef]
  33. Zhao, K. Stability of a nonlinear Langevin system of ML-Type fractional derivative affected by time-varying delays and differential feedback control. Fractal Fract. 2022, 6, 725. [Google Scholar] [CrossRef]
  34. Rizwan, R.; Liu, F.; Zheng, Z.; Park, C.; Paokanta, S. Existence theory and Ulam’s stabilities for switched coupled system of implicit impulsive fractional order Langevin equations. Bound. Value Probl. 2023, 2023, 115. [Google Scholar] [CrossRef]
  35. Priyadharsini, J.; Balasubramaniam, P. Hyers-Ulam stability result for hilfer fractional integrodifferential stochastic equations with fractional noises and non-instantaneous impulses. Evol. Equ. Control. Theory 2024, 13, 173–193. [Google Scholar] [CrossRef]
  36. Thabet, S.; Vivas-Cortez, M.; Kedim, I.; Samei, M.E.; Ayari, M.I. Solvability of a ϱ-Hilfer fractional snap dynamic system on unbounded domains. Fractal Fract. 2023, 7, 607. [Google Scholar] [CrossRef]
  37. Shah, S.; Rizwan, R.; Xia, Y.; Zada, A. Existence, uniqueness, and stability analysis of fractional Langevin equations with anti-periodic boundary conditions. Math. Methods Appl. Sci. 2023, 46, 17941–17961. [Google Scholar] [CrossRef]
  38. Zhao, K. Solvability and GUH-stability of a nonlinear CF-fractional coupled Laplacian equations. AIMS Math. 2023, 8, 13351–13367. [Google Scholar] [CrossRef]
  39. Alam, M.; Khalid, K. Analysis of q-fractional coupled implicit systems involving the nonlocal Riemann–Liouville and Erdelyi-Kober q-fractional integral conditions. Math. Methods Appl. Sci. 2023, 46, 12711–12734. [Google Scholar] [CrossRef]
  40. Khan, T.; Markarian, C.; Fachkha, C. Stability analysis of new generalized mean-square stochastic fractional differential equations and their applications in technology. AIMS Math. 2023, 8, 27840–27856. [Google Scholar] [CrossRef]
  41. Kiskinov, H.; Milev, M.; Cholakov, S.; Zahariev, A. Fundamental matrix, integral representation and stability analysis of the solutions of neutral fractional systems with derivatives in the Riemann–Liouville sense. Fractal Fract. 2024, 8, 195. [Google Scholar] [CrossRef]
  42. Bensassa, K.; Dahmani, Z.; Rakah, M.; Sarikaya, M.Z. Beam deflection coupled systems of fractional differential equations: Existence of solutions, Ulam-Hyers stability and travelling waves. Anal. Math. Phys. 2024, 14, 29. [Google Scholar] [CrossRef]
  43. Zhao, K. Generalized UH-stability of a nonlinear fractional coupling (p1,p2)-Laplacian system concerned with nonsingular Atangana-Baleanu fractional calculus. J. Inequal. Appl. 2023, 2023, 96. [Google Scholar] [CrossRef]
  44. Shah, K.; Sarwar, M.; Abdeljawad, T.; Shafiullah. Sufficient criteria for the existence of solution to nonlinear fractal-fractional order coupled system with coupled integral boundary conditions. J. Appl. Math. Comput. 2024, 70, 1771–1785. [Google Scholar] [CrossRef]
  45. Zhang, L.; Liu, X.; Jia, M.; Yu, Z. Piecewise conformable fractional impulsive differential system with delay: Existence, uniqueness and Ulam stability. J. Appl. Math. Comput. 2024, 70, 1543–1570. [Google Scholar] [CrossRef]
  46. Guo, D.; Lakshmikantham, V. Nonlinear Problems in Abstract Cone; Academic Press: Orlando, FL, USA, 1988. [Google Scholar]
  47. Zhao, K.; Liu, J.; Lv, X. A unified approach to solvability and stability of multipoint BVPs for Langevin and Sturm–Liouville equations with CH–fractional derivatives and impulses via coincidence theory. Fractal Fract. 2024, 8, 111. [Google Scholar] [CrossRef]
  48. Zhao, K. Global stability of a novel nonlinear diffusion online game addiction model with unsustainable control. AIMS Math. 2022, 7, 20752–20766. [Google Scholar] [CrossRef]
Figure 1. Numerical simulation of solutions in Example 1.
Figure 1. Numerical simulation of solutions in Example 1.
Symmetry 16 00774 g001
Figure 2. Numerical simulation of GUH stability with ϵ = 0.01 in Example 1.
Figure 2. Numerical simulation of GUH stability with ϵ = 0.01 in Example 1.
Symmetry 16 00774 g002
Figure 3. Numerical simulation of solutions in Example 2.
Figure 3. Numerical simulation of solutions in Example 2.
Symmetry 16 00774 g003
Figure 4. Numerical simulation of GUH stability with ϵ = 0.01 in Example 2.
Figure 4. Numerical simulation of GUH stability with ϵ = 0.01 in Example 2.
Symmetry 16 00774 g004
Figure 5. Numerical simulation of solutions in Example 3.
Figure 5. Numerical simulation of solutions in Example 3.
Symmetry 16 00774 g005
Figure 6. Numerical simulation of GUH stability with ϵ = 0.01 in Example 3.
Figure 6. Numerical simulation of GUH stability with ϵ = 0.01 in Example 3.
Symmetry 16 00774 g006
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

Lv, X.; Zhao, K.; Xie, H. Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions. Symmetry 2024, 16, 774. https://doi.org/10.3390/sym16060774

AMA Style

Lv X, Zhao K, Xie H. Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions. Symmetry. 2024; 16(6):774. https://doi.org/10.3390/sym16060774

Chicago/Turabian Style

Lv, Xiaojun, Kaihong Zhao, and Haiping Xie. 2024. "Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions" Symmetry 16, no. 6: 774. https://doi.org/10.3390/sym16060774

APA Style

Lv, X., Zhao, K., & Xie, H. (2024). Stability and Numerical Simulation of a Nonlinear Hadamard Fractional Coupling Laplacian System with Symmetric Periodic Boundary Conditions. Symmetry, 16(6), 774. https://doi.org/10.3390/sym16060774

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