Next Article in Journal
Generalized Mann Viscosity Implicit Rules for Solving Systems of Variational Inequalities with Constraints of Variational Inclusions and Fixed Point Problems
Previous Article in Journal
Some Characterizations of Generalized Null Scrolls
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Generalized Solutions of the nth Order Cauchy–Euler Equation

by
Amornrat Sangsuwan
1,
Kamsing Nonlaopon
1,*,
Somsak Orankitjaroen
2 and
Ismail Mirumbe
3
1
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
3
Department of Mathematics, Makerere University, Kampala 7062, Uganda
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(10), 932; https://doi.org/10.3390/math7100932
Submission received: 23 August 2019 / Revised: 25 September 2019 / Accepted: 3 October 2019 / Published: 10 October 2019

Abstract

:
In this paper, we use the Laplace transform technique to examine the generalized solutions of the nth order Cauchy–Euler equations. By interpreting the equations in a distributional way, we found that whether their solution types are classical, weak or distributional solutions relies on the conditions of their coefficients. To illustrate our findings, some examples are exhibited.

1. Introduction

The nth order Cauchy–Euler equations
a n t n y ( n ) ( t ) + a n 1 t n 1 y ( n 1 ) ( t ) + + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where a 0 , a 1 , , a n are real constant coefficients and t R , are often one of the first higher order ordinary homogeneous linear differential equations with variable coefficients introduced in an undergraduate level course. Naturally, we will discuss the second order Cauchy–Euler equations first. The appropriate form for its solution is y = t r where r is a parameter to be resolved. Replacing y with t r in the Cauchy–Euler equations yields the characteristic polynomial whose roots determine the forms of the general solution (e.g., see the textbooks [1,2]). This same technique can be carried over to solve the higher order Cauchy–Euler equations.
In the framework of distribution theory, R. P. Kanwal [3] classifies solution the type of ordinary homogeneous linear differential equations
a n ( t ) y ( n ) ( t ) + a n 1 ( t ) y ( n 1 ) ( t ) + + a 1 ( t ) y ( t ) + a 0 ( t ) y ( t ) = 0 ,
where the coefficient functions a 0 ( t ) , a 1 ( t ) , , a n ( t ) are infinitely differentiable and t R . The type can be explained as follows. The solution is a classical solution if it is at least n times continuously differentiable so that the differentiation in Equation (2) can be achieved in the ordinary sense with an identity result. The solution is a weak solution if it is less than n times continuously differentiable and thus it does not satisfy Equation (2) in the ordinary sense but in the weak or distributional sense. The solution is a distributional solution if it is a singular distribution satisfying Equation (2) in the weak sense. All of these are referred to as generalized solutions. It is widely known that the normal form of Equation (2) does not have weak or distributional, but classical solutions. Of particular interest are the singular distributions appearing as a finite series of the Dirac delta function and its derivatives. They can arise as a distributional solution for certain classes of ordinary differential equations with singular coefficients (see J. Wiener [4] in 1982). The applications of the distribution theory to differential equations have been examined by L. Schwartz [5] and A. H. Zemanian [6]. In 1983, J. Wiener and S. M. Shah [7] provided an overview of research in the distributional field and proposed a unified way in the investigation of both distributional and entire solutions to some classes of linear ordinary differential equations. Many mathematicians have also studied the distributional solutions in the field of theory of distributions, as can be seen in [8,9,10,11,12].
A. Kananthai [13] in 1999 considered certain third order Cauchy–Euler equations
t 3 y ( t ) + t 2 y ( t ) + t y ( t ) + m y ( t ) = 0 ,
where m is an integer and t R . He constructed a formula for m corresponding to each type of generalized solution of Equation (3), which are Laplace transformable. In 2017, A. Liangprom and K. Nonlaopon [14] extended the same study to certain fourth order Cauchy–Euler equations, a natural extension of Equation (3). The result for the general nth order Cauchy–Euler equations in this form was finally established by A. Sangsuwan, K. Nonlaopon and S. Orankitjaroen [15] one year after.
In 2018, P. Jodnok and K. Nonlaopon [16] presented the generalized solutions of the fifth order Cauchy–Euler equations of the form
t 5 y 5 ( t ) + a 4 t 4 y 4 ( t ) + a 3 t 3 y ( t ) + a 2 t 2 y ( t ) + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where a 0 , a 1 , , a 4 are real constants and t R . Depending on the values of a 0 , a 1 , , a 4 , they showed that the solutions of Equation (4) are either the weak solutions or the distributional solutions.
In 2015, S. Nanta [17] studied the distributional solutions of the nth order Cauchy–Euler
a n t n y ( n ) ( t ) + a n 1 t n 1 y ( n 1 ) ( t ) + + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where a i , i = 0 , 1 , , n are real constants using Fourier transform. She found that the type of solutions of Equation (5) depend on the conditions of a i .
Here we aim to seek the generalized solutions of the nth order Cauchy–Euler equations of the form of Equation (5) in the space of right-sided distributions. The solutions are obtained by applying the Laplace transform technique. Our work is an improved version of that of A. Sangsuwan et al. [15].
The present paper is arranged into three sections. In Section 2, we provide related definitions and lemmas necessary to obtain our main results. We then proceed to prove our results together with supported examples in Section 3.

