Next Article in Journal
Information Measures Based on T-Spherical Fuzzy Sets and Their Applications in Decision Making and Pattern Recognition
Previous Article in Journal
Smooth Versions of the Mann–Whitney–Wilcoxon Statistics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Theorems in Solving Families of Improper Integrals

Department of Mathematics, Faculty of Science, Zarqa University, Zarqa 13110, Jordan
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(7), 301; https://doi.org/10.3390/axioms11070301
Submission received: 11 April 2022 / Revised: 22 May 2022 / Accepted: 24 May 2022 / Published: 21 June 2022

Abstract

:
Many improper integrals appear in the classical table of integrals by I. S. Gradshteyn and I. M. Ryzhik. It is a challenge for some researchers to determine the method in which these integrations are formed or solved. In this article, we present some new theorems to solve different families of improper integrals. In addition, we establish new formulas of integrations that cannot be solved by mathematical software such as Mathematica or Maple. In this article, we present three main theorems that are essential in generating new formulas for solving improper integrals. To show the efficiency and the simplicity of the presented techniques, we present some applications and examples on integrations that cannot be solved by regular methods. Furthermore, we acquire new results for integrations and compare them to that obtained in the classical table of integrations. Some previous results, become special cases of our outcomes or generalizations to acquire new integrals.

1. Introduction

In the history of mathematics, there are many improper integrals that appeared in the fields of science, physics, engineering and etc., [1,2,3,4,5]. Some of these integrals can be solved easily and others need some techniques to be solved. some of these integrals can’t be solved up to now, some of them are solved numerically, and there still some integrations that cannot be determined exactly or need much effort to be solved. The importance of computing improper integrals has come from the wide usage in applied math, physics, engineering and etc., [6,7,8,9,10,11,12,13,14,15].
Evaluation of Improper integrals is not typically subjected to general rules or assignable methods at first glance. Various ways, techniques and methods have been studied and presented by scientists and researchers for providing a closed-form expression for definite integrals. Series methods, double integrals technique, Cauchy residue theorem, differentiating under the integral sign, and plenty of other techniques are used in solving complicated improper integrals exactly or approximately, [16,17,18,19,20]. It is nearly always that one blindly looks for a way that leads to the answer, without certainty on the success.
Within these theorems, we present the closed form of the integral that can be found only by identifying the function that formed the improper integral. These theorems are considered a solid tool for unraveling some families of improper integrals and creating an infinite number of challenging integrals.
In 1826. A.L. Cauchy came up with the Residue Theorem, which is one of the most important achievements in complex analysis. Nevertheless, applications of the residue theorem to solve integrals over real line require rigorous conditions that must be met to solve the integrals, such as determining the appropriate closed contour, finding the poles, and also that the process of finding residues may be difficult and exhausting in some integrals. In his published memoir, Cauchy had come up with mathematical formulas using the residue theorem [21,22,23,24,25,26,27], the formulas that he achieved are considered very simple comparing to the results that will be presented; considering that the theorems in this article don’t depend on the residue theorem.
Another important achievement in the world of improper integrals, is Ramanujan’s master theorem, which provides an explicit expression for the Mellin transform of a function in terms of the analytic continuation of its Taylor coefficients [16]. It was widely used by Ramanujan as a tool in computing definite integrals and infinite series. The theorems presented in this article are just as powerful and efficient as Ramanujan’s master theorem in solving some families of improper integrals, for more about integral transforms, see [5,6,7,8].
This article is organized as follows, in the second section, we present some basic definitions and theorems that are needed in our research, Main results and theorems are presented in section three. Some applications and examples are presented in section four. Finally, some mathematical remarks are presented in section five.

2. Basic Definitions and Theorems

Before going into our results, we need to introduce some basic definitions and theorems, to be used in the construction of the new results.
Definition 1
([18]). Let f be analytic function in an open set Ω which is a subset of the complex plane, and D is a disc centered at z 0   subset of D ; whose closure is contained in the region Ω, then f has a power series expansion at z 0   as,
f ( z ) = n = 0 a n ( z z 0 ) n ,   n N * ,     a n   is   a   real   coefficient   of   the   series .    
Definition 2
([19]). Let  f be a real analytic function, that is, infinitely differentiable such that Taylor series at any point x 0 in it is domain is,
T ( x ) = n = 0 f ( n ) ( x 0 ) n !   ( x x 0 ) n ,
converges to f ( x ) in a neighborhood of x 0 pointwisely.
Definition 3
([20]). Cauchy principal value of a finite integral of a function f   about a point c , with a c b   is given by,
P V a b f ( x ) d x lim ε 0 + [ a c ε f ( x ) d x + c + ε b f ( x ) d x   ] .
Definition 4
((Gamma function) [21]). Euler integral of the second kind, is defined by the integral,
Γ ( z ) = 0 x z 1 e x d x ,   where     R e ( z ) > 0 ,  
and Euler’s reflection formula for the gamma function is,
  Γ ( z ) Γ ( 1 z ) = π sin ( π z ) ,       z
Definition 5
((Beta function) [21]). It is also called Euler integral of the first kind that is defined by the integral,
β ( p , y ) = 0 1 t p 1 ( 1 t ) y 1   d t ,
for complex numbers inputs p   a n d   y such that R e   p > 0   a n d   R e   y > 0 .
Integral (5) can be rewritten as,
β ( p , y ) = 0 t p 1 ( 1 + t ) p + y d t                   ; R e ( p ) ,   R e ( y ) > 0 .
There is a relation between the gamma function and the beta function as follows,
β ( x , y ) = Γ ( x ) Γ ( y ) Γ ( x + y ) .
In the following arguments, we present some formulas of integrations, that are essential in our work. Herein, we mention that these integrals could be evaluated by many methods, beta function, residue theorem and other techniques.
Lemma 1.
For any real number x and a constant p , we have the following integrals,
(i)
I 1 = 0 x p x 2 + 1   d x = π 2 cos ( π p 2 ) ,   where   1 < Re ( p ) < 1.
(ii)
I 2 = P V 0 x p 1 x 2 1 d x = π cot (   π p 2 ) 2 ,   where   0 < R e ( p ) < 2.
Proof of (i).
This is a direct proof that can be obtained from the relations between the gamma function and the beta function. To prove the first integral, let x 2 = t ,   then I 1   using the fact in Equation (6) becomes,
          I 1 = 1 2 β ( p + 1 2 , 1 p 2 ) .  
