Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Novel Hybrid MSA-CSA Algorithm for Cloud Computing Task Scheduling Problems
Next Article in Special Issue
Kneser-Type Oscillation Criteria for Half-Linear Delay Differential Equations of Third Order
Previous Article in Journal
The λ-Fold Spectrum Problem for the Orientations of the Eight-Cycle
Previous Article in Special Issue
Solutions of Magnetohydrodynamics Equation through Symmetries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Geometric Interpretations of Euler’s Substitutions

by
Jan L. Cieśliński
* and
Maciej Jurgielewicz
Faculty of Physics, University of Bialystok, 15-245 Białystok, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(10), 1932; https://doi.org/10.3390/sym15101932
Submission received: 27 September 2023 / Revised: 12 October 2023 / Accepted: 16 October 2023 / Published: 18 October 2023

Abstract

:
We consider a classical case of integrals containing an irrational integrand in the form of a square root of a quadratic polynomial. It is known that such “irrational integrals” can be expressed in terms of elementary functions by one of three of Euler’s substitutions. It is less well known that the Euler substitutions have a geometric interpretation. In the framework of this interpretation, one can see that the number 3 is not the most suitable. We show that it is natural to introduce a fourth Euler substitution. In his original treatise, Leonhard Euler used two substitutions which are sufficient to cover all cases.

1. Introduction

Integrals of rational functions can be expressed in terms of elementary functions. Therefore, a natural method of integration consists of using suitable substitutions and integration by parts to reduce our problem to integration of rational functions.
In this paper, we consider irrational integrals containing the quadratic root of a quadratic polynomial, i.e., integrals of the form
R ( x , y ) d x ,
where R is a rational fuction (a quotient of two polynomials) of x and y, and
y = a x 2 + b x + c .
The subject is, in principle, known. A standard method to deal with such integrals consists of using one of the so-called Euler’s substitutions [1,2,3]. However, there are some details which need to be clarified. We will describe in detail a geometric approach to this problem and explain how many Euler substitutions actually exist.
In fact, to the best of our knowledge, all sources and textbooks mention exactly three types of substitutions in this context. It is not clear who was the first to introduce such classification. Leonhard Euler himself used only two of these substitutions (which is sufficient to cover all cases). Three Euler substitutions are usually introduced and discussed in Russian sources; see, e.g., [4,5,6] (Leonhard Euler, although of Swiss origin, lived and worked in Saint Petersburg for many years). Surprisingly enough, the three substitutions appeared in an old textbook, published in 1892 by a Harvard professor, William E. Byerly [7], without any reference to Euler.
In our paper, we present a clear geometric intepretation of this problem, shortly mentioned in some sources, mainly of Russian origin [2,8]. The textbook [8] is not translated into English. Another book by the same author [4] does not mention this geometric approach in the section on Euler’s substitutions.
The main novelty of this paper is the introduction of the fourth Euler substitution, which is a natural consequence of the geometric approach discussed in our paper.

2. Three Classical Euler’s Substitutions

The main idea of Euler’s substitutions consists of expressing a x 2 + b x + c as a linear function of x and a new parameter t in such a way that the resulting equation is linear with respect to x. In this paper, we use the most common numbering of these three substitutions, compare [1,2,4,5,7]. In some sources, a different order is used; see [6,9,10].

2.1. First Euler Substitution

This substitution can be done only in the case a > 0 :
a x 2 + b x + c = ± x a + t .
Squaring both sides we get:
a x 2 + b x + c = a x 2 ± 2 x t a + t 2 .
Terms quadratic in x cancel out and the resulting equation is linear in x. Computing x, we get a rational dependence on t:
x = t 2 c b 2 t a .
Then, from (2) and (3), we get
y = t 2 a + t b c a b 2 t a .

2.2. Second Euler Substitution

This substitution can be done only in the case c > 0 :
a x 2 + b x + c = x t ± c .
Squaring both sides we get:
a x 2 + b x + c = x 2 t 2 ± 2 x t c + c .
The constant c cancels out and dividing both sides by x we again derive an equation linear in x. Hence, similarly as in the previous case,
x = b 2 t c t 2 a , y = b t ( t 2 + a ) c t 2 a .

2.3. Third Euler Substitution

This substitution can be done only in the case Δ > 0 , where
Δ b 2 4 a c
is the discriminant of the quadratic polynomial. Then the polynomial has two distinct real roots x 1 and x 2 , and the third Euler substitution is given by:
a x 2 + b x + c = ( x x 1 ) t .
Squaring both sides we get:
a ( x x 1 ) ( x x 2 ) = ( x x 1 ) 2 t 2 a ( x x 2 ) = ( x x 1 ) t 2 .
Computing x from the resulting equation and then using (10) and (2) we obtain
x = t 2 x 1 a x 2 t 2 a , y = ( x 1 x 2 ) a t t 2 a ,
where, of course,
x 1 , 2 = b ± b 2 4 a c 2 a .