2. Preliminaries

The space D (the space of distributions) is the dual space of D , the space of testing functions.
Definition 1.
A distribution T D is a continuous linear functional on D . The value of T acting on a testing function ϕ ( t ) is written as T , ϕ or T , ϕ ( t ) , and T , ϕ ( t ) C , where C is the set of complex numbers.
Distributions that are most useful are those generated by locally integrable functions. In fact, every locally integrable function f ( t ) generates a distribution, which is defined by
f , ϕ = R f ( t ) ϕ ( t ) d t .
Definition 2.
The kth order derivative of a distribution T is defined by
T ( k ) , ϕ ( t ) = ( 1 ) k T , ϕ ( t )
for all ϕ ( t ) D .
Definition 3.
Let f ( t ) be a locally integrable function which satisfies the following conditions:
(i) 
f ( t ) = 0 for all t < 0 ;
(ii) 
there exists a real number c such that e c t f ( t ) is absolutely integrable over R .
Then the Laplace transform of f ( t ) is
F ( s ) = L { f ( t ) } = f ( t ) , e s t ,
where s is a complex variable.
Furthermore, if f is continuous, then its Laplace transform F ( s ) is analytic on the half-plane R e ( s ) > σ a , where σ a is an abscissa of absolute convergence for L { f ( t ) } .
Definition 4.
Let f ( t ) be a function satisfying the conditions in Definition 3, and L { f ( t ) } = F ( s ) . The inverse Laplace transform of F ( s ) is defined by
f ( t ) = L 1 { F ( s ) } = 1 2 π i lim ω c i ω c + i ω F ( s ) e s t d s ,
where R e ( s ) > σ a .
Recall that the Laplace transform G ( s ) of a locally integrable function g ( t ) that satisfies the conditions of Definition 3 is
G ( s ) = L { g ( t ) } = g ( t ) , e s t ,
where R e ( s ) > σ a .
Definition 5.
Let f ( t ) be a distribution satisfying the following properties:
(i) 
f is a right-sided distribution, that is, f D R .
(ii) 
There exists a real number c for which e c t f ( t ) is a tempered distribution.
The Laplace transform of a right-sided distribution f ( t ) satisfying (ii) is defined by
F ( s ) = L { f ( t ) } = e s t f ( t ) , X ( t ) e ( s c ) t ,
where X ( t ) is an infinitely differentiable function with support bounded on the left, which equals to 1 over the neighborhood of the support of f ( t ) .
For R e ( s ) > c , X ( t ) e ( s c ) t is a testing function in the space S of testing functions of rapid descent and e c t f ( t ) is in the space S of tempered distributions. Equation (9) can be deduced to
F ( s ) = L { f ( t ) } = f ( t ) , e s t ,
then Equation (10) posses the sense given by the right-hand side of Equation (9). Now, F ( s ) is a function of s defined over the right half-plane R e ( s ) > c . A. H. Zemanian [6] proved that F ( s ) is an analytic function in the region of convergence R e ( s ) > σ 1 , where σ 1 is the abscissa of convergence for which e c t f ( t ) S for some real c > σ 1 . For more details about the Laplace transform of distributions, see [18,19] and the references therein.
Example 1.
Let H ( t ) be the Heaviside function, δ ( t ) be the Dirac delta function and f ( t ) be a Laplace-transformable distribution in D R . If k is a positive integer, then
(i) 
L { ( t k 1 H ( t ) ) / ( k 1 ) ! } = 1 / s k , R e ( s ) > 0 .
(ii) 
L δ ( t ) = 1 , < R e ( s ) < .
(iii) 
L δ ( k ) ( t ) = s k , < R e ( s ) < .
(iv) 
L t k f ( t ) = ( 1 ) k F ( k ) ( s ) , R e ( s ) > σ 1 .
(v) 
L f ( k ) ( t ) = s k F ( s ) , R e ( s ) > σ 1 .
Lemma 1.
Let ψ ( t ) be an infinitely differentiable function. Then
ψ ( t ) δ ( m ) ( t ) = ( 1 ) m ψ ( m ) ( 0 ) δ ( t ) + ( 1 ) m 1 m ψ ( m 1 ) ( 0 ) δ ( t ) + ( 1 ) m 2 m ( m 1 ) 2 ! ψ ( m 2 ) ( 0 ) δ ( t ) + + ψ ( 0 ) δ ( m ) ( t ) .
and
[ ψ ( t ) H ( t ) ] ( m ) = ψ ( m ) ( t ) H ( t ) + ψ ( m 1 ) ( 0 ) δ ( t ) + ψ ( m 2 ) ( 0 ) δ ( t ) + + ψ ( 0 ) δ ( m 1 ) ( t ) .
We refer the reader to [3] for a proof of Lemma 1.
A useful formula that follows from Equation (11), for any monomial ψ ( t ) = t n , is
t n δ ( m ) ( t ) = 0 , for m < n ; ( 1 ) n m ! ( m n ) ! δ ( m n ) ( t ) , for m n .
Lemma 2.
If the equation
i = 0 n a i ( t ) t i y ( i ) ( t ) = 0
with coefficients a i ( t ) C n and a n ( 0 ) 0 has a solution
y ( t ) = i = 0 k b i δ ( i ) ( t ) , b k 0 ,
of order k (order of distribution Equation (15)), then we have
i = 0 n ( 1 ) i a i ( 0 ) ( k + i ) ! = 0 .
Conversely, if k is the smallest non-negative integer root of Equation (16), there exists a kth order solution of Equation (15) at t = 0 .
We refer the reader to [4] for a proof of Lemma 2.