Using the relation between the beta and the gamma function in Equation (7), we acquire,
I 1 = 1 2 Γ ( p + 1 2 ) Γ ( 1 p 2 ) .
Now, using the reflection formula (4) of the the gamma function we have,
I 1 = π 2 cos ( π p 2 ) .
 □
Proof of (ii).
By factorizing the denominator of I 2 , and decompose the integral into partial fractions, we have,
  I 2 = P V 0 x p 1 2 ( x 1 ) x p 1 2 ( x + 1 ) d x = P V ( 0 x p 1 2 ( x 1 ) d x 0 x p 1 2 ( x + 1 ) d x ) .
Let M = 0 x p 1 x + 1 d x ,   and = P V 0 x p 1 ( x 1 ) .
M can be solved easily, by using Equation (6), followed by the relation between beta and gamma function in Equation (7), and then using the reflection formula (4) we acquire,
      M = 0 x p 1 x + 1 d x = π sin ( π p ) .      
The value of N can be obtained as follows, using the value in Equation (10), we acquire,
M 2 = 0 0 x p 1 y p 1 ( x + 1 ) ( y + 1 ) d x   d y .
Let y = z x d y = d z x , thus after simple computations and using the idea of interchanging the order of integration, we have,
M 2 = P V 0 x p 1 x + 1 0 x 1 p   z p 1 ( 1 + z x ) x d z     d x = P V 0 0 z p 1 ( x + 1 ) ( z + x ) d x   d z .
After simple computations, one can obtain,
M 2 = P V 0 z p 1 z 1 ln ( z ) d z ,
which can be re written as,
                                            M 2 = d d p 0 z p 1 z 1 ln ( z ) d z .
Thus,
M 2 d p = P V 0 z p 1 z 1 d z
Now, using the value of M in Equation (10), we have,
      N = P V 0 z p 1 z 1 d z = 0 π 2 sin 2 ( π p ) d p = π cot ( p   π ) .    
Thus, we the proof is complete and we have,
I 2 = π 2 cot ( p   π ) π 2 sin ( π p ) = π cot ( p π 2 ) 2 .
 □

3. Main Theorems

In this section, we introduce new attractive results for integrations, that can be used to generate improper integrals and solve difficult applications. To achieve our goal, we need to present some relations about analytic functions that can be found in [20,21,22].
Let f ( z ) = f ( β e θ x ) be analytic function around 0, then according to Taylor’s series where θ and β denote positive or negative real quantities, we have,
f ( z ) = k = 0 f ( k ) ( 0 ) k ! ( z ) k ,   and   f ( β e θ x ) = k = 0 f ( k ) ( 0 ) k ! β k e θ k x
Using the formula e i θ x + e i θ x = 2 cos ( θ x ) and e i θ x e i θ x = 2 i sin ( θ x ) , one can acquire for any p ,
e i π p 2 f ( β e i θ x ) + e i π p 2 f ( β e i θ x ) =   k = 0 f ( k ) ( 0 ) k ! β k ( e i π p 2 e i θ k x + e i π p 2 e i θ k x ) = 2 k = 0 f ( k ) ( 0 ) k ! β k cos ( k θ x π p 2 ) .          
Similarly,
1 i ( e i π p 2 f ( β e i θ x ) e i π p 2 f ( β e i θ x ) ) = 1 i   k = 1 f ( k ) ( 0 ) k ! β k ( e i π p 2 e i θ k x e i π p 2 e i θ k x ) = 2 k = 0 f ( k ) ( 0 ) k ! β k sin ( k θ x π p 2 ) .  
The results obtained in Equations (13) and (14) can be extended in a modified form as follows.
Lemma 2.
Let g ( z ) beanalytic function that can beexpanded in a series form,
  g ( z ) = k = 0 M k e k z ,  
whether z   be real or imaginary   ,   a n d   k = 0 M k   be absolutely convergent, then
Equations (13) and (14) can be re written in a modified form as follows,
1   2 ( e i π p 2 g ( i θ x ) + e i π p 2 g ( i θ x ) ) = 1 2 k = 0 M k ( e i π p 2 e i k θ x + e i π p 2 e i k θ x )   = k = 0 M k cos ( k θ x π p 2 )  
and,
1   2 i ( e i π p 2 g ( i θ x ) e i π p 2 g ( i θ x ) ) = 1 2 i k = 0 M k ( e i π p 2 e i k θ x e i π p 2 e i k θ x )   = k = 0 M k sin ( k θ x π p 2 ) ,      
where θ > 0   and x is any real number.
The following theorems are new results that can be used for solving some difficult improper integrals. Herein, it is worth mentioning that some of Cauchy’s results in [26], can be obtained as corollaries from the following theorems as will be stated.
Theorem 1.
Let f ( z ) be an analytic function around 0, then we have the following result,
    I = 0 x p   ( e i π p 2 f ( β e i θ x ) + e i π p 2 f ( β e i θ x ) ) 1 + x 2 d x = π   f ( β e θ ) ,  