2.4. Original Euler’s Approach

It is interesting that Leonhard Euler himself, in his famous monograph, used only two of these substitutions, see [11]. He considered two cases: Δ > 0 and Δ < 0 . In the first case ( Δ > 0 ) he proposed the substitution (6), while in the second case ( Δ < 0 ) he proposed the substitution (3) in a slightly modified form:
a x 2 + b x + c = x a t c .
Obviously, the case Δ = 0 is not included because then the quadratic polynomial is a square of the linear function in x and y is linear is x as well. Hence, the integrand in (1) is rational in x from the very beginning.

3. Geometric Interpretation

It is convenient to square both sides of (2) resulting in the equation of a quadratic curve
y 2 = a x 2 + b x + c .
We will denote this curve (a conic section) by Q a b c , i.e., ( x , y ) Q a b c .

3.1. Elliptic Case: a < 0

The canonical form of the quadratic polynomial yields:
y 2 + | a | x b 2 | a | 2 = c b 2 4 a .
We can distinguish three cases, depending on the sign of the discriminant Δ :
Δ < 0 Q a b c = ( empty set ) ,
Δ = 0 Q a b c = b 2 | a | , 0 ( sin gle point ) ,
Δ > 0 Q a b c is an ellipse .
Only in the last case, we get a non-degenerated quadratic curve.

3.2. Parabolic Case: a = 0

For a = 0 (and b 0 ) the conic Q a b c is a parabola with the symmetry axis y = 0 .

3.3. Hyperbolic Case: a > 0

The canonical form of the quadratic polynomial yields:
y 2 a x + b 2 a 2 = c b 2 4 a .
We can distinguish three cases, depending on the sign of the discriminant Δ :
Δ < 0 Q a b c is a hyperbola with vertices at the line x = b 2 a ,
Δ = 0 Q a b c is a pair of intersection lines ,
Δ > 0 Q a b c is a hyperbola with vertices at x axis .
Therefore, for Δ 0 we get a non-degenerated quadratic curve.

3.4. Rational Parameterization: Standard Approach

The key idea leading to a rational parameterization consists of fixing an arbitrary point P 0 = ( x 0 , y 0 ) on the conic Q a b c and assigning to any other point P = ( x , y ) of this conic the line P 0 P . Taking as a parameter t the slope of this line, we obtain a rational parameterization of the conic Q a b c [2,8]. Thus, we have the system of three equations:
y 2 = a x 2 + b x + c , y 0 2 = a x 0 2 + b x 0 + c , y y 0 = t ( x x 0 ) .
The points ( x , y ) and ( x 0 , y 0 ) belong to the conic Q a b c and t is the slope of the straight line passing through ( x , y ) and ( x 0 , y 0 ) . Subtracting the second equation from the first one we get:
( y y 0 ) ( y + y 0 ) = ( x x 0 ) ( a ( x + x 0 ) + b ) , y 0 2 = a x 0 2 + b x 0 + c , y y 0 = t ( x x 0 ) .
Substituting the last equation into the first one we obtain:
t ( y + y 0 ) a ( x + x 0 ) b ( x x 0 ) = 0 , y 0 2 = a x 0 2 + b x 0 + c , y y 0 = t ( x x 0 ) .
Assuming x x 0 , we get
t ( y + y 0 ) = a ( x + x 0 ) + b , y 0 2 = a x 0 2 + b x 0 + c , y y 0 = t ( x x 0 ) .
Now, the first and the last equation form a system of two linear equations for two variables x , y , which can be solved in the standard way. As a result, we obtain:
x = x 0 t 2 2 y 0 t + a x 0 + b t 2 a , y = y 0 t 2 + ( 2 a x 0 + b ) t a y 0 t 2 a ,
which means that we expressed x and y as rational functions of the parameter t.
Corollary 1.
There are many Euler-like substitutions. Each of them is determined by the choice of x 0 , provided that a x 0 2 + b x 0 + c 0 . Then the point P 0 ( x 0 , y 0 ) is given by:
P 0 = x 0 , ± a x 0 2 + b x 0 + c ,
and other points P = ( x , y ) Q a b c are parameterized by (28).
In particular, the second Euler substitution corresponds to x 0 = 0 (provided that the graph of the quadric Q a b c intersects the axis y), see Figure 1 and Figure 2. The third Euler substitution corresponds to x 0 being a root of the polynomial a x 2 + b x + c (provided that the graph of Q a b c intersects the axis x), see Figure 3 and Figure 4.
The first Euler substitution apparently does not fit this picture. However, its geometric interpretation is even simpler and more evident. The Formula (3) describes the family of lines parallel to asymptotes of the corresponding hyperbola, see Figure 5. We may treat it as a special case of (28) when the point ( x 0 , y 0 ) lies at a very large number. Note that points ( x 0 , ± x 0 a ) belong to the conic (15) in the limit for x 0 .

