Abstract
We discuss the onesided Green’s function, associated with an initial value problem and the twosided Green’s function related to a boundary value problem. We present a specific calculation associated with a differential equation with constant coefficients. For both problems, we also present the Laplace integral transform as another methodology to calculate these Green’s functions and conclude which is the most convenient one. An incursion in the socalled fractional Green’s function is also presented. As an example, we discuss the isotropic harmonic oscillator.
1 Introduction
There are several methods to discuss a secondorder linear partial differential equation. Among them we mention the simplest one, the method of separation of variables, and the method of integral transforms, particularly the Laplace transform, which is in many cases most convenient [1].
On the other hand, after the method of separation of variables, we get the following general secondorder linear ordinary differential equation:
on the interval
where
the socalled selfadjoint form known also as an ordinary differential equation in
the SturmLiouville form. In this equation,
The methodology of the Laplace integral transform is adequate to discuss the onesided Green’s function^{a} because the initial conditions, in general, are given in terms of the own function and the first derivative. A simple question arises when we discuss the twosided Green’s function associated with a problem involving boundary conditions, i.e., is the Laplace transform methodology convenient to discuss this problem? The answer depends on the sort of problem we are studying, as we will see in the following sections.
On the other hand, the fractional harmonic oscillator was discussed in a series of papers by Narahari et al.[25] where they presented the dynamic of the fractional harmonic oscillator, including also the damping, and by Tofighi [6] who discusses the intrinsic damping.
In this paper we discuss Eq. (1) associated with an initial value problem and a boundary value problem. In both cases, we present two methodologies, the Laplace integral transform and the Green’s function methodology. After that we conclude which methodology is the most convenient one. We sum the paper up presenting the corresponding fractional case where we discuss the Green’s function associated with the fractional harmonic oscillator. Finally, we present our concluding remarks.
2 Onesided Green’s function
To solve the selfadjoint differential equation, we introduce the socalled onesided
Green’s function, also called influence function, a twoparameter function, denoted
by
where
For a fixed value of ξ, the onesided Green’s function is the solution of the corresponding homogeneous
initial value problem, i.e., for
and the initial conditions
2.1 Constant coefficients
As a particular case, we consider the problem involving an ordinary differential equation with constant coefficients, which can represent a problem related to the classical harmonic oscillator and transmission lines, for example,
where a and b are positive constants, and the homogeneous initial conditions are given by
First of all, we write the ordinary differential equation in the corresponding selfadjoint ordinary differential equation
where we identify
Thus, the onesided Green’s function, denoted by
and the initial conditions
Two linearly independent solutions of the homogeneous ordinary differential equation are given by
which furnishes for the general solution (the Green’s function)
where
Solving the system above, we obtain
which furnish for the onesided Green’s function, for
and, as we know, satisfies the property
Finally, the solution of our initial value problem is given by
2.2 Laplace transform
Another way to calculate this onesided Green’s function is by means of the methodology
of Laplace transform. Multiplying Eq. (3) by the kernel of Laplace transform,
where
Using the initial conditions, we obtain an algebraic equation whose solution is given by
By means of the convolution product and the relation involving the inverse Laplace transform
we get
which is exactly Eq. (5).
3 Twosided Green’s function
In the case that we have a twopoint boundary value problem, i.e., when the boundary conditions are fixed on the extremes of the interval, we have a twosided Green’s function, also called Green’s function, only.
We consider the secondorder ordinary differential equation in the SturmLiouville form
for
The twosided Green’s function, denoted by
for
This Green’s function must satisfy also the continuity
and a jump discontinuity on the first derivative [1]
The solution
where
3.1 Constant coefficients
As an example, we discuss the following ordinary differential equation:
with the homogeneous boundary conditions
First, we construct the corresponding Green’s function, as we have seen before. Two linearly independent solutions of the homogeneous ordinary differential equation are
Imposing that
On the other hand, to satisfy
where
Using the continuity at
Thus, the Green’s function is given by the following expression:
Using this Green’s function, the solution of the boundary value problem can be written as follows:
3.2 Laplace transform
As we have already said, the Laplace transform converts the differential equation with constant coefficients into an algebraic equation whose solution is
where
Using the inverse Laplace transform and the convolution product, we get
Thus, substituting the boundary conditions, we obtain two algebraic equations, a
system involving
which can be rewritten as follows:
where the Green’s function is given by
which is the same expression as that obtained in Section 3.1.
At this point we conclude that for a problem involving initial conditions, the Laplace
integral transform is more convenient since
4 Fractional Green’s function
Fractional calculus is one of the most accurate tools to refine the description of natural phenomena. The usual way to use this tool is to replace the integerorder derivatives of the partial differential equation that describes one specific phenomenon by a derivative of noninteger order. For many expected reasons, the solution of a fractional partial differential equation used to be much more complicated than the solution of the corresponding integerorder partial differential equation.
On the other hand, many important results and generalizations were obtained using this procedure in several areas such as fluid flow, diffuse transport, electrical networks, probability, biomathematics and others [712]. Here, as a generalization to the integer case, we present a calculation associated with the socalled fractional onesided and twosided Green’s function relative to the fractional differential equation with constant coefficients, i.e., we obtain the fractional Green’s function for the fractional differential equation whose coefficients are constants. We discuss the problem by means of the Laplace integral transform, and as an application, we present explicitly the Green’s function associated with the fractional harmonic oscillator.
4.1 Fractional onesided Green’s function
Let a, b and c be real constants. We present the solution of the fractional differential equation
where
Introducing the Laplace transform and using the initial conditions, we obtain an algebraic equation whose solution can be written as follows:
where
which is valid for
Using the Laplace transform of the generalized MittagLeffler function and its corresponding inverse [13],
with
which is the solution of the fractional differential equation.
Thus, the onesided fractional Green’s function can be written as follows:
Taking
which is the fractional onesided Green’s function associated with the fractional
harmonic oscillator. The onesided Green’s function associated with the classical
harmonic oscillator is recovered by introducing
which can also be written as follows:
which is the same expression as that obtained in Eq. (4) in the case
To conclude this section, we plot, in Figure 1, a graphic
Figure 1. Graphics for
4.2 Fractional twosided Green’s function
Let a, b and c be real constants. We present the solution of the fractional differential equation,
Eq. (6), with the homogeneous boundary conditions,
By means of the Laplace transform, we can write
where
From the inverse Laplace transform and the convolution theorem, we get
which can be rewritten in the following way:
As we have already said in Section 3, to explicitly calculate the solution, we must
substitute the homogeneous boundary conditions in the last equation and determine
the
5 Concluding remarks
We have presented and discussed the socalled onesided and twosided Green’s function to study, respectively, an initial value problem and a boundary value problem. Besides that, we studied the same problems by means of the Laplace transform methodology in order to conclude which methodology was most accurate for each problem. We also obtained the fractional generalization of the onesided and twosided Green’s function in terms of the generalized MittagLeffler function.
We conclude that for the initial value problem, the Laplace integral transform methodology is more convenient; on the other hand, for the boundary value problem, the twosided Green’s function provides a much more simple calculation. It is important to note that in the present manuscript we did not consider the problem involving physical dimensions. This problem has been discussed in a recent paper by Inizan [16].
A natural continuation of this work would be to study the problems involving partial differential equations with nonconstant coefficients and their fractional versions, which could provide a better description of the phenomena related to those equations. Study in this direction is upcoming.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
RFC had the idea of the work. ECO had the idea to consider the fractional version. AOC had the fractional harmonic oscillator, idea. All authors have made the calculations involved. All authors read and approved the final manuscript.
Acknowledgements
We are grateful to Prof. J. Vaz Jr., Dr. J. Emílio Maiorino and Dr. E. Contharteze Grigoletto for several useful discussions. Besides that, we are thankful to the referees for several important suggestions which improved this article a lot.
End notes