where 1 < Re ( p ) < 1 ,   and   θ > 0 .
Proof. 
Since f is analytic around 0, thus, using the fact in Equation (13), we have,
e i π p 2 f ( β e i θ x ) + e i π p 2 f ( β e i θ x ) = 2 ( f ( 0 ) cos ( π p 2 ) + k = 1 f ( k ) ( 0 ) k !     β k cos ( k θ x π p 2 ) ) .      
Substituting Equation (19) in Equation (18) to acquire,
I = 0 x p   ( 2 ( f ( 0 ) cos ( π p 2 ) + k = 1 f ( k ) ( 0 ) k !     β k cos ( k θ x π p 2 ) ) ) 1 + x 2 d x .       = 2 f ( 0 ) cos ( π p 2 ) 0 x p 1 + x 2 d x + 2 k = 1 f ( k ) ( 0 ) k ! β k 0 x p   cos ( k θ x π p 2 ) 1 + x 2 d x .
Use the result in Equation (8) to evaluate 0 x p 1 + x 2 d x , thus we obtain,
I = π f ( 0 ) + 2 k = 1 f ( k ) ( 0 ) k ! β k 0 x p   cos ( k θ x π p 2 ) 1 + x 2 d x .  
The integral 0 x p   cos ( k θ x π p 2 ) 1 + x 2 d x , can be solved using the Laplace transform as follows.
Define,
    j ( t ) = 0 x p ( cos ( π p 2 ) cos ( θ t k x ) + sin ( π p 2 ) sin ( θ t k x ) ) 1 + x 2 d x ,  
where θ , t > 0 , 1 < Re ( p ) < 1 ,     k ,     k 0 , and our desired integral is j ( 1 ) .
Thus, by taking Laplace transform for both sides of Equation (21), we have,
J ( s ) = 0 0 x p   e s t [ cos ( π p 2 ) cos ( θ t k x ) + sin ( π p 2 ) sin ( θ t k x ) ] 1 + x 2 d x   d t .
Using Fubini’s theorem, the double integral yields a finite answer when the integral is replaced by its absolute value, i.e., converges in the Riemann sense. Thus, we can interchange the order of the integrations, to acquire
J ( s ) = 0 x p 1 + x 2 ( 0 e s t [ cos ( π p 2 ) cos ( θ t k x ) + sin ( π p 2 ) sin ( θ t k x ) ]   d t )   d x . = 0 x p 1 + x 2 { cos ( π p 2 ) cos ( θ t k x ) + sin ( π p 2 ) sin ( θ t k x ) }   d x .    
Running the Laplace transform on both sides of Equation (22), we acquire,
J ( s ) = 0 x p   s   cos ( p π 2 ) ( s 2 + k 2 x 2 θ 2 ) ( 1 + x 2 ) + x p + 1 k θ   sin ( p π 2 ) ( 1 + x 2 ) ( s 2 + k 2 x 2 θ 2 )   d x .  
Now, use the partial fraction decomposition as,
1 ( 1 + x 2 ) ( s 2 + k 2 x 2 θ 2 ) = 1 ( s 2 k 2 θ 2 ) ( 1 x 2 + 1   ( k θ ) 2   s 2 + k 2 x 2 θ 2 ) .
Substituting Equation (24) in Equation (23) to acquire,
J ( s ) = s cos ( π p 2 ) s 2 ( θ k ) 2 0 ( x p x 2 + 1   θ 2 k 2 x p θ 2 k 2 x 2 + s 2 )   d x + θ k sin ( π p 2 ) s 2 ( θ k ) 2 0 ( x p + 1 x 2 + 1   θ 2 k 2 x p + 1 θ 2 k 2 x 2 + s 2 )   d x .
Using the results in Equation (8), one can obtain,
J ( s ) = s cos ( π p 2 ) s 2 ( θ k ) 2 ( π 2 cos ( π p 2 ) ( s θ k ) p 1 π 2 cos ( π p 2 ) ) + θ k sin ( π p 2 ) s 2 ( θ k ) 2 ( π 2 sin ( π p 2 ) + π ( s θ k ) p 2 sin ( π p 2 ) ) .    
After simple computations, Equation (25) becomes,
  J ( s ) = π 2 ( s + θ k ) .  
Taking the inverse Laplace transform for both sides of Equation (24), we acquire,
j ( t ) = π 2 e θ k t .  
Thus, j ( 1 ) = π 2 e θ k , which is the desired integral.
Substituting the value of j(1) in Equation (20), we have,
I = π f ( 0 ) + π k = 1 f ( k ) ( 0 ) k ! β k   e θ k . I = π k = 0 f ( k ) ( 0 ) k ! β k   e θ k .    
Using the fact in Equation (12), we acquire,
  I = π f ( β e θ ) .  
This completes the proof of Theorem 1. □
In the following arguments, we present some results and corollaries from Theorem (1) that are useful in solving integrals and acquire new applications.
Remark 1.
Taking the simple case p = 0 ,   and β = 1 in Theorem 1, we obtain,
        0   ( f ( e i θ x ) + f ( e i θ x ) ) 1 + x 2 d x = π   f ( e θ ) ,  
which is identical to Cauchy’s theorem obtained in [26] (P.62 Theorem 8), and in the classical table of integrals [27] (3.037 Theorem 1).
Corollary 1.
Let g ( z ) be analytic function that has the series representation in Equation (15), then,
      0 x p   ( e i π p 2 g ( i θ x ) + e i π p 2 g ( i θ x ) ) 1 + x 2 d x = π   g ( θ ) ,  
where 1 < Re ( p ) < 1 , and   θ > 0 .
Corollary 2.
Let f ( z ) be analytic function around 0, then we have the following result,
    0 x i p   ( e π p 2 f ( β e i θ x ) + e π p 2 f ( β e i θ x ) ) 1 + x 2 d x = π   f ( β e θ ) ,  