4. New Insights from the Geometric Interpretation

The description given in the previous section is more or less known (see, e.g., [2,8]), although we are not aware of any reference containing all these details. We are going to derive from this geometric picture more quite interesting consequences.
First of all, we identify characteristic points on the graph of a quadratic curve which can be chosen as P 0 in the most natural way: vertices ( M 1 , M 2 , R 1 , R 2 ) and intersections with coordinate axes ( R 1 , R 2 , V 1 , V 2 ); see Figure 6 and Figure 7.
In particular, in the case of the second Euler substitution, P 0 = V 2 (see Figure 1 and Figure 2) or P 0 = V 1 , while in the case of the third Euler substitution P 0 = R 1 (see Figure 3) or P 0 = R 2 (see Figure 4). The first Euler substitution is related to P 0 .

4.1. Fourth Euler’s Substitution

The geometric approach presented above includes all three classical Euler’s substitutions, but it is still missing vertices M 1 and M 2 . Therefore, it is natural to introduce another (fourth) Euler’s substitution, geometrically related to missing vertices: P 0 = M 1 (see Figure 8 and Figure 9) or P 0 = M 2 .
The algebraic description of the fourth Euler substitution is based on the canonical form of the quadratic polynomial:
y = a ( x p ) 2 + q
where
p = b 2 a , q = Δ 4 a .
The fourth Euler substitution is defined by:
y = q + ( x p ) t .
Squaring both sides we get:
a ( x p ) 2 + q = q + 2 ( x p ) t q + ( x p ) 2 t 2 .
The constant q cancels out and dividing both sides by x p , we obtain
a ( x p ) = 2 t q + ( x p ) t 2 ,
which is linear in x. Hence
x p = 2 t q a t 2 ,
and using (32) we get
y = a + t 2 a t 2 q .
Thus we have a rational dependence of x and y on the parameter t. Moreover,
d x d t = 2 ( a + t 2 ) q ( a t 2 ) 2 ,
and we can easily transform the irrational integral function (1) into an integral function rational with respect to t.

4.2. Simplifying Euler’s First Substitution

A geometric approach suggests some modifications or new variants of the existing rational parameterizations. Introducing a new parameter τ
τ = b 2 t a t = ( τ b ) 2 a ,
and substituting it into (4) and (5), we obtain the following simplification of the first Euler substitution:
x = 1 4 a τ + Δ τ 2 b , y = 1 4 a τ Δ τ .
Geometrically, the parameter t is related to intersections with the y axis (compare Figure 5), while the parameter τ is related to intersections with the vertical symmetry axis (i.e., the line x = p ). Indeed, the parameter τ = 0 corresponds to the line passing through the point ( p , 0 ) and this is one of two asymptotes (that is why x and y for τ 0 ).

4.3. Euler’s First Substitution as a Limit of the Generic Case

We are going to show that the first Euler substitution can be derived from the generic case (28) by taking a suitable limit x 0 and y 0 . We consider the pencil of lines y y 0 = t ( x x 0 ) (compare (24)) but as a parameter we take the ordinate τ of the intersection of the line P 0 P with the y-axis (i.e., y = τ for x = 0 ). Hence
t = y 0 τ x 0 .
This change of variable, replacing t with τ , works for any x 0 0 . Therefore, the second Euler substitution (related to the case x 0 = 0 ) is excluded.
Substituting (40) into (28) we obtain
x = ( τ 2 c ) x 0 b x 0 + c + τ 2 2 τ y 0 , y = τ ( b x 0 + 2 c ) ( c + τ 2 ) y 0 b x 0 + c + τ 2 2 τ y 0 .
Assuming a > 0 and taking into account y 0 2 = a x 0 2 + b x 0 + c we take the limit | x 0 | . Hence y 0 ± x 0 a and
x | x 0 | τ 2 c b 2 τ a , y | x 0 | τ b ( c + τ 2 ) a b 2 τ a .
Comparing (42) with (4) and (5) we easily see that both solutions are identical, provided that we identify τ with t. Note, of course, that the t parameter given by (40) is different from the t parameter used in Section 2.1.

