Generalized Solutions of the Third-Order Cauchy-Euler Equation in the Space of Right-Sided Distributions via Laplace Transform

Using the Laplace transform technique, we investigate the generalized solutions of the third-order Cauchy-Euler equation of the form t3y′′′(t) + at2y′′(t) + by′(t) + cy(t) = 0, where a, b, and c ∈ Z and t ∈ R. We find that the types of solutions in the space of right-sided distributions, either distributional solutions or weak solutions, depend on the values of a, b, and c. At the end of the paper, we give some examples showing the types of solutions. Our work improves the result of Kananthai (Distribution solutions of the third order Euler equation. Southeast Asian Bull. Math. 1999, 23, 627–631).


Introduction
Differential equations are used to construct models of reality.Sometimes, the reality we are modeling suggests that some solutions of the differential equations with singularity in the coefficients need not be differentiable in the classical sense.This is where the concept of distributions or generalized functions comes from, which contains the continuous function as a subset.Kanwal [1] was the first who classified the types of solution of the linear differential equation of order n of the form: where a m (t) ∈ C ∞ (R) with a n (t) = 0 and f (t) is a given distribution.Solution types of the differential Equation ( 1) can be classified into 3 groups.The first one is the classical solution; it is the one with n-order continuously differentiable and satisfies (1) in the ordinary sense.The second one is the weak solution; it is the one with less smoothness but still locally integrable, and satisfies (1) in the distributional sense.The third one is the distributional solution; it is in fact a singular distribution and satisfies (1) in the distributional sense.All these solutions are called generalized solutions.
It is well known that the linear homogeneous ordinary differential equations with infinitely-smooth coefficients have no generalized solutions in the sense of the distribution; while the ordinary differential equations with polynomial coefficients such as the Cauchy-Euler equation defined by Equation (1), where a m (t) = c m t m with c m ∈ R, c n = 1 and f (t) = 0, may have a classical solution or a generalized solution in the sense of the distribution; see [2][3][4] for more details.The method for solving the classical solution of the Cauchy-Euler equation was explained in [5][6][7][8][9].
The generalized solutions in the sense of the distribution of ordinary differential equations can be derived from the theory of distributions.Research in these areas, still developing continually, has opened up many aspects and properties in the theory of differential and functional differential equations.
In case of the weak solutions of certain differential equations, many of them have already been studied.For example, Kananthai and Nonlaopon [10] studied the weak solutions of the compound ultra-hyperbolic equation while Sarikaya and Yildirim [11] studied those of the compound Bessel ultra-hyperbolic equation.More details related to the weak solutions of certain differential equations in the field of theory of distributions are referred to [12][13][14].
Regarding the distributional solutions, specifically as a series of Dirac delta function and its derivatives, they have been used in several areas of applied mathematics such as the theory of partial differential equations, operational calculus, and functional analysis; in Physics such as quantum electrodynamics.We refer the readers to the papers [15][16][17][18][19][20] for more details.
For the generalized solutions, Nonlaopon et al. [21] used the Laplace transform technique to study those satisfying the differential equation where m, n ∈ Z with n ≥ 2 and t ∈ R. Opio et al. [22] studied those of the differential equation with polynomial coefficients of the form where n, p, m ∈ N, n ≥ 2, and t ∈ R. Now, we consider the third-order Cauchy-Euler equation of the form where a, b, and c ∈ Z and t ∈ R. The classical solutions of the above equation are in the space C(R) of continuous functions of the form t λ , where λ is a real or a complex number; see [8] for more details.
The classical solutions of this equation using the Laplace transform technique were also studied by Kim [23].Moreover, Kananthai [24] studied the generalized solutions of (2), where a, b = 1, c is some integer, and t ∈ R using the Laplace transform technique.He found that one distributional solution and one weak solution depend on the values of m.Next, Sacorn et al. [25] studied the generalized solutions of (2), where a, b, c are real constants and t ∈ R using the Laplace transform technique.They found that one distributional solution and one weak solution of (2) depend on the values of a, b, and c.For this work our goal is to investigate the generalized solutions in the space of right-sided distributions of (2) for the case of one distributional solution, two distributional solutions, three distributional solutions, containing both the distributional solution and weak solution, the linear combination of two distributional solutions and one weak solution, and the linear combination of one distributional solution and two weak solutions.This paper is organized as follows.Next section we discuss a fundamental idea of Laplace transform.We proceed to Section 3 with the use of Laplace transform technique to our equation with a further investigation.Many examples are then presented to support our main results.Finally, we summarize our work in Section 4.

