Advances in Mathematical Physics

Volume 2018, Article ID 7283518, 9 pages

https://doi.org/10.1155/2018/7283518

## An Exact Solution of the Second-Order Differential Equation with the Fractional/Generalised Boundary Conditions

^{1}Institute of Computer and Information Sciences, Czestochowa University of Technology, Czestochowa, Poland^{2}Institute of Mathematics, Czestochowa University of Technology, Czestochowa, Poland

Correspondence should be addressed to Tomasz Blaszczyk; lp.zcp.mi@kyzczsalb.zsamot

Received 30 December 2017; Accepted 26 February 2018; Published 11 April 2018

Academic Editor: Giorgio Kaniadakis

Copyright © 2018 Mariusz Ciesielski and Tomasz Blaszczyk. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We analysed the initial/boundary value problem for the second-order homogeneous differential equation with constant coefficients in this paper. The second-order differential equation with respect to the fractional/generalised boundary conditions is studied. We presented particular solutions to the considered problem. Finally, a few illustrative examples are shown.

#### 1. Introduction

The second-order differential equations provide an important mathematical tool for modelling the phenomena occurring in dynamical systems. Examples of linear or nonlinear equations appear in almost all of the natural and engineering sciences and arise in many fields of physics.

Many scientists have studied various aspects of these problems, such as physical systems described by the Duffing equation [1], noncommutative harmonic oscillators [2], oscillators in quantum physics [3], the dynamic properties of biological oscillators [4], the analysis of single and coupled low-noise microwave oscillators [5], the Mathieu oscillator [6], the relativistic oscillator [7], or the Schrodinger type oscillator [8].

Classical differential equations are defined by using the integer order derivatives. In recent years, the class of differential equations containing fractional derivatives (known as fractional differential equations) have become an important topic. There are two approaches to obtain these types of equations. The first one is to replace the integer order derivative in classical differential equations with a fractional derivative (see, e.g., [9–14]).

The second approach is a generalisation of a method known in classical and quantum mechanics, where the differential equations are obtained from conservative Lagrangian or Hamiltonian functions. These equations are known in the literature as fractional Euler-Lagrange equations, and they contain both the left and right fractional derivatives. New mechanics models for nonconservative systems, in terms of fractional derivatives, were developed by Riewe in [15, 16] and extended by Klimek [17, 18] and Agrawal [19, 20]. Since then, many authors have studied the fractional differential equations of the variational type (see [21–26]).

In contrast to the above-mentioned references, where the authors analysed the integer order differential equations with classical boundary/initial conditions or fractional differential equations with the Dirichlet or natural boundary conditions, in this paper we consider the second-order differential equation with the fractional/generalised boundary conditions.

#### 2. Statement of the Problem

In this paper, we solve the second-order differential equationwith respect to the following fractional/generalised boundary conditions:where , , , and the operators denote the left and right Riemann-Liouville derivatives, defined, respectively, by [27]and .

Let us consider two particular cases of (1).

*Case i. *Here, we study the following differential equation:It is well known that (4) has a general solution, for , given byThe solution (5) contains two arbitrary independent constants of integration and . A particular solution can be derived from the general solution by applying the set of initial or boundary conditions.

The fractional differentiation of general solution (5) (using the left-side fractional operator) gives usand differentiation by using the right-side operator leads toNow, we formulate the following properties for Riemann-Liouville derivatives of the sine and cosine functions.

*Property 1 (the left-sided Riemann-Liouville fractional derivatives of the sine and cosine functions). *Let and . Then, the following relations hold:

*Property 2 (the right-sided Riemann-Liouville fractional derivatives of the sine and cosine functions). *Let and . Then, the following relations hold:

*Proof (Properties 1 and 2). *We use Taylor’s series expansions of sine and cosine functions [28]and properties of the left and right-sided fractional derivatives of power functions [27] Also we apply the known fundamental trigonometric identitiesThenwhereand in a similar way we obtainFinally, putting (17)–(20) into (13)–(16), we obtain the formulas in Properties 1 and 2.

*Remark 3. *Note that the infinite series included in formulas (17)–(20) can be expressed by formulas containing the Mittag-Leffler function. This observation leads us to the following expressions:where denotes the Mittag-Leffler function defined in [27]The above-mentioned notations can be useful in case when one uses the built-in function (the Mittag-Leffler function) in a mathematical software.

*Case ii. *The second problem has the following form:In this case the general solution of (23), for , is given byThe fractional differentiation of solution (24) (using the left-side operator) givesand for the right-side derivative we haveNext, we formulate the following properties for Riemann-Liouville derivatives of the hyperbolic sine and hyperbolic cosine functions.

*Property 4 (the left-sided Riemann-Liouville fractional derivatives of the hyperbolic sine and hyperbolic cosine functions). *Let and . Then, the following relations hold:

*Property 5 (the right-sided Riemann-Liouville fractional derivatives of the hyperbolic sine and hyperbolic cosine functions). *Let and . Then, the following relations hold:

*Proof (Properties 4 and 5). *Here we apply Taylor’s series expansions of the hyperbolic sine and cosine functions [28]and the trigonometric identitiesThenwhereFinally, putting (32) into (31), we obtain the formulas in Properties 4 and 5.

*Remark 6. *In formulas (32), the infinite series can be also expressed by using the Mittag-Leffler function, and we obtain

#### 3. Examples of the Determination of Particular Solutions

The boundary conditions, written in the general form (2), can be used in many combinations. Now, we show three selected examples. Other combinations of the particular boundary conditions can be easily adopted by the reader (in a similar way).

*Example 7. *Equation (4) with the following boundary conditions given on both sides of the domain:We substitute the general solution (5) into (34) and we have The independent constants of integration and can be determined from the solution of the following system of equations:The analytical solution of (36) is of the form

*Example 8. *Equation (23) with the conditions given on the left side of the domain (): this case corresponds to the initial value problemWhen we put the general solution (24) into (38), then we obtain the following linear system of equations: that can be also written in the matrix form

*Example 9. *Equation (4) with the following set of boundary conditions: This case corresponds to the generalisation of the Robin boundary condition given at . Also here we put general solution (5) into (41) and we obtain orfrom which we can easily determine the constants of integration and .

*Remark 10. *One can note that the fractional boundary conditions and for integral values of parameters *α* and *β* take the classical forms of boundary conditions; this means and . In particular, it should be noted that the difference occurs, among others, in the boundary condition and this form should be taken into account.

#### 4. Example of Solutions

On the basis of the proposed method, to find the particular solutions to the considered equations (4) and (23), we calculated constants of integration and occurring in the general solutions that satisfy the sets of the given initial or boundary conditions (various combinations). In Figures 1–3, numerous examples of solutions have been presented.