5. Euler’s Substitutions versus Trigonometric Substitutions

Another popular method for computing irrational integrals (1) consists of making a suitable trigonometric or hyperbolic substitution. We use the canonical form of the quadratic curve (compare (30)):
y 2 = a ( x p ) 2 + q .
Assuming q 0 (otherwise y depends linearly on x) we introduce new variables ξ , η as folows:
η = y | q | , ξ = ( x p ) | a | | q | .
Then (43) becomes
η 2 = ( sgn a ) ξ 2 + sgn q ,
because a / | a | = sgn a , etc.
Thus we have three separate cases (in the fourth case –both signs negative– there are no real solutions), where trigonometric or hyperbolic substitutions are well known:
η = ξ 2 1 ξ = cosh ϑ , η = sinh ϑ , η = 1 ξ 2 ξ = cos ϑ , η = sin ϑ , η = ξ 2 + 1 ξ = sinh ϑ , η = cosh ϑ .
Is it better than Euler’s substitutions? This is a matter of taste. Perhaps it is easier to memorize, however, one has to remember that integrals of trigonometric or hyperbolic functions have to be converted into integrals of rational functions by another substitution:
t = tan θ 2 or t = tanh θ 2 .

6. Conclusions

We presented and discussed a geometric approach to Euler substitutions. One consequence of this thorough discussion was the introduction of a fourth Euler substitution, in addition to three traditionally mentioned Euler substitutions. In fact, we can say that more (one parameter family) Euler-like substitutions exist and can be further modified or simplified by suitable linear or fractional linear transformations.
Surprisingly, the subject of constructing rational parametrization of algebraic curves (rationalizing roots) has recently become important in the context of Feynman integrals and computations in high energy particle physics [12,13]. Furthermore, Euler’s substitutions were applied for reducing square roots in some mathematical finance calculations [14]. It would be interesting to appply, in those fields, some geometric ideas presented in this paper.

Author Contributions

Conceptualization, J.L.C.; methodology, J.L.C.; formal analysis, J.L.C. and M.J.; investigation, J.L.C. and M.J.; visualization, M.J.; writing—original draft preparation, J.L.C.; writing—review and editing, J.L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Euler Substitution. Wikipedia. Available online: https://en.wikipedia.org/wiki/Euler_substitution (accessed on 8 August 2023).
  2. Euler Substitutions. Encyclopedia of Mathematics. Available online: https://encyclopediaofmath.org/wiki/Euler_substitutions (accessed on 18 July 2023).
  3. Radożycki, T. Solving Problems in Mathematical Analysis; Part I, Section 14.5; Springer: Cham, Switzerland, 2020. [Google Scholar]
  4. Fikhtengol’ts, G.M. The Fundamentals of Mathematical Analysis; Chapter 10, Section 170; Pergamon Press: Oxford, UK, 1965; Volume 1. [Google Scholar]
  5. Piskunov, N. Differential and Integral Calculus; Mir Publishers: Moscow, Russia, 1969. [Google Scholar]
  6. Ryzhik, I.M.; Gradshteyn, I.S. Table of Integrals, Sums, Series, and Products; Gosudarstvennoe Izdatel’stvo Tehniko-Teoretičeskoj Literatury: Moscow, Russia, 1951. (In Russian) [Google Scholar]
  7. Byerly, W.E. Elements of the Integral Calculus; Ginn and Company: Boston, MA, USA, 1892. [Google Scholar]
  8. Fikhtengol’ts, G.M. Kurs Differentsial’nogo i Integral’nogo Ischisleniia; Chapter 8, Section 282; Fizmatgiz: Moscow, Russia, 1959; Volume 2. (In Russian) [Google Scholar]
  9. Lucht, P. The Euler Substitutions. Technical Report. Available online: https://www.researchgate.net/publication/309907150_The_ Eu\protect\discretionary{\char\hyphenchar\font}{}{}ler_Substitutions (accessed on 18 July 2023).
  10. Boyadzhiev, K.N. Special Techniques for Solving Integrals: Examples and Problems; World Scientific: Singapore, 2021. [Google Scholar]
  11. Euler, L. Institutiones Calculi Integralis; Part 1, Section 1, Chapter 2; Bruce, I., Translator; Imperial Academy of Sciences: St. Petersburg, Russia, 1768; Volume 1, Available online: http://www.17centurymaths.com/ (accessed on 8 August 2023).
  12. Besier, M.; van Straten, D.; Weinzierl, S. Rationalizing roots: An algorithmic approach. Commun. Number Theory Phys. 2019, 13, 253–297. [Google Scholar] [CrossRef]
  13. Besier, M.; Festi, D. Rationalizability of square roots. J. Symb. Comput. 2021, 106, 48–67. [Google Scholar] [CrossRef]
  14. Trinh, Y.T.; Hanzon, B. Option Pricing and CVA Calculations using the Monte Carlo-Tree (MC-Tree) Method. arXiv 2022, arXiv:2202.00785. [Google Scholar]