Preliminaries
Before we proceed to our main results, the following definitions and concepts are required.Definition 1.Let D be the space consisting of all real-valued functions ϕ(t) with continuous derivatives of all orders and compact support.The support of ϕ(t) is the closure of the set of all elements t ∈ R such that ϕ(t) = 0.Then, ϕ(t) is called a test function.Definition 2. A distribution T is a continuous linear functional on the space D of the real-valued functions with infinitely-differentiable and bounded support.The space of all such distributions is denoted by D .
For every T ∈ D and ϕ(t) ∈ D, the value that T has on ϕ(t) is denoted by T, ϕ(t) .Note that T, ϕ(t) ∈ R.

Example 1.
(i) The locally-integrable function f (t) is a distribution generated by the locally-integrable function f (t).
A distribution T generated by a locally-integrable function is called a regular distribution; otherwise, it is called a singular distribution.Definition 3. The kth-order derivative of a distribution T, denoted by T (k) , is defined by T Definition 4. Let α(t) be an infinitely-differentiable function.We define the product of α(t) with any distribution T in D by α(t)T, ϕ(t) = T, α(t)ϕ(t) for all ϕ(t) ∈ D.
Definition 5. Let M ∈ R and f (t) be a locally-integrable function satisfying the following conditions: (i) f (t) = 0 for all t < M; (ii) There exists a real number c such that e −ct f (t) is absolutely integrable over R.
The Laplace transform of f (t) is defined by: where s is a complex variable.
It is well known that if f (t) is continuous, then F(s) is an analytic function on the half-plane (s) > σ a , where σ a is an abscissa of absolute convergence for L { f (t)}.Definition 6.Let f (t) be a function satisfying the same conditions as in Definition 5 and L { f (t)} = F(s).The inverse Laplace transform of F(s) is defined by: where (s) > σ a .
Recall that the Laplace transform G(s) of a locally-integrable function g(t) satisfying the conditions of Definition 5, that is, where (s) > σ a , can be written in the form G(s) = g(t), e −st .
Definition 7. Let S be the space of test functions of rapid decay containing the complex-valued functions φ(t) having the following properties: , as well as its derivatives of all orders vanish at infinity faster than the reciprocal of any polynomial, which is expressed by the inequality: where C pk is a constant depending on p, k, and φ(t).Then, φ(t) is called a test function in the space S.
Definition 8.A distribution of slow growth or tempered distribution T is a continuous linear functional over the space S of the test function of rapid decay containing the complex-valued functions, i.e., a complex number T, φ(t) assigned with the properties: We shall let S denote the set of all distributions of slow growth.Definition 9. Let f (t) be a distribution satisfying the following properties: (ii) There exists a real number c such that e −ct f (t) is a tempered distribution.
The Laplace transform of a right-sided distribution f (t) satisfying (ii) is defined by: where X(t) is an infinitely-differentiable function with support bounded on the left, which equals one over a neighborhood of the support of f (t).
For (s) > c, the function X(t)e −(s−c)t is a testing function in the space S and e −ct f (t) is in the space S .Equation ( 6) can be reduced to: Now, F(s) is a function of s defined over the right half-plane (s) > c.Zemanian [26] proved that F(s) is an analytic function in the region of convergence (s) > σ 1 , where σ 1 is the abscissa of convergence and e −ct f (t) ∈ S for some real number c > σ 1 .
Example 3. Let δ(t) be the Dirac delta function, H(t) be the Heaviside function, and f (t) be a Laplace-transformable distribution in D R .If k is a positive integer, then the following hold: with infinitely-differentiable coefficients a i (t) and a n (0) = 0 has a solution: of order p, then: Conversely, if p is the smallest nonnegative integer root of (10), then there exists a pth-order solution of (9) at t = 0.
The proof of this Lemma is given in [18].Lemma 2. Let ψ(t) be an infinitely-differentiable function.Then, and: The proof of Lemma 2 is given in [1].
A useful formula that follows from (11), for any monomial ψ(t) = t n , is that:

