International Journal of Engineering Mathematics
Volume 2015 (2015), Article ID 523043, 4 pages
http://dx.doi.org/10.1155/2015/523043

Research Article

An Efficient Method to Find Solutions for Transcendental Equations with Several Roots

Rogelio Luck,¹ Gregory J. Zdaniuk,² and Heejin Cho¹

¹Department of Mechanical Engineering, Mississippi State University, P.O. Box 9552, Mississippi State, MS 39762, USA
²Cofely UK, Kings Yard, 1 Waterden Road, Queen Elizabeth Olympic Park, London E15 2GP, UK

Received 17 July 2015; Revised 6 October 2015; Accepted 7 October 2015

Academic Editor: Song Cen

Copyright © 2015 Rogelio Luck et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents a method for obtaining a solution for all the roots of a transcendental equation within a bounded region by finding a polynomial equation with the same roots as the transcendental equation. The proposed method is developed using Cauchy’s integral theorem for complex variables and transforms the problem of finding the roots of a transcendental equation into an equivalent problem of finding roots of a polynomial equation with exactly the same roots. The interesting result is that the coefficients of the polynomial form a vector which lies in the null space of a Hankel matrix made up of the Fourier series coefficients of the inverse of the original transcendental equation. Then the explicit solution can be readily obtained using the complex fast Fourier transform. To conclude, the authors present an example by solving for the first three eigenvalues of the 1D transient heat conduction problem.

1. Introduction

The determination of roots of transcendental functions is a problem commonly encountered in a broad spectrum of engineering applications, such as heat transfer (e.g., [1–4]), dynamics and control (e.g., [5]), sound and vibration (e.g., [6, 7]), and quantum mechanics (e.g., [8]). There exist a wide variety of numerical methods that are useful for approximating the solution to any desired degree of accuracy. General descriptions of common root finding techniques are available in most textbooks in the area of numerical analysis, for example, Chapra and Canale [9].

Explicit methods have also been developed to find roots of transcendental equations. Explicit solutions can provide more insight into the problem under consideration. They are also useful, for example, in the development of analytical derivatives for uncertainty analyses and sensitivity studies or for checking convergence of approximate root finding techniques [3]. Leathers and McCormick [10] addressed a methodology for obtaining explicit solutions to several transcendental problems arising in heat transfer. The general approach used was based on methods of Muskhelishvili [11] and was developed by Burniston and Siewert [12]. This approach “depends on formulating an appropriate Riemann problem of complex variable theory and then expressing the solution of the transcendental equation in terms of a canonical solution of that problem” [10]. Luck and Stevens [13] formulated an explicit expression for a single root of any analytic transcendental function. Their method is based on Cauchy’s integral theorem and uses only basic concepts of complex integration to determine the root of a function by locating the singularity of the reciprocal of the function. It has now been found that this idea was initially proposed by Jackson [14, 15]. The current paper extends their idea presenting a method for finding several roots of transcendental equations within a bounded region by transforming the problem into an equivalent problem of finding roots of a polynomial equation with exactly the same roots. It follows that an explicit solution is readily available through the complex fast Fourier transform (cfft) when dealing with up to four roots at a time. For five or more roots the solution is implicit according to the Abel-Ruffini theorem [16]. The implicitness of finding roots of fifth- and higher-order polynomials is not really an issue as algorithms designed for this task are readily available (McNamee [17–20], McNamee and Pan [21]) or larger problems can be divided into smaller problems by analyzing nonoverlapping bounded regions with four roots or less in each region.

2. Method Development

Suppose is a transcendental function with roots within an interval of width and center such that . Let be an th-order polynomial with roots located at exactly the same locations as the roots of in the interval . The ratio is an analytic function in the interval because the poles of in the interval are cancelled out by the roots of . The next step in the analysis is to normalize the variable of interest. Let

such that . Functions and can be recast as

where

and is the number of roots in the interval.

Now consider the function in the interval . If the analysis is extended into the complex plane (by substituting the real variable with the complex variable ), integration around a circular path of radius 1 yields

which is a consequence of Cauchy’s theorem regarding the line integral of a function that is analytic on and inside of the contour defined by the closed path of integration.

Equation (4) can be used to find the coefficients of polynomial as follows:

Note that, given a positive integer , it is also true that

because the function being integrated remains analytic when multiplied by . Following the approach presented in (5)

