“…The Caputo derivative has some advantages over the Riemann-Liouville derivative [6]. Some of these can be summarized as follows: The Caputo derivative used in this study is very useful when dealing with a real-world problem because it allows conventional initial and boundary conditions to be included in the formulation of the problem, and in addition the derivative of a constant is zero [7]. Fractional-order differential equations (FODEs), being a form of interest to scientists in mathematical modeling in recent years, are used in the fields such as physics [8], thermodynamics [9], viscoelastic [10], electrical theory [11], mechatronics [12], medicine [13], chemistry [14,15], chaos theory [16], finance and economics [17], biological systems [18,19], etc.…”