Discrete Dynamics in Nature and Society
Volume 2018 (2018), Article ID 9142124, 9 pages
https://doi.org/10.1155/2018/9142124

Research Article

Optimized Direct Padé and HPM for Solving Equation of Oxygen Diffusion in a Spherical Cell

M. A. Sandoval-Hernandez,¹ H. Vazquez-Leal,² A. Sarmiento-Reyes,¹ U. Filobello-Nino,² F. Castro-Gonzalez,¹ A. L. Herrera-May,3,4 R. A. Callejas-Molina,⁵ and R. Ruiz-Gomez¹

¹National Institute for Astrophysics, Optics and Electronics, Luis Enrique Erro No. 1, Sta. María Tonantzintla, 72840, Puebla, Mexico
²Facultad de Instrumentación Electrónica, Universidad Veracruzana, Cto. Gonzalo Aguirre Beltrán S/N, 91000, Xalapa, Veracruz, Mexico
³Centro de Investigación de Micro y Nanotecnología, Universidad Veracruzana, Calzada Ruíz Cortines 455, Boca del Río, 94292, Veracruz, Mexico
⁴Maestría en Ingeniería Aplicada, Facultad de Ingeniería de la Construcción y el Hábitat, Universidad Veracruzana, Calzada Ruíz Cortines 455, 94294, Boca del Río, Veracruz, Mexico
⁵Instituto Tecnológico de Celaya, Antonio García Cubas Pte No. 600 Esq. Av. Tecnológico, Celaya, 38010, Guanajuato, Mexico

Correspondence should be addressed to M. A. Sandoval-Hernandez; m.sandoval@inaoep.mx

Received 25 January 2018; Accepted 11 July 2018; Published 2 September 2018

Academic Editor: Guang Zhang

Copyright © 2018 M. A. Sandoval-Hernandez 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 work presents the application of homotopy perturbation method (HPM) and Optimized Direct Padé (ODP) to obtain a handy and easily computable approximate solution of the nonlinear differential equation to model the oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics. On one hand, the obtained HPM solution is fully symbolic in terms of the coefficients of the equation, allowing us to use the same solution for different values of the maximum reaction rate, the Michaelis constant, and the permeability of the cell membrane. On the other hand, the numerical experiments show the high accuracy of the proposed ODP solution, resulting in as the lowest absolute relative error (A.R.E.) for a set of coefficients. In addition, a novel technique is proposed to reduce the number of algebraic operations during the process of application of ODP method through the use of the Taylor series, which help to simplify the algebraic expressions used. The powerful process to obtain the solution shows that the Optimized Direct Padé and homotopy perturbation method are suitable methods to use.

1. Introduction

Michaelis-Menten kinetics describes the rate of enzymatic reactions. This model is valid when the concentration of a substrate is higher than the concentration of the enzymes, and for steady state conditions, that is, when the concentration of the complex enzyme-substrate is constant [1]. There are several works about the Michaelis-Menten oxygen uptake kinetics. For example, in [2], the relation between Michaelis-Menten direct and inverse kinetics, chemical kinetics of approximation, and second-order kinetics is presented. Moreover, in [3], the subject of transport phenomena is approached by means of the study of transport of quantity of motion (viscous flow), transport of energy (heat conduction, convection, and radiation), and matter transport (Diffusion) [4, 5]. Consequently, transport phenomena are employed for solving different problems in the area of sciences such as chemistry, biochemistry, soil physics, meteorology, biology, and semiconductors, disciplines in which the use of Bessel functions, differential equations, and Laplace transform is required [6–10]. In [11], the behaviour of dopamine released from a small iontophoresis electrode and its voltammetric detection by a carbon fiber sensor 100pm away is presented as a basis for developing a new paradigm for measuring dopamine kinetics in intact rat neostriatum. In [11], the presented model was derived from the work of diffusion in a punctual iontophoretic source where the lineal term of absorption is replaced by a nonlinear expression that describes a Michaelis-Menten kinetic given by a constant , the Michaelis-Menten constant .