3. Main Results

Equipped with the Laplace transform technique, we are now ready to prove our main results.
Theorem 1.
Consider the nth order Cauchy–Euler equations of the form
a n t n y ( n ) ( t ) + a n 1 t n 1 y ( n 1 ) ( t ) + + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where a i , i = 0 , 1 , 2 , , n are real constants, a n 0 , n is any integers with n 2 and t R . The types of Laplace transformable solutions in D R of Equation (17) depend on the value of a i , and are given by the following cases:
(i) 
If there exists a non-negative integer k such that
i = 0 n ( 1 ) i ( k + i ) ! a i = 0 ,
then there exists a distributional solution of Equation (17), which is a singular distribution of the Dirac delta function and its derivatives.
(ii) 
If there exists a non-negative integer k less than or equal to n such that
i = 0 k a i ( k i ) ! = 0 ,
then there exists a weak solution of Equation (17). Moreover, the solution is continuous if k is greater than or equal to 1.
(iii) 
If there exists a positive integer k such that
i = 0 n a i ( n + k i ) ! = 0 ,
then there exists a classical solution of Equation (17).
Proof. 
We rewrite Equation (17) in brief as
i = 0 n a i t i y ( i ) ( t ) = 0 ,
where a n = 1 . Applying Laplace transform to Equation (21) with a notation L { y } = Y ( s ) , we now refer to properties (iv) and (v) in Example 1 to get
i = 0 n ( 1 ) i a i d i d s i s i Y ( s ) = 0 ,
that is,
i = 0 n ( 1 ) i a i m = 0 i i m 2 m ! s i m Y ( i m ) ( s ) = 0 ,
or
i = 0 n ( 1 ) i a i i ! j = 0 i i j 1 j ! s j Y ( j ) ( s ) = 0 .
Consider a solution of Equation (23) in a simple form Y ( s ) = s r , where r is a real constant that must be determined. Replacing Y ( i ) ( s ) for i = 1 , 2 , 3 , , n in Equation (23) gives
i = 0 n ( 1 ) i a i i ! j = 0 i i j r j s r = 0 .
The identity
j = 0 i i j r j = r + i i
and s r 0 imply that
i = 0 n ( 1 ) i i ! r + i i a i = 0 .
We now examine the proposed three cases of the value r.
Case (i). If r is a non-negative integer, then substituting r = k for k N { 0 } into Equation (24), we obtain Equation (18). Thus, if the condition of Equation (18) holds, then the solution of Equation (21) is Y ( s ) = s k . Obviously Y ( s ) is analytic over the whole s-plane. Taking inverse Laplace transform to Y ( s ) and applying property (iii) in Example 1, we obtain the distributional solutions of Equation (17) in the form
y ( t ) = δ ( k ) ( t ) .
Case (ii). If r is a negative integer which is no less than ( n + 1 ) , then substituting r = ( k + 1 ) for k { 0 , 1 , 2 , , n } into Equation (24), we obtain
i = 0 n ( 1 ) i i ! ( k + 1 ) + i i a i = i = 0 k ( 1 ) i i ! ( k + 1 ) + i i a i = i = 0 k ( 1 ) i i ! ( ( k + 1 ) + i ) ( ( k + 1 ) + i 1 ) ( ( k + 1 ) + i i + 1 ) i ! a i = i = 0 k k ! ( k i ) ! a i = 0 .
Thus, if the condition of Equation (19) holds, then the solution of Equation (21) is Y ( s ) = s ( k + 1 ) . Now we take the inverse Laplace transform to Y ( s ) , applying property (i) in Example 1, and we obtain the weak solutions of Equation (17), since k n in the form
y ( t ) = H ( t ) t k k ! .
Observe that the solution is continuous for k 1 .
Case (iii). If r is a negative integer less than ( n + 1 ) , then a substitution of r = ( n + k + 1 ) for k N into Equation (24), gives
i = 0 n ( 1 ) i i ! ( n + k + 1 ) + i i a i = i = 0 n ( 1 ) i i ! ( ( n + k + 1 ) + i ) ( ( n + k + 1 ) + i 1 ) ( ( n + k + 1 ) + i i + 1 ) i ! a i = i = 0 k ( n + k ) ! ( n + k i ) ! a i = 0 .
Thus, if the condition of Equation (20) holds, then the solution of Equation (21) is Y ( s ) = s ( n + k + 1 ) . Now we take the inverse Laplace transform to Y ( s ) and, applying property (i) in Example 1, we obtain the classical solutions of Equation (17) because k 1 and the solutions are
y ( t ) = H ( t ) t n + k ( n + k ) ! .
 □