Following is the procedure for obtaining coefficients . First, the integrals are evaluated by substituting

Let

It follows that

Another simplification can be implemented by defining

Note that can be interpreted as the th complex Fourier series coefficient of . Finally, substituting into (7) and cancelling yield

Equation (12) is the equation used to find coefficients . Note that (12) is true for any positive integer value of . This means that one can “create” as many equations as there are unknown coefficients . Varying from 0 to , the following system of equations is obtained:

Because (13) represents a system of linearly dependent equations, that is, having fewer equations than unknowns, additional equations are necessary to obtain a solution. One possible equation consists of setting equal to any real nonzero number, which allows moving the last column of the matrix to the right hand side.

Equation (13) is worthy of comment. Matrix is composed solely of Fourier series coefficients, implying that the equations do not depend on normalization factors used to compute the Fourier series coefficients. In practice, the coefficients of are approximated through the complex fast Fourier transform (cfft). The cfft approximation is less accurate for the last column of the matrix because this column corresponds to high-frequency terms. It follows that moving the last column of the matrix to the right hand side of the equation results in a slightly more robust solution. Furthermore, the system of equations constitutes a Hankel matrix for which an inverse can be easily obtained [22]. Therefore, (13) can be quickly implemented and solved for the coefficients . Finally, the symmetry or “structure” of (13) leads the authors to believe there is a deeper meaning to this system of equations.

3. Example

To illustrate the current method, a simple yet helpful example is described below. The solution of the unsteady heat conduction equation yields the characteristic equation or eigenfunction of the following form [23]:

Equation (14) is a transcendental equation with an infinite number of roots called characteristic values or eigenvalues. The characteristic equation in this case is implicit. Arbitrarily assuming , (14) can be rewritten as

Let

The plot of is shown on Figure 1.

Figure 1: Plot of .

The method described in the previous section is now applied to solve for the first three positive roots of (16). Figure 1 illustrates that the three roots occur in the interval , suggesting , , and . Note that one simple way of approximating coefficients is to discretize for and to use the complex fast Fourier transform (cfft) readily available in many commercial mathematical software packages (e.g., MATLAB and Mathcad), rather than performing the integral of (11). With and , (13) becomes

Discretizing into 256 equidistant points, the cfft yields

Solving this system of equations gives

which translates into

Figure 2 plots the original transcendental function and the polynomial . The plot clearly shows that the polynomial has the same roots as the original function ; that is, , , and . The roots of a cubic polynomial can be found using the explicit algebraic expressions given by many authors, for example, Abramowitz and Stegun [24], Press et al. [25], and Weisstein [26]. Finally, undoing the variable normalization yields , , and as the roots of or eigenvalues of the original problem.

Figure 2: Plot of and with locations of the roots.

4. Conclusion

The current paper presents a method for finding simultaneously several roots of a transcendental equation within a bounded region, by transforming the problem into an equivalent problem of finding roots of a polynomial equation with exactly the same roots. The method requires that the ratio be analytic in the interval under study. The concept shown here is that the correlation between the roots of a transcendental function and the Fourier series coefficients of its inverse can be expressed with a polynomial and proven with Cauchy’s integral theorem.

The advantage of the current method is the possibility of solving for several roots in one shot and is most evident when dealing with an implicit problem as shown in the previous example. In this case, an implicit problem was transformed into a simpler explicit one. The current method remains explicit when dealing with up to four roots. However, the implicitness of finding roots of higher-order polynomials is not really an issue as algorithms designed for this task are readily available or larger problems can be divided into smaller problems with four roots each by selecting appropriate intervals as illustrated in Section 3. In any case, the current method provides a relatively simple solution procedure for a group of problems encountered in many engineering applications.

Notation

: | Biot number
: | th coefficient of the polynomial
: | Function of
: | Function of
: | Function of
: | th Fourier series coefficient of function
: | Center of a circle in the complex plane
: | Number of roots
: | Polynomial, function of
: | Polynomial, function of
: | Positive integer
: | Radius of a circle in the complex plane
: | Independent variable
: | Normalized variable
: | Complex variable
: | Angle in the complex plane.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

J. S. Phelps III and M. J. Funcheon, “Solutions of a transcendental equation basic to the study of composite media heat transfer,” Journal of Computational and Applied Mathematics, vol. 10, no. 2, pp. 135–139, 1984.