In [12], the equation that models oxygen diffusion in a spherical cell, including nonlinear absorption kinetics, is solved by transforming the Lane-Emden equation into its equivalent Volterra integral form, then Adomian decomposition was employed to solve the nonlinear Fredholm-Volterra integral. Furthermore, in [13] a fully symbolic solution was proposed to solve the aforementioned equation by means of a modified Taylor series obtaining an accuracy of as the lowest mean square error without the use of complicated integrals.

There are several methods of solutions for solving nonlinear differential equations such as inverse scattering transformation [14, 15], Darboux transformation [16, 17], bilinear method [18, 19], the tanh-function method [20–22], the variable separation approach [23–25], the symmetry method [26, 27], sine-cosine method [28–30], Adomian decomposition method (ADM) [12, 31–33], and homotopy perturbation method (HPM) [34–36]. HPM is based on the use of power series of , which transforms a differential equation into a set of linear differential equations.

In [37] was presented a procedure to apply Padé method to find approximate solutions for nonlinear differential equations, which consist in that the solution of a differential equation can be directly expressed as a rational power series of the independent variable as a Padé approximant. From (6) ODP employs a polynomial-like rational expression as the proposal of approximation of the nonlinear differential equation to be solved. In general terms, it works by means of substituting the rational expression in the differential equation and then regroups the powers of the independent variable. It is important to note that due to the rational expression a large amount of algebraic operations is generated. However, this work proposes a Direct Padé (DP) modification oriented to reduce this algebraic operation by means of the Taylor expansion of the rational function.

This paper is organized as follows. In Section 2, a brief description of Michaelis-Menten kinetics and its equation is given. In Section 3, a brief description of Lane-Emden equation is presented. The equation to be solved and its boundary conditions are introduced in Section 4. Section 5 describes in detail Optimized Direct Padé. In Section 6, the HPM solution is presented. Section 7 presents the analysis of solution for the case study by HPM and ODP. In Section 8, its numerical simulations, comparisons, and discussion are presented. Lastly, in Section 9 a brief conclusion of this work is presented.

2. Michaelis-Menten Equation

General principles of kinetics in chemical reactions are applicable to catalysed reactions by enzymes in living things. However, this shows a characteristic side which is not observed in nonenzymatic catalyst, the substrate saturation, in terms of enzyme molecules active sites occupation [1]. The study of the effect that substrate concentration has over enzyme activity is not simple task; logically thinking, substrate concentration lowers as the reaction increases. A simplification on kinetics experiments consists in measuring initial velocity . If time is short enough, substrate lowering will be minimal and this could be considered, thus, constant. This behaviour is characteristic of most enzymes and was studied by Michaelis and Menten in 1913 [38]. Figure 1 presents the three phases of enzyme kinetics:

(i)For a low substrate concentration, velocity reaction is directly proportional to substrate concentration (linear relation), first-order kinetics.(ii)For a high substrate concentration, velocity reaction is practically constant and independent of substrate concentration; kinetics is considered zero order.(iii)For medium concentrations of substrate, velocity of the process becomes nonlinear and this phase is called mixed kinetics.

Figure 1: Enzyme kinetics: Michaelis-Menten model.

The curve which expresses the relationship between substrate concentration and initial velocity has the same shape for most enzymes, a rectangular hyperbola, whose algebraic expression is given by Michaelis-Menten equation (1), which describes the relationship between initial velocity and substrate concentration, [1]:

, maximum velocity, and , Michaelis constant, are two characteristic kinetics parameters of every enzyme that can be experimentally determined.

