Next Article in Journal
Trapezium-Type Inequalities for an Extension of Riemann–Liouville Fractional Integrals Using Raina’s Special Function and Generalized Coordinate Convex Functions
Next Article in Special Issue
Inverse Problem for a Mixed Type Integro-Differential Equation with Fractional Order Caputo Operators and Spectral Parameters
Previous Article in Journal
Strong Convergence of Extragradient-Type Method to Solve Pseudomonotone Variational Inequalities Problems
Previous Article in Special Issue
Asymptotic Properties of Neutral Differential Equations with Variable Coefficients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Series of Distributional Solutions for Certain Linear Differential Equations

by
Nipon Waiyaworn
1,
Kamsing Nonlaopon
1,* and
Somsak Orankitjaroen
2
1
Department of Mathematics, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand
*
Author to whom correspondence should be addressed.
Axioms 2020, 9(4), 116; https://doi.org/10.3390/axioms9040116
Submission received: 31 August 2020 / Revised: 30 September 2020 / Accepted: 2 October 2020 / Published: 13 October 2020

Abstract

:
In this paper, we present the distributional solutions of the modified spherical Bessel differential equations t 2 y ( t ) + 2 t y ( t ) [ t 2 + ν ( ν + 1 ) ] y ( t ) = 0 and the linear differential equations of the forms t 2 y ( t ) + 3 t y ( t ) ( t 2 + ν 2 1 ) y ( t ) = 0 , where ν N { 0 } and t R . We find that the distributional solutions, in the form of a finite series of the Dirac delta function and its derivatives, depend on the values of ν . The results of several examples are also presented.

1. Introduction

It is well known that the linear differential equation of the form
n = 0 m a n ( t ) y ( n ) ( t ) = 0 , a m ( t ) 0 ,
where a n ( t ) is an infinitely smooth coefficient for each n, and has no distributional solutions other than the classical ones. However, if the leading coefficient a m ( t ) has a zero, the classical solution of (1) may cease to exist in a neighborhood of that zero. In that case, (1) may have a distributional solution. It was not until 1982 that Wiener [1] proposed necessary and sufficient conditions for the existence of an Nth-order distributional solution to the differential equation (1). The Nth-order distributional solution that Wiener proposed is a finite sum of Dirac delta function and its derivatives:
y ( t ) = n = 0 N b n δ ( n ) ( t ) , b N 0 .
It can be easily verified by (10) that δ ( t ) is a zero order distributional solution of the equation
t y ( t ) + 2 y ( t ) + t y ( t ) = 0 ;
the Bessel equation
t 2 y ( t ) + t y ( t ) + ( t 2 1 ) y ( t ) = 0 ;
the confluent hypergeometric equation
t y ( t ) + ( 2 t ) y ( t ) y ( t ) = 0 ;
and the second order Cauchy–Euler equation
t 2 y ( t ) + 3 t y ( t ) + y ( t ) = 0 .
The distributional solutions with higher order of Cauchy–Euler equations were studied by many researchers; see [2,3,4,5,6,7,8] for more details.
The infinite order distributional solution of the form
y ( t ) = n = 0 b n δ ( n ) ( t )
to various differential equations in a normal form with singular coefficients was studied by many researchers [9,10,11,12,13]. Furthermore, a brief introduction to these concepts is presented by Kanwal [14].
In 1984, Cooke and Wiener [15] presented the existence theorems for distributional and analytic solutions of functional differential equations. In 1987, Littlejohn and Kanwal [16] studied the distributional solutions of the hypergeometric differential equation, whose solutions are in the form of (3). In 1990, Wiener and Cooke [17] presented the necessary and sufficient conditions for the simultaneous existence of solutions to linear ordinary differential equations in the forms of rational functions and (2).
As mentioned in abstract, we propose the distributional solutions of the modified spherical Bessel differential equations
t 2 y ( t ) + 2 t y ( t ) [ t 2 + ν ( ν + 1 ) ] y ( t ) = 0
and the linear differential equations of the forms
t 2 y ( t ) + 3 t y ( t ) ( t 2 + ν 2 1 ) y ( t ) = 0 ,
where ν N { 0 } and t R . The modified spherical Bessel differential equation is just the spherical Bessel equation with a negative separation constant. The spherical Bessel equation occurs when dealing with the Helmholtz equation in spherical coordinates of various problems in physics such as a scattering problem [18].
We use the simple method, consisting of Laplace transforms of right-sided distributions and power series solution, for searching the distributional solutions of these equations. We find that the solutions are in the forms of finite linear combinations of the Dirac delta function and its derivatives depending on the values of ν .