where p ,   θ > 0
Proof. 
The proof comes directly form Theorem 1. Since, we have 0 < Re ( p ) < 1 , thus if we set p = a + i b   where a and b , 1 < a < 1 . Then i p = i a b , and Re ( i p ) = Im ( p ) . Thus, Equation (31) fulfills for p . □
Corollary 3.
Let g ( z ) be analytic function that has the series representation in Equation (15), then, Corollary 2 can be rewritten as,
        0 x i p   ( e π p 2 g ( i θ x ) + e π p 2 g ( i θ x ) ) 1 + x 2 d x = π   g ( θ ) ,    
where p ,   θ > 0 .
Corollary 4.
Let f ( z ) be analytic function around 0, then we have the following result,
        0 x p 1   ( e i π p 2 f ( β e i θ x ) e i π p 2 f ( β e i θ x ) ) i   ( 1 + x 2 ) d x = π f ( β e θ ) ,    
where 0 < Re ( p ) < 2 , and θ > 0 .
Proof. 
This is a direct result of Theorem 1, that can be deduced by setting p = p 1 in Equation (18). □
Corollary 5.
Let g ( z ) be analytic function that has the series represention in Equation (15) and using the fact in Equation (17) then, we have the following result:
  0 x p 1   ( e i π p 2 g ( i θ x ) e i π p 2 g ( i θ x ) ) i   ( 1 + x 2 ) d x = π g ( θ ) ,  
where 0 < Re   ( p ) < 2 , and θ > 0 .
Proof. 
The proof comes directly from Theorem 1 and Corollary 4. □
Corollary 6.
Let f ( z ) be analytic function around 0, then we have the following result:
0 x i p 1   ( e π p 2 f ( β e i θ x ) e π p 2 f ( β e i θ x ) ) i ( 1 + x 2 ) d x = π f ( β e θ ) ,  
where 2 < I m ( p ) < 0 , and   θ > 0 .
Proof. 
The proof comes directly from Corollary 4. Since, we have 0 < Re ( p ) < 2 , thus if we set p = a + i b   where a and b ,   0 < a < 2 . Then i p = i a 1 , and Re ( i p ) = Im ( p ) . Thus, Equation (31) fulfills for 2 < Im ( p ) < 0 . □
Corollary 7.
Let g ( z ) be an analytic function that has the series representation in Equation (15) and using the fact in Equation (17), then Corollary 6. can be rewritten as,
0 x i p 1   ( e π p 2 g ( i θ x ) e π p 2 g ( i θ x ) ) 1 + x 2 d x = π   g ( θ ) ,  
where 2 < I m ( p ) 0   a n d   θ 0 .
Theorem 2.
Let f ( z ) be analytic function around 0, then we have the following result,
  P V 0 x p 1   ( e i π p 2 f ( β e i θ x ) e i π p 2 f ( β e i θ x ) ) i ( x 2 1 ) d x = π 2 ( e i π p 2 f ( β   e i θ ) + e   i π p 2 f ( β e i   θ ) ) ,  
where θ > 0 ,   and   0 < R e ( p ) < 2.
Proof. 
Since f is analytic around 0, thus, using the facts in Equations (12) and (14), we acquire,
    1 i ( e i π p 2 f ( β e i θ x ) e i π p 2 f ( β e i θ x ) ) = 2 ( k = 1 f ( k ) ( 0 ) k ! β k sin ( k θ x π p 2 ) f ( 0 ) sin ( π p 2 ) )
Substituting Equation (37) in Equation (36) to acquire,
I = P V 0 x p 1   ( 2 ( k = 1 f ( k ) ( 0 ) k ! β k sin ( k θ x π p 2 ) f ( 0 ) sin ( π p 2 ) ) ) x 2 1 d x . = 2 f ( 0 ) sin ( π p 2 ) 0 x p 1 x 2 1 d x + 2 k = 1 f ( k ) ( 0 ) k ! β k   ( P V 0 x p 1   sin ( k θ x π p 2 ) x 2 1 d x ) .
Using the result in Equation (9) to evaluate P V 0 x p 1 x 2 1 d x ,thus we obtain,
  I = π f ( 0 ) cos ( π p 2 ) + 2 k = 1 f ( k ) ( 0 ) k ! β k ( P V 0 x p 1   sin ( k θ x π p 2 ) x 2 1 d x   ) .
