Rational biparameter homotopy perturbation method and Laplace-Padé coupled version. (English) Zbl 1268.65094

Vazquez-Leal, Hector; Sarmiento-Reyes, Arturo; Khan, Yasir; Filobello-Nino, Uriel; Diaz-Sanchez, Alejandro

J. Appl. Math. 2012, Article ID 923975, 21 p. (2012).

Summary: The fact that most of the physical phenomena are modelled by nonlinear differential equations underlines the importance of having reliable methods for solving them. This work presents the rational biparameter homotopy perturbation method (RBHPM) as a novel tool with the potential to find approximate solutions for nonlinear differential equations. The method generates the solutions in the form of a quotient of two power series of different homotopy parameters. Besides, in order to improve accuracy, we propose the Laplace-Padé rational biparameter homotopy perturbation method (LPRBHPM), when the solution is expressed as the quotient of two truncated power series. The usage of the method is illustrated with two case studies. On one side, a Riccati nonlinear differential equation is solved and a comparison with the homotopy perturbation method (HPM) is presented. On the other side, a nonforced Van der Pol oscillator is analysed and we compare results obtained with RBHPM, LPRBHPM, and HPM in order to conclude that the LPRBHPM and RBHPM methods generate the most accurate approximated solutions.

MSC:

65L05	Numerical methods for initial value problems involving ordinary differential equations
34A34	Nonlinear ordinary differential equations and systems
34C10	Oscillation theory, zeros, disconjugacy and comparison theory for ordinary differential equations

Keywords:

numerical examples; rational biparameter homotopy perturbation method; truncated power series; nonlinear differential equation; Van der Pol oscillator

References:

[1]	J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2221-2230, 2009. · Zbl 1189.65246
[2]	J. Biazar and H. Ghazvini, “Convergence of the homotopy perturbation method for partial differential equations,” Nonlinear Analysis. Real World Applications, vol. 10, no. 5, pp. 2633-2640, 2009. · Zbl 1173.35395
[3]	D. A and T. Özi\cs, “An analytical study for Fisher type equations by using homotopy perturbation method,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 602-609, 2010. · Zbl 1201.65187
[4]	T. Özi\cs and D. A, “He’s homotopy perturbation method for solving heat-like and wave-like equations with variable coefficients,” Physics Letters A, vol. 372, no. 38, pp. 5944-5950, 2008. · Zbl 1223.35294
[5]	A. Yildirim and T. Özi\cs, “Solutions of singular IVPs of Lane-Emden type by homotopy perturbation method,” Physics Letters A, vol. 369, no. 1-2, pp. 70-76, 2007. · Zbl 1209.65120
[6]	T. Özi\cs and A. Yildirim, “A note on He’s homotopy perturbation method for van der Pol oscillator with very strong nonlinearity,” Chaos, Solitons and Fractals, vol. 34, no. 3, pp. 989-991, 2007.
[7]	H. Vázquez-Leal, “Rational homotopy perturbation method,” Journal of Applied Mathematics, vol. 2012, Article ID 490342, 14 pages, 2012. · Zbl 1251.65119
[8]	A. Barari, M. Omidvar, A. R. Ghotbi, and D. D. Ganji, “Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations,” Acta Applicandae Mathematicae, vol. 104, no. 2, pp. 161-171, 2008. · Zbl 1197.34015
[9]	J. H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527-539, 2004. · Zbl 1062.65074
[10]	J. H. He, “An elementary introduction to the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 410-412, 2009. · Zbl 1165.65374
[11]	J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257-262, 1999. · Zbl 0956.70017
[12]	J. H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287-292, 2004. · Zbl 1039.65052
[13]	J. H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73-79, 2003. · Zbl 1030.34013
[14]	J. H. He, “Coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37-43, 2000. · Zbl 1068.74618
[15]	J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695-700, 2005. · Zbl 1072.35502
[16]	J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87-88, 2006. · Zbl 1195.65207
[17]	Y. Khan, H. Vazquez-Leal, and Q. Wu, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications. In press.
[18]	J. Biazar and M. Eslami, “A new homotopy perturbation method for solving systems of partial differential equations,” Computers and Mathematics with Applications, vol. 62, no. 1, pp. 225-234, 2011. · Zbl 1228.65199
[19]	J. Biazar and B. Ghanbari, “The homotopy perturbation method for solving neutral functional-differential equations with proportional delays,” Journal of King Saud University, vol. 24, no. 1, pp. 33-37, 2012.
[20]	D. D. Ganji, H. Tari, and M. B. Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 1018-1027, 2007. · Zbl 1141.65384
[21]	M. Sheikholeslami, H. R. Ashorynejad, D. D. Ganji, and A. Yildirim, “Homotopy perturbation method for three-dimensional problem of condensation film on inclined rotating disk,” Scientia Iranica, vol. 19, no. 3, pp. 437-442, 2012.
[22]	D. D. Ganji and A. Sadighi, “Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24-34, 2007. · Zbl 1120.65108
[23]	D. D. Ganji, “The application of He’s homotopy perturbation method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 355, no. 4-5, pp. 337-341, 2006. · Zbl 1255.80026
[24]	M. Rafei, D. D. Ganji, and H. Daniali, “Solution of the epidemic model by homotopy perturbation method,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 1056-1062, 2007. · Zbl 1112.92054
[25]	D. D. Ganji and A. Rajabi, “Assessment of homotopy-perturbation and perturbation methods in heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 33, no. 3, pp. 391-400, 2006.
[26]	N. Faraz and Y. Khan, “Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface,” Ain Shams Engineering Journal, vol. 2, no. 34, pp. 221-226, 2011.
[27]	Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers and Mathematics with Applications, vol. 61, no. 11, pp. 3391-3399, 2011. · Zbl 1222.76014
[28]	M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1937-1945, 2011. · Zbl 1219.65121
[29]	H. Aminikhah, “The combined laplace transform and new homotopy perturbation methods for stiff systems of odes,” Applied Mathematical Modelling, vol. 36, no. 8, pp. 3638-3644, 2012. · Zbl 1252.34018
[30]	F. I. Compean, D. Olvera, F. J. Campa, L. N. Lopez de Lacalle, A. Elias-Zuniga, and C. A. Rodriguez, “Characterization and stability analysis of a multivariable milling tool by the enhanced multistage homotopy perturbation method,” International Journal of Machine Tools and Manufacture, vol. 57, pp. 27-33, 2012.
[31]	A. M. A. El-Sayed, A. Elsaid, I. L. El-Kalla, and D. Hammad, “A homotopy perturbation technique for solving partial differential equations of fractional order in finite domains,” Applied Mathematics and Computation, vol. 218, no. 17, pp. 8329-8340, 2012. · Zbl 1245.65141
[32]	Y. Khan, H. Vazquez-Leal, and N. Faraz, “An auxiliary parameter method using adomian polynomials and laplace transformation for nonlinear differential equations,” Applied Mathematical Modelling. In press. · Zbl 1352.65172
[33]	Y. Khan, H. Vazquez-Leal, and L. Hernandez-Martinez, “Removal of noise oscillation term appearing in the nonlinear equation solution,” Journal of Applied Mathematics, vol. 2012, Article ID 387365, 9 pages, 2012. · Zbl 1251.34059
[34]	S. T. Mohyud-Din, A. Yildirim, and M. Inc, “Coupling of homotopy perturbation and modified Lindstedt-Poincaré methods for traveling wave solutions of the nonlinear Klein-Gordon equation,” Journal of King Saud University, vol. 24, no. 2, pp. 187-191, 2012.
[35]	H. Vazquez-Leal, R. Castaneda-Sheissa, U. Filobello-Nino, A. Sarmiento-Reyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution related integrals,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012. · Zbl 06173038
[36]	H. Vázquez-Leal, U. Filobello-Niño, R. Castañeda-Sheissa, L. Hernández-Martínez, and A. Sarmiento-Reyes, “Modified HPMS inspired by homotopy continuation methods,” Mathematical Problems in Engineering, vol. 2012, Article ID 309123, 19 pages, 2012. · Zbl 1264.65131
[37]	U. Filobello-Nino, H. Vazquez-Leal, R. Castaneda-Sheissa et al., “An approximate solution of blasius equation by using HPM method,” Asian Journal of Mathematics and Statistics, vol. 5, pp. 50-59, 2012.
[38]	J. Biazar, F. Badpeima, and F. Azimi, “Application of the homotopy perturbation method to Zakharov-Kuznetsov equations,” Computers and Mathematics with Applications, vol. 58, no. 11-12, pp. 2391-2394, 2009. · Zbl 1189.65244
[39]	J. Biazar and H. Ghazvini, “Homotopy perturbation method for solving hyperbolic partial differential equations,” Computers and Mathematics with Applications, vol. 56, no. 2, pp. 453-458, 2008. · Zbl 1155.65395
[40]	J. Biazar and H. Ghazvini, “Exact solutions for non-linear Schrödinger equations by He’s homotopy perturbation method,” Physics Letters A, vol. 366, no. 1-2, pp. 79-84, 2007. · Zbl 1203.65207
[41]	J. Biazar, H. Ghazvini, and M. Eslami, “He’s homotopy perturbation method for systems of integro-differential equations,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1253-1258, 2009. · Zbl 1197.65106
[42]	J. Biazar, M. Eslami, and H. Aminikhah, “Application of homotopy perturbation method for systems of Volterra integral equations of the first kind,” Chaos, Solitons and Fractals, vol. 42, no. 5, pp. 3020-3026, 2009. · Zbl 1198.65254
[43]	J. Biazar and H. Ghazvini, “He’s homotopy perturbation method for solving systems of Volterra integral equations of the second kind,” Chaos, Solitons and Fractals, vol. 39, no. 2, pp. 770-777, 2009. · Zbl 1197.65219
[44]	J. Biazar and H. Ghazvini, “Numerical solution for special non-linear Fredholm integral equation by HPM,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 681-687, 2008. · Zbl 1132.65115
[45]	U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “Hpm applied to solve nonlinear circuits: a study case,” Applied Mathematical Sciences, vol. 6, no. 85-88, pp. 4331-4344, 2012. · Zbl 1263.37081
[46]	Y. C. Jiao, Y. Yamamoto, C. Dang, and Y. Hao, “An aftertreatment technique for improving the accuracy of Adomian’s decomposition method,” Computers and Mathematics with Applications, vol. 43, no. 6-7, pp. 783-798, 2002. · Zbl 1005.34006
[47]	M. A. El-Tawil, A. A. Bahnasawi, and A. Abdel-Naby, “Solving Riccati differential equation using Adomian’s decomposition method,” Applied Mathematics and Computation, vol. 157, no. 2, pp. 503-514, 2004. · Zbl 1054.65071
[48]	S. H. Behiry, H. Hashish, I. L. El-Kalla, and A. Elsaid, “A new algorithm for the decomposition solution of nonlinear differential equations,” Computers and Mathematics with Applications, vol. 54, no. 4, pp. 459-466, 2007. · Zbl 1126.65059
[49]	Y.-G. Wang, W.-H. Lin, and N. Liu, “A homotopy perturbation-based method for large deflection of a cantilever beam under a terminal follower force,” International Journal for Computational Methods in Engineering Science and Mechanics, vol. 13, pp. 197-201, 2012.
[50]	J. H. He, “An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineering,” International Journal of Modern Physics B, vol. 22, no. 21, pp. 3487-3578, 2008. · Zbl 1149.76607
[51]	M. Merdan, A. Gökdo\vgan, and A. Yildirim, “On the numerical solution of the model for HIV infection of CD4+ T cells,” Computers and Mathematics with Applications, vol. 62, no. 1, pp. 118-123, 2011. · Zbl 1228.65137
[52]	G. A. Baker, Essentials of Padé Approximations, Academic Express, London, UK, 1975. · Zbl 0315.41014
[53]	N. H. Sweilam and M. M. Khader, “Exact solutions of some coupled nonlinear partial differential equations using the homotopy perturbation method,” Computers and Mathematics with Applications, vol. 58, no. 11-12, pp. 2134-2141, 2009. · Zbl 1189.65259
[54]	S. Momani, G. H. Erjaee, and M. H. Alnasr, “The modified homotopy perturbation method for solving strongly nonlinear oscillators,” Computers and Mathematics with Applications, vol. 58, no. 11-12, pp. 2209-2220, 2009. · Zbl 1189.65168
[55]	Y. Khan and N. Faraz, “Application of modified Laplace decomposition method for solving boundary layer equation,” Journal of King Saud University, vol. 23, no. 1, pp. 115-119, 2011.
[56]	D. Bahuguna, A. Ujlayan, and D. N. Pandey, “A comparative study of numerical methods for solving an integro-differential equation,” Computers and Mathematics with Applications, vol. 57, no. 9, pp. 1485-1493, 2009. · Zbl 1186.65158
[57]	S. Momani and V. S. Ertürk, “Solutions of non-linear oscillators by the modified differential transform method,” Computers and Mathematics with Applications, vol. 55, no. 4, pp. 833-842, 2008. · Zbl 1142.65058
[58]	A. Gökdoêan, M. Merdan, and A. Yildirim, “The modified algorithm for the differential transform method to solution of Genesio systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 45-51, 2012.
[59]	P. Y. Tsai and C. K. Chen, “An approximate analytic solution of the nonlinear Riccati differential equation,” Journal of the Franklin Institute, vol. 347, no. 10, pp. 1850-1862, 2010. · Zbl 1210.34016
[60]	A. E. Ebaid, “A reliable aftertreatment for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 528-536, 2011. · Zbl 1221.34009
[61]	W. H. Enright, K. R. Jackson, S. P. Norsett, and P. G. Thomsen, “Interpolants for runge-kutta formulas,” ACM Transactions on Mathematical Software, vol. 12, no. 3, pp. 193-218, 1986. · Zbl 0617.65068
[62]	E. Fehlberg, “Klassische runge-kutta-formeln vierter und niedrigerer ordnung mit schrittweiten-kontrolle und ihre anwendung auf waermeleitungsprobleme,” Computing, vol. 6, pp. 61-71, 1970. · Zbl 0217.53001
[63]	K. Reif, K. Weinzierl, A. Zell, and R. Unbehauen, “A homotopy approach for nonlinear control synthesis,” IEEE Transactions on Automatic Control, vol. 43, no. 9, pp. 1311-1318, 1998. · Zbl 0957.93014
[64]	N. Aslam and A. K. Sunol, “Homotopy continuation based prediction of azeotropy in binary and multicomponent mixtures through equations of state,” Physical Chemistry Chemical Physics, vol. 6, no. 9, pp. 2320-2326, 2004.
[65]	F. Kubler and K. Schmedders, “Computing equilibria in stochastic finance economies,” Computational Economics, vol. 15, no. 1-2, pp. 145-172, 2000. · Zbl 0969.91008
[66]	R. C. Melville, L. Trajkovic, S. C. Fang, and L. T. Watson, “Artificial parameter homotopy methods for the DC operating point problem,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 12, no. 6, pp. 861-877, 1993. · Zbl 05448987
[67]	K. Yamamura, T. Sekiguchi, and Y. Inoue, “A fixed-point homotopy method for solving modified nodal equations,” IEEE Transactions on Circuits and Systems I, vol. 46, no. 6, pp. 654-665, 1999. · Zbl 0952.94016
[68]	R. Geoghegan, J. C. Lagarias, and R. C. Melville, “Threading homotopies and DC operating points of nonlinear circuits,” SIAM Journal on Optimization, vol. 9, no. 1, pp. 159-178, 1999. · Zbl 0994.65058
[69]	J. Lee and H. D. Chiang, “Constructive homotopy methods for finding all or multiple dc operating points of nonlinear circuits and systems,” IEEE Transactions on Circuits and Systems I, vol. 48, no. 1, pp. 51-66, 2001. · Zbl 0997.65082
[70]	A. Ushida, Y. Yamagami, Y. Nishio, I. Kinouchi, and Y. Inoue, “An efficient algorithm for finding multiple DC solutions based on the SPICE-oriented Newton homotopy method,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 21, no. 3, pp. 337-348, 2002. · Zbl 05449956
[71]	H. Vázquez-Leal, L. Hernández-Martínez, and A. Sarmiento-Reyes, “Double-bounded homotopy for analysing nonlinear resistive circuits,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS ’05), pp. 3203-3206, May 2005.
[72]	M. Van Barel, K. D. Ikramov, and A. A. Chesnokov, “A continuation method for solving symmetric Toeplitz systems,” Computational Mathematics and Mathematical Physics, vol. 48, no. 12, pp. 2126-2139, 2008.
[73]	H. Vázquez-Leal, L. Hernández-Martínez, A. Sarmiento-Reyes, and R. Castañeda-Sheissa, “Numerical continuation scheme for tracing the double bounded homotopy for analysing nonlinear circuits,” in Proceedings of the International Conference on Communications, Circuits and Systems, pp. 1122-1126, May 2005.
[74]	L. Trajkovic, R. C. Melville, and S. C. Fang, “Passivity and no-gain properties establish global convergence of a homotopy method for DC operating points,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 914-917, May 1990.
[75]	M. M. Green, “Efficient continuation method for use in globally convergent DC circuit simulation,” in Proceedings of the 3rd URSI International Symposium on Signals, Systems and Electronics (ISSSE ’95), pp. 497-500, October 1995.
[76]	L. B. Goldgeisser and M. M. Green, “A method for automatically finding multiple operating points in nonlinear circuits,” IEEE Transactions on Circuits and Systems I, vol. 52, no. 4, pp. 776-784, 2005. · Zbl 1374.94899
[77]	D. M. Wolf and S. R. Sanders, “Multiparameter homotopy methods for finding dc operating points of nonlinear circuits,” IEEE Transactions on Circuits and Systems I, vol. 43, no. 10, pp. 824-838, 1996.
[78]	L. O. Chua and A. Ushida, “A switching-parameter algorithm for finding multiple solutions of nonlinear resistive circuits,” International Journal of Circuit Theory and Applications, vol. 4, no. 3, pp. 215-239, 1976. · Zbl 0326.65034
[79]	K. S. Chao, D. K. Liu, and C. T. Pan, “A systematic search method for obtaining multiple solutions of simultaneous nonlinear equations,” IEEE Transactions on Circuits and Systems, vol. 22, no. 9, pp. 748-752, 1975.
[80]	F. H. Branin Jr., “Widely convergent method for finding multiple solutions of simultaneous nonlinear equations,” IBM Journal of Research and Development, vol. 15, no. 9, pp. 2779-2780, 1973.
[81]	L. Vandenberghe and J. Vandewalle, “Variable dimension algorithms for solving resistive circuits,” International Journal of Circuit Theory and Applications, vol. 18, no. 5, pp. 443-474, 1990. · Zbl 0705.94021
[82]	M. Kojima and Y. Yamamoto, “Variable dimension algorithms: basic theory, interpretations and extensions of some existing methods,” Mathematical Programming, vol. 24, no. 2, pp. 177-215, 1982. · Zbl 0509.90070
[83]	Y. Yamamoto, “A variable dimension fixed point algorithm and the orientation of simplices,” Mathematical Programming, vol. 30, no. 3, pp. 301-312, 1984. · Zbl 0575.65043
[84]	H. Vázquez-Leal, L. Hernández-Martínez, A. Sarmiento-Reyes, and R. S. Murphy-Arteaga, “Improving multi-parameter homotopy via symbolic analysis techniques for circuit simulation,” in Proceedings of the European Conference on Circuit Theory and Design, vol. 2, pp. 402-405, 2003.
[85]	J. S. Roychowdhury and R. C. Melville, “Homotopy techniques for obtaining a DC solution of large-scale MOS circuits,” in Proceedings of the 33rd Annual Design Automation Conference, pp. 286-291, June 1996.
[86]	J. Roychowdhury and R. Melville, “Delivering global DC convergence for large mixed-signal circuits via homotopy/continuation methods,” IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 25, no. 1, pp. 66-78, 2006. · Zbl 05449547
[87]	L. Trajkovic, R. C. Melville, and S. C. Fang, “Finding DC operating points of transistor circuits using homotopy methods,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 758-761, June 1991.
[88]	L. Trajkovic and A. N. Willson Jr., “Theory of dc operating points of transistor networks,” Archiv fur Elektronik und Ubertragungstechnik, vol. 46, no. 4, pp. 228-241, 1992.
[89]	H. Vazquez-Leal, R. Castaneda-Sheissa, A. Yıldırım et al., “Biparameter homotopybased direct current simulation of multistable circuits,” British Journal of Mathematics and Computer Science, vol. 2, no. 3, pp. 137-150, 2012.
[90]	Y. Inoue, S. Kusanobu, K. Yamamura, and T. Takahashi, “Newton-fixed-point homotopy method for finding dc operating points of nonlinear circuits,” in Proceedings of the International Technical Conference on Circuits/Systems, Computer and Communications (ITC-CSCC ’01), vol. 1, pp. 370-373, Tokushima, Japan, July 2001.
[91]	L. Trajković and W. Mathis, “Parameter embedding methods for finding dc operating points: formulation and implementation,” in Proceedings of the NOLTA, pp. 1159-1164, Dicember 1995.
[92]	L. Trajković, R. C. Melville, and S. C. Fang, “Passivity and no-gain properties establish global convergence of a homotopy method for DC operating points,” in Proceedings of the IEEE International Symposium on Circuits and Systems, pp. 914-917, May 1990.
[93]	H. Vazquez-Leal, L. Hernandez-Martinez, A. Sarmiento-Reyes, R. Castaneda-Sheissa, and A. Gallardo-Del-Angel, “Homotopy method with a formal stop criterion applied to circuit simulation,” IEICE Electronics Express, vol. 8, no. 21, pp. 1808-1815, 2011.