A recent historical review about George Green can be found in ref. [17].
References

Capelas de Oliveira, E: Special Functions and Applications, Editora Livraria da Física, São Paulo (2012) (in Portughese)

Clarke, T, Narahari Achar, BN, Hanneken, JW: MittagLeffler functions and transmission lines. J. Mol. Liq.. 114, 159–163 (2004). Publisher Full Text

Narahari Achar, BN, Hanneken, JW, Clarke, T: Response characteristics of a fractional oscillator. Physica A. 309, 275–288 (2002). Publisher Full Text

Narahari Achar, BN, Hanneken, JW, Clarke, T: Damping characteristic of a fractional oscillator. Physica A. 339, 311–319 (2004). Publisher Full Text

Narahari Achar, BN, Hanneken, JW: Dynamic response of the fractional relaxoroscillator to a harmonic driving force. In: Sabatier J (ed.) Advances in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering, Springer, Heidelberg (2007)

Tofighi, A: The intrinsic damping of the fractional oscillator. Physica A. 329, 29–34 (2003). Publisher Full Text

Capelas de Oliveira, E, Mainardi, F, Vaz, J Jr..: Models based on MittagLeffler functions for anomalous relaxation in dielectrics. Eur. Phys. J. Spec. Top.. 193, 161–171 (2011). Publisher Full Text

Capelas de Oliveira, E, Vaz, J Jr..: Tunneling in fractional quantum mechanics. J. Phys. A, Math. Theor.. 44, Article ID 185303 (2011)

Figueiredo Camargo, R, Charnet, R, Capelas de Oliveira, E: On the fractional Green function. J. Math. Phys.. 50, Article ID 043514 (2009)

Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations, Elsevier, Amsterdam (2006) Edited by Jan Van Mill

Klafter J, Lim SC, Metzler R (eds.): Fractional Dynamics, Recent Advances, World Scientific, Singapore (2011)

De Paoli, AL, Rocca, MC: The fractionary Schrödinger equation, Green functions and ultradistribution. Phys. A, Stat. Mech. Appl.. 392, 111–122 (2013). Publisher Full Text

Figueiredo Camargo, R, Chiacchio, AO, Capelas de Oliveira, E: Differentiation to fractional orders and the fractional telegraph equation. J. Math. Phys.. 49, Article ID 033505 (2008)

Mainardi, F: Fractional Calculus and Waves in Linear Viscoelasticity, Imperial College Press, London (2010)

Gorenflo, R, Mainardi, F: Fractional calculus, integral and differential equations of fractional order. In: Carpinteri A, Mainardi F (eds.) Fractals and Fractional Calculus in Continuum Mechanics, pp. 223–276. Springer, Wien (1997) (Eprint arxiv:0805.3823), in particular, pp. 247252

Inizan, P: Homogeneous fractional embeddings. J. Math. Phys.. 49, Article ID 082901 (2008)

Challis, L, Sheard, F: The Green of Green function. Phys. Today. 56, 41–46 (2003)