American Journal of Computational and Applied Mathematics

p-ISSN: 2165-8935    e-ISSN: 2165-8943

2011;  1(1): 15-19

doi: 10.5923/j.ajcam.20110101.04

Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials

K. Maleknejad , E. Hashemizadeh

Department of Mathematics, Karaj Branch, Islamic Azad University, Karaj, Iran

Correspondence to: K. Maleknejad , Department of Mathematics, Karaj Branch, Islamic Azad University, Karaj, Iran.

Email:

Copyright © 2012 Scientific & Academic Publishing. All Rights Reserved.

Abstract

This paper proposed a numerical method for nonlinear singular ordinary differential equations, that arises in biology and some diseases. We solved these nonlinear problems by a new method based on shifted Legendre polynomials. Operational matrices of derivatives for this function are presented to reduce the nonlinear singular boundary value problems to a system of nonlinear algebraic equations. The method is computationally very simple and attractive, and applications are demonstrated through illustrative examples. The results obtained are compared by the known results.

Keywords: Nonlinear Singular Boundary Value Problem, Legendre Polynomials, Operational Matrix of Derivative, Collocation Method, Biology

Cite this paper: K. Maleknejad , E. Hashemizadeh , "Numerical Solution of Nonlinear Singular Ordinary Differential Equations Arising in Biology Via Operational Matrix of Shifted Legendre Polynomials", American Journal of Computational and Applied Mathematics , Vol. 1 No. 1, 2011, pp. 15-19. doi: 10.5923/j.ajcam.20110101.04.

1. Introduction

The aim of this paper is to introduce a new method for the numerical solution of the following class of singular boundary value problems
(1)
(2)
(3)
which arising in biology and physiology problems. We assume that is continuous, exists and is continuous and also . The boundary value problem (1)-(3) with and arise in the study of various tumor growth problems, see ([1-6]), with linear and with nonlinear of the form
(4)
A mathematical model of tumor growth is a mathematical expression of the dependence of tumor size on time.
And when in the study of oxygen diffusion problem in a spherical cell with Michaelis-Menten Kinetics, see ([7-9]). A similar problem arise with and in modelling of heat conduction in human head, see[10-13], with of the form
(5)
Existence-uniqueness results for such problems have been established by several researchers[14-16]. In recent years, finding numerical solutions of singular differential equations, particularly those arising in physiology, has been the focus of a number of authors, which you can see some of them in[17-20].
The purpose of this paper is to introduce a novel method based on operational matrices of derivatives of shifted Legandre polynomials that have been introduced recently in Saadatmandi and Dehghan work's[21] for the numerical solution of the class of singular second-order boundary value problems given in the (1-3) that arise in physiology. In this work by use of shifted Legendre polynomials as basis and operational matrices of derivatives of them we convert these kinds of equations to algebraic equations. The advantage of this method analogy to other existed method for these problems is its trusty and simply in implementation, we compared our results with some existed results to prove this claim.
This paper is organized as follows: Section 2 represents preliminaries, in this section we introduced shifted Legendre polynomials, and some properties of them, specially the operational matrices of derivatives, in Section 3 we implemented them on physiology problems. In Section 4, a number of applied models in physiology are discussed to show the efficiency and accuracy of the proposed method, the results obtained are compared by the known results. Finally, Section 5 includes a conclusion for the paper.

2. Definitions and Properties of Shifted Legendre Polynomials

Consider the Legendre polynomials on the interval
the set in Hilbert space is a complete orthogonal set[22,23]. In order to use these polynomials on the interval we define the so-called shifted Legendre polynomials by introducing the change of variable . Let the shifted Legendre polynomials be denoted by . Then can be obtained as follows:
(6)
where and . The analytic form of the shifted Legendre polynomials of degree given by
(7)
Note that and . The orthogonality condition is
Any function can be expanded in terms of shifted Legendre polynomials as
where the coefficients are given by
In practice, only the first -terms shifted Legendre polynomials are considered. Then we have
(9)
where the shifted Legendre coefficient vector and the shifted Legendre vector are given by:
(10)
(11)
The derivative of the vector can be expressed by
(12)
where is the operational matrix of derivative given by[21]
as if is odd and if is even . For example for even we have
By using Eq. (12), it is clear that
(13)
where and the superscript, in , denote matrix powers. Thus
(14)

3. Implementation of Shifted Legendre Polynomials Method on Physiology Problems

In this section we solve nonlinear singular boundary value problem of the form Eq.(1) with the mixed conditions (2) and (3) by using shifted Legendre polynomials.
From Eq. (9) we can approximate our unknown as
(15)
where and are defined in Eqs.(10) and (11). By using Eqs.(12) and (13) we have
(16)
and
(17)
By substituting Eqs.(15), (16) and (17) in Eq. (1) we have
(18)
Also by using Eqs.(2), (3), (15) and (16) we have
(19)
(20)
Eqs.(19) and (20) give two linear equations. Since the total unknowns for vector in Eq.(15) is , we collocate Eq.(18) in points in the interval that are roots of shifted Legendre polynomial , then we have,
(21)
for . Now the resulting Eqs. (19), (20) and (21) generate a system of nonlinear equations which can be solved using Newton's iterative method[24,25]. We used the Mathematica 7 software to solve this nonlinear system.

4. Illustrative Examples and Applied Models

To show the efficiency of the proposed numerical method, we implement it on three nonlinear singular boundary problems that arise in real physiology applications. Our results are compared with result in Refs.[17-20]. The austerity of our method in implementation in analogy to other existed methods and its trusty answers is considerable.

4.1. Example 1

Consider the following oxygen diffusion problem
with the boundary conditions:
Table 1 shows the numerical results for various number of meshes, and present method solutions are compared with results in Refs.[17] and[18].
Table 1. Approximate solutions for Example 1.


Present methodwith

Method in[17]With

Method in[18]with

0.00.828483290359680.828483294813550.82848327295802
0.10.829706092433800.829706096887900.82970607521884
0.20.833374733591000.833374738043080.83337471691089
0.30.839489913953700.839489918339860.83948989814383
0.40.848052784996060.848052788760510.84805277036165
0.50.859064927169230.859064927530320.85906491397434
0.60.872528319958280.872528315698550.87252830841853
0.70.888445305623190.888445299497020.88844529589927
0.80.906818548066800.906818541799650.90681854026297
0.90.927650988365580.927650983052560.92765098252660
1.00.950945798496480.950945794805230.95094579461056
Table 2. Numerical errors for Example 2.


Present methodwith

Present methodwith

Approach II[17]with

0.0





0.1





0.2





0.3





0.4





0.5





0.6





0.7





0.8





0.9





1.0





4.2. Example 2

Consider the following singular two point boundary value problem:
with the exact solution
,
where . Table 2 shows numerical errors of this example in analogy to errors for this example in[17].

4.3. Example 3

Consider this problem that is coincide by heat conduction model of the human head,
we consider the solution of this problem with conditions as follows:
Table 3 illustrates results for this example by proposed method alongside numerical solutions for this example that have been given in Refs[19-20].
Table 3. Approximate solutions for Example 3.


Present methodwith

Method in[19] with forth-orderMethod in[20]
0.00.36751681510.36751810740.3675169710
0.10.36636232920.36636375610.3663623697
0.20.36289406610.36289593780.3628941066
0.30.35709754570.35709914290.3570975842
0.40.34894842060.34894999030.3489484612
0.50.33841214870.33841365810.3384121893
0.60.32544352240.32544500190.3254435631
0.70.30998604020.30998785670.3099860810
0.80.29197110300.29197896540.2919711440
0.90.27131701010.27131856370.2713170512
1.00.24792772330.24792928370.2479277646
Table 4. The maximum absolute errors in solution of Example 4 for
     .


Case (i)

Case (i)

Case (ii)

Case (ii)

15







20







25







30







Table 5. The maximum absolute errors in solution of Example 4 for
     .


Case (i)

Case (i)

Case (ii)

Case (ii)

15







20







25







30







4.4. Example 4

Consider the following singular two point boundary value problem:
for the following two cases:
with the exact solution
Maximum absolute errors for this problem have been displayed for in Table 4 and for in Table 5, which show the accuracy of proposed method and these results in analogy to exhibited results for this example in[19-20] show advantage of this method.

5. Conclusions

This paper present a new approach, based on shifted Legendre polynomials for the numerical solution of a class of singular boundary value problems arising in biology and physiology problems. By use of shifted Legendre polynomials as basis and operational matrices of derivatives of these functions we convert such problems to an algebraic system. The implementation of current approach in analogy to existed methods is more convenient and the accuracy is high and we can execute this method in a computer speedy with minimum CPU time used. The numerical applied models that have been presented in the paper and the compared results support our claim.

ACKNOWLEDGEMENTS

The authors would like to thank Islamic Azad University of Karaj Branch for partially financially supporting this re-search and providing facilities and encouraging this work.

References