Figure 1. Geometric interpretation of the second Euler substitution in the case a < 0 and c > 0 . The point P is parameterized by the slope t of the line P 0 P .
Figure 1. Geometric interpretation of the second Euler substitution in the case a < 0 and c > 0 . The point P is parameterized by the slope t of the line P 0 P .
Symmetry 15 01932 g001
Figure 2. Geometric interpretation of the second Euler substitution in the case a > 0 and c > 0 . The point P is parameterized by the slope t of the line P 0 P .
Figure 2. Geometric interpretation of the second Euler substitution in the case a > 0 and c > 0 . The point P is parameterized by the slope t of the line P 0 P .
Symmetry 15 01932 g002
Figure 3. Geometric interpretation of the third Euler substitution in the case a < 0 and Δ > 0 . The point P is parameterized by the slope t of the line P 0 P .
Figure 3. Geometric interpretation of the third Euler substitution in the case a < 0 and Δ > 0 . The point P is parameterized by the slope t of the line P 0 P .
Symmetry 15 01932 g003
Figure 4. Geometric interpretation of the third Euler substitution in the case a > 0 . The point P is parameterized by the slope t of the line P 0 P .
Figure 4. Geometric interpretation of the third Euler substitution in the case a > 0 . The point P is parameterized by the slope t of the line P 0 P .
Symmetry 15 01932 g004
Figure 5. Geometric interpretation of the first Euler substitution. The points P and P 1 are parameterized by intersections t and t 1 , respectively, of the y-axis with the line parallel to one of the asymptotes of the hyperbola y 2 = a x 2 + b x + c .
Figure 5. Geometric interpretation of the first Euler substitution. The points P and P 1 are parameterized by intersections t and t 1 , respectively, of the y-axis with the line parallel to one of the asymptotes of the hyperbola y 2 = a x 2 + b x + c .
Symmetry 15 01932 g005
Figure 6. Characteristic points on the graph of an ellipse: intersections with the coordinate axes (provided that they exist) and extremes (minimum M 1 and maximum M 2 ).
Figure 6. Characteristic points on the graph of an ellipse: intersections with the coordinate axes (provided that they exist) and extremes (minimum M 1 and maximum M 2 ).
Symmetry 15 01932 g006
Figure 7. Characteristic points on the graphs of hyperbolas (two hyperbolas with the same | a | are presented): intersections with the coordinate axes ( V 1 , V 2 , R 1 , R 2 ) and extremes ( M 1 , M 2 ).
Figure 7. Characteristic points on the graphs of hyperbolas (two hyperbolas with the same | a | are presented): intersections with the coordinate axes ( V 1 , V 2 , R 1 , R 2 ) and extremes ( M 1 , M 2 ).
Symmetry 15 01932 g007
Figure 8. Geometric interpretation of the fourth Euler substitution in the case a > 0 . The point P is parameterized by the slope t of the line P 0 P , where P 0 = M 1 .
Figure 8. Geometric interpretation of the fourth Euler substitution in the case a > 0 . The point P is parameterized by the slope t of the line P 0 P , where P 0 = M 1 .
Symmetry 15 01932 g008
Figure 9. Geometric interpretation of the fourth Euler substitution in the case a < 0 . The point P is parameterized by the slope t of the line P 0 P , where P 0 = M 1 .
Figure 9. Geometric interpretation of the fourth Euler substitution in the case a < 0 . The point P is parameterized by the slope t of the line P 0 P , where P 0 = M 1 .
Symmetry 15 01932 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cieśliński, J.L.; Jurgielewicz, M. On Geometric Interpretations of Euler’s Substitutions. Symmetry 2023, 15, 1932. https://doi.org/10.3390/sym15101932

AMA Style

Cieśliński JL, Jurgielewicz M. On Geometric Interpretations of Euler’s Substitutions. Symmetry. 2023; 15(10):1932. https://doi.org/10.3390/sym15101932

Chicago/Turabian Style

Cieśliński, Jan L., and Maciej Jurgielewicz. 2023. "On Geometric Interpretations of Euler’s Substitutions" Symmetry 15, no. 10: 1932. https://doi.org/10.3390/sym15101932

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