Research article

Fractional calculus in beam deflection: Analyzing nonlinear systems with Caputo and conformable derivatives

  • In this paper, our study is divided into two parts. The first part involves analyzing a coupled system of beam deflection type that involves nonlinear equations with sequential Caputo derivatives. The also system incorporates the Caputo derivatives in the initial conditions, which adds a layer of complexity and realism to the problem. We focus on proving the existence of a unique solution for this system, and highlighting the robustness and applicability of fractional derivatives in modeling complex physical phenomena. In the second part of the paper, we employ conformable fractional derivatives, as defined by Khalil, to examine another system consisting of two coupled evolution equations. By the Tanh method, we derive new progressive waves. The connection between these two parts lies in the use of fractional calculus to extend and enhance classical problems.

    Citation: Abdelkader Lamamri, Iqbal Jebril, Zoubir Dahmani, Ahmed Anber, Mahdi Rakah, Shawkat Alkhazaleh. Fractional calculus in beam deflection: Analyzing nonlinear systems with Caputo and conformable derivatives[J]. AIMS Mathematics, 2024, 9(8): 21609-21627. doi: 10.3934/math.20241050

    Related Papers:

    [1] Kamsing Nonlaopon, Muhammad Uzair Awan, Sadia Talib, Hüseyin Budak . Parametric generalized (p,q)-integral inequalities and applications. AIMS Mathematics, 2022, 7(7): 12437-12457. doi: 10.3934/math.2022690
    [2] Syed Ghoos Ali Shah, Shahbaz Khan, Saqib Hussain, Maslina Darus . q-Noor integral operator associated with starlike functions and q-conic domains. AIMS Mathematics, 2022, 7(6): 10842-10859. doi: 10.3934/math.2022606
    [3] Aisha M. Alqahtani, Rashid Murtaza, Saba Akmal, Adnan, Ilyas Khan . Generalized q-convex functions characterized by q-calculus. AIMS Mathematics, 2023, 8(4): 9385-9399. doi: 10.3934/math.2023472
    [4] Shuhai Li, Lina Ma, Huo Tang . Meromorphic harmonic univalent functions related with generalized (p, q)-post quantum calculus operators. AIMS Mathematics, 2021, 6(1): 223-234. doi: 10.3934/math.2021015
    [5] Xhevat Z. Krasniqi . Approximation of functions in a certain Banach space by some generalized singular integrals. AIMS Mathematics, 2024, 9(2): 3386-3398. doi: 10.3934/math.2024166
    [6] Bo Wang, Rekha Srivastava, Jin-Lin Liu . Certain properties of multivalent analytic functions defined by q-difference operator involving the Janowski function. AIMS Mathematics, 2021, 6(8): 8497-8508. doi: 10.3934/math.2021493
    [7] Maimoona Karim, Aliya Fahmi, Shahid Qaisar, Zafar Ullah, Ather Qayyum . New developments in fractional integral inequalities via convexity with applications. AIMS Mathematics, 2023, 8(7): 15950-15968. doi: 10.3934/math.2023814
    [8] Li Xu, Lu Chen, Ti-Ren Huang . Monotonicity, convexity and inequalities involving zero-balanced Gaussian hypergeometric function. AIMS Mathematics, 2022, 7(7): 12471-12482. doi: 10.3934/math.2022692
    [9] Ebrahim Amini, Mojtaba Fardi, Shrideh Al-Omari, Rania Saadeh . Certain differential subordination results for univalent functions associated with q-Salagean operators. AIMS Mathematics, 2023, 8(7): 15892-15906. doi: 10.3934/math.2023811
    [10] Ye-Cong Han, Chuan-Yu Cai, Ti-Ren Huang . Monotonicity, convexity properties and inequalities involving Gaussian hypergeometric functions with applications. AIMS Mathematics, 2022, 7(4): 4974-4991. doi: 10.3934/math.2022277
  • In this paper, our study is divided into two parts. The first part involves analyzing a coupled system of beam deflection type that involves nonlinear equations with sequential Caputo derivatives. The also system incorporates the Caputo derivatives in the initial conditions, which adds a layer of complexity and realism to the problem. We focus on proving the existence of a unique solution for this system, and highlighting the robustness and applicability of fractional derivatives in modeling complex physical phenomena. In the second part of the paper, we employ conformable fractional derivatives, as defined by Khalil, to examine another system consisting of two coupled evolution equations. By the Tanh method, we derive new progressive waves. The connection between these two parts lies in the use of fractional calculus to extend and enhance classical problems.



    In this paper, we primarily focus on investigating different properties of generalized (p,q)-elliptic integrals and the generalized (p,q)-Hersch-Pfluger distortion function. In recent years, mathematicians have made significant progress in studying inequalities and various properties related to complete elliptic integrals, especially Legendre elliptic integrals and generalized complete elliptic integrals of the first and second types [3,4,5,10,11,23,26].

    We first introduce some necessary notation. For complex numbers a,b,c with c0,1,2,, and x(1,1), the Gaussian hypergeometric function [9] is defined as follows:

    F(a,b;c;x)=2F1(a,b;c;x)=n=0(a,n)(b,n)(c,n)xnn!, (1.1)

    where (a,n)a(a+1)(a+n1) is the shifted factorial function for nN+, and (a,0)=1 for a0. As we all know, 2F1(a,b;c;x) has many important applications in the theory of geometric functions and several other contexts [19]. Many special functions in mathematical physics are special or limit cases of this function [1]. Bhayo studied a new form of the generalized (p,q)-complete elliptic integrals as an application of generalized (p,q)-trigonometric functions [6]. In recent years, the generalization of classical trigonometric functions has attracted significant interest [8,20]. For this, we need the generalized arcsine function arcsinp,q(x) and the generalized πp,q. For p,q(1,), set

    arcsinp,q(x)x0dt(1tq)1/p,x[0,1],

    and the generalized πp,q is the number defined by

    πp,q=2arcsinp,q(1)210dt(1tq)1/p=2qB(11p,1q),

    where B is the beta function. For Rex>0 and Rey>0, the classical gamma function Γ(x) and beta function B(x,y) are respectively defined as

    Γ(x)=0ettx1dt,B(x,y)=Γ(x)Γ(y)Γ(x+y).

    Clearly, arcsinp,q(x) is an increasing homeomorphism from [0,1] onto [0,πp,q/2], and its inverse function is the generalized (p,q)-sine function sinp,q is defined on the interval [0,πp,q/2]. Moreover, the function sinp,q can be extended to the interval [0,πp,q] by

    sinp,q(x)=sinp,q(πp,qx),x[πp,q/2,π].

    sinp,q can be also extended to the whole R, and the generalized (p,q)-sine function reduces to the classical sine function for p=q=2.

    Applying the definitions of sinp,q(x) and πp,q, we can define the generalized (p,q)-elliptic integrals of the first kind Kp,q and of the second kind Ep,q by

    Kp,q(r)=πp,q/20dt(1rqsinqp,q(t))11/p=10dt(1tq)1/p(1rqtq)11/p

    and

    Ep,q(r)=πp,q/20(1rqsinqp,q(t))1/pdt=10(1rqtq1tq)1/pdt

    respectively, for p,q(1,), r(0,1).

    As a special case of the Gaussian hypergeometric function, these generalized (p,q)-elliptic integrals can be represented by Gaussian hypergeometric functions [15] as

    {Kp,q=Kp,q(r)=πp,q22F1(11p,1q;11p+1q;rq)Kp,q=Kp,q(r)=Kp,q(r)Kp,q(0)=πp,q2,Kp,q(1)= (1.2)

    and

    {Ep,q=Ep,q(r)=πp,q22F1(1p,1q;11p+1q;rq)Ep,q=Ep,q(r)=Ep,q(r)Ep,q(0)=πp,q2,Ep,q(1)=1, (1.3)

    where p,q(1,), r(0,1), r=(1rq)1/q. If p=q=2, we can derive the classical complete elliptic integrals K and E, which are well-known complete elliptic integrals of the first kind and second kind, respectively. These complete elliptic integrals play an important role in many branches of quasiconformal mapping, complex analysis, and physics.

    In 2017, Yang et al. [24] showed that the ratio K(r)/ln(c/r) is strictly concave if and only if c=e4/3 on (0,1), and K(r)/ln(1+4/r) is strictly convex on (0,1).

    In 2019, Wang et al. [22] presented the convexity of the function (Ear2Ka)/r2 and some properties of the function αKa(rα) with respect to the parameter α.

    In 2020, Huang et al. [13] established monotonicity properties for certain functions involving the complete p-elliptic integrals of the first and second kinds. They also presented the inequality π/2log2+log(1+1/r)+α(1r)<Kp(r)<π/2log2+log(1+1/r)+β(1r) which holds for all r(0,1) with the best possible constants α and β. Moreover, these generalized elliptic integrals have significant applications in the theory of geometric functions and in the theory of mean values. More properties and applications of these integrals are given in [7,13,17,21,22,24,25].

    The generalized (p,q)-elliptic integrals of the first kind Kp,q and of the second kind Ep,q satisfy the following Legendre relation:

    Kp,q(r)Ep,q(r)+Kp,q(r)Ep,q(r)Kp,q(r)Kp,q(r)=πp,q2. (1.4)

    This relation has important applications in many areas of mathematics and physics, including celestial mechanics, quantum mechanics, and statistical mechanics. when p=q=2, the equation reduces to the classical Legendre relation.

    Inspired by the papers [22], [24] and [13], we can consider extending the results of K(r), Ka(r) and Kp(r) to the generalized (p,q)-elliptic integrals of the first kind Kp,q.

    A generalized (p,q)-modular equation of degree h>0 is

    2F1(a,b;c;1sq)2F1(a,b;c;sq)=h2F1(a,b;c;1rq)2F1(a,b;c;rq),r(0,1), (1.5)

    which a,b,c>0 with a+bc. Using the decreasing homeomorphism μp,q:(0,1)(0,) defined by

    μp,q(r)=πp,q2Kp,q(r)Kp,q(r),

    for p,q(1,). The function μp,q is called the generalized (p,q)-Gr¨otzsch ring function, we can rewrite (1.5) as

    μp,q(s)=hμp,q(r),r(0,1). (1.6)

    The solution of (1.6) is given by

    s=φp,qK(r)=μ1p,q(μp,q(r)/K). (1.7)

    For p,q(1,), r(0,1), K(0,), we have

    φp,qK(r)q+φp,q1/K(r)q=1. (1.8)

    The function φp,qK(r) is referred to as the generalized (p,q)-Hersch-Pfluger distortion function with degree K=1/h. For p=q=2, the functions μp,q(r) and φp,qK(r) reduces to well-known special cases that Gr¨otzsch ring function μ(r) and Hersch-Pfluger distortion function φK(r), respectively, which play important role in the theory of plane quasiconformal mappings.

    In 2015, Alzer et al. [2] studied the monotonicity, convexity, and concavity properties of the function μ(rα)/α, and established various Gr¨otzsch ring functional inequalities based on these properties.

    In [3], the authors focused on studying the properties of the generalized Gr¨otzsch ring function μa(r) and the generalized Hersch-Pfluger distortion function φaK(r). They derived various inequalities involving μa(r) and r1/KφaK(r) by utilizing the monotonicity, concavity, and convexity of the functions μa(r), r1/KφaK(r), log(1/φaK(ex)), and φaK(rx)/φaK(x).

    In 2022, Lin et al. [16] explored the monotonicity and convexity properties of the function μp,q(r) and obtained sharp functional inequalities that sharpen and extend some existing results on the modulus of μ(r).

    Inspired by the papers [3], [12], and [16], our motivation is to extend the existing results of the functions φK(r) and φaK(r) to the generalized (p,q)-Hersch-Pfluger distortion function. Our goal is to gain the properties of the generalized (p,q)-Hersch-Pfluger distortion function and derive sharp functional inequalities for this function.

    Our main objective of this paper is to investigate various properties of the generalized (p,q)-elliptic integrals and the generalized (p,q)-Hersch-Pfluger distortion function. We specifically focus on establishing complete monotonicity, logarithmic, geometric concavity, and convexity properties of certain functions involving these generalized integrals and arcsine functions. Additionally, they derive several sharp functional inequalities for the generalized (p,q)-Hersch-Pfluger distortion function, which improve upon and generalize existing results. Apart from the introduction, this paper consists of three additional sections. Section 2 contains some preliminaries as well as several formulas and lemmas. In Section 3, we present some of the major results regarding the generalized (p, q)-elliptic integrals and provide their proof. In Section 4, we study the generalized (p, q)-Hersch-Pfluger distortion function, and present some of the major results and provide their proof.

    In this section, we present several formulas and lemmas that have been extensively utilized in the paper. These formulas and lemmas play a crucial role in the analysis and proofs of the major results. Throughout this paper, we denote p,q(1,), r(0,1) and r=(1rq)1/q.

    Lemma 2.1. [16] Derivative formulas:

    (1)dKp,qdr=Ep,qrqKp,qrrq,(2)dEp,qdr=q(Ep,qKp,q)pr,(3)dμp,q(r)dr=π2p,q4rrqK2p,q,(4)d(Ep,qrqKp,q)dr=(pq)(Kp,qEp,q)+p(q1)rqKp,qpr,(5)d(Kp,qEp,q)dr=(pqrq)Ep,q+(qp)rqKp,qprrq.

    Based on the derivative formula Lemma 2.1, we derive the derivative formulas of the function φp,qK(r) in the following lemma.

    Lemma 2.2. Let p,q(1,), then

    φp,qK(r)r=ssqKp,q(s)Kp,q(s)rrqKp,q(r)Kp,q(r)=1KssqKp,q(s)2rrqKp,q(r)2=KssqKp,q(s)2rrqKp,q(r)2

    for r(0,1).

    Proof. By the definitions of s=φp,qK(r), μp,q(s)=μp,q(r)/K, and the derivative formulas of Kp,q(r), we can derive the following equation

    π2p,q4ssqKp,q(s)2sr=1Kπ2p,q4rrqKp,q(r)2,

    then the derivative formulas of φp,qK(r) is following.

    Lemma 2.3. [4] For p[0,), let I=[0,p), and suppose that f,g:I[0,) are functions such that f(x)/g(x) is decreasing on I{0} and g(0)=0, g(x)>0 for x>0. Then

    f(x+y)(g(x)+g(y))g(x+y)(f(x)+f(y)),

    for x,y,x+yI. Moreover, if the monotoneity of f(x)/g(x) is strict, then the above inequality is also strict on I{0}.

    The following result is a monotone form of L'Hˆopital's Rule [4] and will be useful in deriving monotoneity properties and obtaining inequalities.

    Lemma 2.4. [4] For <a<b<, let f,g:[a,b]R be continuous on [a,b], and be differentiable on (a,b). Let g(x)0 on (a,b). If f(x)/g(x) is increasing (decreasing) on (a,b), then so are

    [f(x)f(a)]/[g(x)g(a)]and[f(x)f(b)]/[g(x)g(b)].

    If f(x)/g(x) is strictly monotone, then the monotonicity in the conclusion is also strict.

    The following lemma presents some known results of generalized (p,q)-elliptic integrals, which can be utilized to prove the main results of this paper.

    Lemma 2.5. [14] For p,q(1,), r(0,1), a=11/p,b=a+1/q, then the functions

    (1) h1(r)=Ep,q(r)rqKp,q(r)rq is strictly increasing and convex from (0,1) onto (aπp,q/(2b),1).

    (2) h2(r)=Ep,q(r)rqKp,q(r)rqKp,q(r) is strictly decreasing from (0,1) onto (0,a/b).

    (3) h3(r)=rcKp,q(r) is decreasing(increasing) on (0,1) iff ca/b (c0 respectively) with h3((0,1))=(0,πp,q/2) if ca/b.

    (4) h4(r)=Kp,q(r)Ep,q(r)rqKp,q(r) is increasing from (0,1) onto (1/(qb),1).

    (5) h5(r)=rqKp,q(r)/Ep,q(r) is strictly decreasing from (0,1) onto itself.

    Lemma 2.6. For r(0,1), K,p,q(1,), let s=φp,qK(r),t=φp,q1/K(r).

    (1) The function f(r)=Kp,q(s)/Kp,q(r) is increasing from (0,1) onto (1,K).

    (2) For q>3p4p1, the function g(r)=sq/2Kp,q(s)2rq/2Kp,q(r)2 is decreasing from (0,1) onto (0,1).

    Proof. (1) According to the Lemma 2.1, we have

    K2p,q(r)f(r)=Kp,q(s)rrqKp,q(r){Kp,q(s)[Ep,q(s)sqKp,q(s)]Kp,q(r)[Ep,q(r)rqKp,q(r)]}.

    Denote

    f1(r)=Kp,q(r)[Ep,q(r)rqKp,q(r)]=rqKp,q(r)[Ep,q(r)rqKp,q(r)rq].

    By applying Lemma 2.5(1)(3) and considering s>r, we can conclude that f1(r) is increasing. Hence, we can determine that f(r) is positive. Therefore, we can deduce that f(r) is an increasing function. For the limiting values, we have limr0+f(r)=1 and limr1f(r)=K.

    (2) Let

    g1(r)=Kp,q(r)[qrqKp,q(r)4(Ep,q(r)rqKp,q(r))]=rqKp,q(r)Kp,q(r)[q4Ep,q(r)rqKp,q(r)rqKp,q(r)].

    By differentiation, we have

    [rq/2Kp,q(r)2]2g(r)=sq/2Kp,q(s)2Kp,q(r)2rrq/2Kp,q(r)(g1(s)g1(r)).

    If q>3p4p1, g1(r) is increasing by Lemma 2.5(2) and (3), thus g(r) is negative for s>r. Hence g(r) is decreasing, the limiting values follow from the definitions (1.2) and (1.7)

    limr0+g(r)=1,limr1g(r)=0.

    Lemma 2.7. For r(0,1) and p,q(1,), the inequality

    rqKp,q(r)2N(r)>1q

    holds.

    Proof. Let M(r)=rqKp,q(r)2/N(r), M1(r)=N(r)/(rqKp,q(r)), we have

    M1(r)=1rqKp,q(r){qpKp,q(r)(Ep,q(r)rqKp,q(r))+Ep,q(r)(Kp,q(r)Ep,q(r))(qpq+1)rqKp,q(r)Ep,q(r)}=qpEp,q(r)rqKp,q(r)rq+Ep,q(r)Kp,q(r)Ep,q(r)rqKp,q(r)(qpq+1)Ep,q(r).

    According to Lemma 2.5, we have

    M1(r)<qp+q(11p)Ep,q(r)<qπp,q2qKp,q(r),

    and

    M(r)=Kp,q(r)/M1(r)>1/q.

    Lemma 2.8. For r(0,1), p,q(1,), the function

    h(r)=1log(1/r)Ep,q(r)rqKp,q(r)rqKp,q(r)

    is strictly increasing from (0,1) onto (0,).

    Proof. By differentiation, we have

    h(r)=1r(log(1/r))2N(r)r(rqKp,q(r))2=N(r)rq+1(log(1/r))2Kp,q(r)2[rqKp,q(r)2N(r)(rq/2rqlog(1r))2],

    where N(r)=qpKp,q(r)(Ep,q(r)rqKp,q(r))+Ep,q(r)(Kp,q(r)Ep,q(r))(qpq+1)rqKp,q(r)Ep,q(r).

    Let h2(r)=rq/2log(1/r), h3(r)=rq, then h1(r)=h2(r)/h3(r), and h2(1)=h3(1)=0.

    h2(r)h3(r)=1qrq/2[q2log(1/r)1].

    By Lemma 2.4, the function h1(r) is strictly increasing from (0,1) onto (0,1/q). According to Lemma 2.7, we conclude that rqKp,q(r)2/N(r)>1/q. Therefore, it is easy to check that h(r) is increasing. For the limiting values, limr0+h(r)=0. Since

    limr1rqKp,q(r)=0,limr1rqlog(1/r)=q,limr1Kp,q(r)=,

    we have

    limr1h(r)=limr11rqKp,q(r){rqKp,q(r)log(1/r)(Ep,q(r)rqKp,q(r))}=.

    In this section, we present some of the main results regarding the generalized (p,q)-elliptic integrals.

    Theorem 3.1. For p,q(1,), the function Fp,q(r)=(Ep,q(r)rqKp,q(r))/rq is concave on (0,r0) and convex on (r0,1) for some point r0(0,1).

    Proof. Let F(r)=(Ep,q(r)rqKp,q(r))/rq, by the definitions (1.2) and (1.3), which can be expressed as

    F(r)=Ep,q(r)rqKp,q(r)rq=πp,q2rq[F(1p,1q;11p+1q;rq)rqF(11p,1q;11p+1q;rq)]=πp,q2rq[n=0((1p,n)(1q,n)(11p+1q,n)n!(11p,n)(1q,n)(11p+1q,n)n!)rqn+n=0(11p,n)(1q,n)(11p+1q,n)n!rq(n+1)]=πp,q2n=0(11p,n)(1q,n)(11p+1q,n+1)(n+1)!(11p)(n+1)rqn=aπp,q2(a+b)2F1(a,b;a+b+1;rq),

    with a=11/p,b=1/q. This implies that

    Fp,q(r)=Ep,q(r)rqKp,q(r)rq=aπp,q2(a+b)2F1(a,b;a+b+1;rq).

    By differentiation, we have

    Fp,q(r)=a2bqπp,q2(a+b)(a+b+1)rq12F1(a+1,b+1;a+b+2;rq).

    Hence,

    2(a+b)(a+b+1)a2bqπp,qFp,q(r)=(q1)rq22F1(a+1,b+1;a+b+2;rq)q(a+1)(b+1)a+b+2r2q22F1(a+2,b+2;a+b+3;rq).

    According to [4, Theorem 1.19(10)],

    2F1(a,b;c;x)=(1x)cab2F1(ca,cb;c;x)(a+b>c,a,b,c>0),

    we have

    2(a+b)(a+b+1)a2bqπp,qFp,q(r)=(q1)rq22F1(a+1,b+1;a+b+2;rq)q(a+1)(b+1)a+b+2rq22F1(a+1,b+1;a+b+3;rq)=(q1)rq22F1(a+1,b+1;a+b+3;rq)[2F1(a+1,b+1;a+b+2;rq)2F1(a+1,b+1;a+b+3;rq)qq1(a+1)(b+1)a+b+2]. (3.1)

    Let

    F1(x)=2F1(a+1,b+1;a+b+2;x)2F1(a+1,b+1;a+b+3;x)=n=0anxnn=0bnxn,

    it is easy to obtain that F1(x) is strictly increasing from (0,1) onto (1,), since

    anbn=a+b+2+na+b+2>1.

    Similarly, we can deduce that

    r2F1(a+1,b+1;a+b+2;rq)2F1(a+1,b+1;a+b+3;rq)qq1(a+1)(b+1)a+b+2

    is strictly decreasing from (0,1) onto (1b(2a+b+2)(1b)(a+b+2),). Therefore, the sign of Fp,q(r) changes from negative to positive on (0,1) by (3.1), we know that there exists r0(0,1) such that Fp,q(r) is concave on (0,r0) and convex on (r0,1).

    Theorem 3.2. For p,q(1,), r(0,1), α>0, Let Hp,q(α)=αKp,q(rα), Gp,q(α)=Kp,q(rα)/α. Then

    (1) The function αHp,q(α) is strictly increasing and log-concave on (0,);

    (2) The function α1/Hp,q(α) is strictly convex on (0,);

    (3) The function αGp,q(α) is strictly decreasing and log-convex on (0,).

    Proof. (1) Let t=rα, then dt/dα=tlogr<0, and

    dHp,q(α)dα=Kp,q(t)+αEp,q(t)tqKp,q(t)tqlogr=Kp,q(t)Ep,q(t)tqKp,q(t)tqlog(1t)=Kp,q(t)log(1t)[1log(1/t)Ep,q(t)tqKp,q(t)tqKp,q(t)].

    Hence, the monotonicity of Hp,q(α) follows from Lemma 2.8.

    By logarithmic differentiation,

    d(logHp,q(α))dα=1αEp,q(t)tqKp,q(t)tqKp,q(t)log(1r)=log(1r)[1log(1/t)Ep,q(t)tqKp,q(t)tqKp,q(t)].

    It is not difficult to verify that d(logHp,q(α))/dα is strictly increasing with respect to t by Lemma 2.8, and is strictly decreasing with respect to α. Thus the function αHp,q(α) is strictly increasing and log-concave on (0,).

    (2) Since t=rα, then α=log(1/t)/log(1/r). Differentiating 1/Hp,q(α) yields

    ddα(1Hp,q(α))=1Hp,q(α)2Kp,q(t)log(1t)[1log(1/t)Ep,q(t)tqKp,q(t)tqKp,q(t)]=1α2Kp,q(t)2Kp,q(t)log(1t)h(t)=(log1r)2h(t)log(1/t)Kp,q(t), (3.2)

    where h(t) is defined in Lemma 2.8. Let f(r)=log(1/r)Kp,q(r), f1(r)=log(1/r), f2(r)=1/Kp,q(r), We clearly see that f1(1)=f2(1)=0, then

    f1(r)f2(r)=rqKp,q(r)2Ep,q(r)rqKp,q(r)=rqEp,q(r)rqKp,q(r)rqKp,q(r)2rq,

    which is decreasing follows from Lemma 2.5(1), (3). Hence the function f(r) is decreasing from (0,1) onto (0,) by Lemma 2.4. Therefore, it follows from Lemma 2.8 and (3.2) that the function α1/Hp,q(α) is strictly convex on (0,).

    (3) Since t=rα, and dt/dα=tlogr<0, simple computations yields

    dGp,q(α)dα=1α[Ep,q(t)tqKp,q(t)tqlogrKp,q(t)α]=1α[Ep,q(t)tqKp,q(t)tqlogrKp,q(t)logrlogt]=Kp,q(t)logrα[Ep,q(t)tqKp,q(t)tqKp,q(t)+1log(1/t)].

    Let g(r)=Ep,q(r)rqKp,q(r)rqKp,q(r)+1log(1/r), we have

    g(r)=N(r)rq+1(log(1/r))2Kp,q(r)2[rqKp,q(r)2N(r)+(rq/2rqlog(1r))2],

    is positive by Lemma 2.8, where the defition of N(r) is in Lemma 2.8. Hence the function g(r) is increasing from (0,1) onto (0,). Since logr<0, the monotonicity of the function Gp,q(α) follows immediately.

    Since logGp,q(α)=logKp,q(t)logα, by differentiation, we obtain that

    d(logGp,q(α))dα=Ep,q(t)tqKp,q(t)tqKp,q(t)logr1α=Ep,q(t)tqKp,q(t)tqKp,q(t)logrlogrlogt=logr[Ep,q(t)tqKp,q(t)1+1log(1/t)]. (3.3)

    It follows from (3.3) and Lemma 2.5(5) that d(logGp,q(α))/dα is strictly decreasing with respect to t. Therefore, Gp,q(α) is log-convex on (0,) with respect to α.

    Next, we apply Theorem 3.2 to obtain the inequality involving the generalized (p,q)-elliptic integrals Kp,q.

    Corollary 3.1. For p,q(1,).

    (1) Let α,β be positive numbers with α>β>0. The double inequality

    1<Kp,q(rβ)Kp,q(rα)<αβ

    holds for all r(0,1).

    (2) Inequality

    Kp,q(xy)2log(1/x)log(1/y)log[1/(xy)]Kp,q(x)Kp,q(y)

    holds with equality if and only if x=y for all x,y(0,1).

    (3) Inequality

    4log[1/(xy)]Kp,q(xy)1log(1/x)Kp,q(x)+1log(1/y)Kp,q(y)

    holds with equality if and only if x=y for all x,y(0,1).

    (4) Let α,β be positive numbers with α>β>0. The double inequality

    βα<Kp,q(rβ)Kp,q(rα)

    holds for all r(0,1).

    (5) Inequality

    Kp,q(xy)12log(1xy)Kp,q(x)Kp,q(y)log(1/x)log(1/y)

    holds with equality if and only if x=y for all x,y(0,1).

    Proof. (1) By utilizing the monotonicity of the function Hp,q(α) stated in Theorem 3.2, along with the monotonicity of the function Kp,q(r), we can establish that αKp,q(rα)>βKp,q(rβ). Consequently, the double inequality holds.

    (2) Since the function αHp,q(α) is strictly log-concave on (0,), we can deduce that

    logHp,q(α+β2)12(logHp,q(α)+logHp,q(β))Hp,q(α+β2)Hp,q(α)Hp,q(β)

    with equality if and only if α=β for α,β>0. For x,y(0,1) and set

    α=log(1/x)log(1/r),β=log(1/y)log(1/r).

    Simple computations yields

    Hp,q(α)=αKp,q(rα)=log(1/x)log(1/r)Kp,q(x),Hp,q(β)=βKp,q(rβ)=log(1/y)log(1/r)Kp,q(y), (3.4)
    Hp,q(α+β2)=12log[1/(xy)]log(1/r)Kp,q(xy). (3.5)

    Hence, the inequality

    Kp,q(xy)2log(1/x)log(1/y)log[1/(xy)]Kp,q(x)Kp,q(y)

    hold with equality if and only if x=y.

    (3) Since the function α1/Hp,q(α) is strictly convex on (0,), we get

    1Hp,q(α+β2)12(1Hp,q(α)+1Hp,q(β)) (3.6)

    with equality if and only if α=β for α,β>0. Set

    α=log(1/x)log(1/r),β=log(1/y)log(1/r).

    From (3.4)–(3.6), we conclude that the inequality hold with equality if and only if x=y.

    (4) Since the function αGp,q(α) is strictly decreasing, and the monotonicity of the function Kp,q(r), we have

    Kp,q(rα)α<Kp,q(rβ)β.

    (5) Since the function αGp,q(α) is log-convex on (0,),

    logGp,q(α+β2)12(logGp,q(α)+logGp,q(β))Gp,q(α+β2)Gp,q(α)Gp,q(β)

    with equality if and only if α=β for α,β>0. Set

    α=log(1/x)log(1/r),β=log(1/y)log(1/r).

    Simple computations yields

    Gp,q(α)=Kp,q(rα)α=log(1/r)log(1/x)Kp,q(x),Gp,q(β)=Kp,q(rβ)β=log(1/r)log(1/y)Kp,q(y),
    Gp,q(α+β2)=2log(1/r)log[1/(xy)]Kp,q(xy).

    Hence, the inequality

    Kp,q(xy)12log(1xy)Kp,q(x)Kp,q(y)log(1/x)log(1/y)

    hold with equality if and only if x=y.

    Remark 3.1. In [22], Wang et al. provided the proof for the convexity of the function (Ear2Ka)/r2 and presented certain properties of the functions αKa(rα) and 1/αKa(rα) with respect to the parameter α. It is worth noting that Theorem 3.1 and Theorem 3.2(1), (2) can be reduced to [22, Theorem 1.1, Theorem 1.3] if p=q=1/a.

    In this section, we study the complete monotonicity, logarithmic, geometric concavity and convexity of the generalized (p,q)-Hersch-Pfluger distortion function, and present some of the main results about φp,qK.

    Theorem 4.1. For K,p,q(1,), and q>3p4p1, let a=11/p,b=1/q, f,g be defined on (0,1] by

    f(r)=r1/Kφp,qK(r),g(r)=rKφp,q1/K(r).

    Then f is decreaseing from (0,1] onto [1,ebR(a,b)(11/K)), and g is increasing from (0,1] onto (ebR(a,b)(1K),1].

    Proof. Let s=φp,qK(r), then f(r)=sr1/K, we have

    (r1/K)2f(r)=sKr1/K1[(sq/2Kp,q(s)rq/2Kp,q(r))21].

    Hence,

    f(r)f(r)=1Kr[(sq/2Kp,q(s)rq/2Kp,q(r))21].

    According to Lemma 2.6(2), f(r) is negative. Combine with f(1)=1 and by [18, Theorem 2],

    limr0+log(r1/Ks)=limr0+[(μp,q(s)+logs)1K(μp,q(r)+logr)]=bR(a,b)(11K).

    Let t=φp,q1/K(r), thus r=φp,qK(t) and

    g(r)=φp,qK(t)Kt=(t1/Kφp,qK(t))K=f(t)K.

    According to the monotonicity of f(r), g(r) is increasing on (0,1]. The limiting values can also be derived from [18, Theorem 2].

    Next, we utilize Theorem 4.1 to derive the inequality concerning the generalized (p,q)-Hersch-Pfluger distortion function φp,qK.

    Corollary 4.1. For K,p,q(1,), and q>3p4p1, let a=11/p,b=1/q, then

    (1) The double inequality

    |φp,qK(r)φp,qK(s)|φp,qK(|rs|)ebR(a,b)(11/K)|rs|1/K (4.1)

    hold with equality if and only if r=s.

    (2) The double inequality

    |φp,q1/K(r)φp,q1/K(s)|φp,q1/K(|rs|)ebR(a,b)(1K)|rs|K (4.2)

    hold with equality if and only if r=s.

    Proof. (1) According to Lemma 2.3 and the monotonicity of the function f(r)=r1/Kφp,qK(r), we can conclude that

    φp,qK(x+y)φp,qK(x)+φp,qK(y)

    for x,y(0,1). Set r=x+y and s=y, we get

    |φp,qK(r)φp,qK(s)|φp,qK(|rs|).

    According to f(r) is decreaseing from (0,1] onto [1,ebR(a,b)(11/K)), we obtain that

    φp,qK(|rs|)ebR(a,b)(11/K)|rs|1/K

    with equality if and only if r=s.

    (2) By Lemma 2.3 and the monotonicity of g(r), we have

    φp,q1/K(x)+φp,q1/K(y)φp,q1/K(x+y)

    for x,y(0,1). Set r=x+y and s=y, we obtain that

    |φp,q1/K(r)φp,q1/K(s)|φp,q1/K(|rs|)

    with equality if and only if r=s.

    According to g(r) increasing from (0,1] onto (ebR(a,b)(1K),1], we get

    φp,q1/K(|rs|)ebR(a,b)(1K)|rs|K

    with equality if and only if r=s. Therefore, the double inequality (4.2) hold.

    Theorem 4.2. For K,p,q(1,), q>3p4p1, the function f(x)=log(1/φp,qK(ex)) is increasing and convex on (0,), g(x)=log(1/φp,q1/K(ex)) is increasing and concave on (0,), and

    φp,qK(r)φp,qK(t)(φp,qK(rt))2,φp,q1/K(r)φp,q1/K(t)(φp,q1/K(rt))2,

    with equality if and only if K=1 for each r,t(0,1).

    Proof. Let r=ex,s=φp,qK(r), according to the Lemma 2.6, we have

    f(x)=1K(sq/2Kp,q(s)rq/2Kp,q(r))2

    is positive and increasing with respect to x. Thus f is increasing and convex. Therefore,

    f(x+y2)12(f(x)+f(y)),

    and putting r=ex,t=ey, we obtain

    φp,qK(r)φp,qK(t)(φp,qK(rt))2,

    with equality if and only if K=1 for each r,t(0,1). The proof for g(x) follows a similar approach.

    Theorem 4.3. For K,p,q(1,), q>3p4p1, r(0,1), the function f(x)=φp,qK(rx)/φp,qK(x) is increasing from (0,1) onto (r1/K,φp,qK(r)), while function g(x)=φp,q1/K(rx)/φp,q1/K(x) is decreasing from (0,1) onto (φp,q1/K(r),rK). In particular,

    φp,qK(rt)φp,qK(r)φp,qK(t),φp,q1/K(rt)φp,q1/K(r)φp,q1/K(t),

    with equality if and only if K=1 for each r,t(0,1).

    Proof. Let t=rx,u=φp,qK(t),s=φp,qK(x),

    f(x)=uKsx[(uq/2Kp,q(u)tq/2Kp,q(t))2(sq/2Kp,q(s)xq/2Kp,q(x))2],

    then

    f(x)f(x)=1Kx[(uq/2Kp,q(u)tq/2Kp,q(t))2(sq/2Kp,q(s)xq/2Kp,q(x))2].

    Since t<x and sq/2Kp,q(s)xq/2Kp,q(x) is decreasing with respect to r by Lemma 2.6, f(x) is positive on (0,1). For the limiting values, by using L'Hˆopital's Rule, we get

    limr0+f(r)=r1/K,limr1h(r)=φp,qK(r).

    Since the monotonicity of the function f(x), along with the definition of the function φp,qK(r), we obtain

    φp,qK(rt)φp,qK(r)φp,qK(t),

    with equality if and only if K=1 for each r,t(0,1). The proof for g(x) follows a similar approach. As a result, we will omit the detailed proof.

    Theorem 4.4. For K,p(1,), r(0,1), q>3 and p>q(q3)q22q1, the function f(r) defined by

    f(r)=arcsin(φp,qK(r))arcsin(r1/K)

    is strictly decreasing from (0,1] into [1,e(11/K)bR(a,b)), the function g(r) defined by

    g(r)=arcsin(φp,q1/K(r))arcsin(rK)

    is strictly increasing from (0,1] into (e(1K)bR(a,b),1], where a=11/p,b=a+1/q.

    Proof. Let s=φp,qK(r), f1(r)=arcsin(s), f2(r)=arcsin(r1/K), f(r)=f1(r)/f2(r). Since f1(0)=f2(0)=0, we have

    f1(r)f2(r)=sr1/K(1r2/K1r2)1/2sq1Kp,q(s)2rq1Kp,q(r)2.

    Let f3(r)=sq1Kp,q(s)2rq1Kp,q(r)2, according to the Lemma 2.1, we have

    [rq1Kp,q(r)2]2f3(r)=sq1Kp,q(s)2Kp,q(r)rrKp,q(r)(f4(s)f4(r)),

    where

    f4(r)=Kp,q(r)[(q1)rqKp,q(r)2(Ep,q(r)rqKp,q(r))]=rqKp,q(r)Kp,q(r)[(q1)2Ep,q(r)rqKp,q(r)rqKp,q(r)].

    Using Lemma 2.5(2) and (3), we can observe that f4(r) is positive. Consequently, we can deduce that f3(r) is decreasing. As a result, we obtain that f(r) is decreasing on the interval (0,1] by Theorem 4.1. Furthermore, it can be deduced that the monotonicity of g(r) is similar to that of f(r).

    Remark 4.1. Theorems 4.1 to 4.4 can be seen as variations and extensions of the results presented in [3, Theorem 1.14, Theorem 1.15, Theorem 6.7, Theorem 6.13]. When p=q=1/a, the results obtained in Theorems 4.1 to 4.4 can be reduced to those obtained in [3].

    In this paper, we investigate the properties of the generalized (p,q)-elliptic integrals and the generalized (p,q)-Hersch-Pfluger distortion function. Through our analysis, we have established complete monotonicity, logarithmic, geometric concavity, and convexity properties for certain functions involving these integrals and arcsine functions. These properties provide valuable insights into the behavior of these functions. Furthermore, we have derived several sharp functional inequalities for the generalized (p,q)-elliptic integrals and the generalized (p,q)-Hersch-Pfluger distortion function. These inequalities not only improve upon existing results but also generalize them.

    The authors declare that they have not used Artificial Intelligence tools in the creation of this article.

    The research was supported by the Natural Science Foundation of China (Grant Nos. 11601485, 11401531).

    The authors declare that they have no conflicts of interest.



    [1] T. Abdeljawad, On conformable fractional calculus, J. Comput. Appl. Math., 279 (2015), 57–66. https://doi.org/10.1016/j.cam.2014.10.016 doi: 10.1016/j.cam.2014.10.016
    [2] A. R. Aftabizadeh, Existence and uniqueness theorems for fourth-order boundary value problems, J. Math. Anal. Appl., 116 (1986), 415–426. https://doi.org/10.1016/S0022-247X(86)80006-3 doi: 10.1016/S0022-247X(86)80006-3
    [3] R. Agarwal, On fourth-order boundary value problems arising in beam analysis, Differ. Integr. Equ., 2 (1989), 91–110.
    [4] B. Ahmad, S. K. Ntouyas, Existence results for a coupled system of Caputo type sequential fractional differential equations with nonlocal integral boundary conditions, Appl. Math. Comput., 266 (2015), 615–622. https://doi.org/10.1016/j.amc.2015.05.116 doi: 10.1016/j.amc.2015.05.116
    [5] A. Alsaedi, M. Alnahdi, B. Ahmad, S. K. Ntouyas, On a nonlinear coupled Caputo-type fractional differential system with coupled closed boundary conditions, AIMS Math., 8 (2023), 17981–17995. https://doi.org/10.3934/math.2023914 doi: 10.3934/math.2023914
    [6] E. Alvarez, H. Cabrales, T. Castro, Optimal control theory for a system of partial differential equations associated with stratified fluids, Mathematics, 9 (2021), 1–23. https://doi.org/10.3390/math9212672 doi: 10.3390/math9212672
    [7] A. Anber, Z. Dahmani, The SGEM method for solving some time and Space-Conformable fractional evolution problems, Int. J. Open Prob. Comput. Math., 16 (2023), 33–44.
    [8] A. Anber, I. Jebril, Z. Dahmani, N. Bedjaoui, A. Lamamri, The Tanh method and the (G'/G)-expansion method for solving the space-time conformable FZK and FZZ evolution equations, Int. J. Innov. Comput. Inf. Contr., 20 (2024), 557–573. https://doi.org/10.24507/ijicic.20.02.557 doi: 10.24507/ijicic.20.02.557
    [9] I. M. Batiha, S. Alshorm, I. H. Jebril, M. A. Hammad, A brief review about fractional calculus, J. Open Prob. Comput. Math., 15 (2022), 39–56.
    [10] I. M. Batiha, S. A. Njadat, R. M. Batyha, A. Zraiqat, A. Dababneh, S. Momani, Design fractional-order PID controllers for Single-Joint robot arm model, Int. J. Adv. Soft Comput. Appl., 14 (2022), 96–114. https://doi.org/10.15849/IJASCA.220720.07 doi: 10.15849/IJASCA.220720.07
    [11] K. Bensaassa, R. Wael Ibrahim, Z. Dahmani, Existence, uniqueness and numerical simulation for solutions of a class of fractional differential problems, Submitted, 2023.
    [12] K. Bensaassa, Z. Dahmani, M. Rakah, M. Z. Sarikaya, Beam deflection coupled systems of fractional differential equations: Existence of solutions, Ulam-Hyers stability and travelling waves, Anal. Math. Phys., 14 (2024). https://doi.org/10.1007/s13324-024-00890-6
    [13] A. Carpinteri, F. Mainardi, Fractional calculus in continuum mechanics, 2 Eds., New York: Academic Press, 1997. https://doi.org/10.1007/978-3-7091-2664-6
    [14] Z. Dahmani, A. Anber, Y. Gouari, M. Kaid, I. Jebril, Extension of a method for solving nonlinear evolution equations via conformable fractional approach, Int. Conf. Infor. Tech., 2021, 38–42. http://doi.org/10.1109/ICIT52682.2021.9491735
    [15] Z. Dahmani, A. Anber, I. Jebril, Solving conformable evolution equations by an extended numerical method, Jordan J. Math. Stat., 15 (2022), 363–380. https://doi.org/10.47013/15.2.14 doi: 10.47013/15.2.14
    [16] M. A. Del Pino, R. F. Manasevich, Existence for a fourth-order boundary value problem under a two-parameter nonresonance condition, Proc. American Math. Soc., 112 (1991), 81–86. https://doi.org/10.2307/2048482 doi: 10.2307/2048482
    [17] B. M. Dia, M. S. Goudiaby, O. Dorn, Boundary feedback stabilization of Two-Dimensional shallow water equations with viscosity term, Mathematics, 10 (2022), 4036,132–143. https://doi.org/10.3390/math10214036 doi: 10.3390/math10214036
    [18] Y. Gouari, Z. Dahmani, I. Jebril, Application of fractional calculus on a new differential problem of duffing type, Adv. Math. Sci. J., 9 (2020), 10989–11002. https://doi.org/10.37418/amsj.9.12.82 doi: 10.37418/amsj.9.12.82
    [19] Y. Gouari, Z. Dahmani, S. E. Farooq, F. Ahmad, Fractional singular differential systems of Lane-Emden type: Existence and uniqueness of solutions, Axioms, 9 (2020), 95. https://doi.org/10.3390/axioms9030095 doi: 10.3390/axioms9030095
    [20] Y. Gouari, Z. Dahmani, Stability of solutions for two classes of fractional differential equations of Lane-Emden type, J. Int. Math., 24 (2021), 2087–2099. http://doi.org/10.1080/09720502.2020.1856343 doi: 10.1080/09720502.2020.1856343
    [21] R. Khalil, M. Al Horani, A. Yousef, M. Sababheh, A new definition of fractional derivative, J. Comput. Appl. Math., 264 (2014), 65–70. https://doi.org/10.1016/j.cam.2014.01.002 doi: 10.1016/j.cam.2014.01.002
    [22] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Amsterdam: Elsevier B.V., 2006.
    [23] P. Li, Y. Lu, C. Xu, J. Ren, Dynamic exploration and control of bifurcation in a fractional-order Lengyel-Epstein model owing time delay, MATCH Commun. Math. Comput. Chem., 92 (2024), 437–482.
    [24] P. Li, C. Xu, M. Farman, A. Akgül, Y. Pang, Qualitative and stability analysis with lyapunov function of emotion panic spreading model insight of fractional operator, Fractals, 32 (2024), 2440011. http://doi.org/10.1142/S0218348X24400115 doi: 10.1142/S0218348X24400115
    [25] W. Malfliet, W. Hereman, The Tanh method: I. Exact solutions of nonlinear evolution and wave equations, Phys. Scripta, 54 (1996), 563–568. http://doi.org/10.1088/0031-8949/54/6/003 doi: 10.1088/0031-8949/54/6/003
    [26] M. Marin, A. Öchsner, M. M. Bhatti, Some results in Moore-Gibson-Thompson thermoelasticity of dipolar bodies, ZAMM J. Appl. Math. Mech., 100 (2020), e202000090. https://doi.org/10.1002/zamm.202000090 doi: 10.1002/zamm.202000090
    [27] M. Marin, A. Hobiny, I. Abbas, Finite element analysis of nonlinear bioheat model in skin tissue due to external thermal sources, Mathematics, 9 (2021), 1459. https://doi.org/10.3390/math9131459 doi: 10.3390/math9131459
    [28] M. Rakah, Z. Dahmani, A. Senouci, New uniqueness results for fractional differential equations with a Caputo and khalil derivatives, Appl. Math. Inf. Sci., 16 (2022), 943–952. http://dx.doi.org/10.18576/amis/160611 doi: 10.18576/amis/160611
    [29] M. Rakah, Y. Gouari, R. Ibrahim, Z. Dahmani, H. Kahtan, Unique solutions, stability and travelling waves for some generalized fractional differential problems, Appl. Math. Sci. Engineer., 23 (2023). https://doi.org/10.1080/27690911.2023.2232092
    [30] U. Sadiya, M. Inc, M. A. Arefin, M. H. Uddin, Consistent travelling waves solutions to the non-linear time fractional Klein-Gordon and Sine-Gordon equations through extended tanh-function approach, J. Taibah Univ. Sci., 16 (2022), 594–607. https://doi.org/10.1080/16583655.2022.2089396 doi: 10.1080/16583655.2022.2089396
    [31] A. Tudorache, R. Luca, On a system of sequential Caputo fractional differential equations with nonlocal boundary conditions, Frac. Fract., 7 (2023), 1–23. https://doi.org/10.3390/fractalfract7020181 doi: 10.3390/fractalfract7020181
    [32] Q. Wang, L. Yang, Positive solution for a nonlinear system of fourth-order ordinary differential equations, Electr. J. Differ. Equat., 2020 (2020), 1–15.
    [33] A. M. Wazwaz, The Tanh method for compact and non compact solutions for variants of the KdV-Burger and the K(n, n)-Burger equations, Phys. Nonlinear Phen., 213 (2006), 147–151. https://doi.org/10.1016/j.physd.2005.09.018 doi: 10.1016/j.physd.2005.09.018
    [34] M. Xia, X. Zhang, D. Kang, C. Liu, Existence and concentration of nontrivial solutions for an elastic beam equation with local nonlinearity, AIMS Math., 7 (2021), 579–605. https://doi.org/10.3934/math.2022037 doi: 10.3934/math.2022037
    [35] Y. Yang, Q. Qi, J. Hu, J. Dai, C. Yang, Adaptive Fault-Tolerant control for consensus of nonlinear fractional order Multi-Agent systems with diffusion, Frac. Fract., 7 (2023), 760. https://doi.org/10.3390/fractalfract7100760 doi: 10.3390/fractalfract7100760
    [36] M. Younis, Soliton solutions of fractional order KdV-Burger's equation, J. Adv. Phys., 3 (2014), 325–328.
    [37] J. L. Zhou, Y. B. He, S. Q. Zhang, H. Y. Deng, X. Y. Lin, Existence and stability results for nonlinear fractional integrodifferential coupled systems, Boundary Value Pro., 10 (2023), 1–14. https://doi.org/10.1186/s13661-023-01698-2 doi: 10.1186/s13661-023-01698-2
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(986) PDF downloads(218) Cited by(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog