Research article Special Issues

Some properties for certain class of bi-univalent functions defined by q-Cătaş operator with bounded boundary rotation

  • Received: 01 August 2021 Accepted: 27 September 2021 Published: 18 October 2021
  • MSC : 30C45

  • Throughout the paper, we introduce a new subclass Hn,q,λ,lα,μ,ρ,m,β f(z) by using the Bazilevič functions with the idea of bounded boundary rotation and q-analogue Cătaş operator. Also we find the estimate of the coefficients for functions in this class. Finally, in the concluding section, we have chosen to reiterate the well-demonstrated fact that any attempt to produce the rather straightforward (p,q)-variations of the results, which we have presented in this article, will be a rather trivial and inconsequential exercise, simply because the additional parameter p is obviously redundant.

    Citation: S. M. Madian. Some properties for certain class of bi-univalent functions defined by q-Cătaş operator with bounded boundary rotation[J]. AIMS Mathematics, 2022, 7(1): 903-914. doi: 10.3934/math.2022053

    Related Papers:

    [1] Aissa Boukarou, Kaddour Guerbati, Khaled Zennir, Mohammad Alnegga . Gevrey regularity for the generalized Kadomtsev-Petviashvili I (gKP-I) equation. AIMS Mathematics, 2021, 6(9): 10037-10054. doi: 10.3934/math.2021583
    [2] Xiaochun Sun, Yulian Wu, Gaoting Xu . Global well-posedness for the 3D rotating Boussinesq equations in variable exponent Fourier-Besov spaces. AIMS Mathematics, 2023, 8(11): 27065-27079. doi: 10.3934/math.20231385
    [3] Muhammad Imran Liaqat, Fahim Ud Din, Wedad Albalawi, Kottakkaran Sooppy Nisar, Abdel-Haleem Abdel-Aty . Analysis of stochastic delay differential equations in the framework of conformable fractional derivatives. AIMS Mathematics, 2024, 9(5): 11194-11211. doi: 10.3934/math.2024549
    [4] William O. Bray, Ellen Hunter . Wave equations & energy. AIMS Mathematics, 2019, 4(3): 463-481. doi: 10.3934/math.2019.3.463
    [5] Phuong Nguyen Duc, Ahmet Ocak Akdemir, Van Tien Nguyen, Anh Tuan Nguyen . Remarks on parabolic equation with the conformable variable derivative in Hilbert scales. AIMS Mathematics, 2022, 7(11): 20020-20042. doi: 10.3934/math.20221095
    [6] Kedong Wang, Xianguo Geng, Mingming Chen, Ruomeng Li . Long-time asymptotics for the generalized Sasa-Satsuma equation. AIMS Mathematics, 2020, 5(6): 7413-7437. doi: 10.3934/math.2020475
    [7] Said Mesloub, Faten Aldosari . Well posedness for a singular two dimensional fractional initial boundary value problem with Bessel operator involving boundary integral conditions. AIMS Mathematics, 2021, 6(9): 9786-9812. doi: 10.3934/math.2021569
    [8] Zhongdi Cen, Jian Huang, Aimin Xu . A posteriori mesh method for a system of singularly perturbed initial value problems. AIMS Mathematics, 2022, 7(9): 16719-16732. doi: 10.3934/math.2022917
    [9] 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
    [10] Dumitru Baleanu, Babak Shiri . Generalized fractional differential equations for past dynamic. AIMS Mathematics, 2022, 7(8): 14394-14418. doi: 10.3934/math.2022793
  • Throughout the paper, we introduce a new subclass Hn,q,λ,lα,μ,ρ,m,β f(z) by using the Bazilevič functions with the idea of bounded boundary rotation and q-analogue Cătaş operator. Also we find the estimate of the coefficients for functions in this class. Finally, in the concluding section, we have chosen to reiterate the well-demonstrated fact that any attempt to produce the rather straightforward (p,q)-variations of the results, which we have presented in this article, will be a rather trivial and inconsequential exercise, simply because the additional parameter p is obviously redundant.



    Fractional calculus (FC) is a study of derivatives and integration of arbitrary order. At the initial phases of its development, it was considered an abstract mathematical idea with nearly no applications. But, now the situation is entirely different with FC and is considered to be the most useful and important topic among the scientific community. The various important phenomena of science and engineering have been well described with the use of fractional derivatives, including, partial bed-load transport, diffusion model, dynamics of earthquakes, viscoelastic systems, biological systems, chaos, and wave propagation, (see [1,2,3,7,8,9,10,41,42]) and references therein. The nonlocal property of the fractional operators makes them more efficient for modeling the various problems of physics, fluid dynamics and their related disciplines, (see [17,18,43,44]).

    Most of the classical (integer order) and FDEs are not possible to solve analytically, therefore, some reliable numerical methods and semi-analytical methods have been developed to solve them. Some of them are presented as follows: In [34], the authors extended the study of [45] to solve numerically the Caputo type tempered fractional integro-differential equations. In [36], the authors solved the various kind of Riemann problems by developing a mixed-method, namely, the characteristic Adomian decomposition method (CADM) which is the improved form of ADM. In [38], the authors constructed the power series solution of the time-fractional Majda-Biello system by applying the Lie group analysis. For more information on Lie group analysis and its applications, (see [35]). In [39], the authors investigated the behaviour of cavitation models together with thermodynamic effects by developing a compressible, two-phase, one-fluid Navier-Stokes solver. In [40], the authors presented the analytical and computational study of the drift-flux two-phase flow model. The analytical solution of the model is based on the analytical solution of the Riemann problem which consists of two different nonlinear algebraic equations in velocity equilibrium and pressure nonequilibrium between the two-phase systems. The numerical solution is obtained by using the finite volume methods. In [37], the authors implemented a void ratio transport-equation model in a one-fluid two-phase compressible software by proposing a various numerical techniques which are based on a finite volume approach, including, HLLC, Rusanov, VF Roe, Jameson-Schmidt-Turkel and AUSM. For more detail on transport-equation and its applications, (see [33]). In [43], the authors proposed the fully discretized finite element approximations to solve the time-fractional diffusion equations of variable-order.

    Orthogonal polynomials coupled with operational matrices of fractional derivatives operators defined with singular or nonsingular kernels have played a significant role in the development of the spectral methods. The most commonly used methods are Spectral tau methods, and Spectral collocation methods. These methods have been successfully used to find the approximate solution of a large class of multi-order ordinary and partial FDEs. Some of them are presented as follows: In [13], the authors solved the multi-order (MOFDEs) by combining the technique of operational matrices with Spectral tau and Spectral collocation methods. In [14], the authors solved the MOFDEs by using the Chebyshev operational matrix (COM) together with Spectral tau and Spectral collocation methods. In [15], the authors solved the MOFDEs by reducing them into a system of FDEs, then combining the COM with Spectral tau and Spectral collocation methods. In [19], the authors derived the derivative operational matrix in Caputo sense of FLFVs which then applied with Spectral tau method to solve MOFDEs. In [20], the authors derived the integral COM in Riemann-Liouville sense which then applied with Spectral tau method to solve numerically the MOFDEs. In [22], the author derived the Jacobi integral operational matrix in Riemann-Liouville sense which then applied with Spectral tau method to solve MOFDEs. However, these methods have certain numerical difficulties in applying: in spectral tau methods, the residual function must be expanded as a series of orthogonal polynomials and then the initial or boundary conditions are applied as constraints; whereas in spectral collocation methods, the requirement of the choice of the suitable collocation points are necessary and FDEs must be satisfied exactly at these points.

    In this article, we propose a numerical method which is completely based on the fractional derivative and integral operational matrices of FLFVs. The proposed method is independent of the choice of the suitable collocation points and expansion of the residual function as a series of orthogonal polynomials. Consequently, the proposed method produces high efficient numerical results as compared to the Spectral tau method [13], Bessel collocation method [6], Taylor matrix method [5], function approximation theory [21], and stochastic technique [4]. The other novel aspect of our article is the development of new integral and derivative operational matrices in Riemann-Liouville and Caputo senses respectively. As an applications of the proposed method, we solve various MOFDEs corresponding to initial conditions. Our approach has ability to reduce the MOFDEs into a system of Sylvester types matrix equations which can be solved using MATLAB built in function lyap(.). The rest of the article is organized as follows:

    In Section 2, some preliminaries definitions of FC are recalled. In Section 3, the analytical expressions and some useful properties of CSLPs are reviewed. Moreover, in the same section some important properties of FLFVs and its analytical expression are recalled. In Section 4, the new generalized integral operational matrix and derivative operational matrix of FLFVs in the senses of Riemann-Liouville and Caputo are derived. In Section 5, based on the results of Section 4, a numerical method is proposed which is an excellent tool to solve numerically the MOFDEs corresponding to initial conditions. In Section 6, the validity and the efficiency of the proposed numerical method is tested by solving various examples and comparing the results with other methods in the literature. Some concluding remarks and future directions are given in Section 7.

    Some important definitions and results are recalled in this section which are indispensable to construct the proposed numerical algorithm.

    Definition 1. The (left-hand side) Riemann-Liouville integral of order Ω>0 is given as (see [12]):

    RLJΩu(x)=1Γ(Ω)x0(xt)Ω1u(t)dt,x>0. (2.1)

    Definition 2. The (left-hand side) Caputo derivative of order Ω>0 is given as (see [12]):

    CDΩu(x)={dnu(x)dxn,Ω=nN,RLJnΩu(n)(x),n1<Ω<n ,nN, (2.2)

    where n is an integer, x>0, and u(x)Cn[0,1].

    For the Caputo derivative, we have the following observations:

    CDΩK=0,Kis a constant.  (2.3)
    CDΩxμ={0,forμR+,andμ<Ω,Γ(μ+1)Γ(μ+1Ω)xμΩ,forμR+,andμΩ. (2.4)

    Many approximate techniques are available in the literature to solve numerically the FDEs but the most frequent used approach is the approximation of the solutions via series expansion of the type, Mk=0dktkΩ. Thus the choice of the appropriate orthogonal function of the type Ψ(t)=Mk=0dktkΩ is very important to get the good approximate results.

    In recent, the authors in [23] worked on the development of the fractional representation of the classical Legendre polynomials using the well-known Rodrigues formula. However, it is difficult to solve FDEs with this proposed fractional extension due to complexity involved in these functions, so in order to obtain the solution of FDEs, in a much simpler and efficient manner, Kazem and coauthors generated the orthogonal FLFVs based on the orthogonal CSLPs to find the approximate solutions of linear and nonlinear FDEs (see [19]). In this section, the definition and some useful properties of orthogonal CSLPs and its fractional extension are reviewed.

    The following recurrence formulae is used to evaluate the orthogonal Legendre polynomials (LPs) on the interval of orthogonality [1,1] (see [11]):

    Rj+1(x)=2j+1j+1xRj(x)jj+1Rj1(x),j=1,2,,whereR0(x)=1,R1(x)=x. (3.1)

    Now by setting, x=2t1 in (3.1), the CSLPs on the interval of orthogonality [0,1] can be expressed as following, (see [11])

    Rj(t)=(2j+1)(2t1)j+1Rj(t)jj+1Rj1(t),j=1,2,,whereR0(t)=1,R1(t)=2t1. (3.2)

    The analytical expression of (3.2) is as following

    Rj(t)=jr=0ar,jtr,wherear,j=(1)r+jΓ(1+r+j)Γ(1r+j)(Γ(1+r))2. (3.3)

    The orthogonality conditions of CSLPs can be expressed as

    10Rj(t)Ri(t)={12j+1,forj=i,0,forji. (3.4)

    The orthogonal FLFVs is defined in the following relation by changing the variable, t with zΩ,ΩR+ in orthogonal CSLPs defined in (3.2). For our convenience, we use the notation FRΩj(z) to represent the FLFVs (Rj(zΩ)).

    FRΩj(z)=(2j+1)(2zΩ1)j+1FRΩj(z)jj+1FRΩj1(z),j=1,2,,whereFRΩ0(z)=1,FRΩ1(z)=2zΩ1. (3.5)

    The analytical representation of (3.5) is as following

    FRΩj(z)=jq=0aq,jzqΩ,whereaq,j=(1)q+jΓ(1+q+j)Γ(1q+j)(Γ(1+q))2. (3.6)

    The orthogonality conditions of FLFVs can be expressed as

    10FRΩj(z)FRΩi(z)w(z)dz={1Ω(2j+1),forj=i,0,forji, (3.7)

    where w(z)=zΩ1 is a weight function.

    If u(z)L2[0,1], then it can be demonstrated as a series expansion of orthogonal FLFVs, as

    u(z)=j=0ejFRΩj(z). (3.8)

    Using (3.7), the coefficients ej can be determined as

    ej=Ω(2j+1)10u(z)FRΩj(z)w(z)dz,j=0,1,2,. (3.9)

    For the practicality, considering the first n+1-terms of (3.8), we have

    u(z)nj=0ejFRΩj(z),=(Λ)TΨ(z), (3.10)

    where

    (Λ)T=[e0,e1,,en],

    and

    Ψ(z)=[FRΩ0(z),FRΩ1(z),FRΩ2(z),,FRΩn(z)]T. (3.11)

    For n=4,Ω=1/2, we have

    Ψ(z)=(12z16z6z+112z30z+20z32190z+70z220z140z32+1),

    and for u(z)=1+z, we have

    (Λ)T=(43121600).

    Operational matrices approach together with Spectral Tau approach is among of the approaches that has been frequently used to solve numerically the ordinary and partial FDEs (see [13,14,15,19,20,21,22,24,25,26]) and references therein. Under this approach, the numerical solutions are obtained by reducing the under study FDEs both ordinary and partial to an equivalent system of algebraic equations. The solution of the algebraic equations leads to the solution of the original FDEs. In recent, the authors in [19] extended the study of the generalized derivative operational matrix of CSLPs [13] to the generalized derivative operational matrix of orthogonal FLFVs in Caputo sense to solve numerically the multi-order FDEs. In this section, based on the study in [13,19], we develop the new generalized integral operational matrix of orthogonal FLFVs in Riemann-Liouville sense.

    Lemma 4.1. The Caputo derivative of order β>0 of FLFVs can be expressed as

    CDβFRΩj(z)=jq=0a(q,j)Γ(qΩ+1)Γ(qΩβ+1)zqΩβ, (4.1)

    where a(q,j)=0 when qΩR+ and qΩ<β. For other cases, a(q,j)=a(q,j).

    Proof. For the proof (see [19]).

    Lemma 4.2. The Riemann Liouville integral of order α>0 of FLFVs can be expressed as

    RLJαFRΩj(z)=jq=0a(q,j)Γ(qΩ+1)Γ(qΩ+α+1)zqΩ+α. (4.2)

    Proof. Using linearity property and Eq (3.5), we have

    RLJαFRΩj(z)=jq=0a(q,j)RLJαzqΩ (4.3)

    Now the lemma can be proved easily using ([27], Eq (22)).

    Lemma 4.3. Suppose, ΩR+, and α>0, then

    zα+qΩ=ni=0c(r,i)FRΩi(z),wherec(r,i)=Ω(2i+1)ir=0(1)(i+r)Γ(i+r+1)Γ(ir+1)(Γ(r+1))2Ω(q+1+r)+α,i=0,1,,n. (4.4)

    Proof. Using (n+1) orthogonal FLFVs, we may approximate zβ+qΩ, as

    zα+qΩni=0c(r,i)FRΩi(z). (4.5)

    Using Eq (3.9), and after doing some lengthy calculations, we have

    c(r,i)=Ω(2i+1)ir=0(1)(i+r)Γ(i+r+1)Γ(ir+1)(Γ(r+1))2Ω(q+1+r)+α. (4.6)

    Now using the Eqs (4.5) and (4.6), we can obtain the required result.

    Lemma 4.4. Suppose, ΩR+, and β>0, then

    zqΩβ=ni=0c(r,i)FRΩi(z),wherec(i)=Ω(2i+1)ir=0(1)(i+r)Γ(i+r+1)Γ(ir+1)(Γ(r+1))2Ω(q+1+r)β,i=0,1,,n. (4.7)

    Proof. Using (n+1) orthogonal FLFVs, we may approximate zqΩβ, as

    zqΩβni=0c(r,i)FRΩi(z). (4.8)

    Using Eq (3.9), and after doing some lengthy calculations, we have

    c(i)=Ω(2i+1)ir=0(1)(i+r)Γ(i+r+1)Γ(ir+1)(Γ(r+1))2Ω(q+1+r)β. (4.9)

    Now using the Eqs (4.8) and (4.9), we can obtain the required result.

    In this section, the generalized operational matrices of fractional integral and fractional derivative for FLFVs in Riemann Liouville and Caputo senses are studied respectively.

    Theorem 4.5. Suppose Ψ(z) be the orthogonal FLFVs as defined in (3.11), and α>0 be the order of the integral consider in Riemann-Liouville sense, then

    RLJαΨ(z)Dα(n+1,n+1)Ψ(z), (4.10)

    where Dα is the (n+1)×(n+1) operational matrix which can be computed as

    Dα(n+1,n+1)=jq=0a(q,j)c(r,i)Γ(qΩ+1)(qΩ+α+1),j=0,1,,n. (4.11)

    Proof. Using (3.11) and linearity of Riemann-Liouville fractional integral, we have

    RLJαΨ(z)=jq=0a(q,j)RLJα(zqΩ) (4.12)

    Using lemma (4.2), Eq (4.11) can be expressed as

    RLJαΨ(z)=jq=0a(q,j)Γ(qΩ+1)Γ(qΩ+α+1)zα+qΩ. (4.13)

    Now Eq (4.13) and Lemma 4.3 proves the required result.

    Remark 4.6. The operational matrix of integration for CSLPs studied by Khan and Khalil [24] is a special case of Theorem 4.5 for Ω=1.

    Theorem 4.7. Suppose Ψ(z) be the orthogonal FLFVs as defined in (3.11), and β>0 be the order of the derivative consider in Caputo sense, then

    CDβΨ(z)Qβ(n+1,n+1)Ψ(z), (4.14)

    where Qβ is the (n+1)×(n+1) operational matrix which can be computed as

    Qβ(n+1,n+1)=jq=βΩa(q,j)c(r,i)Γ(qΩ+1)(qΩβ+1),j=βΩ,βΩ+1,. (4.15)

    Proof. Using (3.11) and linearity of Caputo fractional derivative, we have

    CDβΨ(z)=jq=0a(q,j)CDβ(zqΩ) (4.16)

    Using lemma (4.1), Eq (4.16) can be expressed as

    CDβΨ(z)=jq=βΩa(q,j)Γ(qΩ+1)Γ(qΩβ+1)zqΩβ. (4.17)

    Now Eq (4.17) and lemma (4.4) proves the required result.

    Remark 4.8. The operational matrix of derivatives for CSLPs studied by Saadatmandi and Dehghan [13] is a special case of Theorem 4.7 for Ω=1.

    In this section, we are interested in the development of the computational algorithm to solve numerically the following generalized multi-order FDEs

    CDαu(z)=g(z,u(z),CDβ1u(z),CDβ2u(z),,CDβnu(z)),u(s)(0)=hs,s=0,1,,n. (5.1)

    where n1<αn, 0<β1<β2<<βn<α, and CDα is a Caputo derivative of order α defined in (2.2).

    We are interested to approximate the solution in terms of orthogonal FLFVs, such that

    CDαu(z)=ΛTΨ(z). (5.2)

    Now integrating (5.2) in Riemann-Liouville sense and using the initial conditions u(s)(0)=hs,s=0,1,,n, we have

    u(z)=ΛTJαΨ(z)+ns=0hszs. (5.3)

    Using Theorem 4.5, Eq (5.3) can also be written as

    u(z)=ΛTDα(n+1,n+1)Ψ(z)+ns=0hszs. (5.4)

    The terms in ns=0hszs can also be approximated using FLFVs. So, Eq (5.4) can also be written as

    u(z)=ΛTDα(n+1,n+1)Ψ(z)+HTΨ(z). (5.5)

    Now using Theorem 4.7 and Eq (5.5), the right hand side terms of Eq (5.1) can be approximated as

    {CDβ0u(z)=ΛTDα(n+1,n+1)Qβ0(n+1,n+1)Ψ(z)+HTQβ0(n+1,n+1)Ψ(z),CDβ1u(z)=ΛTDα(n+1,n+1)Qβ1(n+1,n+1)Ψ(z)+HTQβ1(n+1,n+1)Ψ(z),CDβnu(z)=ΛTDα(n+1,n+1)Qβn(n+1,n+1)Ψ(z)+HTQβn(n+1,n+1)Ψ(z),andg(z)=GTΨ(z). (5.6)

    Now inserting (5.2) and (5.6) in (5.1), we have

    ΛTΨ(z)=ΛTDα(n+1,n+1)(Qβ0(n+1,n+1)+Qβ1(n+1,n+1)++Qβn(n+1,n+1))Ψ(z)+HT(Qβ0(n+1,n+1)+Qβ1(n+1,n+1)++Qβn(n+1,n+1))Ψ(z)+GTΨ(z).

    After simplification, we have

    ΛTΛTDα(n+1,n+1)ˆQHTˆQGT, (5.7)

    where ˆQ=(Qβ0(n+1,n+1)+Qβ1(n+1,n+1)++Qβn(n+1,n+1)).

    Equation (5.7) is a matrix equation of Sylvester type that can be easily solved for the unknown ΛT using any computational software. By inserting the values of ΛT in Eq (5.5), we get the approximate solution of (5.1).

    In this section, the applicability of the method is studied by solving various examples and comparing their analytical solutions with the approximate solutions obtained using our developed computational numerical algorithm.

    Example 1. We consider the following non-homogeneous Bagley-Torvik equation [13,19,28]

    CD2u(z)+CD32u(z)+u(z)=1+z, (6.1)

    subject to the following initial conditions

    u(0)=1=u(0). (6.2)

    The exact solution is given as under

    u(z)=1+z. (6.3)

    Here, we have

    GT=[32,12,0],HT=[32,12,0],ΛT=[0,0,0]. (6.4)

    and

    Q32=[0000009.02705.41621.2896],D2=[0.16670.25000.08330.08330.100000.016700.0238]. (6.5)

    Using (6.4) and (6.5) in (5.5), we get the exact solution of the problem (6.1)-(6.2).

    Remark 6.1. The results of Example 1 at various values of n obtained by using our method are compared with the results obtained by using the method [4] at n=10 in Tables 1-4 and Figure 5. We observe that at various values of n (the number of terms of FLFVs), the exact solution and the approximate solution obtained by using our method coincide, (see Table 2 and Figure 5). However, the results obtained by using the method in [4] yield less accurate results, (see Table 2-4 and Figure 5).

    Table 1.  Comparison of the approximate solution (AS) for Example 1 obtained by using our proposed method (PM) with the AS obtained by using the method ([4], Example 2).
    z Exact Solution AS using PM at n=2 AS in [4] at n=10 (number of neurons)
    0 1.00 1.00 1.024862
    0.1 1.10 1.10 1.121206
    0.2 1.20 1.20 1.220821
    0.3 1.30 1.30 1.323041
    0.4 1.40 1.40 1.426952
    0.5 1.50 1.50 1.531330
    0.6 1.60 1.60 1.634569
    0.7 1.70 1.70 1.734591
    0.8 1.80 1.80 1.828738
    0.9 1.90 1.90 1.913640
    1.0 2.00 2.00 1.985057

     | Show Table
    DownLoad: CSV
    Table 2.  Comparison of the Absolute errors for Example 1 obtained by using our proposed method (PM) at n=2 (number of terms of FLFVs) with the absolute errors obtained by using the method ([4], Example 2) at n=10 (number of neurons).
    z Absolute error using [4] Absolute error using PM
    0 2.30e2 0
    0.1 2.69e2 0
    0.2 3.13e2 0
    0.3 3.45e2 0
    0.4 3.45e2 0
    0.5 2.87e2 0
    0.6 1.36e2 0
    0.7 1.49e2 0
    0.8 2.30e2 0
    0.9 2.69e2 0
    1.0 3.13e2 0

     | Show Table
    DownLoad: CSV
    Table 3.  Comparison of the Absolute errors for Example 1 obtained by using our proposed method (PM) at n=8 (number of terms of FLFVs) with the absolute errors obtained by using the method ([4], Example 2) at n=10 (number of neurons).
    z Absolute error using [4] Absolute error using PM
    0 2.30e2 0
    0.1 2.69e2 0
    0.2 3.13e2 0
    0.3 3.45e2 0
    0.4 3.45e2 0
    0.5 2.87e2 0
    0.6 1.36e2 0
    0.7 1.49e2 0
    0.8 2.30e2 0
    0.9 2.69e2 0
    1.0 3.13e2 0

     | Show Table
    DownLoad: CSV
    Table 4.  Comparison of the Absolute errors for Example 1 obtained by using our proposed method (PM) at n=10 (number of terms of FLFVs) with the absolute errors obtained by using the method ([4], Example 2) at n=10 (number of neurons).
    z Absolute error using [4] Absolute error using PM
    0 2.30e2 0
    0.1 2.69e2 0
    0.2 3.13e2 0
    0.3 3.45e2 0
    0.4 3.45e2 0
    0.5 2.87e2 0
    0.6 1.36e2 0
    0.7 1.49e2 0
    0.8 2.30e2 0
    0.9 2.69e2 0
    1.0 3.13e2 0

     | Show Table
    DownLoad: CSV
    Figure 1.  FLFVs are plotted for Ω=1/2 and for various choices of n.
    Figure 2.  FLFVs are plotted taking n=4 and for various choices of Ω.
    Figure 3.  Exact and approximated graphs of the function u(z)=1+z using FLFVs at Ω=1/2 , and n=4.
    Figure 4.  Exact and approximated graphs of the function u(z)=sin(z)+1+z2 using FLFVs at Ω=1/2, and n=4.
    Figure 5.  Comparison of the approximate solution of Example 1 obtained using the proposed method at n=2 with the approximate solution obtained using the method of [4] at n=10.

    Remark 6.2. In Tables 5 and 6, the results of Example 1 obtained by using our method are compared with the results obtained by using the methods in [5] and [6]. Our proposed method provides the exact solution at n=2 (first three terms of FLFVs). However, the methods in [5] and [6] provide the exact solution at n=6 (first seven terms of Taylor series, and first seven terms of Bessel series). This shows that our proposed method is numerically more efficient.

    Table 5.  Comparison of the approximate solution (AS) for Example 1 obtained by using our proposed method (PM) with the AS obtained by using the method ([5], Example 1).
    z Exact solution AS using PM at n=2 AS in [5] at n=6 (number of terms of Taylor series)
    0 1.00 1.00 1.00
    0.1 1.10 1.10 1.10
    0.2 1.20 1.20 1.20
    0.3 1.30 1.30 1.30
    0.4 1.40 1.40 1.40
    0.5 1.50 1.50 1.50
    0.6 1.60 1.60 1.60
    0.7 1.70 1.70 1.70
    0.8 1.80 1.80 1.80
    0.9 1.90 1.90 1.90
    1.0 2.00 2.00 2.00

     | Show Table
    DownLoad: CSV
    Table 6.  Comparison of the approximate solution (AS) for Example 1 obtained by using our proposed method (PM) with the AS obtained by using the method ([6], Example 1).
    z Exact solution AS using PM at n=2 AS in [6] at n=6 (number of terms of Bessel series)
    0 1.00 1.00 1.00
    0.1 1.10 1.10 1.10
    0.2 1.20 1.20 1.20
    0.3 1.30 1.30 1.30
    0.4 1.40 1.40 1.40
    0.5 1.50 1.50 1.50
    0.6 1.60 1.60 1.60
    0.7 1.70 1.70 1.70
    0.8 1.80 1.80 1.80
    0.9 1.90 1.90 1.90
    1.0 2.00 2.00 2.00

     | Show Table
    DownLoad: CSV

    Example 2. Consider the following non-homogeneous multi-order fractional problem [21,32]

    CDαu(z)=aCDβ0u(z)+bCDβ1u(z)+cCDβ2u(z)+dCDβ3u(z)+G(z),z[0,1],0<α<2, (6.6)

    subject to the following initial conditions

    u(0)=0,u(0)=0. (6.7)

    The source term is as under

    G(z)=4zz26776z324503+42z514z6+z7+1516z13256292.

    The exact solution corresponding to α=2,a=c=1,b=2,d=0,β0=0,β1=1,β2=12 is given below

    u(z)=z7z2.

    We study the Example 2 by comparing the approximate results obtained using our algorithm with the approximate results achieved using the methods in [21,32]. We observe that our method produces better approximate results and higher precision in the approximate solution, see Table 7, Figure 6-9.

    Table 7.  Comparison of absolute error of Example 2 computed using our proposed method (PM) and the absolute error computed with [21] at n=4,6,7.
    z n=4 [21] n=4 PM n=6 [21] n=6 PM n=7 [21] n=7 PM
    0.2 0.0844 0.0132 0.0044 0.00003 2.81025203108243e-15 0
    0.4 0.3501 0.0105 0.0079 0.00035 6.63358257213531e-15 0
    0.6 0.6734 0.0349 0.0143 0.00034 3.27515792264421e-15 0
    0.8 1.0234 0.0387 0.0214 0.00069 4.25770529943748e-14 0
    1 1.6700 0.0780 0.0280 0.00103 2.43819897540083e-13 0

     | Show Table
    DownLoad: CSV
    Figure 6.  Graphs of exact solution and approximate solution of Example 2 at Ω=1,α=2 and n=4.
    Figure 7.  Graphs of exact solution and approximate solution of Example 2 at Ω=1,α=2 and n=5.
    Figure 8.  Graphs of exact solution and approximate solution of Example 2 at Ω=1,α=2 and n=6.
    Figure 9.  Comparison of the error plots of Example 2 with the results of [21] at n=4,6,7.

    The applicability of the numerical algorithm is tested at various scale levels. We observe that with the increase in scale level, the approximate solution and the exact solution have great resemblance, see Figure 7. Also, at n=7, the approximate solution is equal to the exact solution, u(z)=z7z2 of the problem (6.6)-(6.7), see Table (7, Column 7]. We also calculate the amount of the absolute errors at scale level, n=4,6,7. we analyze that with the increase in scale level, the amount of the absolute errors decreases significantly, see Table 7 and Figure 9. It is worth to mention that by taking few terms of orthogonal FLFVs, the good match is obtained between the approximate solution and the exact solution. For example compare the results of ([21], Figures 4-5 and Table 1) and ([32], Figure 2) with the results of our paper (Figures 6-8 and Table (7, Column 7).

    Example 3. Consider the following non-homogeneous multi-order fractional problem [21]

    CDαu(z)=aCDβ0u(z)+bCDβ1u(z)+cCDβ2u(z)+dCDβ3u(z)+G(z),z[0,1],0<α<2, (6.8)

    subject to the following initial conditions

    u(0)=0,u(0)=0. (6.9)

    The source term is as under

    G(z)=6z+z312Γ(73)z43+6Γ(103)z73.

    The exact solution corresponding to α=2,a=c=1,b=0,d=2,β0=0,β2=23(0,1), and β3=53(1,2) is given below

    u(z)=z3.

    Here for n=3, we have

    {GT=[0.42150.68670.02160.1984],HT=[0000],ΛT=[2.99913.00000.00000.0000]. (6.10)

    and

    {Q23=[00001.67960.71990.23910.12920.71993.02361.38440.50391.43980.66453.95941.9994],Q53=[0000000010.07864.31941.43980.77537.199030.235813.84425.0393],D2=[0.16670.25000.083300.08330.100000.01670.016700.0238000.007100.0111]. (6.11)

    Using (6.10) and (6.11) in (5.5), we get the exact solution, u(z)=z3 of the problem (6.8) and (6.9).

    We study the Example 3 by comparing the approximate results obtained using our algorithm with the approximate results achieved using the methods in [21]. We observe that our method produces better approximate results and higher precision in the approximate solution, see Table 8, Figures 10 and 11.

    Table 8.  Comparison of absolute error of Example 3 computed using our proposed method (PM) and the absolute error computed with [21] at n=2,3.
    z Exact solution n=2 [21] n=2 PM n=3 [21] n=3 PM
    0.2 0.0080 0.0614 0.0291 2.55004350968591e-16 0
    0.4 0.0640 0.1541 0.0885 2.77555756156289e-16 0
    0.6 0.2160 0.3339 0.1239 3.33066907387547e-16 0
    0.8 0.5120 0.6488 0.0875 4.44089209850063e-16 0
    1 1.0000 1.1468 0.0689 6.66133814775094e-16 0

     | Show Table
    DownLoad: CSV
    Figure 10.  Graphs of exact solution and approximate solution of Example 3 at Ω=1,α=2 and n=3.
    Figure 11.  Comparison of the error plots of Example 3 with the results of [21] at n=2,3.

    We observe that by increasing the scale level, the absolute error amount decreases significantly, see Table 8. Also at low scale level, n=3, the exact solution of the problem, (6.8)-(6.9) is obtained which shows that our proposed algorithm is more efficient as compared to [21], see Figures 10 and 11 and Table 8.

    Example 4. Consider the following fractional problem [19,30]

    CDαu(z)=u(z)+G(z),z[0,1],0<α1, (6.12)

    subject to the following initial condition

    u(0)=0. (6.13)

    The source term is as under

    G(z)=2z2αΓ(3α)z1αΓ(2α)+z2z.

    The exact solution of the problem (6.12)-(6.13) is as under

    u(z)=z2z.

    We study the Example 4 to investigate the applicability of orthogonal FLFVs at various fractional values of Ω, and α, see Figures 12-14. We observe that at various choices of α=0.5,0.7,0.9, and Ω=12,13, the approximate results have a great compatibility with the exact solution. Also at n=2(r+1), the exact solution of the problem (6.12)-(6.13) can be obtained for Ω=1r+1>α2,rN. For instance, by inserting (6.14) in (5.5), we can obtain the exact solution, u(z)=z2z of (6.12)-(6.13), see Figure 14.

    Figure 12.  Graphs of exact solution and approximate solution of Example 4 at Ω=13,n=4, and various values of α.
    Figure 13.  Graphs of exact solution and approximate solution of Example 4 at Ω=12,n=4 , and various values of α.
    Figure 14.  Graphs of exact solution and approximate solution of Example 4 at Ω=α=12,n=4.

    Here for n=4,Ω=12, and α=12, we have

    {Ψ(z)=[12z16z6z+112z30z+20z32190z+70z220z140z32+1],D1/2=[0.56420.56420000.02660.32200.2954000.07990.06340.24210.225700.002490.08490.07130.20010.18990.00800.00830.07770.06940.1752],GT=[47997861414693660616975379209941493415679873386578213105753065538344692221904871836977163678129868313597467731564262115061310111169354070572322450563288195426556510315139133741760548032570170],[ΛT=0.188060.11280.37610.07520.0000]. (6.14)

    Example 5. Consider the following fractional problem [19,31]

    CDαu(z)u(z)=1,z[0,1],0<α1, (6.15)

    subject to the following initial condition

    u(0)=0. (6.16)

    The exact solution of the problem (6.15)-(6.16) is as under

    u(z)=j=1zαjΓ(αj+1).

    We study the Example 5 to investigate the applicability of orthogonal FLFVs at various fractional values of Ω, and α, see Figures 15 and 16. We observe that for various choices of Ω=0.4,0.6,0.8,1, and α=0.9, there is a good compatibility between the exact solution and the approximate solution. Also, we get improve approximate results with the increase in scale level, see Figure 16. We also determine the amount of the absolute error at scale level, n=10 and for various choices of Ω to investigate the efficiency and accuracy of the proposed algorithm for fractional choices of Ω, see Table 9. We also analyze the effect of the fractional parameter, α by comparing the approximate solution with the exact solution, see Figure 17.

    Figure 15.  Graphs of exact solution and approximate solution of Example 5 at α=0.9,n=3 and various choices of Ω.
    Figure 16.  Graphs of exact solution and approximate solution of Example 5 at α=0.9,n=4 and various choices of Ω.
    Table 9.  Absolute error of Example 5 at n=10 and various choices of Ω=α.
    z Ω=α=34 Ω=α=12 Ω=α=23
    0.2 4.273327e-11 5.870078e-07 9.236112e-10
    0.4 8.335039e-09 2.999352e-05 1.401714e-07
    0.6 2.455070e-07 3.010034e-04 2.873981e-06
    0.8 2.731391e-06 1.568348e-03 2.478263e-05
    1 1.782081e-05 5.703710e-03 1.328454e-04

     | Show Table
    DownLoad: CSV
    Figure 17.  Graphs of exact solution and approximate solution of Example 5 at n=7,Ω=1 and various choices of α.

    Example 6. Consider the following fractional problem [13]

    CDαu(z)=u(z),z[0,1],0<α2, (6.17)

    subject to the following initial conditions

    u(0)=1,u(0)=0. (6.18)

    The condition, u(0)=0 is applicable only when α>1.

    The exact solution of the problem (6.17)-(6.18) is as under

    u(z)=j=0(zα)jΓ(αj+1).

    We study Example 6 to check the applicability of our proposed technique when 0<α<1, and 1<α<2, see Figures 18-23, Tables 10 and 11. We also investigate the accuracy and validity of the FLFVs for various choices of α=Ω, see Tables 10 and 11. We also observe the effect of the fractional parameter α by comparing the exact solution with the approximate solution, see Figures 21 and 22. We observe that with the increase in scale level n and α2, the approximate solution of the problem (6.17)-(6.18) approaches to the solution of classical order differential equation, see Figures 21 and 22. We also observe that when 0<α<1 then at very low scale level n=2,3, the approximate solution shows great resemblance with the exact solution, see Figures 18-19. We also determine the absolute errors at scale level n=10, and for various choices of α=Ω. We observe that our technique is more efficient as compared to the method used in [13], see Table 11 and Figure 23.

    Figure 18.  Graphs of exact solution and approximate solution of Example 6 at Ω=1,α=1 and n=2.
    Figure 19.  Graphs of exact solution and approximate solution of Example 6 at Ω=1,α=1 and n=3.
    Figure 20.  Error plot of Example 6 at Ω=1,α=1 and n=3.
    Figure 21.  Graphs of exact solution and approximate solution of Example 6 at Ω=1,n=3, and various values of α.
    Figure 22.  Graphs of exact solution and approximate solution of Example 6 at Ω=0.95,n=4, and various values of α.
    Figure 23.  Comparison of the error plots of Example 6 with the results of [13] at n=10, 0<α<2 and Ω=1.
    Table 10.  Absolute error of Example 6 at n=10 and various choices of Ω=α.
    z Ω=α=12 Ω=α=34 Ω=α=23 Ω=α=54 Ω=α=74 Ω=α=32
    0.2 4.212086e-07 2.509733e-11 6.980854e-10 8.742474e-17 2.671553e-18 1.176754e-17
    0.4 1.792713e-05 6.879759e-09 1.070398e-07 4.165262e-16 9.565113e-19 1.463810e-16
    0.6 1.593986e-04 1.891806e-07 2.018518e-06 2.007704e-14 3.257643e-18 3.281789e-16
    0.8 7.478411e-04 1.975211e-06 1.612752e-05 1.013113e-12 4.161064e-18 7.977363e-16
    1 2.473455e-03 1.213707e-05 8.053908e-05 2.159538e-11 2.954844e-17 1.218018e-14

     | Show Table
    DownLoad: CSV
    Table 11.  Comparison of absolute error of Example 6 computed using our proposed method (PM) with the absolute error computed using [13] at n=10,Ω=1 and at various choices of α.
    α z=0.1 z=0.1PM z=0.3 z=0.3PM z=0.5 z=0.5PM z=0.7 z=0.7PM z=0.9 z=0.9PM
    1.2 3.1e-03 2.2e-05 2.8e-03 1.9e-05 4.5e-03 1.2e-05 3.6e-03 1.6e-05 1.8e-03 2.0e-05
    1.4 1.0e-03 1.6e-05 7.0e-04 1.4e-05 1.3e-03 7.9e-06 1.1e-03 1.1e-05 2.4e-04 1.4e-05
    1.6 3.0e-04 7.5e-06 1.3e-04 6.3e-06 3.1e-04 3.4e-06 3.0e-04 4.9e-06 6.2e-07 6.1e-06
    1.8 6.1e-05 2.2e-06 1.4e-05 1.8e-06 4.9e-05 9.2e-07 5.3e-05 1.4e-06 8.8e-06 1.7e-06
    0.2 2.9e-01 2.3e-03 4.5e-01 1.9e-02 7.4e-01 5.3e-02 3.7e-01 1.0e-01 2.0e-01 1.8e-01
    0.4 3.9e-01 7.0e-04 5.1e-01 1.0e-03 7.3e-01 1.5e-03 3.3e-01 2.4e-03 2.2e-01 8.2e-03
    0.6 6.7e-03 3.8e-04 2.0e-05 4.4e-04 5.2e-03 3.32e-04 4.4e-03 3.42e-04 4.6e-03 3.41e-04
    0.8 1.1e-03 1.1e-04 2.1e-04 1.14e-04 8.4e-04 7.8e-05 8.7e-04 9.5e-05 5.8e-04 1.2e-04

     | Show Table
    DownLoad: CSV

    We proposed a numerical method which is completely based on the fractional derivative and integral operational matrices of FLFVs. The proposed method is independent of the choice of the suitable collocation points and expansion of the residual function as a series of orthogonal polynomials. Consequently, the proposed method produced high efficient numerical results as compared to the Spectral tau method [13], Bessel collocation method [6], Taylor matrix method [5], function approximation theory [21], and stochastic technique [4]. The other novel aspect of our article is the development of new integral and derivative operational matrices in Riemann-Liouville and Caputo senses respectively. As an applications of the proposed method, we solved various MOFDEs corresponding to initial conditions. Our approach has ability to reduce the MOFDEs into a system of Sylvester types matrix equations which can be solved using MATLAB built in function lyap(.). The proposed results are not sufficient to solve numerically the fractional boundary value problems (BVPs). We are interested in extending the developed numerical results for fractional BVPs as well to partial FDEs by developing some more operational matrices.

    We are grateful to the reviewers for their comments and suggestions which improve the quality of article. The research is supported by Cankaya University, Department of Mathematics, Turkey.

    The authors declare no conflict of interest.



    [1] F. Al-Aboudi, On univalent functions defined by a generalized Sălăgean operator, Int. J. Math. Math. Sci., 27 (2004), 481–494. doi: 10.1155/S0161171204108090. doi: 10.1155/S0161171204108090
    [2] B. Ahmad, M. G. Khan, B. A. Frasin, M. K. Aouf, T. Abdeljawad, W. K. Mashwani, et al., On q-analogue of meromorphic multivalent functions in lemniscate of Bernoulli domain, AIMS Math., 6 (2021), 3037–3052. doi: 10.3934/math.2021185. doi: 10.3934/math.2021185
    [3] Q. Z. Ahmad, N. Khan, M. Raza, M. Tahir, B. Khan, Certain q-difference operators and their applications to the subclass of meromorphic q-starlike functions, Filomat, 33 (2019), 3385–3397. doi: 10.2298/FIL1911385A. doi: 10.2298/FIL1911385A
    [4] M. K. Aouf, A generalization of functions with real part bounded in the mean on the unit disc, Math. Japonica., 33 (1988), 175–182.
    [5] M. K. Aouf, Generalization of certain subclasses of multivalent functions with negative coefficients defined by using a differential operator, Math. Comput. Model., 50 (2009), 1367–1378. doi: 10.1016/j.mcm.2008.08.026. doi: 10.1016/j.mcm.2008.08.026
    [6] M. K. Aouf, On certain multivalent functions with negative coefficients defined by using a differential operator, Mat. Vesn., 62 (2010), 23–35. doi: 0025-51651001023A.
    [7] M. K. Aouf, S. M. Madian, Inclusion and properties neighborhood for certain p-valent functions associated with complex order and q-p-valent Cătaş operator, J. Taibah Univ. Sci., 14 (2020), 1226–1232. doi: 10.1080/16583655.2020.1812923. doi: 10.1080/16583655.2020.1812923
    [8] M. K. Aouf, S. M. Madian, Coefficient bounds for bi-univalent classes defined by Bazilevič functions and convolution, Bol. Soc. Mat. Mex., 26 (2020), 1045–1062. doi: 10.1007/s40590-020-00304-0. doi: 10.1007/s40590-020-00304-0
    [9] M. K. Aouf, S. M. Madian, Certain classes of analytic functions associated with q-analogue of p-valent Cătaş operator, Moroccan J. Pure Appl. Anal., 7 (2021), 430–447. doi: 10.2478/mjpaa-2021-0029. doi: 10.2478/mjpaa-2021-0029
    [10] M. K. Aouf, S. M. Madian, Subordination factor sequence results for starlike and convex classes defined by q-Cătaş operator, Afr. Mat., 2021, 1–13. doi: 10.1007/s13370-021-00896-4. doi: 10.1007/s13370-021-00896-4
    [11] M. K. Aouf, S. M. Madian, A. O. Mostafa, Bi-univalent properties for certain class of Bazilevič functions defined by convolution and with bounded boundary rotation, J. Egypt. Math. Soc., 27 (2019), 11. doi: 10.1186/s42787-019-0012-2. doi: 10.1186/s42787-019-0012-2
    [12] I. E. Bazilevič, On a case of integrability in quadratures of the Lowner-Kufarev equation, Mat. Sb., 37 (1955), 471–476.
    [13] M. Çaǧlar, E. Deniz, H. M. Srivastava, Second Hankel determinant for certain subclasses of bi-univalent functions, Turk. J. Math., 41 (2017), 694–706. doi: 10.3906/mat-1602-25
    [14] A. Cătaş, On certain classes of p-valent functions defined by multiplier transformations, Proceedings of the international symposium on geometric function theory and applications: GFTA, 2007,241–250.
    [15] P. L. Duren, Univalent functions, New York: Springer Verlag, 1983.
    [16] P. Goswami, B. S. Alkahtani, T. Bulboaca, Estimate for initial Maclaurin coefficients of certain subclasses of bi-univalent functions, arXiv. Available from: https://arXiv.org/abs/1503.04644.
    [17] Q. Hu, H. M. Srivastava, B. Ahmad, N. Khan, M. G. Khan, W. K. Mashwani, et al., A subclass of multivalent Janowski type q-starlike functions and its consequences, Symmetry, 13 (2021), 1275. doi: 10.3390/sym13071275. doi: 10.3390/sym13071275
    [18] Q. H. Xu, H. G. Xiao, H. M. Srivastava, A certain general subclass of analytic and bi-univalent functions and associated coefficient estimate problems, Appl. Math. Comput., 218 (2012), 11461–11465. doi: 10.1016/j.amc.2012.05.034. doi: 10.1016/j.amc.2012.05.034
    [19] Q. Khan, M. Arif, M. Raza, G. Srivastava, H. Tang, S. ur Rehman, Some applications of a new integral operator in q-analog for multivalent functions, Mathematics, 7 (2019), 1178. doi: 10.3390/math7121178. doi: 10.3390/math7121178
    [20] B. Khan, Z. G. Liu, H. M. Srivastava, N. Khan, M. Darus, M. Tahir, A study of some families of multivalent q-starlike functions involving higher-order q-derivatives, Mathematics, 8 (2020), 1470. doi: 10.3390/math8091470. doi: 10.3390/math8091470
    [21] B. Khan, Z. G. Liu, H. M. Srivastava, N. Khan, M. Tahir, Applications of higher-order derivatives to subclasses of multivalent q-starlike functions, Maejo Int. J. Sci. Technol., 15 (2021), 61–72.
    [22] N. Khan, M. Shafiq, M. Darus, B. Khan, Q. Z. Ahmad, Upper bound of the third Hankel determinant for a subclass of q-starlike functions associated with Lemniscate of Bernoulli, J. Math. Inequal., 14 (2020), 51–63. doi: 10.7153/jmi-2020-14-05. doi: 10.7153/jmi-2020-14-05
    [23] B. Khan, H. M. Srivastava, N. Khan, M. Darus, Q. Z. Ahmad, M. Tahir, Applications of certain conic domains to a subclass of q-starlike functions associated with the Janowski functions, Symmetry, 13 (2021), 574. doi: 10.3390/sym13040574. doi: 10.3390/sym13040574
    [24] B. Khan, H. M. Srivastava, N. Khan, M. Darus, M. Tahir, Q. Z. Ahmad, Coefficient estimates for a subclass of analytic functions associated with a certain leaf-like domain, Mathematics, 8 (2020), 1334. doi: 10.3390/math8081334. doi: 10.3390/math8081334
    [25] B. Khan, H. M. Srivastava, M. Tahir, M. Darus, Q. Z. Ahmad, N. Khan, Applications of a certain q-integral operator to the subclasses of analytic and bi-univalent functions, AIMS Math., 6 (2020), 1024–1039. doi: 10.3934/math.2021061. doi: 10.3934/math.2021061
    [26] W. Ma, D. Minda, A unified treatment of some special classes of univalent functions, In: Z. Li, F. Ren, L. Lang, S. Zhang, Proceedings of the conference on complex analysis, International Press Inc., 1 (1994), 157–169.
    [27] N. Magesh, T. Rosy, S. Varma, Coefficient estimate problem for a new subclass of bi-univalent functions, J. Compl. Anal., 2013 (2013), 474231. doi: 10.1155/2013/474231. doi: 10.1155/2013/474231
    [28] H. Orhan, M. Kamali, Neighborhoods of a class of analytic functions with negative coefficients, Acta Math. Acad. Paedagog. Nyhazi., 21 (2005), 55–61.
    [29] H. Orhan, N. Magesh, V. K. Balaji, Certain classes of bi-univalent functions with bounded boundary variation, Tbilisi Math. J., 10 (2017), 17–27. doi: 10.1515/tmj-2017-0042. doi: 10.1515/tmj-2017-0042
    [30] H. Orhan, N. Magesh, V. K. Balaji, Initial coefficient bounds for a general class of bi-univalent functions, Filomat, 29 (2015), 1259–1267. doi: 10.2298/FIL1506259O. doi: 10.2298/FIL1506259O
    [31] K. S. Padmanabhan, R. Paravatham, Properties of a class of functions with bounded boundary rotation, Ann. Pol. Math., 31 (1976), 311–323. doi: 10.4064/ap-31-3-311-323
    [32] B. Pinchuk, Functions of bounded boundary rotation, Israel J. Math., 10 (1971), 6–16. doi: 10.1007/BF02771515. doi: 10.1007/BF02771515
    [33] S. Prema, B. S. Keerthi, Coefficient bounds for certain subclasses of analytic functions, J. Math. Anal., 4 (2013), 22–27.
    [34] M. S. Robertson, Variational formulas for several classes of analytic functions, Math. Z., 118 (1976), 311–319. doi: 10.1007/BF01109867. doi: 10.1007/BF01109867
    [35] G. S. Sălăgean, Subclasses of univalent functions, In: C. A. Cazacu, N. Boboc, M. Jurchescu, I. Suciu, Complex analysis-fifth Romanian-finnish seminar, Lecture Notes in Mathematics, Springer-Verlag, 1013 (1983), 362–372. doi: 10.1007/BFb0066543.
    [36] R. Singh, On Bazilevič functions, Proc. Am. Math. Soc., 38 (1973), 261–271. doi: 10.1090/S0002-9939-1973-0311887-9. doi: 10.1090/S0002-9939-1973-0311887-9
    [37] H. M. Srivastava, Operators of basic (or q-) calculus and fractional q-calculus and their applications in geometric function theory of complex analysis, Iran. J. Sci. Technol. Trans. Sci., 44 (2020), 327–344. doi: 10.1007/s40995-019-00815-0. doi: 10.1007/s40995-019-00815-0
    [38] H. M. Srivastava, Univalent functions, fractional calculus, and associated generalized hypergeometric functions, In: Univalent functions, fractional calculus and their applications, Halsted Press, 1989,329–354.
    [39] H. M. Srivastava, M. Arif, M. Raza, Convolution properties of meromorphically harmonic functions defined by a generalized convolution q-derivative operator, AIMS Math., 6 (2021), 5869–5885. doi: 10.3934/math.2021347. doi: 10.3934/math.2021347
    [40] H. M. Srivastava, D. Bansal, Coefficient estimates for a subclass of analytic and bi-univalent functions, J. Egypt. Math. Soc., 23 (2015), 242–246. doi: 10.1016/j.joems.2014.04.002. doi: 10.1016/j.joems.2014.04.002
    [41] H. M. Srivastava, S. S. Eker, R. M. Ali, Coefficient estimates for a certain class of analytic and bi-univalent functions, Filomat, 29 (2015), 1839–1845. doi: 10.2298/FIL1508839S. doi: 10.2298/FIL1508839S
    [42] H. M. Srivastava, S. Gaboury, F. Ghanim, Initial coefficient estimates for some subclasses of m-fold symmetric bi-univalent functions, Acta Math. Sci., 36 (2016), 863–871. doi: 10.1016/S0252-9602(16)30045-5. doi: 10.1016/S0252-9602(16)30045-5
    [43] H. M. Srivastava, S. Gaboury, F. Ghanim, Coefficient estimates for some general subclasses of analytic and bi-univalent functions, Afr. Mat., 28 (2017), 693–706. doi: 10.1007/s13370-016-0478-0. doi: 10.1007/s13370-016-0478-0
    [44] H. M. Srivastava, S. B. Joshi, S. S. Joshi, H. Pawar, Coefficient estimates for certain subclasses of meromorphically bi-univalent functions, Palest. J. Math., 5 (2016), 250–258.
    [45] H. M. Srivastava, A. K. Mishra, P. Gochhayat, Certain subclasses of analytic and bi-univalent functions, Appl. Math. Lett., 23 (2010), 1188–1192. doi: 10.1016/j.aml.2010.05.009. doi: 10.1016/j.aml.2010.05.009
    [46] H. M. Srivastava, A. O. Mostafa, M. K. Aouf, H. M. Zayed, Basic and fractional q-calculus and associated Fekete-Sezgö problem for p-valent q-starlike functions and p-valent q-convex functions of complex order, Miskolc Math. Notes, 20 (2019), 489–509. doi: 10.18514/MMN.2019.2405. doi: 10.18514/MMN.2019.2405
    [47] H. M. Srivastava, M. Tahir, B. Khan, Q. Z. Ahmad, N. Khan, Some general families of q-starlike functions associated with the Janowski functions, Filomat, 33 (2019), 2613–2626. doi: 10.2298/FIL1909613S. doi: 10.2298/FIL1909613S
    [48] H. M. Srivastava, M. Tahir, B. Khan, Q. Z. Ahmad, N. Khan, Some general classes of q-starlike functions Aassociated with the Janowski functions, Symmetry, 11 (2019), 292. doi: 10.3390/sym11020292. doi: 10.3390/sym11020292
    [49] H. M. Srivastava, B. Khan, N. Khan, M. Tahir, S. Ahmad, N. Khan, Upper bound of the third Hankel determinant for a subclass of q-starlike functions associated with the q-exponential function, Bull. Sci. Math., 167 (2021), 102942. doi: 10.1016/j.bulsci.2020.102942. doi: 10.1016/j.bulsci.2020.102942
    [50] H. Tang, H. M. Srivastava, S. Sivasubramanian, P. Gurusamy, The Fekete-Szegö functional problems for some subclasses of m -fold symmetric bi-univalent functions, J. Math. Ineq., 10 (2016), 1063–1092. doi: 10.7153/jmi-10-85. doi: 10.7153/jmi-10-85
    [51] D. K. Thomas, On Bazilevič functions, Am. Math. Soc., 132 (1968), 353–361.
    [52] D. K. Thomas, On the coefficients of gamma-starlike functions, J. Korean Math. Soc., 55 (2018), 175–184. doi: 10.4134/JKMS.j170105. doi: 10.4134/JKMS.j170105
  • This article has been cited by:

    1. Jacques Demongeot, Cécile Fougère, mRNA COVID-19 Vaccines—Facts and Hypotheses on Fragmentation and Encapsulation, 2022, 11, 2076-393X, 40, 10.3390/vaccines11010040
    2. George Engelhard, Functional Approaches for Modeling Unfolding Data, 2023, 0013-1644, 001316442211434, 10.1177/00131644221143474
    3. Kayode Oshinubi, Sana S. Buhamra, Noriah M. Al-Kandari, Jules Waku, Mustapha Rachdi, Jacques Demongeot, Age Dependent Epidemic Modeling of COVID-19 Outbreak in Kuwait, France, and Cameroon, 2022, 10, 2227-9032, 482, 10.3390/healthcare10030482
    4. Jong-Min Kim, Copula Dynamic Conditional Correlation and Functional Principal Component Analysis of COVID-19 Mortality in the United States, 2022, 11, 2075-1680, 619, 10.3390/axioms11110619
    5. Paula Vicente, Abdul Suleman, COVID-19 in Europe: from outbreak to vaccination, 2022, 22, 1471-2458, 10.1186/s12889-022-14454-5
    6. Ibrahim M. Almanjahie, Zoulikha Kaid, Ali Laksaci, Mustapha Rachdi, Estimating the Conditional Density in Scalar-On-Function Regression Structure: k-N-N Local Linear Approach, 2022, 10, 2227-7390, 902, 10.3390/math10060902
    7. Tobias Wistuba, Andreas Mayr, Christian Staerk, Estimating the course of the COVID-19 pandemic in Germany via spline-based hierarchical modelling of death counts, 2022, 12, 2045-2322, 10.1038/s41598-022-13723-y
    8. Antonio Glaría, Rodrigo Salas, Stéren Chabert, Pablo Roncagliolo, Alexis Arriola, Gonzalo Tapia, Matías Salinas, Herman Zepeda, Carla Taramasco, Kayode Oshinubi, Jacques Demongeot, A Step Forward to Formalize Tailored to Problem Specificity Mathematical Transforms, 2022, 8, 2297-4687, 10.3389/fams.2022.855862
    9. Adam Melnyk, Lena Kozarov, Sebastian Wachsmann-Hogiu, A deconvolution approach to modelling surges in COVID-19 cases and deaths, 2023, 13, 2045-2322, 10.1038/s41598-023-29198-4
    10. Guido Tapia-Riera, Lenin Riera-Segura, Christian Calle-Cárdenas, Isidro R. Amaro, Saba Infante, 2022, Chapter 11, 978-3-031-18271-6, 152, 10.1007/978-3-031-18272-3_11
    11. Syazwana Aziz, Muhammad Shahimi Ariffin, Kok Beng Gan, Normurniyati Abd Shattar, 2022, Smoothing of Hip Angle Kinematics Data During Parachute Landing Using Functional Data Analysis Approach, 978-1-6654-7488-7, 34, 10.1109/SCOReD57082.2022.9973921
    12. Carey Yu-Fan Ling, Piau Phang, Siaw-Hong Liew, Vivek Jason Jayaraj, Benchawan Wiwatanapataphee, Exploration of COVID-19 data in Malaysia through mapper graph, 2024, 13, 2192-6670, 10.1007/s13721-024-00472-3
    13. Jacques Demongeot, Pierre Magal, Data-driven mathematical modeling approaches for COVID-19: A survey, 2024, 50, 15710645, 166, 10.1016/j.plrev.2024.08.004
    14. Wan Anis Farhah Wan Amir, Md Yushalify Misro, Mohd Hafiz Mohd, 2024, 3203, 0094-243X, 100004, 10.1063/5.0224405
    15. Janet O. Agbaje, Oluwatosin Babasola, Kabiru Michael Adeyemo, Abraham Baba Zhiri, Aanuoluwapo Joshua Adigun, Samuel Adefisoye Lawal, Oluwole Adegoke Nuga, Roseline Toyin Abah, Umar Muhammad Adam, Kayode Oshinubi, Modeling COVID-19 Disease with Deterministic and Data-Driven Models Using Daily Empirical Data in the United Kingdom, 2024, 4, 2673-8112, 289, 10.3390/covid4020020
    16. Jacques Demongeot, Pierre Magal, Kayode Oshinubi, Forecasting the changes between endemic and epidemic phases of a contagious disease, with the example of COVID-19, 2024, 1477-8599, 10.1093/imammb/dqae012
    17. Manuel Ribeiro, Leonardo Azevedo, André Peralta Santos, Pedro Pinto Leite, Maria João Pereira, Delfina Fernandes Hlashwayo, Understanding spatiotemporal patterns of COVID-19 incidence in Portugal: A functional data analysis from August 2020 to March 2022, 2024, 19, 1932-6203, e0297772, 10.1371/journal.pone.0297772
    18. Jules Waku, Kayode Oshinubi, Umar Muhammad Adam, Jacques Demongeot, Forecasting the Endemic/Epidemic Transition in COVID-19 in Some Countries: Influence of the Vaccination, 2023, 11, 2079-9721, 135, 10.3390/diseases11040135
    19. Adejimi Adeniji, Emmanuel Addai, Shatalov Michael, Kekana Malesela, Joshua Kiddy K. Asamoah, Kayode Oshinubi, Analysis of school bullying menace incorporating family education: a mathematical modeling approach, 2025, 10, 2297-4687, 10.3389/fams.2024.1502500
    20. Jacques Demongeot, Jean-Gabriel Minonzio, A signal-processing tool adapted to the periodic biphasic phenomena: the Dynalet transform, 2024, 1477-8599, 10.1093/imammb/dqae025
  • Reader Comments
  • © 2022 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(2553) PDF downloads(75) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog