The importance of convex and non-convex functions in the study of optimization is widely established. The concept of convexity also plays a key part in the subject of inequalities due to the behavior of its definition. The principles of convexity and symmetry are inextricably linked. Because of the considerable association that has emerged between the two in recent years, we may apply what we learn from one to the other. In this study, first, Hermite-Hadamard type inequalities for LR-p-convex interval-valued functions (LR-p-convex-I-V-F) are constructed in this study. Second, for the product of p-convex various Hermite-Hadamard (HH) type integral inequalities are established. Similarly, we also obtain Hermite-Hadamard-Fejér (HH-Fejér) type integral inequality for LR-p-convex-I-V-F. Finally, for LR-p-convex-I-V-F, various discrete Schur's and Jensen's type inequalities are presented. Moreover, the results presented in this study are verified by useful nontrivial examples. Some of the results reported here for be LR-p-convex-I-V-F are generalizations of prior results for convex and harmonically convex functions, as well as p-convex functions.
Citation: Jorge E. Macías-Díaz, Muhammad Bilal Khan, Muhammad Aslam Noor, Abd Allah A. Mousa, Safar M Alghamdi. Hermite-Hadamard inequalities for generalized convex functions in interval-valued calculus[J]. AIMS Mathematics, 2022, 7(3): 4266-4292. doi: 10.3934/math.2022236
[1] | Apassara Suechoei, Parinya Sa Ngiamsunthorn . Extremal solutions of φ−Caputo fractional evolution equations involving integral kernels. AIMS Mathematics, 2021, 6(5): 4734-4757. doi: 10.3934/math.2021278 |
[2] | Adisorn Kittisopaporn, Pattrawut Chansangiam . Approximate solutions of the 2D space-time fractional diffusion equation via a gradient-descent iterative algorithm with Grünwald-Letnikov approximation. AIMS Mathematics, 2022, 7(5): 8471-8490. doi: 10.3934/math.2022472 |
[3] | Nataliya Gol’dman . Nonlinear boundary value problems for a parabolic equation with an unknown source function. AIMS Mathematics, 2019, 4(5): 1508-1522. doi: 10.3934/math.2019.5.1508 |
[4] | Feng Hu, Yingxiang Xu . Optimized Schwarz waveform relaxation for heterogeneous Cattaneo-Vernotte non-Fourier heat transfer. AIMS Mathematics, 2025, 10(3): 7370-7395. doi: 10.3934/math.2025338 |
[5] | Jundong Zhou, Yawei Chu . The complex Hessian quotient flow on compact Hermitian manifolds. AIMS Mathematics, 2022, 7(5): 7441-7461. doi: 10.3934/math.2022416 |
[6] | Murat A. Sultanov, Vladimir E. Misilov, Makhmud A. Sadybekov . Numerical method for solving the subdiffusion differential equation with nonlocal boundary conditions. AIMS Mathematics, 2024, 9(12): 36385-36404. doi: 10.3934/math.20241726 |
[7] | Quincy Stévène Nkombo, Fengquan Li, Christian Tathy . Stability properties of Radon measure-valued solutions for a class of nonlinear parabolic equations under Neumann boundary conditions. AIMS Mathematics, 2021, 6(11): 12182-12224. doi: 10.3934/math.2021707 |
[8] | Abdulghani R. Alharbi . Numerical solutions to two-dimensional fourth order parabolic thin film equations using the Parabolic Monge-Ampere method. AIMS Mathematics, 2023, 8(7): 16463-16478. doi: 10.3934/math.2023841 |
[9] | Zhanwei Gou, Jincheng Shi . Blow-up phenomena and global existence for nonlinear parabolic problems under nonlinear boundary conditions. AIMS Mathematics, 2023, 8(5): 11822-11836. doi: 10.3934/math.2023598 |
[10] | Costică Moroşanu . Modeling of the continuous casting process of steel via phase-field transition system. Fractional steps method. AIMS Mathematics, 2019, 4(3): 648-662. doi: 10.3934/math.2019.3.648 |
The importance of convex and non-convex functions in the study of optimization is widely established. The concept of convexity also plays a key part in the subject of inequalities due to the behavior of its definition. The principles of convexity and symmetry are inextricably linked. Because of the considerable association that has emerged between the two in recent years, we may apply what we learn from one to the other. In this study, first, Hermite-Hadamard type inequalities for LR-p-convex interval-valued functions (LR-p-convex-I-V-F) are constructed in this study. Second, for the product of p-convex various Hermite-Hadamard (HH) type integral inequalities are established. Similarly, we also obtain Hermite-Hadamard-Fejér (HH-Fejér) type integral inequality for LR-p-convex-I-V-F. Finally, for LR-p-convex-I-V-F, various discrete Schur's and Jensen's type inequalities are presented. Moreover, the results presented in this study are verified by useful nontrivial examples. Some of the results reported here for be LR-p-convex-I-V-F are generalizations of prior results for convex and harmonically convex functions, as well as p-convex functions.
Nonlinear difference equations have long interested both mathematics and other sciences. Since these equations play a key role in many applications such as the natural model of a discrete process, they appear in many disciplines such as population biology, optics, economics, probability theory, genetics, psychology. See e.g., [1,2,3,4,5,6] and the references therein. For the last two decades, there has been interest in studying the global attractivity, the boundedness character and the periodic nature of nonlinear difference equations. For some recent results see, for example, [7,8,9,10,11,12,13,14,15,16]. However, for the last decade, some researchers have focused on the solvability of nonlinear difference equations and their systems. For some recent results see, for example, [17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33].
In this paper, we consider the following system of second-order nonlinear difference equations
xn+1=a1+ynxn−1, yn+1=b1+xnyn−1, n∈N0, | (1.1) |
where the parameters a, b, and the initial values x−1, x0, y−1, y0 are arbitrary real numbers such that the solution {( xn,yn)}n≥−1 exists. System (1.1) can be obtained systematically as follows. First, we consider the following difference equation
xn+1=a+bxnc+dxn, ad≠bc, d≠0, n∈N0, | (1.2) |
where the parameters a, b, c, d, and the initial value x0 are arbitrary such that xn are defined. Equation (1.2) is called a first order linear fractional difference equation. Equation (1.2) is solvable by virtue of several changes of variables. The most common method is to transform Eq (1.2) into the second order linear equation by using the change of variables c+dxn=zn/zn−1. For a detailed background on Eq (1.2), see e.g, [34,35]. Also, for other equations related to Eq (1.2), see [36,37,38,39]. A different case occurs when we get b=0. This case yields the following difference equation
xn+1=ac+dxn, acd≠0, n∈N0. | (1.3) |
Equation (1.3) can also be transformed into the second order linear equation by using the change of variables xn=zn−1/zn and so is solvable. Some generalizations of Eq (1.3) can inherit its solvability property. For example, the following difference equation
xn+1=ac+dxnxn−1, acd≠0, n∈N0, | (1.4) |
where the parameters a, c, d, and the initial values x−1, x0 are arbitrary such that xn are defined, is also solvable by using the change of variables xn=zn−1/zn. Hence, the general solutions of (1.2)–(1.4) follow from the general solutions of the associated linear equations and the corresponding changes of variables. Note that both Eq (1.2) and Eq (1.3) can be reduced equations with one parameter or two parameters. If we choose xn=cdun, adc2=α, and xn=√cdvn, ac√dc=β, then they are reduced equations with one parameter in un and vn, respectively. Therefore, we can take c=d=1 under favorable conditions.
Based on the above considerations, we investigate a two-dimensional generalization that maintains the solvability characteristic of Eq (1.4). So, we get a further generalization of (1.4), that is, the system given in (1.1). System (1.1) can also be transformed into a system of third-order linear equations by using the changes of variables xn=un−1/vn, yn=vn−1/un, and so can be solved. But, we will use a more practical method introduced firstly in [40] to solve the system.
We need to the following two results in the sequel of our study.
Lemma 1.1. [41] Consider the cubic equation
P(z)=z3−αz2−βz−γ=0. | (1.5) |
Equation (1.5) has the discriminant
Δ=−α2β2−4β3+4α3γ+27γ2+18αβγ. | (1.6) |
Then, the following statements are true:
(i) If Δ<0, then the polynomial P has three distinct real zeros ρ1, ρ2, ρ3.
(ii) If Δ=0, then there are two subcases:
(a) if β=−α23 and γ=α327 , then the polynomial P has the triple root ρ=α3,
(b) if β≠−α23 or γ≠α327, then the polynomial P has the double root r and the simple root ρ.
(iii) If Δ>0, then the polynomial P has one real root p and two complex roots re±iθ,θ∈(0,π).
Theorem 1.1 (Kronecker's theorem). [42] If θ is irrational, the set of points un=nθ−[nθ] is dense in the interval (0,1).
In the above theorem, n is an integer and [nθ] is greatest integer function of nθ.
This section, which contains our main results, is examined in three subsections.
In this subsection, by using an interesting and practical method, we solve system (1.1). If we take a=0 in system (1.1), then we have xn=0 for every n≥1 and yn=b for every n≥2. If we take b=0 in system (1.1), then we have xn=a for every n≥2 and yn=0 for every n≥1. So, to enable the use of the method, we suppose ab≠0 in the sequel of our study.
We start by writing system (1.1) in the following
1x2n+1=1a+y2nx2n−1a, | (2.1) |
1x2n+2=1a+y2n+1x2na, | (2.2) |
1y2n+1=1b+x2ny2n−1b, | (2.3) |
1y2n+2=1b+x2n+1y2nb | (2.4) |
for every n≥0. By multiplying (2.1), (2.2), (2.3) and (2.4) by
1n∏k=0x2k−1n∏k=0y2k, | (2.5) |
1n∏k=0x2kn+1∏k=0y2k−1, | (2.6) |
1n∏k=0x2kn∏k=0y2k−1, | (2.7) |
1n+1∏k=0x2k−1n∏k=0y2k, | (2.8) |
we have the followings
1n+1∏k=0x2k−1n∏k=0y2k=1an∏k=0x2k−1n∏k=0y2k+1an−1∏k=0x2k−1n−1∏k=0y2k, | (2.9) |
1n+1∏k=0x2kn+1∏k=0y2k−1=1an∏k=0x2kn+1∏k=0y2k−1+1an−1∏k=0x2kn∏k=0y2k−1, | (2.10) |
1n∏k=0x2kn+1∏k=0y2k−1=1bn∏k=0x2kn∏k=0y2k−1+1bn−1∏k=0x2kn−1∏k=0y2k−1 | (2.11) |
1n+1∏k=0x2k−1n+1∏k=0y2k=1bn+1∏k=0x2k−1n∏k=0y2k+1bn∏k=0x2k−1n−1∏k=0y2k | (2.12) |
for every n≥0, respectively. In fact, the equalities (2.9)–(2.12) constitute a linear system with respect to (2.5)–(2.8). Hence, we should solve (2.9)–(2.12). By using (2.9) in (2.12), (2.12) in (2.9) and similarly by using (2.10) in (2.11), (2.11) in (2.10), we have the following statements
1n+1∏k=0x2k−1n+1∏k=0y2k=1abn∏k=0x2k−1n∏k=0y2k+2abn−1∏k=0x2k−1n−1∏k=0y2k+1abn−2∏k=0x2k−1n−2∏k=0y2k, | (2.13) |
1n+1∏k=0x2k−1n∏k=0y2k=1abn∏k=0x2k−1n−1∏k=0y2k+2abn−1∏k=0x2k−1n−2∏k=0y2k+1abn−2∏k=0x2k−1n−3∏k=0y2k, | (2.14) |
1n∏k=0x2kn+1∏k=0y2k−1=1abn−1∏k=0x2kn∏k=0y2k−1+2abn−2∏k=0x2kn−1∏k=0y2k−1+1abn−3∏k=0x2kn−2∏k=0y2k−1, | (2.15) |
and
1n+1∏k=0x2kn+1∏k=0y2k−1=1abn∏k=0x2kn∏k=0y2k−1+2abn−1∏k=0x2kn−1∏k=0y2k−1+1abn−2∏k=0x2kn−2∏k=0y2k−1 | (2.16) |
for every n≥2. Note that the equations in (2.13)–(2.16) are linear with respect to (2.5), (2.8), (2.6) and (2.7), respectively, and they can be represented by the following third-order difference equation
zn+1−1abzn−2abzn−1−1abzn−2=0, n≥2, | (2.17) |
whose characteristic equation is the following equation
P(λ)=λ3−1abλ2−2abλ−1ab=0. | (2.18) |
By Lemma 1.1, we see that there are three cases to be considered.
In this case P has three real distinct zeros denoted by ρ1, ρ2, ρ3, respectively. Hence, from (2.13) and (2.17), we can write
zn=C1ρn1+C2ρn2+C3ρn3=1n∏k=0x2k−1n∏k=0y2k | (2.19) |
from which it follows that
1x2n−1y2n=C1ρn1+C2ρn2+C3ρn3C1ρn−11+C2ρn−12+C3ρn−13, | (2.20) |
where C1, C2, C3 are arbitrary real constants given by
C1(x−1,y0)=ρ2ρ3x1y2x3y4−(ρ2+ρ3)x3y4+1(ρ1−ρ2)(ρ1−ρ3)x−1y0x1y2x3y4,C2(x−1,y0)=ρ1ρ3x1y2x3y4−(ρ1+ρ3)x3y4+1(ρ2−ρ1)(ρ2−ρ3)x−1y0x1y2x3y4,C3(x−1,y0)=ρ1ρ2x1y2x3y4−(ρ1+ρ2)x3y4+1(ρ3−ρ1)(ρ3−ρ2)x−1y0x1y2x3y4, |
for every n≥0. By using (2.20) in the first equation of system (1.1), we have
x2n+1=a(C1ρn1+C2ρn2+C3ρn3)C1(ρ1+1)ρn−11+C2(ρ2+1)ρn−12+C3(ρ3+1)ρn−13 | (2.21) |
for every n≥−1. On the other hand, the first equation of system (1.1) can be written as follows
y2n=a−x2n+1x2n+1x2n−1 | (2.22) |
for every n≥0. By using (2.21) and its backward shifted one from n to n−1 in (2.22), we have
y2n=C1(ρ1+1)ρn−21+C2(ρ2+1)ρn−22+C3(ρ3+1)ρn−23a(C1ρn1+C2ρn2+C3ρn3) | (2.23) |
for every n≥0. Now, we consider Eq (2.16), which is linear with respect to (2.7) and (2.17). Hence, we have
zn=C′1ρn1+C′2ρn2+C′3ρn3=1n∏k=0x2kn∏k=0y2k−1 | (2.24) |
from which it follows that
1x2ny2n−1=C′1ρn1+C′2ρn2+C′3ρn3C′1ρn−11+C′2ρn−12+C′3ρn−13, | (2.25) |
where C′1, C′2, C′3 are arbitrary real constants given by
C′1=C1(y−1,x0), C′2=C2(y−1,x0), C′3=C3(y−1,x0), |
for every n≥0. By using (2.25) in the second equation of system (1.1), we have
y2n+1=b(C′1ρn1+C′2ρn2+C′3ρn3)C′1(ρ1+1)ρn−11+C′2(ρ2+1)ρn−12+C′3(ρ3+1)ρn−13 | (2.26) |
for every n≥−1. On the other hand, the second equation of system (1.1) can be written as follows
x2n=b−y2n+1y2n+1y2n−1 | (2.27) |
for every n≥0. By using (2.26) and its backward shifted one from n to n−1 in (2.27), we have
x2n=C′1(ρ1+1)ρn−21+C′2(ρ2+1)ρn−22+C′3(ρ3+1)ρn−23b(C′1ρn1+C′2ρn2+C′3ρn3) | (2.28) |
for every n≥0. Consequently, in the case 4ab<−27, the representation forms of the general solution of system (1.1) are given by (2.21), (2.23), (2.26) and (2.28).
In this case P(λ) has the simple root ρ and the double root r. Moreover, since ab=−427, we have ρ=−3/4 and r=−3. Hence, from (2.13) and (2.17), we have
zn=C1ρn+rn(C2+C3n)=1n∏k=0x2k−1n∏k=0y2k | (2.29) |
from which it follows that
1x2n−1y2n=C1ρn+rn(C2+C3n)C1ρn−1+rn−1(C2+C3(n−1)), | (2.30) |
where C1, C2, C3 are arbitrary real constants given by
C1(x−1,y0)=r2x1y2x3y4−2rx3y4+1(ρ−r)2x−1y0x1y2x3y4,C2(x−1,y0)=ρ(ρ−2r)x1y2x3y4+2rx3y4−1(ρ−r)2x−1y0x1y2x3y4,C3(x−1,y0)=ρrx1y2x3y4−(ρ+r)x3y4+1(r−ρ)rx−1y0x1y2x3y4, |
for every n≥0. By substituting (2.30) in the first equation of system (1.1), we have
x2n+1=a(C1ρn+rn(C2+C3n))C1(ρ+1)ρn−1+rn−1(C2(r+1)+C3(nr+n−1)) | (2.31) |
for every n≥−1. On the other hand, by using (2.31) and its backward shifted one from n to n−1 in (2.22), we have
y2n=C1(ρ+1)ρn−2+rn−2(C2(r+1)+C3((n−1)r+n−2))a(C1ρn+rn(C2+C3n)) | (2.32) |
for every n≥0. Now, by considering Eqs (2.16) and (2.17). we have
zn=C′1ρn+rn(C′2+C′3n)=1n∏k=0x2kn∏k=0y2k−1 | (2.33) |
from which it follows that
1x2ny2n−1=C′1ρn+rn(C′2+C′3n)C′1ρn−1+rn−1(C′2+C′3(n−1)), | (2.34) |
where C′1, C′2, C′3 are arbitrary real constants given by
C′1=C1(y−1,x0), C′2=C2(y−1,x0), C′3=C3(y−1,x0), |
for every n≥0. By substituting (2.34) in the second equation of system (1.1), we have
y2n+1=b(C′1ρn+rn(C′2+C′3n))C′1(ρ+1)ρn−1+rn−1(C′2(r+1)+C′3(nr+n−1)) | (2.35) |
for every n≥−1. On the other hand, by using (2.35) and its backward shifted one from n to n−1 in (2.27), it follows that
x2n=C′1(ρ+1)ρn−2+rn−2(C′2(r+1)+C′3((n−1)r+n−2))b(C′1ρn+rn(C′2+C′3n)) | (2.36) |
for every n≥0. Consequently, in the case 4ab=−27, the representation forms of the general solution of system (1.1) are given by (2.31), (2.32), (2.35) and (2.36).
In this case P(λ) has one real root and two complex roots denoted by ρ and re±iθ,θ∈(0,π), respectively. Hence, from (2.13) and (2.17), we have
zn=C1ρn+rn(C2cosnθ+C3sinnθ)=1n∏k=0x2k−1n∏k=0y2k | (2.37) |
from which it follows that
1x2n−1y2n=C1ρn+rn(C2cosnθ+C3sinnθ)C1ρn−1+rn−1(C2cos(n−1)θ+C3sin(n−1)θ), | (2.38) |
where C1, C2, C3 are arbitrary real constants given by
C1(x−1,y0)=r2x1y2x3y4−2rcosθx3y4+1(ρ2−2ρrcosθ+r2)x−1y0x1y2x3y4,C2(x−1,y0)=ρ(ρ−2rcosθ)x1y2x3y4+2rcosθx3y4−1(ρ2−2ρrcosθ+r2)x−1y0x1y2x3y4,C3(x−1,y0)=ρr(rcos2θ−ρcosθ)x1y2x3y4+(ρ2−r2cos2θ)x3y4+rcosθ−ρrsinθ(ρ2−2ρrcosθ+r2)x−1y0x1y2x3y4, |
for every n≥0. By using (2.38) in the first equation of system (1.1), we have
x2n+1=a(C1ρn+rn(C2cosnθ+C3sinnθ))C1(ρ+1)ρn−1+rn−1(C4cosnθ+C5sinnθ), | (2.39) |
where C4=C2(r+cosθ)−C3sinθ, C5=C3(r+cosθ)+C2sinθ, for every n≥−1. On the other hand, by using (2.39) and its backward shifted one from n to n−1 in (2.22), we have
y2n=C1(ρ+1)ρn−2+rn−2(C4cos(n−1)θ+C5sin(n−1)θ)a(C1ρn+rn(C2cosnθ+C3sinnθ)) | (2.40) |
for every n≥0. Now, by considering Eqs (2.16) and (2.17), we have
zn=C′1ρn+rn(C′2cosnθ+C′3sinnθ)=1n∏k=0x2kn∏k=0y2k−1 | (2.41) |
from which it follows that
1x2ny2n−1=C′1ρn+rn(C′2cosnθ+C′3sinnθ)C′1ρn−1+rn−1(C′2cos(n−1)θ+C′3sin(n−1)θ), | (2.42) |
where C′1, C′2, C′3 are arbitrary real constants given by
C′1=C1(y−1,x0), C′2=C2(y−1,x0), C′3=C3(y−1,x0), |
for every n≥0. By using (2.38) in the second equation of system (1.1), we have
y2n+1=b(C′1ρn+rn(C′2cosnθ+C′3sinnθ))C′1(ρ+1)ρn−1+rn−1(C′4cosnθ+C′5sinnθ) | (2.43) |
where C′4=C′2(r+cosθ)−C′3sinθ, C′5=C′3(r+cosθ)+C′2sinθ, for every n≥−1. On the other hand, by using (2.43) and its backward shifted one from n to n−1 in (2.27), we have
x2n=C′1(ρ+1)ρn−2+rn−2(C′4cos(n−1)θ+C′5sin(n−1)θ)b(C′1ρn+rn(C′2cosnθ+C′3sinnθ)) | (2.44) |
for every n≥0. Consequently, in the case 4ab>−27, the representation forms of the general solution of system (1.1) are given by (2.39), (2.40), (2.43) and (2.44).
The representation forms given in the previous subsection are valid where the denominators are not zero. That is, we can obtain the set of initial values that make the solutions of the system undefined from the forms by equating their denominators to zero. This operation enables us to obtain a set of initial values that produce the well-defined solutions of system (1.1). In the following we give a theorem that helps us characterize such solutions.
Theorem 2.1. Consider system (1.1). Then, the following statements are true:
(a) If 4ab<−27, then the forbidden set of system (1.1) is given by
F={(x−1,x0,y−1,y0):α′n=0 or βn=0 or αn=0 or β′n=0}, |
where
α′n=C′1ρn1+C′2ρn2+C′3ρn3, n≥0,βn=C1(ρ1+1)ρn−11+C2(ρ2+1)ρn−12+C3(ρ3+1)ρn−13, n≥−1,αn=C1ρn1+C2ρn2+C3ρn3,n≥0,β′n=C′1(ρ1+1)ρn−11+C′2(ρ2+1)ρn−12+C′3(ρ3+1)ρn−13, n≥−1. |
(b) If 4ab=−27, then the forbidden set of system (1.1) is given by
F={(x−1,x0,y−1,y0):α′n=0 or βn=0 or αn=0 or β′n=0}, |
where
α′n=C′1ρn+rn(C′2+C′3n), n≥0,βn=C1(ρ+1)ρn−1+rn−1(C2(r+1)+C3(nr+n−1)), n≥−1,αn=C1ρn+rn(C2+C3n),n≥0,β′n=C′1(ρ+1)ρn−1+rn−1(C′2(r+1)+C′3(nr+n−1)), n≥−1, |
and ρ=−3/4, r=−3.
(c) If 4ab>−27, then the forbidden set of system (1.1) is given by
F={(x−1,x0,y−1,y0):α′n=0 or βn=0 or αn=0 or β′n=0}, |
where
α′n=C′1ρn+rn(C′2cosnθ+C′3sinnθ), n≥0βn=C1(ρ+1)ρn−1+rn−1(C4cosnθ+C5sinnθ), n≥−1αn=C1ρn+rn(C2cosnθ+C3sinnθ),n≥0β′n=C′1(ρ+1)ρn−1+rn−1(C′4cosnθ+C′5sinnθ), n≥−1. |
Proof. The proof is simple and follows by equalizing denominators of the representation forms obtained in the previous section to zero.
By considering this theorem, we say that a well-defined solution of system (1.1) is a solution {(xn,yn)}n≥−1 obtained using the initial values such that (x−1,x0,y−1,y0)∈R4∖F.
In this subsection we study the long-term behavior of the solutions of system (1.1) by using the representation forms obtained in the first subsection. We analyze the solutions in the following cases of the parameter ab:
i) Case 4ab<−27: in this case we have ab∈(−427,0).
ii) Case 4ab=−27: in this case we have ab=−427.
iii) Case 4ab>−27: in this case we have ab∈(−∞,−427)∪(0,+∞).
This case yields the following result.
Theorem 2.2. Let {(xn,yn)}n≥−1 be a well-defined solution of system (1.1). Suppose that 4ab<−27. Then, the following statements are true:
(a) If Ci≠0 for i∈{1,2,3} and |ρ|=max{|ρ1|,|ρ2|,|ρ3|}, then x2n+1→aρρ+1 and y2n→bρρ+1 as n→∞.
(b) If C′i≠0 for i∈{1,2,3} and |ρ|=max{|ρ1|,|ρ2|,|ρ3|} , then x2n→aρρ+1 and y2n+1→bρρ+1 as n→∞.
(c) If Ci=0 and CjCk≠0 for i,j,k∈{1,2,3} with i≠j≠k and |ρ|=max{|ρj|,|ρk|}, then x2n+1→aρρ+1 and y2n→bρρ+1 as n→∞.
(d) If C′i=0 and C′jC′k≠0 for i,j,k∈{1,2,3} with i≠j≠k and |ρ|=max{|ρj|,|ρk|}, then x2n→aρρ+1 and y2n+1→bρρ+1 as n→∞.
Proof. (a)–(b) Let us assume without losing generality that |ρ1|=max{|ρ1|,|ρ2|,|ρ3|}. Then, we have the following limits
limn→∞x2n=limn→∞ρn−21ρn1C′1(ρ1+1)+C′2(ρ2+1)(ρ2ρ1)n−2+C′3(ρ3+1)(ρ3ρ1)n−2b(C′1+C′2(ρ2ρ1)n+C′3(ρ3ρ1)n)=ρ1+1bρ21. |
Since ρ1 is a zero of the polynomial P, we have the relation
ρ1+1bρ21=aρ1ρ1+1 |
from (2.18). Hence, we have
limn→∞x2n=aρ1ρ1+1 |
and
limn→∞x2n+1=limn→∞ρn1ρn−11a(C1+C2(ρ2ρ1)n+C3(ρ3ρ1)n)C1(ρ1+1)+C2(ρ2+1)(ρ2ρ1)n−1+C3(ρ3+1)(ρ3ρ1)n−1=aρ1ρ1+1, |
limn→∞y2n=limn→∞ρn−21ρn1C1(ρ1+1)+C2(ρ2+1)(ρ2ρ1)n−2+C3(ρ3+1)(ρ3ρ1)n−2a(C1+C2(ρ2ρ1)n+C3(ρ3ρ1)n)=ρ1+1aρ21=bρ1ρ1+1, |
limn→∞y2n+1=limn→∞ρn1ρn−11b(C′1+C′2(ρ2ρ1)n+C′3(ρ3ρ1)n)C′1(ρ1+1)+C′2(ρ2+1)(ρ2ρ1)n−1+C′3(ρ3+1)(ρ3ρ1)n−1=bρ1ρ1+1. |
The proofs of other cases are similar and so they will be omitted.
Remark 2.1. Note that the cases |C1|+|C2|+|C3|=|Ci| and |C′1|+|C′2|+|C′3|=|Cj|, i,j∈{1,2,3} are impossible. Because, for example, if C1=C2=0, then we need to the common solution of the system
ρ2ρ3x1y2x3y4−(ρ2+ρ3)x3y4+1=0, ρ1ρ3x1y2x3y4−(ρ1+ρ3)x3y4+1=0. |
This case requires that ρ1=ρ2 which is a contradiction.
Corollary 2.1. Suppose that 4ab<−27. Then, every well-defined solution of system (1.1) has a finite limit point.
This case yields the following result.
Theorem 2.3. Let {(xn,yn)}n≥−1 be a well-defined solution of system (1.1). Suppose that 4ab=−27. Then, the following statements are true:
(a) If |C2|+|C3|≠0, then x2n+1→3a2and y2n→3b2 as n→∞.
(b) If |C′2|+|C′3|≠0, then x2n→3a2 and y2n+1→3b2 as n→∞.
(c) If C2=C3=0 and C1≠0, then x2n+1→−3a and y2n→−3b as n→∞.
(d) If C′2=C′3=0 and C′1≠0, then x2n→−3a, and y2n+1→−3b as n→∞.
Proof. (a)–(b) Since ρ=−3/4 and r=−3, we have |ρ|<|r|. So, from (2.31), (2.32), (2.35) and (2.36), we have the following limits
limn→∞x2n=limn→∞rn−2brnC′1(ρ+1)(ρr)n−2+C′2(r+1)+C′3((n−1)r+n−2)C′1(ρr)n+(C′2+C′3n)=r+1br2=−29b=3a2, |
limn→∞x2n+1=limn→∞arnrn−1C1(ρr)n+C2+C3nC1(ρ+1)(ρr)n−1+C2(r+1)+C3(nr+n−1)=arr+1=3a2, |
limn→∞y2n=limn→∞rn−2arnC1(ρ+1)(ρr)n−2+C2(r+1)+C3((n−1)r+n−2)C1(ρr)n+C2+C3n=r+1ar2=−29a=3b2, |
limn→∞y2n+1=limn→∞brnrn−1C′1(ρr)n+C′2+C′3nC′1(ρ+1)(ρr)n−1+C′2(r+1)+C′3(nr+n−1)=brr+1=3b2. |
The proofs of (c) and (d) are clear from the forms in (2.31), (2.32), (2.35) and (2.36).
Corollary 2.2. Suppose that 4ab=−27. Then, every well-defined solution of system (1.1) has a finite limit point.
For this case we first prove the following lemma.
Lemma 2.1. Suppose that 4ab>−27 and the zeros of the polynomial P(λ) are ρ and re±iθ, r>0, θ∈(0,π). Then, the following statements are true:
(a) If ab∈(0,+∞), then r<ρ
(b) If ab∈(−∞,−427), then r>|ρ|
Proof. Let 4ab>−27. Then, since ρ and re±iθ are the zeros of the polynomial P(λ), the relations
ρr2=1ab and ρ3−1ab(ρ+1)2=0 |
are satisfied. We conclude from these relations that abρ>0. This implies that if ab<0, then ρ<0 and if ab>0, then ρ>0. Also, from (2.18), we have
ρ2−1abρρ2−2abρρ−1abρ=ρ2−r2ρ2−2r2ρ−r2=0, |
which implies
r=|ρρ+1|. | (2.45) |
We must consider the following two cases:
(a) If ab∈(0,+∞), then ρ>0 and so, from (2.45), we have
r=ρρ+1<ρ. |
(b) If ab∈(−∞,−427), then we see from (2.45) that ρ<−1. So, from (2.45), we have
r=|ρρ+1|>|ρ|. |
So, the proof is completed.
Remark 2.2. Note that the equality r=|ρ| is impossible. Because, in this case we have 1ab=ρ3 which yields ab= −18∉(−∞,−427).
Theorem 2.4. Let {(xn,yn)}n≥−1 be a well-defined solution of system (1.1). Suppose that 4ab>−27. Then, the following statements are true:
(a) If ab∈(0,+∞) and C1≠0, then x2n+1→aρρ+1 and y2n→bρρ+1 as n→∞.
(b) If ab∈(0,+∞) and C′1≠0, then x2n→aρρ+1 and y2n+1→bρρ+1 as n→∞.
(c) If ab∈(−∞,−427) and |C2|+|C3|≠0, then both x2n+1 and y2n are periodic or converge to a periodic solution or dense in R.
(d) If ab∈(−∞,−427) and |C′2|+|C′3|≠0, then both x2n and y2n+1 are periodic or converge to a periodic solution or dense in R.
Proof. (a)–(b) The proof follows from the formulas given in (2.39), (2.40), (2.43), (2.44) and Lemma 2.1 by taking the limit.
(c)–(d) Since the proof is similar for the forms given in (2.39), (2.40), (2.43), (2.44), we only prove for (2.44). Suppose that C′1=0. Then, from (2.44), we have
x2n=1br2C′4cos(n−1)θ+C′5sin(n−1)θC′2cosnθ+C′3sinnθ, | (2.46) |
where C′4=C′2(r+cosθ)−C′3sinθ, C′5=C′3(r+cosθ)+C′2sinθ. Also, we see that the form given in (2.46) can be written as
x2n=1br2(rcosθ+cos2θ+(rsinθ+sin2θ)tan(nθ−γ)), |
where γ is an arbitrary constant, which corresponds to the arbitrary constants C′2, C′3, and satisfies the equality
cosγC′2=sinγC′3, −π2≤γ≤π2. |
Now, we consider the following two cases:
(i) If θ=pqπ such that p,q are co-prime integers, we have
x2n=1br2(rcosθ+cos2θ+(rsinθ+sin2θ)tan(qθ−γ))=1br2(rcosθ+cos2θ+(rsinθ+sin2θ)tan(pπ−γ))=1br2(rcosθ+cos2θ+(rsinθ+sin2θ)tan(−γ)) |
which implies x2n+q=x2n. Suppose that C′1≠0. Then, since the inequality r>|ρ| holds, from (2.44), we have
x2n=1br2C′1(ρ+1)(ρr)n−2+C′4cos(n−1)θ+C′5sin(n−1)θC′1(ρr)n+C′2cosnθ+C′3sinnθ, |
which leads to
x2n→1br2C′4cos(n−1)θ+C′5sin(n−1)θC′2cosnθ+C′3sinnθ |
for large enough values of n. Hence, the sequence (x2n)n≥0 converges to a periodic solution obtained in the case C′1=0.
(ii) If θ=tπ such that t is irrational, then we have by virtue of the Kronecker's Theorem that the set {(nt−[nt])π:n∈N0 and t is irrational} is dense in the interval (0,π). Hence, we have
tan(nθ−γ)=tan(ntπ−γ)=tan(ntπ−[nt]π−γ) |
which implies the sequence (x2n)n≥0 is dense in (−∞,+∞).
In this paper we conducted a detailed analysis on all solutions of system (1.1). To do this analysis, we obtained the representation forms of general solution of the system by using a practical method. By using these forms, we characterized the well-defined solutions of system (1.1). Finally, we studied the long-term behavior of the well-defined solutions. We can summarize our results as follows:
Consider system (1.1). Then,
(a) If ab∈(−427,0), then every well-defined solution of system (1.1) has a finite limit point.
(b) If ab∈(0,+∞) and C1C′1≠0, then every well-defined solution of system (1.1) has a finite limit point.
(c) If ab∈(−∞,−427) and (|C2|+|C3|)(|C′2|+|C′3|)≠0, then every well-defined solution of system (1.1) is periodic or converges to a periodic solution or dense in R2.
This work was supported by the National Natural Science Foundation of China (No. 71601072), Key Scientific Research Project of Higher Education Institutions in Henan Province of China (No. 20B110006) and the Fundamental Research Funds for the Universities of Henan Province. The authors thank the participants of this study for their valuable contributions.
The authors declare that there are no conflict of interest associated with this publication.
[1] |
M. Z. Sarikaya, A. Saglam, H. Yildirim, On some Hadamard-type inequalities for h-convex functions, J. Math. Inequal., 2 (2008), 335–341. https://doi.org/10.7153/jmi-02-30 doi: 10.7153/jmi-02-30
![]() |
[2] |
M. Bombardelli, S. Varosanec, Properties of h-convex functions related to the Hermite-Hadamard-Fejer inequalities, Comput. Math. Appl., 58 (2009), 1869–1877. https://doi.org/10.1016/j.camwa.2009.07.073 doi: 10.1016/j.camwa.2009.07.073
![]() |
[3] | M. A. Noor, K. I. Noor, M. U. Awan, A new Hermite-Hadamard type inequality for h-convex functions, Comput. Math. Inform., 24 (2015), 191–197. |
[4] | S. S. Dragomir, On the Hadamard's inequlality for convex functions on the co-ordinates in a rectangle from the plane, Taiwan. J. Math., 5 (2001), 775–788. |
[5] |
S. S. Dragomir, B. Mond, Integral inequalities of Hadamard type for log-convex functions, Demonstr. Math., 31 (1998), 354–364. https://doi.org/10.1515/dema-1998-0214 doi: 10.1515/dema-1998-0214
![]() |
[6] | S. S. Dragomir, Refinements of the Hermite-Hadamard integral inequality for log-convex functions, J. Inequal. Pure Appl. Math., 6 (2000), 129–134. |
[7] | S. S. Dragomir, A survey of Jensen type inequalities for log-convex functions of selfadjoint operators in Hilbert spaces, Commun. Math. Anal., 10 (2011), 82–104. |
[8] |
C. P. Niculescu, The Hermite-Hadamard inequality for log-convex functions, Nonlinear Anal-Theor., 75 (2012), 662–669. https://doi.org/10.1016/j.na.2011.08.066 doi: 10.1016/j.na.2011.08.066
![]() |
[9] |
S. S. Dragomir, New inequalities of Hermite-Hadamard type for log-convex functions, Khayyam J. Math., 3 (2017), 98–115. https://doi.org/10.22034/kjm.2017.47458 doi: 10.22034/kjm.2017.47458
![]() |
[10] |
M. A. Noor, F. Qi, M. U. Awan, Some Hermite-Hadamard type inequalities for log-h-convex functions, Analysis, 33 (2013), 367–375. https://doi.org/10.1524/anly.2013.1223 doi: 10.1524/anly.2013.1223
![]() |
[11] |
S. S. Dragomir, On the Hadamard's inequality for convex functions on the co-ordinates in a rectangle from the plane, Taiwan. J. Math., 5 (2001), 775–788. https://doi.org/10.11650/twjm/1500574995 doi: 10.11650/twjm/1500574995
![]() |
[12] |
M. Alomari, M. Darus, On the Hadamard's inequality for log-convex functions on the coordinates, J. Inequal. Appl., 2009 (2009), 283147. https://doi.org/10.1155/2009/283147 doi: 10.1155/2009/283147
![]() |
[13] |
M. Z. Sarikaya, On the Hermite-Hadamard-type inequalities for co-ordinated convex function via fractional integrals, Integr. Transf. Spec. F., 25 (2014), 134–147. https://doi.org/10.1080/10652469.2013.824436 doi: 10.1080/10652469.2013.824436
![]() |
[14] | M. E. Ozdemir, E. Set, M. Z. Sarikaya, Some new Hadamard type inequalities for co-ordinated m-convex and (α, m)-convex functions, Hacet. J. Math. Stat., 40 (2011), 219–229. |
[15] | M. Alomari, M. Darus, Co-ordinated s-convex function in the first sense with some Hadamard-type inequalities, Int. J. Contemp. Math. Sci., 3 (2008), 1557–1567. |
[16] | M. A. Latif, M. Alomari, Hadamard-type inequalities for product two convex functions on the coordinates, Int. Math. Forum., 4 (2009), 2327–2338. |
[17] |
D. Singh, B. A. Dar, D. S. Kim, Sufficiency and duality in non-smooth interval valued programming problems, J. Ind. Manag. Optim., 15 (2019), 647–665. https://doi.org/10.3934/jimo.2018063 doi: 10.3934/jimo.2018063
![]() |
[18] | I. Ahmad, D. Singh, B. A. Dar, Optimality conditions in multiobjective programming problems with interval valued objective functions, Control Cybern., 44 (2015), 19–45. |
[19] |
T. M. Costa, H. Roman-Flores, Y. Chalco-Cano, Opial-type inequalities for interval-valued functions, Fuzzy Set. Syst., 358 (2019), 48–63. https://doi.org/10.1016/j.fss.2018.04.012 doi: 10.1016/j.fss.2018.04.012
![]() |
[20] | Y. Chalco-Cano, A. Flores-Franulic, H. Roman-Flores, Ostrowski type inequalities for interval-valued functions using generalized Hukuhara derivative, Comput. Appl. Math., 31 (2012), 475–472. |
[21] | A. Flores-Franulic, Y. Chalco-Cano, H. Roman-Flores, An Ostrowski type inequality for interval-valued functions, In: IFSA World Congress and NAFIPS Annual Meeting IEEE, 35 (2013), 1459–1462. https://doi.org/10.1109/ifsa-nafips.2013.6608617 |
[22] |
H. Roman-Flores, Y. Chalco-Cano, W. A. Lodwick, Some integral inequalities for interval-valued functions, Comput. Appl. Math., 37 (2016), 1306–1318. https://doi.org/10.1007/s40314-016-0396-7 doi: 10.1007/s40314-016-0396-7
![]() |
[23] |
D. Zhao, T. An, G. Ye, W. Liu, New Jensen and Hermite-Hadamard type inequalities for h-convex interval-valued functions, J. Inequal. Appl., 2018 (2018), 302. https://doi.org/10.1186/s13660-018-1896-3 doi: 10.1186/s13660-018-1896-3
![]() |
[24] |
Y. An, G. Ye, D. Zhao, W. Liu, Hermite-Hadamard type inequalities for interval (h1, h2)-convex functions. Mathematics, 7 (2019), 436. https://doi.org/10.3390/math7050436 doi: 10.3390/math7050436
![]() |
[25] |
H. Budak, T. Tunç, M. Z. Sarikaya, Fractional Hermite-Hadamard-type inequalities for interval-valued functions, Proc. Amer. Math. Soc., 148 (2020), 705–718. https://doi.org/10.1090/proc/14741 doi: 10.1090/proc/14741
![]() |
[26] | H. Budak, T. Tunç, M. Z. Sarikaya, Fractional Hermite-Hadamard-type inequalities for interval-valued functions, Proc. Am. Math. Soc., 148 (2020), 705–718. |
[27] |
L. Chen, M. S. Saleem, M. S. Zahoor, R. Bano, Some inequalities related to interval-valued-convex functions, J. Math., 2021 (2021), 6617074. https://doi.org/10.1155/2021/6617074 doi: 10.1155/2021/6617074
![]() |
[28] |
D. Zhang, C. Guo, D. Chen, G. Wang, Jensen's inequalities for set-valued and fuzzy set-valued functions, Fuzzy Sets Syst., 404 (2021), 178–204. https://doi.org/10.1016/j.fss.2020.06.003 doi: 10.1016/j.fss.2020.06.003
![]() |
[29] |
T. M. Costa, H. Roman-Flores, Some integral inequalities for fuzzy-interval-valued functions, Inform. Sci., 420 (2017), 110–125. https://doi.org/10.1016/j.ins.2017.08.055 doi: 10.1016/j.ins.2017.08.055
![]() |
[30] |
M. B. Khan, M. A. Noor, K. I. Noor, K. S. Nisar, K. A. Ismail, A. Elfasakhany, Some inequalities for LR-(h1, h2)-convex interval-valued functions by means of pseudo order relation, Int. J. Comput. Intell. Syst., 14 (2021), 1–15. https://doi.org/10.1007/s44196-021-00032-x doi: 10.1007/s44196-021-00032-x
![]() |
[31] |
M. B. Khan, P. O. Mohammed, M. A. Noor, D. Baleanu, J. L. G. Guirao, Some new fractional estimates of inequalities for LR-p-convex interval-valued functions by means of pseudo order relation, Axioms, 10 (2021), 1–18. https://doi.org/10.3390/axioms10030175 doi: 10.3390/axioms10030175
![]() |
[32] |
M. B. Khan, M. A. Noor, H. M. Al-Bayatti, K. I. Noor, Some new inequalities for LR-log-h-convex interval-valued functions by means of pseudo order relation, Appl. Math. Inf. Sci, 15 (2021), 459–470. https://doi.org/10.18576/amis/150408 doi: 10.18576/amis/150408
![]() |
[33] |
M. Kunt, İ. İşcan, Hermite-Hadamard-Fejer type inequalities for p-convex functions, Arab J. Math. Sci., 23 (2017), 215–230. https://doi.org/10.1016/j.ajmsc.2016.11.001 doi: 10.1016/j.ajmsc.2016.11.001
![]() |
[34] | L. Fejer, Uber die Fourierreihen Ⅱ, Math. Naturwiss, Anz. Ungar. Akad. Wiss., 24 (1906), 369–390. |
[35] | M. A. Noor, K. I. Noor, M. U. Awan, S. Costache, Some integral inequalities for harmonically h-convex functions, UPB Sci. Bull. Ser. A, Appl. Math. Phys., 77 (2015), 5–16. |
[36] |
V. Lupulescu, Fractional calculus for interval-valued functions, Fuzzy Sets Syst., 265 (2015), 63–85. https://doi.org/10.1016/j.fss.2014.04.005 doi: 10.1016/j.fss.2014.04.005
![]() |
[37] | R. E, Moore, Interval analysis. Prentice Hall, Englewood Cliffs, 1966. |
[38] | J. P. Wan, p-convex functions and their properties, Pure Appl. Math., 23 (2007), 130–133. |
[39] | I. Iscan, Hermite-Hadamard type inequalities for p-convex functions, Int. J. Anal. Appl., 11 (2016), 137–145. |
[40] |
M. Kunt, İ. İşcan, Hermite-Hadamard-Fejer type inequalities for p-convex functions, Arab J. Math. Sci., 23 (2017), 215–230. https://doi.org/10.1016/j.ajmsc.2016.11.001 doi: 10.1016/j.ajmsc.2016.11.001
![]() |