Theorem 2.
The distributional solution of the nth order Cauchy–Euler equations of the form
a n t n y ( n ) ( t ) + a n 1 t n 1 y ( n 1 ) ( t ) + + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where a i , i = 0 , 1 , 2 , , n are real constants, n is any integer and t R , depends on the values of a i , i = 0 , 1 , , n of the form
i = 0 n ( 1 ) i ( k + i ) ! a i = 0 ,
for any k N , being the order of the distribution.
Proof. 
Using Lemma 2 and substituting a i ( 0 ) = a i for i = 0 , 1 , , n 1 and a n ( 0 ) = 1 into Equation (16), we have Equation (28) as required. □
Example 2.
When n = 2 , Equation (17) is just
t 2 y ( t ) + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where t R . Various values of a 1 and a 0 lead to various types of solutions as mentioned in Theorem 1.
If a 1 , a 0 , k are chosen according to Equation (18), for example a 1 = 3 , a 0 = 1 and k = 1 , then Equation (29) becomes
t 2 y ( t ) + 3 t y ( t ) + y ( t ) = 0 ,
and Equation (25) implies that the distributional solution of Equation (30) is y ( t ) = δ ( t ) . Using Equation (13), we can verify easily that Equation (30) is true for y ( t ) = δ ( t ) .
If a 1 , a 0 , k are chosen according to Equation (20), for example a 1 = 7 / 4 , a 0 = 45 / 4 and k = 1 , then Equation (29) becomes
t 2 y ( t ) + 7 4 t y ( t ) 45 4 y ( t ) = 0 ,
and Equation (27) implies that the classical solution of Equation (31) is y ( t ) = t 3 H ( t ) / 3 ! . We can verify easily that Equation (31) is true for y ( t ) = t 3 H ( t ) / 3 ! .
Example 3.
When n = 3 , Equation (17) becomes
t 3 y ( t ) + a 2 t 2 y ( t ) + a 1 t y ( t ) + a 0 y ( t ) = 0 ,
where t R . Various values of a 2 , a 1 and a 0 lead to various types of solutions, as mentioned in Theorem 1.
If a 2 , a 1 , a 0 , k are chosen according to Equation (18), for example a 2 = 1 / 12 , a 1 = 2 / 3 , a 0 = 61 and k = 2 , then Equation (32) becomes
t 3 y ( t ) + 1 12 t 2 y ( t ) + 2 3 t y ( t ) + 61 y ( t ) = 0 ,
and Equation (25) implies that the distributional solution of Equation (33) is y ( t ) = δ ( t ) . Using Equation (13), we can verify easily that Equation (33) is true y ( t ) = δ ( t ) .
If a 2 , a 1 , a 0 , k are chosen according to Equation (19), for example a 2 = 1 , a 1 = 2 / 3 , a 0 = 2 / 3 and k = 2 , then Equation (32) becomes
t 3 y ( t ) + t 2 y ( t ) 2 3 t y ( t ) 2 3 y ( t ) = 0 ,
and Equation (27) implies that the weak solution of Equation (34) is y ( t ) = t 2 H ( t ) / 2 ! . Using Equations (12) and (13), we can verify easily that Equation (34) is true for y ( t ) = t 2 H ( t ) / 2 ! .
Remark 1.
If n = 3 , a 1 = a 2 = a 3 = 1 and a 0 = m , then Theorem 1 is reduced to the condition in [13].
Remark 2.
If n = 5 , then Theorem 1 is reduced to the case of the fifth order Cauchy–Euler equation, as appears in theorem 3.1 of [16].
Remark 3.
If a 1 = a 2 = = a n = 1 and a 0 = m , then Theorem 1 is identical to the condition in [15].