The integral P V 0 x p 1   sin ( k θ x π p 2 ) x 2 1 d x , can be solved using the Laplace transform as follows.
Define,
j ( t ) = P V 0 x p 1 ( cos ( π p 2 ) sin ( θ t k x ) sin ( π p 2 ) cos ( θ t k x ) ) x 2 1 d x ,    
where t , θ > 0 , 0 < R e ( p ) < 2 ,   a n d   k , k 0 , thus our desired integral is j ( 1 ) .
Therefore, by taking Laplace transform for both sides of Equation (39), we acquire,
J ( s ) = P V 0 0 x p 1   e s t [ cos ( π p 2 ) sin ( θ t k x ) sin ( π p 2 ) cos ( θ t k x ) ] x 2 1 d x   d t .
Again, using Fubini’s theorem, we interchange the order of the integration to acquire,
J ( s ) = P V 0 x p 1 x 2 1 ( 0 e s t [ cos ( π p 2 ) sin ( θ t k x ) sin ( π p 2 ) cos ( θ t k x ) ]   d t )   d x . = P V 0 x p 1 x 2 1 { cos ( π p 2 ) sin ( θ t k x ) sin ( π p 2 ) cos ( θ t k x ) }   d x .  
Running the laplace transform on Equation (40), to acquire,
J ( s ) = P V ( 0 x p   θ k   cos ( p π 2 ) ( s 2 + k 2 x 2 θ 2 ) ( x 2 1 ) d x 0 s   x p 1 sin ( p π 2 ) ( x 2 1 ) ( s 2 + k 2 x 2 θ 2 ) d x ) .
Now, using the partial fraction decomposition,
1 ( x 2 1 ) ( s 2 + k 2 x 2 θ 2 ) = 1 ( s 2 + k 2 θ 2 ) ( x 2 1 ) k 2 θ 2 ( s 2 + k 2 θ 2 ) ( s 2 + k 2 x 2 θ 2 ) .
Substituting Equation (42) in Equation (41), we have,
J ( s ) = θ k cos ( π p 2 ) ( s 2 + k 2 θ 2 ) P V 0 ( x p ( x 2 1 ) k 2 θ 2 x p ( s 2 + k 2 x 2 θ 2 ) )   d x s sin ( π p 2 ) ( s 2 + k 2 θ 2 ) P V 0 ( x p 1 ( x 2 1 ) x p 1 ( s 2 + k 2 x 2 θ 2 ) )   d x .
Using the facts in Equations (8) and (9), one can obtain:
J ( s ) = θ k cos ( π p 2 ) ( s 2 + k 2 θ 2 ) ( π 2 tan ( π p 2 ) π ( s θ k ) p 1 2 cos ( π p 2 ) ) s sin ( π p 2 ) ( s 2 + k 2 θ 2 ) ( π 2 cot ( π p 2 ) π ( s θ k ) p 2 2 sin ( π p 2 ) ) .  
After simple computations, Equation (43) becomes,
J ( s ) = π ( s cos ( p π 2 ) + k θ sin ( p π 2 ) ) 2 ( s 2 + k 2 θ 2 ) .  
Taking the inverse Laplace transform for both sides of Equation (44), we acquire,
j ( t ) = π 2   cos ( k t θ p π 2 ) .    
Thus, j ( 1 ) = π 2 cos ( k θ π p 2 ) , which is the desired integral.
Substituting the value of j(1) in Equation (38), we have,
        I = π f ( 0 ) cos ( π p 2 ) + π k = 1 f ( k ) ( 0 ) k ! β k cos ( k θ π p 2 )   I = π k = 0 f ( k ) ( 0 ) k ! β k cos ( k θ π p 2 )      
By rewriting cos ( k θ π p 2 ) in the exponential form, Equation (45) becomes,
  I = π 2 k = 0 f ( k ) ( 0 ) k ! β k ( e i π p 2 e i θ k + 1 2 e i π p 2 e i θ k )
and using the fact in Equation (12), we acquire,
I = π 2 ( e i π p 2 f ( β   e i θ ) + e   i π p 2 f ( β e i   θ ) )
This completes the proof of Theorem 2. □
Now, we introduce some corollaries that comes directly after simple assumptions and computations in the same manner that used in the previous results.
Corollary 8.
Let g ( z ) be analytic function that has the series representation in Equation (15) and using the fact in Equation (17) then, Theorem 2 can be rewritten as,
P V 0 x p 1   ( e i π p 2 g ( i θ x ) e i π p 2 g ( i θ x ) ) i ( x 2 1 ) d x = π 2 ( e i π p 2 g (   i θ ) + e   i π p 2 g ( i θ ) ) .    
Corollary 9.
Let f ( z ) be an analytic function around 0, then we have the following result,
P V 0 x p   ( e i π p 2 f ( β e i θ x ) + e i π p 2 f ( β e i θ x ) ) ( x 2 1 ) d x = π 2 i ( e i π p 2 f ( β   e i θ ) e   i π p 2 f ( β e i   θ ) ) ,    
where θ > 0 , 1 < R e ( p ) < 1.
Remark 2.
Taking the simple case p = 0  and β = 1 in Corollary 9, we have,
P V 0   f ( e i θ x ) + f ( e i θ x ) ( x 2 1 ) d x = π 2 i ( f ( e i θ ) f ( e i θ ) ) ,  
where θ > 0 .
Which is identical to Cauchy’s theorem obtained in [26] (P.62 Theorem 9).
Corollary 10.
Let g ( z ) be an analytic function that has the series representation in Equation (15), and using the fact in Equation (16), then Corollary 9 can be rewritten as,
P V 0 x p   ( e i π p 2 g ( i θ x ) + e i π p 2 g ( i θ x ) ) ( x 2 1 ) d x = π 2 i ( e i π p 2 g ( i θ ) e   i π p 2 g ( i θ ) ) .  
Corollary 11.
Let f ( z ) be an analytic function around 0, then we have the following result,
P V 0 x i p   ( e π p 2 f ( β e i θ x ) + e π p 2 f ( β e i θ x ) ) ( x 2 1 ) d x = π 2 i ( e π p 2 f ( β   e i θ ) e   π p 2 f ( β e i   θ ) ) ,  
where θ > 0   and p .
Corollary 12.
Let g ( z ) be an analytic function that has the series representation in Equation (15) and using the fact in Equation (16), then Corollary 11 can be rewritten as,
  P V 0 x p   ( e π p 2 g ( i θ x ) + e π p 2 g ( i θ x ) ) ( x 2 1 ) d x = π 2 i ( e π p 2 g ( i θ ) e   π p 2 g ( i θ ) ) .  
Corollary 13.
Let f ( z ) be an analytic function around 0, then we have the following result,
P V 0 x i p 1   ( e π p 2 f ( β e i θ x ) e π p 2 f ( β e i θ x ) ) x 2 1 d x = π 2 ( e π p 2 f ( β   e i θ ) + e   π p 2 f ( β e i   θ ) ) ,
where 2 < I m ( p ) < 0 , and   θ > 0 .
Corollary 14.
Let g ( z ) be an analytic function that has the series representation in Equation (15) and using the fact in Equation (17), then Corollary 13 can be written as,
P V 0 x i p 1   ( e π p 2 g ( i θ x ) e π p 2 g ( i θ x ) ) x 2 1 d x = π 2 ( e π p 2 g (   i θ ) + e   π p 2 g ( i   θ   ) ) .  
Theorem 3.
Let f ( z ) be an analytic function around 0, then we have the following result,
  I =   f ( β e i θ x ) + f ( β e i θ x ) x 2 2 x cos ( α ) + 1 d x = π sin ( α   ) ( f ( β e θ sin ( α ) i θ cos ( α ) ) + f ( β e θ sin ( α ) + i cos ( α ) θ ) ) ,  
where 0 < | α | < π   and   θ > 0 .
Proof. 
Since f is analytic around 0, thus Substituting Equation (13) in Equation (54), we acquire,
I = 2 k = 0 f ( k ) ( 0 ) k !     β k cos ( k θ x ) x 2 2 x cos ( α ) + 1 d x .   =   2 k = 0 f ( k ) ( 0 ) k ! β k   cos ( k θ x ) x 2 2 x cos ( α ) + 1 d x .
The integral   cos ( k θ x ) x 2 2 x cos ( α ) + 1 d x can be solved as follows, let,
  J =   cos ( k θ x ) x 2 2 x cos ( α ) + 1 d x ,  
and x = cos ( α ) + y sin ( α )   d x = s i n ( α )   d y .
( x cos ( α ) ) 2 = x 2 2 x cos ( α ) + cos 2 ( α ) = y 2 sin 2 ( α ) .
Thus,
x 2 2 c o s ( α ) x + 1 = sin 2 ( α ) ( 1 + y 2 )
Therefore, Equation (56) becomes,
J = cos ( k θ ( cos ( α ) + y sin ( α ) ) ) sin ( α ) ( 1 + y 2 ) d y = 1 sin ( α ) cos ( k θ cos ( α ) ) cos ( k θ   y s i n ( α ) ) sin ( k θ cos ( α ) ) sin ( k θ   y s i n ( α ) ) ( 1 + y 2 ) d y = 1 sin ( α   ) ( cos ( k θ cos ( α ) ) cos ( k θ   y s i n ( α ) ) ( 1 + y 2 ) d y sin ( k θ cos ( α ) ) sin ( k θ   y s i n ( α ) ) ( 1 + y 2 ) d y ) .    
Since the function sin ( k θ   y s i n ( α ) ) ( 1 + y 2 ) is odd, then,
sin ( k θ   y s i n ( α ) ) ( 1 + y 2 ) d y = 0 .  
Using the fact that can be obtained in [2],
cos ( θ x ) 1 + x 2 d x = π e θ ,
and substituting Equations (58) and (59) in Equation (57), to acquire,
  J = π e θ k   sin ( α ) sin ( α   ) cos ( k θ cos ( α ) ) .  
Substituting Equation (60) in Equation (55), we acquire,
I = 2 k = 0 f ( k ) ( 0 ) k ! β k π e θ k   sin ( α ) sin ( α   ) cos ( k θ cos ( α ) ) .
Rewriting cos ( k θ cos ( α ) ) in the exponential form and after simple computations, we have,
I = π sin ( α   ) k = 0 f ( k ) ( 0 ) k ! β k ( e θ k   sin ( α ) i θ cos ( α ) k + e θ k sin ( α ) + i cos ( α ) θ k ) .
Using the fact in Equation (12), we acquire,
I = π sin ( α   ) ( f ( β e θ sin ( α ) i θ cos ( α ) ) + f ( β e θ sin ( α ) + i cos ( α ) θ ) ) .
This completes the proof of Theorem 3. □
Remark 3.
Taking the simple case β = 1 , α = π 2   in Theorem 3, we obtain Equation (28)which is similar to Cauchy’s theorem obtained in [26] (P.62 Theorem 8) and in the classical table of integrals [27] (3.037 Theorem 1).
Corollary 15.
Let g ( z ) be an analytic function that has the series representation in Equation (15), then we have,
0   ( g ( i θ x ) + g ( i θ x ) ) x 2 2 x cos ( α ) + 1 d x =     = π sin ( α   ) ( g ( θ sin ( α ) + i θ cos ( α ) ) + g ( θ sin ( α ) i cos ( α ) θ ) ) ,  
where 0 < | α | < π   and   θ > 0 .
Corollary 16.
Let h ( z ) be analytic function around 0, then we have the following result,
0 x p   ( e i π p 2 h ( β e i θ x ) + e i π p 2 h ( β e i θ x ) ) x 2 2 x cos ( α ) + 1 d x = π   g ( θ ) ,  
where 1 < Re ( p ) < 1 , 0 < | α | < π ,   and   θ > 0 , where h(z) is analytic function around zero.
Proof. 
The proof comes directly form Theorem 3 by letting,
f ( β e   θ x ) = ( x ) p   h ( β e θ x ) ,
Thus, we have,
f ( β e   i θ x ) + f ( β e i γ x ) = ( i x ) p h ( β e i θ x ) + ( i x ) p h ( β e i   θ x ) = x p ( e i π p 2 h ( β e i θ x ) + e i π p 2   h ( β e i θ x ) ) .
Now, substituting Equation (62) in Theorem 3, we acquire,
0 x p   ( e i π p 2 h ( β e i θ x ) + e i π p 2 h ( β e i θ x ) ) x 2 2 x cos ( α ) + 1 d x = π sin ( α   ) ( ( sin ( α ) + i cos ( α ) ) p h ( β e θ ( sin ( α ) i cos ( α ) ) ) + ( sin ( α ) i cos ( α ) ) p h ( β e θ ( sin ( α ) + i cos ( α ) ) ) ) .
 □

4. Applications and Examples

In this section, we introduce some examples on some complicated integrals, that can’t be solved easily by familiar methods or may take huge efforts to be solved. Herein, we show that using the obtained results in this article, the solution can be determined directly in a simple way.
Example 1.
Solve the following integral,
P V 0   ( e b   tan 1 ( θ x ) e b   tan 1 ( θ x ) ) cos ( b 2 ln ( 1 + θ 2 x 2 ) ) x 2 1 d x ,
Where   θ > 0 ,   a n d   b .
Solution. 
Let g ( z ) = cos ( b ln ( 1 + z ) ) , then we have,
g ( i θ x ) + g ( i θ x ) = cos ( b ln ( 1 i θ x ) ) + cos ( b ln ( 1 + i θ x ) ) = cos ( b 2 ln ( 1 + θ 2 x 2 ) + i b   tan 1 ( θ x ) ) + cos ( b 2 ln ( 1 + θ 2 x 2 ) i b   tan 1 ( θ x ) ) = 2 cosh ( b   tan 1 ( θ x ) ) cos ( b 2 ln ( 1 + θ 2 x 2 ) )
Using Corollary (10) and letting p = 0 , we acquire,
P V 0   ( g ( i θ x ) + g ( i θ x ) ) ( x 2 1 ) d x = π 2 i ( g ( i θ ) g ( i θ ) ) .       P V 0   ( 2 cosh ( b   tan 1 ( θ x ) ) cos ( b 2 ln ( 1 + θ 2 x 2 ) ) ) 1 x 2 d x = π ( sin ( 1 2 b ln ( θ 2 + 1 ) ) sinh ( b t a n 1 ( θ x ) ) ) .
Example 2.
Solve the following integral,
0   x p e cos ( θ x ) cos ( π p 2 s i n θ x ) x 2 + 1 d x ,
where 1 < p < 1.
Solution. 
Let f ( z ) = e z then, we have,
e i π p 2 f ( e i θ x ) + e i π p 2 f ( e i θ x ) = 2 e cos ( θ x ) cos ( π p 2 s i n θ x ) .
Using Theorem 1 and setting β = 1 , we have,
0   x p e cos ( θ x ) cos ( π p 2 s i n θ x ) x 2 + 1 d x = 1 2 π   f ( β e θ ) = π 2   e e θ
Example 3.
Solve the following integral,
0 x p ( α cos ( π p 2 ) cos ( π p 2 + x ) )   ( 1 + x 2 ) ( 1 2 α   cos ( x ) + α 2 )   d x ,
where 1 < R e ( p ) < 1 and α < 1.
Solution. 
Using Theorem 1 and letting β = 1 ,   θ = 1 and f ( z ) = 1 α z , we have,
e i π p 2 f ( e i θ x ) + e i π p 2 f ( e i θ x ) = 2 ( α cos ( p π 2 ) cos ( p π 2 + x θ ) ) ( 1 2 α   cos ( x ) + α 2 ) . 0 x p ( α cos ( π p 2 ) cos ( π p 2 + x ) )   ( 1 + x 2 ) ( 1 2 α   cos ( x ) + α 2 ) d x = π 2   ( 1 α e 1 ) .

5. Conclusions

In this research, we introduced three interesting theorems in complex analysis. The main purpose of these results is to acquire the solutions of some families of improper integrals that take much effort and time to solve, and sometimes they can’t be solved directly even with computer software. The proof of these theorems was obtained in the article using idea of Laplace transform and some techniques depends on some theorems such as Fubini’s theorem. Some corollaries and remarks were discussed to prove the applicability of these theorems and comparisons with previous results are also illustrated. Three applications are presented, to show the reliability and the simplicity of the introduced theorems. These results can be generalized to solve wide families of improper integrals on different domains. Finally, we mention that the restrictions on the domain and the parameters of these results can mitigated, so that we can acquire better results in the future. In the following, we introduce some mathematical remarks about the proposed theorems.
The theorems can be extended for analytic functions around α ,   where α , and we are working on such results.
If θ < 0 , then Theorems 1, Theorems 3 and the corollaries after, can be generalized by switching θ to | θ | .
If θ < 0 , then Theorems 2 and the corollaries after, can be generalized by switching θ to θ .

Author Contributions