Maximum velocity is obtained when velocity of reaction becomes independent of substrate concentration; when this occurs, velocity reaches a maximum value. This value depends on the available amount of enzyme. Michaelis constant indicates substrate concentration in which reaction rate is half of maximum velocity. This parameter is independent of enzyme concentration and it is characteristic of every enzyme according to the employed substrate (or substrates). Constant indicates the affinity that the enzyme has for the substrate, being this higher at a lower value of . For a lower value , the lower the substrate needed to reach half of maximum velocity, thus, the higher the affinity an enzyme has for the substrate.

3. Lane-Emden Equation

Lane-Emden equation (2) is a basic nonlinear equation arising in the study of stellar structure theory and has

where can be chosen to describe a Cartesian geometry, ; cylindrical, ; or spherical, [39]. is a real function and is an analytical function. Initial conditions are

where and are constants.

Equation (2) is used to model physical and astrophysical phenomena such as theory of stellar structure, spherical cloud of gas thermal behaviour, isothermal gas sphere [40], and transport phenomena [5] as it occurs for the case of particle diffusion in biophysical processes [4] as well in physiological processes.

A model of oxygen diffusion and nonlinear absorption in a sphere was originally proposed and resolved in [41]. Later, the same model was proposed, analysed, and solved in [42]. The equation of the model is a steady state diffusion-reaction equation that represents oxygen transport by means of linear diffusion in a sphere.

4. Equation of Oxygen Diffusion in a Spherical Cell

The differential equation that models the oxygen diffusion in a spherical cell with Michaelis-Menten oxygen uptake kinetics (4) [41–46] is obtained by substituting (1) into (2) and given by

with boundary conditions

Variables and represent the oxygen concentration and the radial distance, respectively. Parameters , , and are constants that represent the maximum reaction rate. Michaelis constant is the concentration at half saturation and cell membrane permeability, respectively.

When boundary condition is equal to , oxygen distribution is symmetrical to the center of the sphere. Boundary condition at determines that oxygen flow through the cell membrane is proportional to differential, which is less than the cell membrane normalized concentration [44, 45].

5. Optimized Direct Padé

The method ODP [37] considers that a nonlinear differential equation can be expressed as follows:

with boundary conditions

where and are the linear and nonlinear operators, respectively; is a known analytical function; is a boundary operator; is the boundary of the domain [47].

is the initial condition for (6) that satisfies the boundary condition and is the homotopy perturbation parameter. The solution of (6) can be written as

where and are unknowns determined by the method; are the numerator and denominator order.

There are no systematic criteria to choose the optimal Padé order for a given problem. However, generally a finite number of terms are required to obtain highly precise Padé approximation. The basic procedure for Direct Padé can be described as follows:

(1)Boundary conditions are substituted into (8) to generate an equation for each boundary condition. It is important to note that there exists an algebraic equation for each boundary condition.(2) from (8) is substituted into (6), then regrouping the resulting equation in terms of powers. Afterwards, the regrouping procedure includes eliminating denominator terms arising from Padé approximant (8). In this way, the resulting expression is a power series that represent the residual error from (6).(3)In order to reduce the residual error, from the lowest order, each coefficient from powers in the resulting power series is equalized to zero to obtain an algebraic equation in terms of the unknown coefficients of (8).(4)All the previous steps generate NAES in terms of the unknowns from (8). Finally, NAES are solved to obtain and .

The disadvantage of the method is that it becomes as difficult as the number of terms that the expression has in the numerator as well in the denominator because it is a rational expression. Even more, the substitution into (8) gets more complicated if higher order derivatives are encountered, mainly because the first derivative of a rational expression lies in terms of another rational expression with a squared denominator and if in the nonlinear equation exist several different powers, another rational expression with even more terms is generated. For this reason, a common denominator has to be found, to formulate the system of equations that allows finding and . All this process may be too cumbersome and computationally expensive. In this matter, it is concluded that if the nonlinear equation possesses higher order derivatives of (the independent variable) the number of algebraic operations is notoriously increased when powers are regrouped to form the equation system to find and . Thus, in order to avoid the algebraic complexity that implies the substitution of (8) into (6), in this paper we propose a novel strategy to circumvent such issue by obtaining the Taylor expansion of (8).