Author Contributions

All authors contributed equally to this article. They read and approved the final manuscript.

Funding

This research received no external funding.

Acknowledgments

The second author was financially supported by the National Research Council of Thailand and Faculty of Science, Khon Kaen University 2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boyce, W.E.; DiPrima, R.C. Elementary Differential Equations and Boundary Value Problems, 9th ed.; John Wiley & Sons: New York, NY, USA, 2008. [Google Scholar]
  2. Kohler, W.; Johnson, L. Elementary Differential Equations, 2nd ed.; Addison-Wesley: Boston, MA, USA, 2006. [Google Scholar]
  3. Kanwal, R.P. Generalized Functions: Theory and Technique, 3rd ed.; Springer Science & Business Media: New York, NY, USA, 2004. [Google Scholar]
  4. Wiener, J. Generalized-function solutions of differential and functional differential equations. J. Math. Anal. Appl. 1982, 88, 170–182. [Google Scholar] [CrossRef] [Green Version]
  5. Schwartz, L. Theorie des Distributions; Actualite’s Scientifiques et Industrial: Paris, France, 1959. [Google Scholar]
  6. Zemanian, A.H. Distribution Theory and Transform Analysis; McGraw-Hill: New York, NY, USA, 1965. [Google Scholar]
  7. Wiener, J.; Shah, S.M. Distributional and entire solutions of ordinary differential and functional differential equations. Int. J. Math. Sci. 1983, 6, 243–270. [Google Scholar]
  8. Cooke, K.L.; Wiener, J. Distributional and analytic solutions of functional differential equations. J. Math. Anal. Appl. 1984, 98, 111–129. [Google Scholar] [CrossRef] [Green Version]
  9. Littlejohn, L.L.; Kawal, R.P. Distributional solutions of the hypergeometric differential equation. J. Math. Anal. Appl. 1987, 122, 325–345. [Google Scholar] [CrossRef] [Green Version]
  10. Wiener, J.; Cooke, K.L. Coexistence of analytic and distributional solutions for linear differential equations I. J. Math. Anal. Appl. 1990, 148, 390–421. [Google Scholar] [CrossRef]
  11. Wiener, J.; Cooke, K.L.; Shah, S.M. Coexistence of analytic and distributional solutions for linear differential equations II. J. Math. Anal. Appl. 1991, 159, 271–289. [Google Scholar] [CrossRef]
  12. Hernandezurena, L.G.; Estrada, R. Solution of ordinary differential equations by series of delta functions. J. Math. Anal. Appl. 1995, 191, 40–55. [Google Scholar] [CrossRef]
  13. Kananthai, A. Distribution solutions of the third order Euler equation. Southeast Asian Bull. Math. 1999, 23, 627–631. [Google Scholar]
  14. Liangprom, A.; Nonlaopon, K. On the generalized solutions of a certain fourth order Euler equations. J. Nonlinear Sci. Appl. 2017, 10, 4077–4084. [Google Scholar] [CrossRef] [Green Version]
  15. Sangsuwan, A.; Nonlaopon, K.; Orankitjaroen, S. The generalized solutions of a certain nth order Cauchy–Euler equation. Asian-Eur. J. Math. 2018, 13. [Google Scholar] [CrossRef]
  16. Jodnok, P.; Nonlaopon, K. On the generalized Solutions of the fifth order Euler equations. Far East J. Math. Sci. 2018, 106, 59–74. [Google Scholar] [CrossRef]
  17. Nanta, S. Distribution Solution of the nth Order Euler Equation; Chiang Mai University: Chiang Mai, Thailand, 2015. [Google Scholar]
  18. Bhattacharyya, P.K. Distributions: Generalized Functions with Applications in Sobolev Spaces; Walter de Gruyter: Berlin, Germany, 2012. [Google Scholar]
  19. Liverman, T.P.G. Generalized Functions and Direct Operational Methods; Prentice Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]

Share and Cite

MDPI and ACS Style

Sangsuwan, A.; Nonlaopon, K.; Orankitjaroen, S.; Mirumbe, I. The Generalized Solutions of the nth Order Cauchy–Euler Equation. Mathematics 2019, 7, 932. https://doi.org/10.3390/math7100932

AMA Style

Sangsuwan A, Nonlaopon K, Orankitjaroen S, Mirumbe I. The Generalized Solutions of the nth Order Cauchy–Euler Equation. Mathematics. 2019; 7(10):932. https://doi.org/10.3390/math7100932

Chicago/Turabian Style

Sangsuwan, Amornrat, Kamsing Nonlaopon, Somsak Orankitjaroen, and Ismail Mirumbe. 2019. "The Generalized Solutions of the nth Order Cauchy–Euler Equation" Mathematics 7, no. 10: 932. https://doi.org/10.3390/math7100932

APA Style

Sangsuwan, A., Nonlaopon, K., Orankitjaroen, S., & Mirumbe, I. (2019). The Generalized Solutions of the nth Order Cauchy–Euler Equation. Mathematics, 7(10), 932. https://doi.org/10.3390/math7100932

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop