Solution of Nonlinear Equations in Science through Lagrange’s Inversion Theorem
Pushpa N. Rathie1,, Prabhata K. Swamee2, Luan Carlos de S. M. Ozelim3
1Department of Statistics, University of Brasilia, Brasilia, Brazil
2Department of Civil Engineering, ITM University, Gurgaon, India
3Department of Civil and Environmental Engineering, University of Brasilia, Brasilia, Brazil
Abstract
Nonlinear problems arise in most of the scientific fields. In general, such behavior is represented by a nonlinear equation, whose solution is sought. Analytical and numerical methods have been applied to the solution of this class of equations, notwithstanding, in cases where highly nonlinear phenomena are analyzed, the number of iterations and computational effort necessary to achieve the minimum required accuracy is very high. Lagrange´s Inversion Theorem (LIT) has been applied to solve this kind of problems analytically, giving the solution as an infinite power series. This way, the accuracy can be as high as necessary by taking more terms from the series solution, which is easily computationally implemented. Also, in some cases it is possible to relate the series obtained to the expansion of special and elementary functions, which enables one to exactly solve the desired equation. In the present review paper, a total of eleven applications have been discussed in order to show the role of LIT in various areas of nonlinear sciences.
Keywords: Lagrange´s Inversion Theorem, Civil engineering, statistics, Graph theory, algebraic equations, chemical engineering
Applied Mathematics and Physics, 2013 1 (4),
pp 120-125.
DOI: 10.12691/amp-1-4-4
Received October 22, 2013; Revised November 05, 2013; Accepted November 07, 2013
Copyright © 2013 Science and Education Publishing. All Rights Reserved.Cite this article:
- Rathie, Pushpa N., Prabhata K. Swamee, and Luan Carlos de S. M. Ozelim. "Solution of Nonlinear Equations in Science through Lagrange’s Inversion Theorem." Applied Mathematics and Physics 1.4 (2013): 120-125.
- Rathie, P. N. , Swamee, P. K. , & Ozelim, L. C. D. S. M. (2013). Solution of Nonlinear Equations in Science through Lagrange’s Inversion Theorem. Applied Mathematics and Physics, 1(4), 120-125.
- Rathie, Pushpa N., Prabhata K. Swamee, and Luan Carlos de S. M. Ozelim. "Solution of Nonlinear Equations in Science through Lagrange’s Inversion Theorem." Applied Mathematics and Physics 1, no. 4 (2013): 120-125.
Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |
1. Introduction
Modeling natural phenomena is one of the greatest challenges in every branch of Science. At first, when very little is known about the subject of study, linear approximations and descriptions tend to be the first approach made by scientists. On the other hand, when the phenomena of interest are further investigated, a few peculiarities which, at first, did not seem to be important have to be taken into account for building better models.
In general, when these particularities are considered, nonlinear equations arise. In the process of applying the new relations to the description of phenomena, the success intrinsically relies on the solution of the former. In order to address this issue, analytical and numerical methods are considered.
Since the advent of modern computer sciences, the usage of numerical methods has tremendously increased. This comes from the fact that computers became powerful tools of multiprocessing procedures and routines. On the other hand, as widely known, modeling systems do not depend on the solution of a single equation, but on a great number of them. In special, while considering nonlinear equations, the number of iterations necessary to obtain a desired accuracy is very high. This way, even if the numerical solution of a single equation by standard trial and error methods is worth the computational effort, the solution of large sets surely is not.
Analytical approaches, on the other hand, tend to be applicable only to special cases, in which it is possible to explicit the variable of interest by means of arithmetic manipulation. In order to address more complicated problems, an analytical exact method which does not depend on arithmetical manipulation is the Lagrange´s Inversion Theorem (LIT), which gives the exact analytical solution of a nonlinear equation by means of an infinite series. The theorem can be readily stated as follows [27]: Let y be defined as the following function of constant, function φ, and a parameter δ:
(1) |
then any function ζ(y) is expressed as the following power series in δ:
(2) |
It can be noticed that the right hand side of Eq.(2) contains y through δ defined in Eq.(1).
It is evident that the convergence issues concerning the series in Eq.(2) have to be taken into account for consistency of the solution. Besides, it is worth noticing that the right hand side of Eq.(2) does not depend on y, this way if one takes ζ(y) = y, the once implicit function y is now explicit.
In the presented paper, the usage of both Eqs.(1) and (2) is shown by discussing eleven applications of them in Civil Engineering, Mathematics, Statistics, Graph Theory and Chemical Engineering.
2. Examples of Application in Science
In this section a brief review of where LIT has been applied is shown. In order to better structure the paper, the applications are grouped following the area they are related to.
2.1. Civil EngineeringIn the present subsection, applications in Civil engineering are shown. Since this paper is a review paper, the full theory behind each application is suppressed.
Regarding the alternate and sequent depths applications, one may refer to, for example, [1], [2] and [7].
Regarding the three-parameter infiltration equation application, one may refer to [12] and [11].
2.1.1. Alternate Depths in a Triangular Shaped Channel
In order to obtain the alternate depths for a triangular shaped channel, the following equation has to be taken into account [23]:
(3) |
in which, ec is a reduced energy parameter and z, a reduced depth parameter. Since Eq.(3) has two roots of interest (z1 and z2), one has to invert it accordingly, this way, in order to get the first root, Eq.(3) can be rearranged as:
(4) |
in which . By means of LIT with ζ(v) = v-1/4, the value of z1 is given as:
(5) |
On the other hand, in order to get the value of z2, Eq. (3) may be rearranged as:
(6) |
This way, by means of LIT with ζ(z) = z, the value of z2 is:
(7) |
Alternate Depths are a key parameter in the study of open channel flows.
2.1.2. Alternate Depths in a Parabolic Shaped Channel
In order to obtain the alternate depths for a parabolic shaped channel, the following equation has to be taken into account [23]:
(8) |
in which, ec is a reduced energy parameter and z, a reduced depth parameter. As in the last subsection, Eq.(8) has two roots of interest (z1 and z2) thus one has to invert it accordingly. In order to get the first root, Eq.(8) can be rearranged as:
(9) |
in which . By means of LIT with ζ(h) = h-1/3, the value of z1 is given as:
(10) |
On the other hand, in order to get the value of z2, Eq. (8) may be rearranged as:
(11) |
This way, by means of Eq. (11) and LIT with ζ(z) = z, the value of z2 is:
(12) |
2.1.3. Sequent Depths in a Triangular Shaped Channel
In order to obtain the sequent depths for a triangular shaped channel, the following equation has to be taken into account [24]:
(13) |
in which, mc is a reduced momentum parameter and z, a reduced depth parameter. There are two zeros of Eq. (13) which are of interest (z1 and z2) thus one has to invert it accordingly. In order to get the first root, Eq.(13) can be rearranged as:
(14) |
in which . By means of LIT with ζ(w) = w1/2, the value of z1 is given as:
(15) |
On the other hand, in order to get the value of z2, Eq. (13) may be rearranged as:
(16) |
in which . This way, by means of Eq. (16) and LIT with ζ(z) = r1/3, the value of z2 is:
(17) |
Sequent Depths, as the Alternate Depths, are key parameters in the study of open channel flow.
2.1.4. Sequent Depths in a Parabolic Shaped Channel
In order to obtain the sequent depths for a parabolic shaped channel, the following equation has to be taken into account [24]:
(18) |
in which, mc is a reduced momentum parameter and z, a reduced depth parameter. There are two zeros of Eq. (13) which are of interest (z1 and z2) thus one has to invert it accordingly. In order to get the first root, Eq.(18) can be rearranged as:
(19) |
in which . By means of LIT with ζ(q) = q2/3, the value of z1 is given as:
(20) |
On the other hand, in order to get the value of z2, Eq. (18) may be rearranged as:
(21) |
in which . This way, by means of Eq. (16) and LIT with ζ(s) = s2/5, the value of z2 is:
(22) |
2.1.5. Three-Parameter Infiltration Equation
In a recent paper, Rathie et al. [18] studied the three-parameter infiltration equation, stated as:
(23) |
in which I* is the nondimensional cumulative infiltration and t*, the nondimensional time. Also, α is a transition parameter which pertains the interval [0,1]. An approximate solution has been presented in [11].
In order to provide the nondimensional cumulative infiltration explicitely, Eq. (23) can be rearranged as:
(24) |
in which i = exp((1- α)I*); t = αexp(t*(α – 1)); and a = α – 1. By means of Lagrange´s Inversion Theorem, with ζ(i) = i1/(1-α), the value of interest is given as:
(25) |
Computational programs tend to misinterpret Eq.(25) as the gamma function arguments in the latter became negative integers. It is worth noticing that even if the arguments of individual gamma function in Eq. (25) became problematic, the ratio of the gamma function presented in the reffered equation does not. In order to provide a better way of implementing Eq.(25), let one define the Pochhammer symbol as:
(26) |
This way, by means of Eq. (26) and the multiplication formula for the gamma function, Eq.(25) is easily converted to:
(27) |
More applications of LIT in Civil Engineering have been discussed in [20].
2.2. Algebraic EquationsIn the present subsection, the applications related to obtaining solutions to algebraic equations are discussed. One may refer to [3], [6] and [8].
2.2.1. Solution to Real Degree Equations
Recently, Rathie and Ozelim [17] derived the general solution to two classes of real degree equations. The solutions were obtained by means of LIT and subsequently have been converted to H-functions in order to provide a closed-form representation. Since the aim of the present paper is to review the applications of LIT, only the latter shall be discussed. This way, consider the general real degree equation of the following type:
(28) |
in which α, β, γ and δ are arbitrary complex numbers and θ is a real number. By means of LIT, with ζ(x) = x, Eq. (28) can be solved for x as:
(29) |
On the other hand, consider the real degree equation of following type:
(30) |
in which β* and γ* are arbitrary complex numbers and α* and δ* are real numbers. It can be seen, by comparing Eq.(30) to Eq.(28) that, in general, one is not a special case of the other. Only if , , , and they become explicitly correlated. Consider that Eq.(30) can be rewritten as:
(31) |
in which .By means of LIT, with ζ(y) = yφ/α* and assuming sign(φ) = sign(α*), where φ is an arbitrary constant, Eq. (31) provides:
(32) |
In the present subsection, applications in statistics are discussed. In order to obtain further information about the generalized logistic distribution, one may refer to [15], [21] and [13].
2.3.1. New Skew Generalized Logistic Distribution
Let the probability density function of the Skew generalized Logistic Distribution be defined as follows [16]:
(33) |
where and . For the special case in which c = -1, Eq. (33) provides:
(34) |
The cumulative distribution function is obtained by means of Eq. (34) as:
(35) |
In order to get z as a function o H, Eq. (35) shall be rearranged as:
(36) |
Finally, by applying LIT to Eq. (36), z is explicitly given as:
(37) |
2.3.2. Generalized Loglogistic Distribution
Rathie et al. [19] discussed the cummulative distribution function of the Generalized Loglogistic distribution, defined as:
(38) |
where .
In order to get z as a function of F, Eq. (38) may be rewritten as:
(39) |
By applying LIT to Eq. (39), the following is obtained:
(40) |
In the present subsection, applications in Graph theory are discussed. For further information regarding the theory behind the enumeration of certain maps, one may refer to [10], [25] and [26].
2.4.1. Enumeration of Almost Cubic Maps
While studying the enumeration of almost cubic maps, the following equation has been obtained by Mathai and Rathie [9]:
(41) |
In order to solve Eq. (41) for T, the following rearrangement has to be taken into account:
(42) |
By applying LIT to Eq. (42), T is given as:
(43) |
It is worth noticing that Eq.(43) can be further simplified by means of the binomial expansion. Some applications of LIT in Graph Theory have been discussed by Rathie [14].
2.5. Chemical EngineeringIn the present subsection, applications in chemical engineering are discussed. One may refer to [4] and [5] for further information about the equations solved.
2.5.1. Friction Factor Problems Involving Laminar Flow of Bingham Plastic Fluids
In a recent paper, Swamee et al. [22] analyzed the following implicit equation related to friction factor problems involving laminar flow of Bingham plastic fluids:
(44) |
in which f is the parameter of interest and Re and He are constants. By means of LIT, by taking ζ(f Re) = (f Re), Eq. (44) is solved as:
(45) |
where P = 1 + 6Re/He. One may notice that the series in Eq. (45) is fast converging for cases in which P (or He/Re is small). This way, one has to find another series which fastly converges when He/Re is large. Thus, consider the following rearrangement of Eq. (44):
(46) |
Equation (46) can be readily inverted by means of LIT with ζ (f Re2/He) = f Re2/He as:
(47) |
This way, f is given explicitly. In order to define the applicability of both Eqs. (45) and (47), after an error analysis, Swamee et al. [22] showed that the former is valid for He/Re ≤ 30 and the latter, otherwise.
2.5.2. Diameter Problems Involving Laminar Flow of Bingham Plastic Fluids
Swamee et al. [22] also studied the diameter problems related to the laminar flow of Bingham plastic fluids by means of the following equation:
(48) |
in which T = D* - 4, D* is the parameter of interest and q is a constant. For small q, Eq. (48) can be inverted by means of LIT if one takes y = T2; ζ(y) = y-1/2; χ = π/4q; δ = π/128q and φ(y) = y-1/2(y-1/2 +32/3). This way, the solution process yields:
(49) |
On the other hand, for large values of q, by using LIT with y = T4; ζ(y) = y-1/4; χ = π/128q; δ = π/4q and φ(y) = y1/4/3 + y1/2; the value of D* is given as:
(50) |
This way, the series above give the diameter parameter D* explicitly. Swamee et al. [22] showed that, in order to diminish the error related to the applicability of both Eqs. (49) and (50), the former is valid for q ≤ 1.05 and the latter, otherwise.
3. Conclusions
Nonlinear problems are the core of most of Science’s branches. This way, analytical tools to provide exact solutions to nonlinear equations play a major role in each of these fields. Lagrange´s Inversion Theorem (LIT) shows up as a good alternative in the solution process. A total of eleven applications in four different fields have been discussed, and the importance of the solutions provided by LIT is widely recognized.
Conflict of Interest
The authors declare that there is no conflict of interests regarding the publication of this article.
References
[1] | Bakhmeteff, B. A. Hydraulics of open channels, McGraw-Hill, New York, N.Y., 1932. | ||
In article | |||
[2] | Chow, V. T. Open channel hydraulics, McGraw-Hill, Inc., New York, N.Y., 1959. | ||
In article | PubMed | ||
[3] | Cockle, J. Sketch of a theory of transcendental roots, London Edinburgh Dublin Philos. Mag. J. Sci. 20, pp.145-148, 1860. | ||
In article | |||
[4] | Darby R,. Chemical Engineering Fluid Mechanics., Marcel Dekker, New York City, New York, pp 164-169, 1996. | ||
In article | |||
[5] | Darby, R., Melson, J. How to predict the friction factor for flow of Bingham plastics. Chemical Engineering, 28: pp. 59-61, 1981. | ||
In article | |||
[6] | Glasser, M. L. Hypergeometric functions and the trinomial equation, J. Comput. Appl. Math. 118, pp.169-173, 2000. | ||
In article | CrossRef | ||
[7] | Henderson, F. M. Open channel flow, The MacMillan Co., New York, N.Y., 1966. | ||
In article | |||
[8] | King, R. B. Beyond the Quartic Equation, Birkhäuser, Boston, 1996. | ||
In article | |||
[9] | Mathai, A. M.; Rathie, P. N., Enumeration of Almost Cubic Maps, Journal of Combinatorial Theory, Series B, Vol. 13 (1), p. 83-90, 1972. | ||
In article | CrossRef | ||
[10] | Mullin, R. C. On the average number of trees in certain maps. Canad. J. Math., 18, pp. 33-41, 1966. | ||
In article | CrossRef | ||
[11] | Parlange, J. Y., Barry, D. A., Haverkamp, R. Explicit infiltration equations and the Lambert W-function, Adv. Water Resour., 25 (8-12), pp. 1119-1124, 2002. | ||
In article | |||
[12] | Parlange, J. Y., Lisle, I., Braddock, R. D., Smith, R. E. The three parameter infiltration equation, Soil Sci., 133 (6), pp. 337-341, 1982. | ||
In article | CrossRef | ||
[13] | Rathie, P. N. Normal distribution, Univariate, International Encyclopedia of Statistical Sciences, Springer Verlag, Part 14, pp. 1012-1013, 2011. | ||
In article | |||
[14] | Rathie, P. N., On Some Enumeration Problems in Graph Theory, Annals of Discrete Mathematics, vol. 6, p. 281-288, 1980. | ||
In article | CrossRef | ||
[15] | Rathie, P. N., Swamee, P. K. On a new invertible generalized logistic distribution approximation to normal distribution, Technical Research Report in Statistics, 07/2006, Dept. of Statistics, Univ. of Brasilia, Brasilia, Brazil. 2006. | ||
In article | |||
[16] | Rathie, P. N.; Coutinho, M. A new Skew Generalized Logistic distribution and Approximation to Skew Normal Distribution. Aligarh Journal of Statistics, v.31, p.1-12, 2011. | ||
In article | |||
[17] | Rathie, P. N.; Ozelim, L. C. S. M. . General Solutions to Certain Real Degree Equations and their Applications. Advanced Studies in Contemporary Mathematics (Kyungshang), v. 22, p. 325-341, 2012. | ||
In article | |||
[18] | Rathie, P. N.; Swamee, P. K.; Cavalcante, A. L. B.; Ozelim, L. C. S. M., Lagrange’s Inversion Theorem and Infiltration. World Academy of Science, Engineering and Technology, v. 67, p. 473-478, 2012. | ||
In article | |||
[19] | Rathie, P. N.; Swamee, P. K.; Coutinho, M., Generalized Loglogistic Distribution and Its Applications. In: The 5th International Conference MSAST 2011 of IMBIC, 2011, Kolkata. Proceedings of The 5th International Conference MSAST 2011 of IMBIC, p. 88-99, 2011. | ||
In article | |||
[20] | Rathie, P. N.; Swamee, P. K.; Ozelim, L. C. S. M.. Closed-form exact solutions in hydraulics by means of the H-Function. In: The 9th Annual Conference of the Society for Special Function and Their Application (SSFA), 2010, Gwalior. Proceedings of the 9th Annual Conference of the Society for Special Function and Their Application (SSFA), v. 9. p. 3-14, 2010. | ||
In article | |||
[21] | Swamee, P. K., Rathie, P. N., Invertible Alternatives to Normal and Lognormal Distributions, Journal of Hydrologic Engineering, 12(2), pp. 218-221, 2007. | ||
In article | CrossRef | ||
[22] | Swamee, P. K.; Rathie, P. N.; Aggarwal, N., Exact solution of friction factor and diameter problems involving laminar flow of Bingham plastic fluids. Journal of Petroleum and Gas Exploration Research. v. 2(2), p. 27-32, 2012. | ||
In article | |||
[23] | Swamee, P. K.; Rathie, P. N.; Ozelim, L. C. S. M., Analytical Solutions for Alternate Depths. ISH Journal of Hydraulic Engineering, v. 17, p. 34-42, 2011. | ||
In article | CrossRef | ||
[24] | Swamee, P. K.; Rathie, P. N.; Ozelim, L. C. S. M., Application of Analytical and Numerical Methods to the Sequent Depths Problem in Civil Engineering. Australian Journal of Basic and Applied Sciences, v. 5, p. 38-47, 2011. | ||
In article | |||
[25] | Tutte, W. T. A census of planar maps. Canad. J. Math., 15, pp. 249-271. 1963. | ||
In article | CrossRef | ||
[26] | Tutte, W. T. On the enumeration of planar maps. Bull. Amer. Math. Soc., 74, pp. 64-74. 1968. | ||
In article | CrossRef | ||
[27] | Whittaker E. T., Watson G. N., A Course of Modern Analysis. 4th Edition, Cambridge Univ. Press, Cambridge, U.K., pg. 133, 1990. | ||
In article | |||