where is the Taylor expansion of (8) with respect to and considering , and is the resulting expansion of (8). The replacement of (8) by (9) decreases enormously the number of algebraic steps when it is substituted into the differential equation obtaining a simpler algebraic procedure. Once (9) is replaced into (6), with regrouping in terms of powers, a system of equations is constructed with the aim of finding () and

6. HPM Solution

Multiplying Michaelis-Menten kinetic equation (4) by , the following is obtained:

The proposed homotopy perturbation formulation is

Assuming the solution for (11) is given by a power series of

substituting (12) into (11), and regrouping in terms of ,

Substituting the first two system solutions (13) into (12), we obtain

This is the solution employing HPM, where and are constants to be determined. By substituting (14) into the second boundary condition of (5), i.e., , we obtain

Solving for in (15),

Substituting in (14),

With the end of reducing the residual error, (17) is substituted into (4), at middle point , between boundaries in order to cancel the error at that point, to obtain for

This strategy seeks to reduce the approximation error. Substituting (18) into (17), we obtain the solution in function of , , and .

where and are given for

7. ODP Solution Procedure

The following rational solution is now proposed for (4):

In order to reduce the amount of algebraic steps, it is proposed in this work, and for the first time, to expand the rational expression (21) and employ the Taylor series of order four:

differentiating with respect to and to satisfy the first boundary with yields

With expansion (22) we have the advantage of eliminating the denominator and thus reducing the amount of algebraic terms when substituting (21) into (4). To satisfy the second boundary condition using (21) and (5)

Substituting fourth-order Taylor expansion of (21) into (10),

Substituting the values of , , and solving the simultaneous system of equations given by (23), (24), and (25), the numerical values of , , , , and are obtained to satisfy (21).

8. Numerical Simulations and Discussion

For all the simulations presented in this work, with the purpose of comparison, Maple numerical routines were employed. For example, numerical values of , , were used, with (19) and (21) expressed with the following values:

In order to verify the precision of the solutions (26) and ODP (27) obtained by HPM, their performance is compared with ADM solutions ( and ) reported in [12], as well against the numerical solution given by Maple numerical routines.

and (13) given by MTSM [13] is

Equations (19) and (21) are second-order compact expressions that present a good accuracy at the interval . Figure 2(a) shows, in a graphical manner, that the obtained solutions with ODP and HPM present a better accuracy than the obtained for (28) and (30), (31), for values of , ,. Additionally, Figures 2(b) and 2(c) show the same comparisons with the values for , , and , , .

(a) Case for = 2.5, ,

(b) Case for = 3.5, ,

(c) Case for = 4.5, ,

(d) Case for = 8.5, ,

Figure 2: Comparison of HPM (19) and ODP (21) versus numerical solution, (28), (29), and (30).

As for the previous cases, in Figure 2(d), for higher values of , the obtained solution employing ODP has better accuracy than other solutions. For this case, (32) and (33) describe the algebraic solutions for this value for HPM and ODP, respectively. For the values mentioned before, Figure 2(d) presents the graphics for the solutions employing ODP and HPM compared to the solutions proposed by MTSM and .

In Table 1 is presented and compared the absolute relative error (A.R.E.) for each approximation according to the values of , , , for HPM and ODP solutions, including (28), (29), and (30).

Table 1: Absolute relative error comparison.

Case
ODP
HPM
MTSM [13]
[12]
[12]

, ,

, ,

, ,

, ,

The absolute relative error is defined as the absolute difference of the areas under the curve for the numerical and approximate solutions, divided by the area of the exact solution. The employed integrating method was the Rule of Simpson .

In all of the presented cases, the A.R.E. obtained with ODP is much lower compared to the Adomian and MTSM methods reported in [12, 13] and matches with all the graphics presented in this work.