[1]  J.A. Adam, A simplified mathematical model of tumor growth, Math. Biosci., vol.81, pp.224–229, 1986
[2]  J.A. Adam, A mathematical model of tumor growth II: effect of geometry and spatial non-uniformity on stability, Math. Biosci. vol.86, pp.183–211, 1987
[3]  J.A. Adam, S.A. Maggelakis, Mathematical model of tumor growth IV: effect of necrotic core, Math. Biosci., vol. 97, pp. 121–136, 1989
[4]  A.C. Burton, Rate of growth of solid tumor as a prob-lem of diffusion, Growth, vol.30 , pp.157–176, 1966
[5]  H.P. Greenspan, Models for the growth of solid tumor as a problem by diffusion, Stud. Appl. Math., vol.52, pp.317–340, 1972
[6]  N.S. Asaithambi, J.B. Goodman, Point wise bounds for a class of singular diffusion problems in physiology, Appl. Math. Comput., vol. 30 , pp.215–222, 1989
[7]  H.S. Lin, Oxygen diffusion in a spherical cell with non-linear oxygen uptake kinetics, J. Theor. Biol., vol. 60 , pp.449–457, 1976
[8]  D.L.S. McElwain, A re-examination of oxygen diffusion in a spherical cell with MichaelisMenten oxygen uptake kinetics, J. Theor. Biol., vol. 71, pp.255–263, 1978
[9]  N. Rashevsky, Mathematical Biophysics, vol. 1, Dover, New York, 1960
[10]  U. Flesch, The distribution of heat sources in the hu-man head: a theoretical consideration, J. Theor. Biol., vol. 54, pp. 285–287, 1975
[11]  J.B. Garner, R. Shivaji, Diffusion problems with mixed nonlinear boundary condition, J. Math. Anal. Appl., vol. 148, pp. 422–430, 1990
[12]  B.F. Gray, The distribution of heat sources in the hu-man head: a theoretical consideration, J. Theor. Biol., vol. 82 pp. 473–476, 1980
[13]  R.C. Duggan, A.M. Goodman, Point wise bounds for nonlinear heat conduction model for the human head, Bull. Math. Biol., vol. 48 (2), pp. 229–236, 1989
[14]  R.K. Pandy, On a class of weakly regular singular two point boundary value problems II, J, Differential Equa-tions, vol. 127, pp. 110-123, 1996
[15]  M.M. Chawla, P.N. Shivkumar, On the existence of so-lution of a class of singular two--point nonlinear boun-dary value problems, J. Comput. Appl. Math., vol. 19, pp. 379–388, 1987
[16]  R.D. Russell, L.F. Shampine, Numerical methods for singular boundary value problems, SIAM J. Numer. Anal., vol. 12, pp. 13-36, 1975
[17]  S.A. Khuri, A. Sayfy, A novel approach for the solution of a class of singular boundary value problems arising in physiology, J. Math. Comput. Model., vol. 52, pp. 626–636, 2010
[18]  Hikmet Caglar, Nazan Caglar, Mehmet Ozer, B-spline solution of non-linear singular boundary value problems arising in physiology, Chaos Solitons Fractals, vol. 39, pp. 1232-1237, 2009
[19]  J. Rashidinia, R. Mohammadi, R. Jalilian, The numerical solution of non-linear singular boundary value problems arising in physiology, J. Appl. Math. Comput., vol. 185, pp. 360–367, 2007
[20]  R.K. Pandey, Arvind K. Singh, On the convergence of a finite difference method for a class of singular boundary value problems arising in physiology, J. Comput. Appl. Math., vol. 166, pp. 553–564, 2004
[21]  A. Saadatmandi, M. Dehghan, A new operational ma-trix for solving fractional-order differential equations, J. Comput. Math. Appl., vol. 59, pp. 1326–1336, 2010
[22]  K. Maleknejad, B. Basirat, E. Hashemizadeh, Hybrid Le-gendre polynomials and Block-Pulse functions approach for nonlinear Volterra–Fredholm integro-differential eq-uations, Comput. Math. Appl.,vol. 61, pp. 2821–2828, 2011
[23]  K. Maleknejad, E. Hashemizadeh , A numerical approach for Hammerstein integral equations of mixed type using operational matrices of hybrid functions, Scientific Bulle-tin, Series A: Applied Mathematics and Physics,‎ vol. 73(3), pp. 95-104, 2011
[24]  K. Maleknejad, S. Sohrabi, H. Derili, A new computa-tional method for solution of non-linear Volter-ra-Fredholm integro-differential equations, Int. J. Appl. Math., vol. 87(2), pp. 327–338, 2010
[25]  K. Maleknejad, E. Hashemizadeh, B. Basirat, Computational method based on Bernestein operational matrices for nonlinear Volterra-Fredholm-Hammerstein integral equations, Commun. Nonlinear. Sci. Numer. Simulat., vol. 17, pp.52–61, 2012