Main Results
In this section, we will state our main results and give their proofs.
Theorem 1.Consider the third-order Cauchy-Euler equation of the form: where a, b, and c ∈ Z and t ∈ R. The types of solutions of (14) depend on the values of a, b, and c in the following way: (i) If a = 3m + 6, b = 3m 2 + 9m + 7 and c = (m + 1) 3 for some m ∈ N ∪ {0}, then there exists a distributional solution of the form: (ii) If a = 3m + n + 6, b = 3m 2 + 2mn + 9m + 3n + 7, and c = (m + 1) 2 (m + n + 1) for some m ∈ N ∪ {0} and n ∈ N, then there exist two distributional solutions of the form: (iii) If a = 3m + 2n + 6, b = 3m 2 + 4mn + n 2 + 9m + 6n + 7, and c = (m + 1)(m + n + 1) 2 for some m ∈ N ∪ {0} and n ∈ N, then there exist two distributional solutions of the form: and n ∈ N, then there exist a distributional solution and a weak solution of the form: and c = (m + 1)(n − 1) 2 for some m ∈ N ∪ {0} and n ∈ N, then there exist a distributional solution and a weak solution of the form: Proof.Applying the Laplace transform L{y(t)} = Y(s) to ( 14) and using Example 3(iv),(v), we obtain: Suppose that a solution of ( 20) is of the form Y(s) = s r , where r is a real constant.Substituting Y(s), Y (s), Y (s), and Y (s) into (20), we obtain: Since s r = 0, we have: Case (i): If a = 3m + 6, b = 3m 2 + 9m + 7, and c = (m + 1) 3 for some m ∈ N ∪ {0}, then substituting a, b, and c into (21) yields: and so: (r Thus, we have the repeated real roots of (22), that is r = m.Therefore, the solution of ( 20) is Y(s) = s m .Now, Y(s) is an analytic function over the entire s-plane.Taking the inverse Laplace transform to Y(s) and using Example 3(ii),(iii), we obtain a solution of (14), which is a distributional solution of the form (15). Case (ii): If a = 3m + n + 6, b = 3m 2 + 2mn + 9m + 3n + 7, and c = (m + 1) 2 (m + n + 1) for some m ∈ N ∪ {0}, n ∈ N, then substituting a, b and c into (21) yields: Thus, we have the real roots of (23), which are r = m and r = m + n.Therefore, the solutions of (20) are Y(s) = s m and Y(s) = s m+n .Now, Y(s) are analytic functions over the entire s-plane.Taking the inverse Laplace transform to Y(s) and using Example 3(ii),(iii), we obtain the solutions of ( 14), which are distributional solutions of the form (16).
Case (iv): 2 for some m ∈ N ∪ {0} and n ∈ N, then substituting a, b, and c into (21) yields: and so: Thus, we have the real roots of (25), which are r = m and r = −n.Therefore, the solutions of (20) are Y(s) = s m and Y(s) = s −n .
Clearly Y(s) are analytic functions on the right half-plane (s) > 0. Applying the inverse Laplace transform to Y(s) and using Example 3(i)-(iii), we obtain the solutions of ( 14), which are a distributional solution and a weak solution of the form (18). 2 for some m ∈ N ∪ {0} and n ∈ N, then substituting a, b, and c into (21) yields: Thus, we have the real roots of (26), which are r = m and r = −n.Therefore, the solutions of (20) are Y(s) = s m and Y(s) = s −n .Now, Y(s) are analytic functions over the half-plane (s) > 0. Taking the inverse Laplace transform to Y(s) and using Example 3(i)-(iii), we obtain the solutions of ( 14), which are a distributional solution and a weak solution of the form (19).
Theorem 2. The type of solutions of (14) whose solutions are a linear combination of distributional solutions and weak solutions depends on the values of a, b, and c in the following way: (i) If a = 3m + 2n + p + 6, b = 3m 2 + 9m + n 2 + n(4m + p + 6) + p(2m + 3) + 7, and c = (m + 1)(m + n + 1)(m + n + p + 1) for some n, p ∈ N and m ∈ N ∪ {0}, then all solutions are a linear combination of the distributional solutions of the form: where C 1 , C 2 , and C 3 are arbitrary constants.
for some n, p ∈ N and m ∈ N ∪ {0}, then all solutions are a linear combination of distributional solutions and a weak solution of the form: where C 1 , C 2 , and C 3 are arbitrary constants.
for some n, p ∈ N and m ∈ N ∪ {0}, then all solutions are a linear combination of the distributional solution and weak solutions of the form: where C 1 , C 2 , and C 3 are arbitrary constants.

Conclusions
We used the Laplace transform technique to find the generalized solutions of the third-order Cauchy-Euler equation of the form: t 3 y (t) + at 2 y (t) + bty (t) + cy(t) = 0, where a, b, and c ∈ Z and t ∈ R.Then, we took the inverse Laplace transform to the derived solutions.We found the conditions of a, b, and c, for the case of one distributional solution, two distributional solutions, three distributional solutions, containing both a distributional solution and a weak solution, a linear combination of two distributional solutions and one weak solution, and a linear combination of one distributional solution and two weak solutions.It should be noted here that the inverse Laplace transform of s n , where n ∈ N as shown in Examples 4 and 5, however, for the classical solutions is not mentioned here, which can be found in Greenberg [8].