In the same manner, for the cases () and () for HPM it is observed that the relative absolute error is lower than the ones presented in Adomian and MTSM. Besides, the solution obtained by HPM has the advantage of being simpler because it only consists of two terms, compared against (28), (29), and (30) that possess more terms.

Thus, for Optimized Direct Padé and HPM, the proposed results in this work coincide with [12, 13, 39, 48, 49]; that is, the oxygen concentration is higher at the substrate and diminishes as it approaches the center of spherical cell, which is due to the oxygen diffusing into the cell; it reacts to give rise to the boundary condition at the center. Likewise, results are congruent with the spherical geometry proposed by Lane-Emden (2).

Finally, it is important to remark that although the precision for ODP tends to be the same as DP, it has the advantage of having a simpler algebraic procedure, avoiding the algebraic complexity of substituting a rational expression that generates more rational expressions with higher powers in the denominator o in its case, more numerator terms, making the formulation of the system of equations difficult for finding and from the rational expression of the solution.

9. Conclusions

In this article were proposed two compact solutions for Michaelis-Menten oxygen uptake kinetics equation employing HPM and ODP. ODP presents a quotient of two polynomials, one of them of second order divided by a first-order one. On the other hand, HPM presents only two algebraic terms, the independent term and the second degree term. Likewise, the solutions found in this work were compared to those reported by other authors, showing that ODP has a better accuracy. HPM solution shows a good approximation in comparison to (29) obtained by Adomian; despite the fact that it only has two terms, it concluded that the solution is in good agreement. Furthermore, in this paper, we propose to replace the rational expression (8) by its Taylor series in order to reduce the amount of algebraic steps during the application of DP. The graphs presented in this work show that the solution given by ODP and HPM tends to converge as we approach unlike Adomian which tends to diverge by increasing the absolute error. For the case of , , and with ODP, the A.R.E. was . Also the HPM and ODP methods were compared against the numerical one obtaining very good approximations. In addition, ODP requires fewer iterations to obtain a more compact and accurate expression than the solution proposed by the Adomian method (29).

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

All the authors express their appreciation to the late Professor Jose Antonio Agustin Perez-Sesma whose contribution to this work was of great significance. Additionally, the authors wish to thank M.TN. César Colocia García and Eng. Francisco Martínez Barrios for their support. The first author is currently a postdoctoral research fellow at Institute for Astrophysics, Optics and Electronics and expresses his gratitude for the support received. He also acknowledges the financial support from the Secretary of Public Education of México (SEP) through grant 2703 E476300.0275652.

References

K. Mathews Christopher, K. E. Van Holde, and G. Ahern Kevin, Bioquimica, Pearson Addison Wesley, 2002.

J. Y. Tang, “On the relationships between the Michaelis-Menten kinetics, reverse Michaelis-Menten kinetics, equilibrium chemistry approximation kinetics, and quadratic kinetics,” Geoscientific Model Development, vol. 8, no. 12, pp. 3823–3835, 2015.

V. Galvanauskas, V. Grincas, and R. Simutis, “Modeling of Oxygen Supply Profiles in Stirred-Tank Aggregated Stem Cells Cultivation Process,” International Journal of Biological, Biomolecular, Agricultural, Food and Biotechnological Engineering, vol. 10, no. 9, pp. 584–587, 2016.

S. I. Rubinow, Introduction to Mathematical Biology, John Wiley and Sons, 1975.

R. Byron Bird, E. S. Warren, and N. L. Edwin, “Fenómenos de transporte: un estudio sistemático de los fundamentos del transporte de materia, energía y cantidad de movimiento,” Limusa Noriega Editores, 2006.

L.-M. Warren, C.-S. Julian, and P. Harriot, Operaciones Unitarias En Ingenieria, McGraw Hill Interamericana, 2007.

G. N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press, 1995.

M. S. Simon and K. K. Ng, Physics of Semiconductor Devices, John Wiley & Sons, 2006.

L. Frenzel, Principles of Electronic Communication Systems, McGraw-Hill, Inc, 2007.