All authors contribute equally to this research. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arfken, G.B.; Weber, H.J. Mathematical Methods for Physicists, 5th ed.; Academic Press: Cambridge, MA, USA, 2000. [Google Scholar]
  2. Saadeh, R.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.; Din, U. Application of fractional residual power series algorithm to solve newell-whitehead-segel equation of fractional order. Symmetry 2019, 11, 1431. [Google Scholar] [CrossRef] [Green Version]
  3. Saadeh, R.; Al-Smadi, M.; Gumah, G.; Khalil, H.; Khan, R.A. Numerical investigation for solving two-point fuzzy boundary value problems by reproducing kernel approach. Appl. Math. Inf. Sci. 2016, 10, 1–13. [Google Scholar] [CrossRef]
  4. Gharib, G.; Saadeh, R. Reduction of the Self-dual Yang-Mills Equations to Sinh-Poisson Equation and Exact Solutions. WSEAS Interact. Math. 2021, 20, 540–546. [Google Scholar] [CrossRef]
  5. Burqan, A.; El-Ajou, A.; Saadeh, R.; Al-Smadi, M. A new efficient technique using Laplace transforms and smooth expansions to construct a series solutions to the time-fractional Navier-Stokes equations. Alex. Eng. J. 2021, 61, 1069–1077. [Google Scholar] [CrossRef]
  6. Saadeh, R.; Ghazal, B. A New Approach on Transforms: Formable Integral Transform and Its Applications. Axioms 2021, 10, 332. [Google Scholar] [CrossRef]
  7. Qazza, A.; Hatamleh, R.; Al-hawari, M. Dirichlet Problem in the Simply Connected Domain, Bounded by Unicursal Curve. Int. J. Appl. Math. 2009, 22, 599–614. [Google Scholar]
  8. Qazza, A.; Hatamleh, R. Dirichlet Problem in the simply connected domain, bounded by the nontrivial kind. Adv. Differ. Equ. Control. Processes 2016, 17, 177–188. [Google Scholar] [CrossRef]
  9. Qazza, A.; Hatamleh, R.; Alodat, N. About the Solution Stability of Volterra Integral Equation with Random Kernel. Far East J. Math. Sci. 2016, 100, 671–680. [Google Scholar] [CrossRef]
  10. Qazza, A.; Hatamleh, R. The existence of a solution for semi-linear abstract differential equations with infinite B chains of the characteristic sheaf. Int. J. Appl. Math. 2018, 31, 611. [Google Scholar] [CrossRef]
  11. Saadeh, R. Numerical solutions of fractional convection-diffusion equation using finite-difference and finite-volume schemes. J. Math. Comput. Sci. 2021, 11, 7872–7891. [Google Scholar]
  12. Qazaa, A.; Burqan, A.; Saadeh, R. A New Attractive Method in Solving Families of Fractional Differential Equations by a New Transform. Mathematics 2021, 9, 3039. [Google Scholar] [CrossRef]
  13. Edwan, R.; Saadeh, R.; Hadid, S.; Al-Smadi, M.; Momani, S. Solving Time-Space-Fractional Cauchy Problem with Constant Coefficients by Finite-Difference Method. In Computational Mathematics and Applications; Springer: Berlin/Heidelberg, Germany, 2020; pp. 25–46. [Google Scholar]
  14. Burqan, A.; Saadeh, R.; Qazaa, A. A Novel Numerical Approach in Solving Fractional Neutral Pantograph Equations via the ARA Integral Transform. Symmetry 2022, 14, 50. [Google Scholar] [CrossRef]
  15. Saadeh, R. Numerical algorithm to solve a coupled system of fractional order using a novel reproducing kernel method. Alex. Eng. J. 2021, 60, 4583–4591. [Google Scholar] [CrossRef]
  16. Nahin, P.J. Inside Interesting Integrals: A Collection of Sneaky Tricks, Sly Substitutions, and Numerous Other Stupendously Clever, Awesomely Wicked, and Devilishly Seductive Maneuvers for Computing Nearly 200 Perplexing Definite Integrals from Physics, Engineering, and Mathematics (Plus 60 Challenge Problems with Complete, Detailed Solutions); Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  17. Boros, G.; Moll, V. Irresistible Integrals: Symbolics, Analysis and Experiments in the Evaluation of Integrals; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar] [CrossRef]
  18. Cauchy, A.L. Sur diverses relations qui existent entre les résidus des fonctions et les integrals définies. In Exercices de mathématiques; Reprinted in Oeuvres Complètes d’Augustin Cauchy; Gauthier-Villars: Paris, France, 1903; Volume 6, pp. 124–145. [Google Scholar]
  19. Boas, H.P. Cauchy’s Residue Sore Thumb. Am. Math. Mon. 2018, 125, 16–28. [Google Scholar] [CrossRef] [Green Version]
  20. Amdeberhan, T.; Espinosa, O.; Gonzalez, I.; Harrison, M.; Moll, V.; Straub, A. Ramanujan’s Master Theorem. Ramanujan J. 2012, 29, 103–120. [Google Scholar] [CrossRef]
  21. Stein, E.M.; Shakarchi, R. Complex Analysis; Princeton University Press: Princeton, NJ, USA, 2003. [Google Scholar]
  22. Thomas, G.B.; Finney, J.; Ross, L. Calculus and Analytic Geometry, 9th ed.; Addison Wesley: Boston, MA, USA, 1996; ISBN 0-201-53174-7. [Google Scholar]
  23. Henrici, P. Applied and Computational Complex Analysis; Power Series, Integration, Conformal Mapping, Location of Zeros; John Wiley & Sons: Hoboken, NJ, USA, 1988; Volume 1. [Google Scholar]
  24. Andrews, L. Special Functions of Mathematics for Engineers; SPIE Optical Engineering Press: Bellingham, WA, USA, 1998. [Google Scholar]
  25. Brown, J.W.; Ruel Churchill, V. Complex Variables and Applications; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
  26. M’Emoire sur les Integrales Definies, Prises entre des Limites Imaginaires; reprint of the 1825 original in Oeuvres compl‘etes d’Augustin Cauchy; Gauthier-Villars: Paris, France, 1974; Volume 15, pp. 41–89.
  27. Zwillinger, D. Table of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 2014; pp. 252–253. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abu Ghuwaleh, M.; Saadeh, R.; Burqan, A. New Theorems in Solving Families of Improper Integrals. Axioms 2022, 11, 301. https://doi.org/10.3390/axioms11070301

AMA Style

Abu Ghuwaleh M, Saadeh R, Burqan A. New Theorems in Solving Families of Improper Integrals. Axioms. 2022; 11(7):301. https://doi.org/10.3390/axioms11070301

Chicago/Turabian Style

Abu Ghuwaleh, Mohammad, Rania Saadeh, and Aliaa Burqan. 2022. "New Theorems in Solving Families of Improper Integrals" Axioms 11, no. 7: 301. https://doi.org/10.3390/axioms11070301

APA Style

Abu Ghuwaleh, M., Saadeh, R., & Burqan, A. (2022). New Theorems in Solving Families of Improper Integrals. Axioms, 11(7), 301. https://doi.org/10.3390/axioms11070301

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