2. Preliminaries

In this section, we introduce the basic knowledge and concepts, which are essential for this work.
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 . The space of all such distributions is denoted by D .
For every T D and φ ( t ) D , the value that T acts 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 ) . Then we define f ( t ) , φ ( t ) = Ω f ( t ) φ ( t ) d t , where Ω is the support of φ ( t ) and φ ( t ) D .
(ii) 
The Dirac delta function is a distribution defined by δ ( t ) , φ ( t ) = φ ( 0 ) and the support of δ ( t ) is { 0 } .
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 ( k ) , φ ( t ) = ( 1 ) k T , φ ( k ) ( t ) for all φ ( t ) D .
Example 2.
(i) 
δ ( t ) , φ ( t ) = δ ( t ) , φ ( t ) = φ ( 0 ) ;
(ii) 
δ ( k ) ( t ) , φ ( t ) = ( 1 ) k δ ( t ) , φ ( k ) ( t ) = ( 1 ) k φ ( k ) ( 0 ) .
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.
If y ( t ) is a singular distribution and satisfies the equation
m = 0 n a m ( t ) y ( n ) ( t ) = f ( t ) ,
where a m ( t ) is an infinitely differentiable function and f ( t ) is an arbitrary known distribution, in the sense of distribution, and is called a distributional solution of (4).
Definition 6.
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 c t f ( t ) is absolutely integrable over R .
The Laplace transform of f ( t ) is defined by
F ( s ) = L f ( t ) = M f ( t ) e s t d t ,
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 ) .
Recall that the Laplace transform G ( s ) of a locally integrable function g ( t ) satisfying the conditions of definition 6, that is,
G ( s ) = L g ( t ) = M g ( t ) e s t d t ,
where ( s ) > σ a , can be written in the form G ( s ) = g ( t ) , e s t .
Definition 7.
Let S be the space of test functions of rapid decay containing the complex-valued functions ϕ ( t ) having the following properties:
(i) 
ϕ ( t ) is infinitely differentiable—i.e., ϕ ( t ) C ( R ) ;
(ii) 
ϕ ( t ) , as well as its derivatives of all orders, vanish at infinity faster than the reciprocal of any polynomial which is expressed by the inequality
| t p ϕ ( k ) ( t ) | < C p k ,
where C p k 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 test function of rapid decay and contains the complex-valued functions—i.e., there is assigned a complex number T , ϕ ( t ) with properties:
(i) 
T , c 1 ϕ 1 ( t ) + c 2 ϕ 2 ( t ) = c 1 T , ϕ 1 ( t ) + c 2 T , ϕ 2 ( t ) for ϕ 1 ( t ) , ϕ 2 ( t ) S and constants c 1 , c 2 ;
(ii) 
lim m T , ϕ m ( t ) = 0 for every null sequence { ϕ m ( t ) } S .
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:
(i) 
f ( t ) is a right-sided distribution, that is, f ( t ) D R .
(ii) 
There exists a real number c such that 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 c 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 1 over a neighbourhood 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 c t f ( t ) is in the space S . Equation (7) can be reduced to
F ( s ) = L f ( t ) = f ( t ) , e s t .
Now F ( s ) is a function of s defined over the right half-plane ( s ) > c . Zemanian [19] proved that F ( s ) is an analytic function in the region of convergence ( s ) > σ 1 , where σ 1 is the abscissa of convergence and e c t 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 holds:
(i) 
L { ( t k 1 H ( t ) ) / ( k 1 ) ! } = 1 / s k , ( s ) > 0 ;
(ii) 
L δ ( t ) = 1 , < ( s ) < ;
(iii) 
L δ ( k ) ( t ) = s k , < ( s ) < ;
(iv) 
L t k f ( t ) = ( 1 ) k F ( k ) ( s ) , ( s ) > σ 1 ;
(v) 
L f ( k ) ( t ) = s k F ( s ) , ( s ) > σ 1 .
The proof of following lemma 1 is given in [14].
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 1 ) ( 0 ) δ ( t ) + + ψ ( 0 ) δ ( m ) ( t ) .
A useful formula that follows from (9), for any monomial ψ ( t ) = t n , is that
t n δ ( m ) ( t ) = 0 , if m < n ; ( 1 ) n m ! ( m n ) ! δ ( m n ) ( t ) , if m n .