J. J. D’azzo and C. H. Houpis, Feedback Control Systems: Analysis and Synthesis, McGraw-Hill, 1960.

C. Nicholson, “Interaction between diffusion and Michaelis-Menten uptake of dopamine after iontophoresis in striatum,” Biophysical Journal, vol. 68, no. 5, pp. 1699–1715, 1995.

R. Rach, A.-M. Wazwaz, and J.-S. Duan, “A reliable analysis of oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics,” International Journal of Biomathematics, vol. 7, no. 2, p. 12, 2014.

H. Vazquez-Leal, M. Sandoval-Hernandez, R. Castaneda-Sheissa, U. Filobello-Nino, and A. Sarmiento-Reyes, “Modified Taylor solution of equation of oxygen diffusion in a spherical cell with Michaelis-Menten uptake kinetics,” International Journal of Applied Mathematical Research, vol. 4, no. 2, pp. 253–258, 2015.

M. J. Ablowitz, X.-D. Luo, and Z. H. Musslimani, “Inverse scattering transform for the nonlocal nonlinear Schrödinger equation with nonzero boundary conditions,” Journal of Mathematical Physics, vol. 59, no. 1, 2018.

R. Beals and R. R. Coifman, “Scattering and inverse scattering for first order systems,” Communications on Pure and Applied Mathematics, vol. 37, no. 1, pp. 39–90, 1984.

G. Chaohao, H. Anning, and Z. Zixiang, Darboux Transformations in Integrable Systems: Theory and their Applications to Geometry, vol. 26, Springer Science & Business Media, 2006.

P. Adamopoulou, A. Doikou, and G. Papamikos, “Darboux–Bäcklund transformations, dressing & impurities in multi-component NLS,” Nuclear Physics B, vol. 918, pp. 91–114, 2017.

P. P. Goldstein, “Hints on the Hirota bilinear method,” Acta Physica Polonica A, vol. 112, no. 6, pp. 1171–1184, 2007.

Y. Matsuno, Bilinear Transformation Method, Elsevier, 1984.

W. Malfliet, “The tanh method: a tool for solving certain classes of nonlinear evolution and wave equations,” Journal of Computational and Applied Mathematics, vol. 164-165, pp. 529–541, 2004.

N. H. Abdel-All, M. A. Abdel-Razek, and A. K. Seddeek, “Expanding the Tanh-Function Method for Solving Nonlinear Equations,” Applied Mathematics, vol. 2, no. 9, pp. 1096–1104, 2011.

W. Malfliet and W. Hereman, “The tanh method: I. Exact solutions of nonlinear evolution and wave equations,” Physica Scripta, vol. 54, no. 6, pp. 563–568, 1996.

A. V. Kuznetsov and A. A. Avramenko, “The method of separation of variables for solving equations describing molecular-motor-assisted transport of intracellular particles in a dendrite or axon,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 464, no. 2099, pp. 2867–2886, 2008.

X.-M. Qian, S. Y. Lou, and X.-B. Hu, “Variable separation approach for a differential-difference system: Special Toda equation,” Journal of Physics A: Mathematical and General, vol. 37, no. 6, pp. 2401–2411, 2004.

H.-C. Hu, S.-Y. Lou, and Q.-P. Liu, “Darboux Transformation and Variable Separation Approach: The Nizhnik-Novikov-Veselov Equation,” Chinese Physics Letters, vol. 20, no. 9, pp. 1413–1415, 2003.

P. E. Hydon, “Introduction to symmetry methods in the solution of differential equations that occur in chemistry and chemical biology,” International Journal of Quantum Chemistry, vol. 106, no. 1, pp. 266–277, 2006.

P. E. Hydon and P. E. Hydon, Symmetry Methods for Differential Equations: A Beginner’s Guide, vol. 22, Cambridge University Press, 2000.

A. M. Wazwaz, “A sine-cosine method for handling nonlinear wave equations,” Mathematical and Computer Modelling, vol. 40, no. 5-6, pp. 499–508, 2004.

