Ordinary Differential Equation with Left and Right Fractional Derivatives and Modeling of Oscillatory Systems

: We consider the principle of least action in the context of fractional calculus. Namely, we derive the fractional Euler–Lagrange equation and the general equation of motion with the composition of the left and right fractional derivatives deﬁned on inﬁnite intervals. In addition, we construct an explicit representation of solutions to a model fractional oscillator equation containing the left and right Gerasimov–Caputo fractional derivatives with origins at plus and minus inﬁnity. We derive a representation for the composition of the left and right derivatives with origins at plus and minus inﬁnity in terms of the Riesz potential, and introduce special functions with which we give solutions to the model fractional oscillator equation with a complex coefﬁcient. This approach can be useful for describing dissipative dynamical systems with the property of heredity.


Introduction
The basic framework for Lagrangian and Hamiltonian mechanics is the principle of least action.In problems of mechanics, the principle of stationary action is most commonly used, and is the most important among all extreme principles.The equations of motion obtained using these principles are called the Euler-Lagrange equations.The Euler-Lagrange equation with the Lagrange function can describe any mechanical motion.
In recent years, fractional integro-differentiation has been widely used to model the properties of dynamic systems [1][2][3].Fractional calculus enables one to naturally take into account the heredity, time non-locality and memory effects in dynamic processes [2].
The principle of least action in the context of the fractional calculus is studied in [4][5][6][7][8][9][10][11].In [4,5], it was shown that applying the Lagrangian with fractional derivatives allows one to describe the dynamics of nonconservative systems.In [6], two fractional variational problems were considered.A generalized Euler-Lagrange equation with fractional differentiation was derived.The classical Euler-Lagrange equation is a special case of the latter.In addition, practical examples of the application of the results obtained were considered.It was shown that this approach would make it possible to simulate dynamic systems with greater accuracy.In [8], the Lagrange formalism with fractional derivatives with infinite limits was considered for the first time.In [9], the Hamiltonian and Lagrangian formulations with fractional differentiation were used to describe constrained systems.It was shown that both treatments for systems with linear velocities are equivalent.In [7,12,13] the principle resulted in very exotic differential equations of fractional order, with composition of the Riemann-Liouville or Caputo operators with different origins.These equations and their possible application to modeling dynamic processes can be found in [12][13][14][15][16][17].
We point out the major challenges related to the equations with different origins.First, methods for obtaining solutions to these equations, both analytically and numerically, are not yet well developed.Next, physical interpretation of the fractional Euler-Lagrange equations is currently not well understood.The first-order Euler-Lagrange equation does not contain the limits of integration in the action functional, unlike the fractional Euler-Lagrange equation.Therefore, the following question arises: is it possible to use a fractional Euler-Lagrange equation if the time variable falls outside the interval within which the action functional is considered?Any value can be chosen as a time interval, including a small physical time (for example, the Planck time, which is approximately 5.4 × 10 −44 s).In this case, studying any physical process in this interval is believed to be incorrect.
This ambiguity in a sense can be eliminated by using the Gerasimov-Caputo fractional derivatives (fractional operators defined on an infinite interval [18,19]) in the Euler-Lagrange equation, and we offer this in the present article.
In the present paper, we derive the fractional Euler-Lagrange equation and the general equation of motion with the composition of the left and right fractional derivatives defined on infinite intervals, and construct an explicit representation of solutions to a model fractional oscillator equation containing the left and right Gerasimov-Caputo fractional derivatives with origins at plus and minus infinity.
The rest of the paper is structured as follows: Section 2 contains the definitions of fractional integrals and derivatives; in Section 3, we derive the Euler-Lagrange equation and the equation of motion with fractional derivatives defined on infinite intervals; and, in Section 4, we construct solutions of a model one-dimensional equation for the fractional oscillator; Section 5 presents the main findings of this work.

Euler-Lagrange Equation
Introduce a velocity vector for a material point motion in a certain medium: where r(t) is the radius vector, t is the dimensionless (related to τ) time and τ is some characteristic time for the process.On account of Equation ( 1), the generalized action functional can be written as In Equation ( 2) infinite limits for integration are established.To determine the velocity vector of Equation ( 1), it is necessary to take into account the entire time axis.Assume r 0 (t) is a solution for the extremum problem.Then by Equation ( 2) and the conditions The extremum of the functional satisfies the condition Using Equation ( 3) we obtain the Euler-Lagrange equation It should be noted that the same problem for classical fields was first considered in [8].In our case, Equation (4) defines the action function for a mechanical system.Moreover, the other techniques can be applied to get Equation (4), for example, using the results of [22], obtained by exploiting the Riesz fractional derivative.
Next, assume that the space of the considered motion is homogeneous and isotropic.Then the scalar function L should not be explicitly dependent on time.The function L should undergo power series expansion that involves only even powers v(t).Up to the second term of the expansion, we have Expansion coefficients m and U[r(t), t] are independent of the derivative fractional order α.Therefore, there is no indication that these coefficients are different from the mass and potential energy of an object.Substituting Equation (5) into Equation (4), we obtain the equation of motion where F[r(t), t] = −(∂/∂r) U[r(t), t].The right-hand side of Equation ( 6) does not explicitly depend on α, and, consequently, F[r(t), t] determines the conservative force acting on a material point.The parameter α, according to [23,24] takes into account energy dissipation.Therefore, the equation of motion in the form of Equation ( 6) is apparently obtained for the first time.
As α = 1, by Equation ( 6) we have the classical equation of motion

Solution of Fractional Oscillator Equation
Here, inspired by Equation ( 6), we are going to solve the equation which is a model one-dimensional equation for the fractional oscillator.
Definition 1.By a solution of Equation ( 7), we mean a function x(t) such that x (t) and x(t) satisfies Equation (7).

Riesz Potentials
Now, we will rewrite Equation (7) in terms of the Riesz potential.
Definition 2. Let µ ∈ (0, 1).The operator K µ is called the Riesz potential, and, for an integrable function x(t), it is defined by Proof.Let M = max{s, t} and m = min{s, t}.It is easy to check that Combining this with the relations we get Equation (8).
By Equation ( 8) we have and, consequently, Equation ( 7) can be rewritten in the form

Special Functions
Consider the functions and Proof.The formulas in Equation ( 9) are derived by direct differentiation.In the same way, an easy computation gives This proves Equation (10).

Solution Representation
is a regular solution of Equation (7).

Conclusions
The principle of least action is often used to describe dissipative systems, in particular in the case of linear oscillators with viscous friction [26].In this paper, we used this principle in the context of fractional calculus.In particular, we considered the concept of effective velocity for a moving material point, which was defined in terms of the Gerasimov-Caputo fractional derivative.This allowed us to take into account the heredity and nonlocality in dynamic processes of a material point with respect to the time.Using the fractional Euler-Lagrange equation, we derived the general equation of motion with Gerasimov-Caputo fractional derivatives.The equation of motion of such a form (see Equation ( 6)), has apparently been obtained for the first time in this paper.
In addition, we expressed the composition of the left and right fractional derivatives with origins at plus and minus infinities in terms of the Riesz potential, and constructed an explicit representation of solutions to a model equation of motion for a fractional oscillator.
It is of interest to give a detailed physical interpretation for Equations ( 4) and ( 6), and apply these equations to real physical processes.Undoubtedly, it is relevant to develop algorithms for the numerical solution of equations in the form of Equation (7).These issues should be the subjects of separate works.