3. Main Results

In this section, we will state our main results and give their proofs.
Theorem 1.
Consider the differential equation of the form
t 2 y ( t ) + 2 t y ( t ) [ t 2 + ν ( ν + 1 ) ] y ( t ) = 0 ,
where ν N { 0 } and t R . The distributional solutions of (11) are given by
y ( t ) = P ν ( D ) δ ( t ) ,
where
P ν ( D ) = 1 2 ν k = 0 ν / 2 ( 1 ) k ( 2 ν 2 k ) ! k ! ( ν k ) ! ( ν 2 k ) ! D ν 2 k ,
is a Legendre polynomial of distributional derivative operator D = d / d t .
Proof. 
Applying the Laplace transform to both sides of (11) with L { y ( t ) } = F ( s ) , and using Example 3(iv), (v), we obtain
( 1 s 2 ) F ( s ) 2 s F ( s ) + ν ( ν + 1 ) F ( s ) = 0 .
Suppose that a solution of (13) is of the form F ( s ) = n = 0 a n s n . Differentiating term by term, we obtain
F ( s ) = n = 1 n a n s n 1
and
F ( s ) = n = 2 n ( n 1 ) a n s n 2 .
Substituting these terms into (13), we have
2 a 2 + ν ( ν + 1 ) a 0 + ( 3 · 2 ) a 3 ( 2 ν ( ν + 1 ) a 1 s + n = 2 ( n + 2 ) ( n + 1 ) a n + 2 n ( n 1 ) a n 2 n a n + ν ( ν + 1 ) a n s n = 0 .
Since s n 0 , it follows that
2 a 2 + ν ( ν + 1 ) a 0 = 0 , ( 3 · 2 ) a 3 ( 2 ν ( ν + 1 ) a 1 = 0 , ( n + 2 ) ( n + 1 ) a n + 2 n ( n 1 ) a n 2 n a n + ν ( ν + 1 ) a n = 0 , n 2 ,
which leads to a recurrence relation
a n + 2 = ( ν n ) ( ν + n + 1 ) ( n + 1 ) ( n + 2 ) a n .
Thus, we obtain
a 2 = ν ( ν + 1 ) 2 ! a 0 a 4 = ( 1 ) 2 ν ( ν 2 ) ( ν + 1 ) ( ν + 3 ) 4 ! a 0 a 2 n = ( 1 ) n ν ( ν 2 ) ( ν 2 n + 2 ) ( ν + 1 ) ( ν + 3 ) ( ν + 2 n 1 ) ( 2 n ) ! a 0 .
Similarly,
a 3 = ( ν 1 ) ( ν + 2 ) 2 · 3 a 1 a 5 = ( 1 ) 2 ( ν 1 ) ( ν 3 ) ( ν + 2 ) ( ν + 4 ) 5 ! a 1 a 2 n + 1 = ( 1 ) n ( ν 1 ) ( ν 3 ) ( ν 2 n + 1 ) ( ν + 2 ) ( ν + 4 ) ( ν + 2 n ) ( 2 n + 1 ) ! a 1 .
Letting a 0 = a 1 = 1 , we get the two solutions of (13) in the forms
F e ( s ) = 1 + n = 1 ( 1 ) n ν ( ν 2 ) ( ν 2 n + 2 ) ( ν + 1 ) ( ν + 3 ) ( ν + 2 n 1 ) ( 2 n ) ! s 2 n
and
F o ( s ) = s + n = 1 ( 1 ) n ( ν 1 ) ( ν 3 ) ( ν 2 n + 1 ) ( ν + 2 ) ( ν + 4 ) ( ν + 2 n ) ( 2 n + 1 ) ! s 2 n + 1 .
If ν is even, letting ν = 2 m , m N { 0 } , we note that
ν ( ν 2 ) ( ν 2 n + 2 ) = 2 m ( 2 m 2 ) ( 2 m 2 n + 2 ) = 0 , m = 0 ; 2 n m ! ( m n ) ! , m > 0 , n m ; 0 , n > m > 0 ,
and
( ν + 1 ) ( ν + 3 ) ( ν + 2 n 1 ) = ( 2 m + 1 ) ( 2 m + 3 ) ( 2 m + 2 n 1 ) = ( 2 m + 2 n ) ! m ! 2 n ( 2 m ) ! ( m + n ) ! .
Then, in this case, F e ( s ) only becomes the finite series of the form
F e ( s ) = ( m ! ) 2 ( 2 m ) ! k = 0 m ( 1 ) k ( 2 m + 2 k ) ! s 2 k ( m k ) ! ( m + k ) ! ( 2 k ) ! .
If ν is odd, letting ν = 2 m + 1 , m N { 0 } , we note that
( ν 1 ) ( ν 3 ) ( ν 2 n + 1 ) = 2 m ( 2 m 2 ) ( 2 m 2 n + 2 ) = 0 , m = 0 ; 2 n m ! ( m n ) ! , m > 0 , n m ; 0 , n > m > 0 ,
and
( ν + 2 ) ( ν + 4 ) ( ν + 2 n ) = ( 2 m + 3 ) ( 2 m + 5 ) ( 2 m + 2 n + 1 ) = ( 2 m + 2 n + 1 ) ! m ! 2 n ( 2 m + 1 ) ! ( m + n ) ! .
Then, in this case, F o ( s ) only becomes the finite series of the form
F o ( s ) = ( m ! ) 2 ( 2 m + 1 ) ! k = 0 m ( 1 ) k ( 2 m + 2 k + 1 ) ! s 2 k + 1 ( m k ) ! ( m + k ) ! ( 2 k + 1 ) ! .
For ν = 0 , 1 , 2 , , we have F ν ( s ) , as follows:
F 0 ( s ) = 1 = P 0 ( s ) , F 1 ( s ) = s = P 1 ( s ) , F 2 ( s ) = 1 3 s 2 = 2 P 2 ( s ) , F 3 ( s ) = s ( 5 / 3 ) s 3 = ( 2 / 3 ) P 3 ( s ) , F 4 ( s ) = 1 10 s 2 + ( 35 / 3 ) s 4 = ( 8 / 3 ) P 4 ( s ) , F 5 ( s ) = s ( 14 / 3 ) s 3 + ( 21 / 5 ) s 5 = ( 8 / 15 ) P 5 ( s ) , F 6 ( s ) = 1 21 s 2 + 63 s 4 ( 231 / 5 ) s 6 = ( 16 / 5 ) P 6 ( s ) , F 7 ( s ) = s 9 s 3 + ( 99 / 5 ) s 5 ( 429 / 35 ) s 7 = ( 16 / 35 ) P 7 ( s ) ,
and so on, where P n ( s ) is the Legendre polynomial of s for n = 0 , 1 , 2 , . Since (13) is linear, P ν is also its solution for all non-negative integer ν . Taking the inverse Laplace transform to P ν ( s ) , and using Example 3(ii),(iii), we obtain the solutions of (11),
y ( t ) = 1 2 ν k = 0 ν / 2 ( 1 ) k ( 2 ν 2 k ) ! k ! ( ν k ) ! ( ν 2 k ) ! D ν 2 k δ ( t ) , ( D d d t distribution derivative ) = P ν ( D ) δ ( t ) ,
which are the distributional solutions of the form (12). □
Example 4.
Letting ν = 1 , (11) becomes
t 2 y ( t ) + 2 t y ( t ) ( t 2 + 2 ) y ( t ) = 0 .
From Theorem 1, (18) has a solution
y ( t ) = δ ( t ) .
Letting ν = 4 , (11) becomes
t 2 y ( t ) + 2 t y ( t ) ( t 2 + 20 ) y ( t ) = 0 .
From Theorem 1, (20) has a solution
y ( t ) = 35 8 δ ( 4 ) ( t ) 15 4 δ ( t ) + 3 8 δ ( t ) .
By applying (10), it is not difficult to verify that (19) and (21) satisfy (18) and (20), respectively.
Theorem 2.
Consider the equation of the form
t 2 y ( t ) + 3 t y ( t ) ( t 2 + ν 2 1 ) y ( t ) = 0 ,
where ν N { 0 } and t R . The distributional solutions of (22) are given by
y ( t ) = T ν ( D ) δ ( t ) ,
where
T ν ( D ) = ν 2 k = 0 ν / 2 ( 1 ) k 2 ν 2 k ( n k 1 ) ! k ! ( ν 2 k ) ! D ν 2 k ,
is a Chebyshev polynomial of the first kind of distributional derivative operator D = d / d t .
Proof. 
Applying the Laplace transform L { y ( t ) } = F ( s ) to (22) and using Example 3(iv), (v), we obtain
( 1 s 2 ) F ( s ) s F ( s ) + ν 2 F ( s ) = 0 .
Suppose that a solution of (24) is of the form F ( s ) = n = 0 a n s n . Differentiating F ( s ) term by term, we obtain
F ( s ) = n = 1 n a n s n 1
and
F ( s ) = n = 2 n ( n 1 ) a n s n 2 .
Substituting these terms into (24), we get
2 a 2 + ν 2 a 0 + 6 a 3 a 1 + ν 2 a 1 s + n = 2 ( n + 2 ) ( n + 1 ) a n + 2 [ n ( n 1 ) + n ν 2 ] a n s n = 0 .
Since s n 0 , it follows that
2 a 2 + ν 2 a 0 = 0 , 6 a 3 a 1 + ν 2 a 1 = 0 ,
and for n = 2 , 3 , ,
( n + 2 ) ( n + 1 ) a n + 2 ( n 2 ν 2 ) a n = 0 .
Hence,
a 2 = ν 2 2 a 0 , a 3 = 1 ν 2 3 ! a 1 ,
and for n = 2 , 3 , ,
a n + 2 = n 2 ν 2 ( n + 2 ) ( n + 1 ) a n .
Thus, we obtain
a 2 = ν 2 2 a 0 , a 4 = ( 2 2 ν 2 ) ( ν 2 ) 4 ! a 0 ,
and so on. Similarly,
a 3 = 1 ν 2 3 ! a 1 , a 5 = ( 3 2 ν 2 ) ( 1 ν 2 ) 5 ! a 1 ,
and so on. A pattern clearly emerges:
a 2 n = [ ( 2 n 2 ) 2 ν 2 ] [ ( 2 n 4 ) 2 ν 2 ] ( 2 2 ν 2 ) ( ν 2 ) ( 2 n ) ! a 0
and
a 2 n + 1 = [ ( 2 n 1 ) 2 ν 2 ] [ ( 2 n 3 ) 2 ν 2 ] ( 3 2 ν 2 ) ( 1 ν 2 ) ( 2 n + 1 ) ! a 1 .
Hence, we get two solutions of (24) in the forms
F e ( s ) = a 0 + n = 1 [ ( 2 n 2 ) 2 ν 2 ] [ ( 2 n 4 ) 2 ν 2 ] ( 2 2 ν 2 ) ( ν 2 ) ( 2 n ) ! a 0 s 2 n
and
F o ( s ) = a 1 s + n = 1 [ ( 2 n 1 ) 2 ν 2 ] [ ( 2 n 3 ) 2 ν 2 ] ( 3 2 ν 2 ) ( 1 ν 2 ) ( 2 n + 1 ) ! a 1 s 2 n + 1 .
If ν is even, letting ν = 2 m , m N { 0 } , then a 2 m + 2 = 0 , so that a 2 n + 2 = 0 for n m . Hence,
F e ( s ) = 1 + n = 1 m [ ( 2 n 2 ) 2 ( 2 m ) 2 ] [ ( 2 n 4 ) 2 ( 2 m ) 2 ] [ 2 2 ( 2 m ) 2 ] [ ( 2 m ) 2 ] ( 2 n ) ! s 2 n ,
equivalently to
F e ( s ) = 1 + n = 1 m k = 1 n 4 ( n k + m ) ( n k m ) ( 2 n ) ! s 2 n .
If ν is odd, letting ν = 2 m + 1 , m N { 0 } , then a 2 m + 3 = 0 , so that a 2 n + 3 = 0 for n m . Hence,
F o ( s ) = s + n = 1 m [ ( 2 n 1 ) 2 ( 2 m + 1 ) 2 ] [ ( 2 n 3 ) 2 ( 2 m + 1 ) 2 ] [ 3 2 ( 2 m + 1 ) 2 ) ] [ 1 ( 2 m + 1 ) 2 ] ( 2 n + 1 ) ! s 2 n + 1 ,
equivalently to
F o ( s ) = s + n = 1 m k = 1 n 4 ( n k + m + 1 ) ( n k m ) ( 2 n + 1 ) ! s 2 n + 1 .
For ν = 0 , 1 , 2 , , we have F ν ( s ) , as follows:
F 0 ( s ) = 1 = T 0 ( s ) , F 1 ( s ) = s = T 1 ( s ) , F 2 ( s ) = 1 2 s 2 = T 2 ( s ) , F 3 ( s ) = s ( 4 / 3 ) s 3 = ( 1 / 3 ) T 3 ( s ) , F 4 ( s ) = 1 8 s 2 + 8 s 4 = T 4 ( s ) , F 5 ( s ) = s 4 s 3 + ( 16 / 5 ) s 5 = ( 1 / 5 ) T 5 ( s ) , F 6 ( s ) = 1 18 s 2 + 48 s 4 32 s 6 = T 6 ( s ) , F 7 ( s ) = s 8 s 3 + 16 s 5 ( 64 / 7 ) s 7 = ( 1 / 7 ) T 7 ( s ) ,
and so on, where T n ( s ) is the Chebyshev polynomial of the first kind of s for n = 0 , 1 , 2 , . Since (24) is linear, T n ( s ) is also its solution. Taking the inverse Laplace transform to T n ( s ) , and using Example 3(ii),(iii), we obtain the solutions of (22),
y n ( t ) = ν 2 k = 0 ν / 2 ( 1 ) k 2 ν 2 k ( n k 1 ) ! k ! ( ν 2 k ) ! D ν 2 k δ ( t ) , = T n ( D ) δ ( t ) ,
which are the distributional solutions of the forms (23). □
Example 5.
Letting ν = 1 , (22) becomes
t 2 y ( t ) + 3 t y ( t ) t 2 y ( t ) = 0 .
From Theorem 2, (27) has a solution
y ( t ) = δ ( t ) .
Letting ν = 4 , (22) becomes
t 2 y ( t ) + 2 t y ( t ) ( t 2 + 15 ) y ( t ) = 0 .
From Theorem 2, (29) has a solution
y ( t ) = 8 δ ( 4 ) ( t ) 8 δ ( t ) + δ ( t ) .
By applying (10), it is not difficult to verify that (28) and (30) satisfy (27) and (29), respectively.

4. Conclusions

In this paper, we seek the distributional solutions of the modified spherical Bessel differential Equation (11) and the linear differential equation of the form (22) by using the Laplace transforms of right-sided distributions and the power series solutions. The obtained solutions in the forms of the finite linear combinations of the Dirac delta function and its derivatives depend on the value of ν , to which their coefficients regard the coefficients of Legendre and Chebyshev polynomials (see [20] for more details). However, for solutions of (11) and (22) in the usual sense, not mentioned here, they can be seen in many standard and technical textbooks (see, for example, Ross [21]) but, even more, may appear in models related to equilibrium of membrane structures, steady states of evolutive equations or nonlinear science (see studies [22,23,24,25]).

Author Contributions

All authors contributed equally to this article. They read and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This work was supported by the Research Fund for Supporting Lecturer to Admit High Potential Student to Study and Research on His Expert Program Year 2018 from the Graduate School, Khon Kaen University, Thailand (Grant no. 611T105).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wiener, J. Generalized-function solutions of differential and functional differential equations. J. Math. Anal. Appl. 1982, 88, 170–182. [Google Scholar] [CrossRef] [Green Version]
  2. Kananthai, A. Distribution solutions of the third order Euler equation. Southeast Asian Bull. Math. 1999, 23, 627–631. [Google Scholar]
  3. 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]
  4. Sacorn, N.; Nonlaopon, K.; Kim, H. A note on the generalized solutions of the third-order Cauchy-Euler equations. Commun. Math. Appl. 2018, 9, 661–669. [Google Scholar]
  5. 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]
  6. Mirumbe, G.I.; Mango, J.M. On generalized solutions of locally Fuchsian ordinary differential equations. J. Math. Sci. Adv. Appl. 2018, 51, 99–117. [Google Scholar] [CrossRef]
  7. Jhanthanam, S.; Nonlaopon, K.; Orankitjaroen, S. Generalized solutions of the third-order Cauchy-Euler equation in the space of right-sided distributions via Laplace transform. Mathematics 2019, 7, 376. [Google Scholar] [CrossRef] [Green Version]
  8. Sangsuwan, A.; Nonlaopon, K.; Orankitjaroen, S.; Mirumbe, I. The generalized solutions of the nth order Cauchy-Euler equation. Mathematics 2019, 7, 932. [Google Scholar] [CrossRef] [Green Version]
  9. Wiener, J. Generalized-function solutions of linear systems. J. Differ. Equ. 1980, 38, 301–315. [Google Scholar] [CrossRef] [Green Version]
  10. 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] [Green Version]
  11. Hernandez, L.G.; Estrada, R. Solution of ordinary by series of delta functions. J. Math. Anal. Appl. 1995, 191, 40–55. [Google Scholar] [CrossRef]
  12. Kanwal, R.P. Delta series solutions of differential and integral equations. Integral Transform. Spec. Funct. 1998, 6, 49–62. [Google Scholar] [CrossRef]
  13. Nonlaopon, K.; Nuntigrangjana, T.; Putjuso, S. Distributional solutions of nth-order differential equations of the Bessel equation. RMUTSB Acad. J. 2017, 5, 1–10. [Google Scholar]
  14. Kanwal, R.P. Generalized Functions: Theory and Technique, 3rd ed.; Springer: New York, NY, USA, 2004. [Google Scholar]
  15. 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]
  16. Littlejohn, L.L.; Kanwal, R.P. Distributional solutions of the hypergeometric differential equation. J. Math. Anal. Appl. 1987, 122, 325–345. [Google Scholar] [CrossRef] [Green Version]
  17. 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] [Green Version]
  18. Boyvel, L.P.; Jones, A.R. Electromagnetic Scattering and Its Applications; Applied Science Publishers: London, UK, 1981. [Google Scholar]
  19. Zemanian, A.H. Distribution Theory and Transform Analysis; McGraw-Hill: New York, NY, USA, 1965. [Google Scholar]
  20. Wiener, J. Generalized Solutions of Functional Differential Equations; World Scientific: Singapore, 1993. [Google Scholar]
  21. Ross, S.L. Differential Equations, 3rd ed.; John Wiley & Sons Inc.: Singapore, 1984. [Google Scholar]
  22. Li, T.; Viglialoro, G. Analysis and explicit solvability of degenerate tensorial problems. Bound. Value Probl. 2018, 1, 1–13. [Google Scholar] [CrossRef] [Green Version]
  23. Viglialoro, G.; Woolley T., E. Boundedness in a parabolic-elliptic chemotaxis system with nonlinear diffusion and sensitivity and logistic source. Math. Methods Appl. Sci. 2018, 41, 1809–1824. [Google Scholar] [CrossRef]
  24. Az-Zo´bi, E.A.; Al-Khaled, K.; Darweesh, A. Numeric-analytic solutions for nonlinear oscillators via the modified multi-stage decomposition method. Mathematics 2019, 7, 550. [Google Scholar] [CrossRef] [Green Version]
  25. Gonza´lez-Gaxiola, O.; Santiago, J.A.; Ruiz de Cha´vez, J. Solution for the nonlinear relativistic harmonic oscillator via Laplace-Adomian decomposition method. Int. J. Appl. Comput. Math. 2017, 3, 2627–2638. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Waiyaworn, N.; Nonlaopon, K.; Orankitjaroen, S. Finite Series of Distributional Solutions for Certain Linear Differential Equations. Axioms 2020, 9, 116. https://doi.org/10.3390/axioms9040116

AMA Style

Waiyaworn N, Nonlaopon K, Orankitjaroen S. Finite Series of Distributional Solutions for Certain Linear Differential Equations. Axioms. 2020; 9(4):116. https://doi.org/10.3390/axioms9040116

Chicago/Turabian Style

Waiyaworn, Nipon, Kamsing Nonlaopon, and Somsak Orankitjaroen. 2020. "Finite Series of Distributional Solutions for Certain Linear Differential Equations" Axioms 9, no. 4: 116. https://doi.org/10.3390/axioms9040116

APA Style

Waiyaworn, N., Nonlaopon, K., & Orankitjaroen, S. (2020). Finite Series of Distributional Solutions for Certain Linear Differential Equations. Axioms, 9(4), 116. https://doi.org/10.3390/axioms9040116

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