A. J. M. Jawad, “The sine–cosine function method for the exact solutions of nonlinear partial differential equations,” IOSR Journal of Mathematics, vol. 13, no. 1, pp. 186–191, 2012.

A.-M. Wazwaz, “The tanh method and the sine-cosine method for solving the KP-MEW equation,” International Journal of Computer Mathematics, vol. 82, no. 2, pp. 235–246, 2005.

A.-M. Wazwaz, “A new algorithm for calculating Adomian polynomials for nonlinear operators,” Applied Mathematics and Computation, vol. 111, no. 1, pp. 33–51, 2000.

G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.

M. Al-Mazmumy, A. Al-Mutairi, and K. Al-Zahrani, “An Efficient Decomposition Method for Solving Bratu’s Boundary Value Problem,” American Journal of Computational Mathematics, vol. 07, no. 01, pp. 84–93, 2017.

H. Vazquez-Leal, R. Castaneda-Sheissa, U. Filobello-Nino, A. Sarmiento-Reyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution integral,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012.

H. Jafari, H. Hosseinzadeh, M. R. Gholami, and D. D. Ganji, “Application of homotopy perturbation method for heat and mass transfer in the two-dimensional unsteady flow between parallel plates,” International Journal of Applied and Computational Mathematics, vol. 3, no. 3, pp. 1677–1688, 2017.

M. Rabbani and R. Arab, “Extension of some theorems to find solution of nonlinear integral equation and homotopy perturbation method to solve it,” Mathematical Sciences, vol. 11, no. 2, pp. 87–94, 2017.

H. Vazquez-Leal, B. Benhammouda, U. Filobello-Nino et al., “Direct application of Padé approximant for solving nonlinear differential equations,” SpringerPlus, vol. 3, no. 1, p. 11, 2014.

A. Cornish-Bowden, “One hundred years of Michaelis–Menten kinetics,” Perspectives in Science, vol. 4, pp. 3–9, 2015.

M. J. Simpson and A. J. Ellery, “An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell,” Applied Mathematical Modelling, vol. 36, pp. 3329–3334, 2012.

S. Chandrasekhar, Introduction to the Study of Stellar Structure, Dover, 1967.

S. H. Lin, “Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics,” Journal of Theoretical Biology, vol. 60, no. 2, pp. 449–457, 1976.

D. L. S. McElwain, “A re-examination of oxygen diffusion in a spherical cell with Michaelis-Menten oxygen uptake kinetics,” Journal of Theoretical Biology, vol. 71, no. 2, pp. 255–263, 1978.

P. Hiltmann and P. Lory, “On oxygen diffusion in a spherical cell with Michaelis-Menten oxygen uptake kinetics,” Bulletin of Mathematical Biology, vol. 45, no. 5, pp. 661–664, 1983.

M. J. Simpson and A. J. Ellery, “An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell,” Applied Mathematical Modelling, vol. 36, no. 7, pp. 3329–3334, 2012.

A.-M. Wazwaz, “The variational iteration method for solving nonlinear singular boundary value problems arising in various physical models,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 10, pp. 3881–3886, 2011.

R. Rach, A.-M. Wazwaz, and J.-S. Duan, “A reliable analysis of oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics,” International Journal of Biomathematics, vol. 7, no. 2, pp. 1–12, 2014.

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, no. 3, pp. 197–201, 2012.

C. Nicholson and J. M. Phillips, “Ion diffusion modified by tortuosity and volume fraction in the extracellular microenvironment of the rat cerebellum,” The Journal of Physiology, vol. 321, no. 1, pp. 225–257, 1981.

J. W. Schofield, E. A. Gaffney, R. A. Gatenby, and P. K. Maini, “Tumour angiogenesis: The gap between theory and experiments,” Journal of Theoretical Biology, vol. 274, no. 1, pp. 97–102, 2011.