J. W. Stevens and R. Luck, “Explicit approximations for all eigenvalues of the 1-D transient heat conduction equations,” Heat Transfer Engineering, vol. 20, no. 2, pp. 35–41, 1999.

A. Haji-Sheikh and J. V. Beck, “An efficient method of computing eigenvalues in heat conduction,” Numerical Heat Transfer Part B: Fundamentals, vol. 38, no. 2, pp. 133–156, 2000.

F. de Monte, “Multi-layer transient heat conduction using transition time scales,” International Journal of Thermal Sciences, vol. 45, no. 9, pp. 882–892, 2006.

S. Ruan and J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with two delays,” Dynamics of Continuous, Discrete and Impulsive Systems, Series A: Mathematical Analysis, vol. 10, no. 6, pp. 863–874, 2003.

N. P. Aleshin, V. S. Kamenskii, and L. I. Mogil'ner, “Solution of a transcendental equation encountered in diffraction problems,” Journal of Applied Mathematics and Mechanics, vol. 47, no. 1, pp. 139–141, 1983.

C. M. Ruggiero, “Solution of transcendental and algebraic equations with applications to wave propagation in elastic plates,” NASA STI/Recon Technical Report 81, 1981.

G. Pickett and Y. Millev, “On the analytic inversion of functions, solution of transcendental equations and infinite self-mappings,” Journal of Physics A: Mathematical and General, vol. 35, no. 20, pp. 4485–4494, 2002.

S. C. Chapra and R. P. Canale, Numerical Methods for Engineers, McGraw-Hill, New York, NY, USA, 3rd edition, 1998.

R. A. Leathers and N. J. McCormick, “Closed-form solutions for transcendental equations of heat transfer,” Journal of Heat Transfer, vol. 118, no. 4, pp. 970–973, 1996.

N. I. Muskhelishvili, Singular Integral Equations, Noordhoff, Groningen, The Netherlands, 1953.

E. E. Burniston and C. E. Siewert, “The use of Riemann problems in solving a class of transcendental equations,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 73, no. 1, pp. 111–118, 1973.

R. Luck and J. W. Stevens, “Explicit solutions for transcendental equations,” SIAM Review, vol. 44, no. 2, pp. 227–233, 2002.

D. Jackson, “Non-essential singularities of functions of several complex variables,” The Annals of Mathematics, vol. 17, no. 4, pp. 172–179, 1916.

D. Jackson, “Roots and singular points of semi-analytic functions,” The Annals of Mathematics, vol. 19, no. 2, pp. 142–151, 1917.

E. Dehn, Algebraic Equations: An Introduction to the Theories of Lagrange and Galois, Columbia University Press, New York, NY, USA, 1930.

J. M. McNamee, “A bibliography on roots of polynomials,” Journal of Computational and Applied Mathematics, vol. 47, no. 3, pp. 391–394, 1993.

J. M. McNamee, “A supplementary bibliography on roots of polynomials,” Journal of Computational and Applied Mathematics, vol. 78, no. 1, p. 1, 1997.

J. M. McNamee, “An updated supplementary bibliography on roots of polynomials,” Journal of Computational and Applied Mathematics, vol. 110, no. 2, pp. 305–306, 1999.

J. M. McNamee, Numerical Methods for Roots of Polynomials, Part 1, Studies in Computational Mathematics, vol. 14, Elsevier Science, Cambridge, Mass, USA, 2007.

J. M. McNamee and V. Pan, Numerical Methods for Roots of Polynomials, Part II, vol. 16 of Studies in Computational Mathematics, Elsevier Science, Cambridge, Mass, USA, 2013.

G. H. Golub and C. F. Van Loan, Matrix Computations, Johns Hopkins University Press, Baltimore, Md, USA, 3rd edition, 1996.

Y. A. Cengel, Heat Transfer: A Practical Approach, McGraw-Hill, New York, NY, USA, 3rd edition, 2005.

M. Abramowitz and I. A. Stegun, Eds., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, New York, NY, USA, 9th edition, 1972.

H. Press, S. A. Teukolsky, W. T. Vetterling, and B. R. Flannery, Numerical Recipes in Fortran 77: The Art of Scientific Computing, Cambridge University Press, 1992.

E. W. Weisstein, CRC Concise Encyclopedia of Mathematics, CRC Press, Boca Raton, Fla, USA, 2nd edition, 2002.