﻿ An Integral Collocation Approach Based on Legendre Polynomials for Solving Riccati, Logistic and Delay Differential Equations

Applied Mathematics
Vol.5 No.15(2014), Article ID:48886,10 pages DOI:10.4236/am.2014.515228

An Integral Collocation Approach Based on Legendre Polynomials for Solving Riccati, Logistic and Delay Differential Equations

M. M. Khader1,2, A. M. S. Mahdy3, M. M. Shehata3

1Department of Mathematics and Statistics, College of Science, Al-Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh, KSA

2Department of Mathematics, Faculty of Science, Benha University, Benha, Egypt

3Department of Mathematics, Faculty of Science, Zagazig University, Zagazig, Egypt

Email: mohamedmbd@yahoo.com, amr_mahdy85@yahoo.com, dr.maha_32@hotmail.com

Copyright © 2014 by authors and Scientific Research Publishing Inc.

Received 28 May 2014; revised 4 July 2014; accepted 16 July 2014

ABSTRACT

In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations with variable coefficients. The properties of the Legendre polynomials are used to reduce the proposed problems to the solution of non-linear system of algebraic equations using Newton iteration method. We give numerical results to satisfy the accuracy and the applicability of the proposed schemes.

Keywords:Integral Collocation Formulation, Spectral Method, Riccati, Logistic and Delay Differential Equations

1. Introduction

It is well known that the ordinary differential equations (ODEs) have been the focus of many studies due to their frequent appearance in various applications, such as in fluid mechanics, viscoelasticity, biology, physics and engineering applications, for more details, for example [1] -[5] . Consequently, considerable attention has been given to the efficient numerical solutions of ODEs of physical interest, because it is difficult to find exact solutions. Different numerical methods have been proposed in the literature for solving ODEs [6] -[13] .

The Riccati differential equation (RDE) is named after the Italian Nobleman Count Jacopo Francesco Riccati (1676-1754). The book of Reid [14] contains the fundamental theories of Riccati equation, with applications to random processes, optimal control, and diffusion problems. Besides important engineering science applications that today are considered classical, such as stochastic realization theory, optimal control, robust stabilization, and network synthesis, the newer applications include such areas as financial mathematics [15] . The solution of this equation can be reached using classical numerical methods such as, the forward Euler method and RungeKutta method. Bahnasawi et al. [16] presented the usage of Adomian decomposition method to solve the nonlinear RDE in an analytic form. Tan and Abbasbandy [17] employed the analytic technique called homotopy analysis method to solve the quadratic RDE.

The Logistic model can be obtained by applying the derivative operator on the Logistic equation. The model is initially published by Pierre Verhulst in 1838 [18] . The continuous Logistic model is described by first order ODE. The discrete Logistic model is simple iterative equation that reveals the chaotic property in certain regions [19] . There are many variations of the population modeling [20] . The Verhulst model is the classic example to illustrate the periodic doubling and chaotic behavior in dynamical system [19] . The model describs the population growth may be limited by certain factors like population density [20] .

2. Applications of Logistic Equation

A typical application of the Logistic equation is a common model of population growth. Let represents the population size and represents the time where the constant defines the growth rate.

Another application of Logistic curve is in medicine, where the Logistic differential equation is used to model the growth of tumors. This application can be considered an extension of the above mentioned use in the framework of ecology.

The solution of Logistic equation is explained the constant population growth rate which not includes the limitation on food supply or spread of diseases [21] . The solution curve of the model is increase exponentially from the multiplication factor up to saturation limit which is maximum carrying capacity [21] ,

where is the population with respect to time, is the rate of maximum population growth and is the carrying capacity. The solution of continuous Logistic equation is in the form of constant growth rate as in formula where is the initial population [22] .

A delay differential equation (DDE) is a differential equation in which the derivative of the function at any time depends on the solution at previous time. Introduction of delay in the model enriches its dynamics and allows a precise description of the real life phenomena. DDEs are proved useful in control systems [23] , lasers, traffic models [24] , metal cutting, epidemiology, neuroscience, population dynamics [25] , chemical kinetics [26] , etc. In DDE, one has to provide history of the system over the delay interval as the initial condition. Due to this reason delay systems are infinite dimensional in nature. Because of in infinite dimensionality the DDEs are difficult to analyze analytically and hence the numerical solutions play an important role.

In [27] , Mai-Duy, et al. derived an integral collocation approach based on Chebyshev polynomials and used it for solving numerically the bi-harmonic equations. In [28] , Bhrawy and Alofi introduced a new shifted Chebyshev operational matrix of fractional integration of arbitrary order and applied together with spectral tau method for solving linear fractional differential equations. Khader et al. [29] introduced a new approximate formula of the fractional derivative using Legendre series expansion and used it to solve numerically the fractional diffusion equation.

In this article, we extend the previous work and derive some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations.

Our paper is organized as follows: In Section 2, we derive some integration collocation formulations using Legendre series expansion. In Section 3, we give the integral collocation approach for solving Riccati, Logistic and delay differential equations. In Section 4, the paper ends with a brief conclusion and some remarks.

3. Integration Collocation Formulations

The well known Legendre polynomials are defined on the interval and can be determined with the aid of the following recurrence formula [30]

where and. 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

where and. The analytic form of the shifted Legendre polynomials of degree is given by

(1)

Note that and. The orthogonality condition is

In this work, to consider the differential equation of -th order, we build the integration collocation method using the truncated Legendre series of degree to represent the -th derivative of the unknown function in the following manner

(2)

Using the integration we can obtain the lower-order derivatives and the function itself as follows

(3)

(4)

(5)

(6)

from (1) and (2) we have

(7)

We now collocate Equatuions (2)-(6) at points, as

(8)

where, and are integrated matrices.

4. Integral Collocation Approach for Solving Riccati, Logistic and Delay Differential Equations

In this section, we introduce the integral collocation approach using Legendre expansion for solving the Riccati, Logistic and delay differential equations.

4.1. Model 1: Riccati Differential Equation

(9)

we also assume an initial condition

(10)

The exact solution to this problem at is

The procedure of the implementation is given by the following steps:

1) Approximate the function using Formula (6) and its relevant derivatives with, as follows

(11)

where is defined in (7) as

Then the Riccati differential Equation (9) is transformed to the following approximated form

(12)

We now collocate Equation (12) at points, as

(13)

For suitable collocation points we use the roots of shifted Legendre polynomial.

2) Also, by substituting from the initial condition (10) in (11) we can obtain an equation which gives the value of the constant as follows

(14)

Equations (13) and (14) represent a system of non-linear algebraic equations which contains seven equations for the unknowns and.

3) Solve the resulting system using the Newton iteration method to obtain the unknowns as follows

Therefore, from Formula (12) we can obtain the approximate solution in the form

The numerical results of the proposed problem (9) are given in Figure 1 with different values of in the interval at. From this figure, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.

4.2. Model 2: Logistic Differential Equation

(15)

We also assume an initial condition

(16)

The exact solution to this problem is given by

Figure 1. A comparison between the exact solution and the approximate solution at and.

The existence and the uniqueness of the proposed problem (15) are introduced in details in [31] .

The procedure of the implementation is given by the following steps:

1) Approximate the function and its relevant derivatives with, using Formula (11). Then the Logistic differential Equation (15) is transformed to the following approximated form

(17)

We now collocate Equation (17) at points, as

(18)

For suitable collocation points we use roots of shifted Legendre polynomial which are

(19)

2) Also, by substituting from the initial condition (16) in (11) with we can obtain an equation which gives the value of the constant.

Equation (18) represents a system of non-linear algebraic equations which contains six equations for the unknowns.

3) Solve the resulting system using the Newton iteration method to obtain the unknowns as follows

(20)

Therefore, from Formula (11) we can obtain the approximate solution in the form

The numerical results of the proposed problem (15) are given in Figure 2 with different values of in the interval. From this figure, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.

4.3. Model 3: Delay Differential Equation

Consider the linear delay differential equation of third-order

(21)

with the initial conditions

(22)

The exact solution of this model is.

The procedure of the implementation is given by the following steps:

1) Approximate the function using Formula (6) and its relevant derivatives with, as follows

Figure 2. A comparison between the exact solution and the approximate solution at m = 3 and m = 5.

(23)

where, and are defined in (7) as follows

Then the delay differential Equation (21) is transformed to the following approximated form

(24)

We now collocate Equation (24) at points, as

(25)

For suitable collocation points we use the roots of shifted Legendre polynomial.

2) Also, by substituting from the initial conditions (22) in (23) we can obtain equations which give the values of the constants

(26)

Equations (25) and (26) represent a system of linear algebraic equations which contains nine equations for the unknowns and.

3) Solve the resulting system using the conjugate gradient method to obtain the unknowns as follows

(27)

Therefore, from Formula (23) we can obtain the approximate solution in the form

The numerical results of the proposed problem (21) are given in Figure 3 with different values of in the interval. From this figure, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.

5. Conclusion and Remarks

In this article, an integral collocation approach based on Legendre polynomials is introduced for solving

Figure 3. A comparison between the exact solution and the approximate solution at m = 5 and m = 7.

numerically the Riccati, Logistic and delay differential equations. The properties of the Legendre polynomials are used to reduce the proposed problems to system of algebraic equations which are solved by a suitable numerical method. From the obtained numerical results, we can conclude that this method gives results with an excellent agreement with the exact solution. All numerical results are obtained using Matlab program 8.

References

1. Khader, M.M. (2011) On the Numerical Solutions for the Fractional Diffusion Equation. Communications in Nonlinear Science and Numerical Simulation, 16, 2535-2542. http://dx.doi.org/10.1016/j.cnsns.2010.09.007
2. Khader, M.M. (2013) Numerical Treatment for Solving the Perturbed Fractional PDEs Using Hybrid Techniques. Journal of Computational Physics, 250, 565-573. http://dx.doi.org/10.1016/j.jcp.2013.05.032
3. Sweilam, N.H., Khader, M.M. and Mahdy, A.M.S. (2013) Numerical Study for the Fractional Differential Equations Generated by Optimization Problem Using Chebyshev Collocation Method and FDM. Applied Mathematics and Information Science, 7, 2013-2020.
4. Khader, M.M. and Hendy, A.S. (2013) A Numerical Technique for Solving Fractional Variational Problems. Mathematical Methods in Applied Sciences, 36, 1281-1289. http://dx.doi.org/10.1002/mma.2681
5. Sweilam, N.H., Khader, M.M. and Nagy, A.M. (2011) Numerical Solution of Two-Sided Space-Fractional Wave Equation Using Finite Difference Method. Journal of Computional and Applied Mathematics, 235, 2832-2841.http://dx.doi.org/10.1016/j.cam.2010.12.002
6. Sweilam, N.H., Khader, M.M. and Adel, M. (2012) On the Stability Analysis of Weighted Average Finite Difference Methods for Fractional Wave Equation. Fractional Differential Calculus, 2, 17-75. http://dx.doi.org/10.7153/fdc-02-02
7. Sweilam, N.H. and Khader, M.M. (2010) A Chebyshev Pseudo-Spectral Method for Solving Fractional Order Integro-Differential Equations. ANZIAM, 51, 464-475. http://dx.doi.org/10.1017/S1446181110000830
8. Sweilam, N.H., Khader, M.M. and Al-Bar, R.F. (2008) Homotopy Perturbation Method for Linear and Nonlinear System of Fractional Integro-Differential Equations. International Journal of Computational Mathematics and Numerical Simulation, 1, 73-87.
9. Sweilam, N.H., Khader, M.M. and Mahdy, A.M.S. (2012) On the Numerical Solution for the Linear Fractional Klein-Gordon Equation Using Legendre Pseudospectral Method. International Journal of Mathematics and Computer Applications Research, 2, 1-10.
10. Khader, M.M., EL-Danaf, T.S. and Hendy, A.S. (2013) A Computational Matrix Method for Solving Systems of High Order Fractional Differential Equations. Applied Mathematical Modelling, 37, 4035-4050.http://dx.doi.org/10.1016/j.apm.2012.08.009
11. Khader, M.M., Mahdy, A.M.S. and Shehata, M.M. (2014) Approximate Analytical Solution to the Time-Fractional Biological Population Model Equation. Jokull, 64, 378-394.
12. Sweilam, N.H., Khader, M.M. and Mahdy, A.M.S. (2012) Numerical Studies for Solving Fractional-Order Logistic Equation. International Journal of Pure and Applied Mathematics, 78, 1199-1210.
13. Sweilam, N.H., Khader, M.M. and Mahdy, A.M.S. (2012) Numerical Studies for Fractional-Order Logistic Differential Equation with Two Different Delays. Journal of Applied Mathematics, 2012, Article ID: 764894.
14. Reid, W.T. (1972) Riccati Differential Equations Mathematics in Science and Engineering. Academic Press, New York.
15. Lasiecka, I. and Triggiani, R. (1991) Differential and Algebraic Riccati Equations with Application to Boundary/Point Control Problems: Continuous Theory and Approximation Theory. Lecture Notes in Control and Information Sciences, Springer, Berlin.
16. Bahnasawi, A.A., El-Tawil, M.A. and Abdel-Naby, A. (2004) Solving Riccati Differential Equation Using ADM. Applied Mathematics and Computation, 157, 503-514. http://dx.doi.org/10.1016/j.amc.2003.08.049
17. Tan, Y. and Abbasbandy, S. (2008) Homotopy Analysis Method for Quadratic Riccati Differential Equation. Communications in Nonlinear Science and Numerical Simulation, 13, 539-546. http://dx.doi.org/10.1016/j.cnsns.2006.06.006
18. Cushing, J.M. (1998) An Introduction to Structured Population Dynamics, Society for Industrial and Applied Mathematics. http://dx.doi.org/10.1137/1.9781611970005
19. Alligood, K.T., Sauer, T.D. and Yorke, J.A. (1996) An Introduction to Dynamical Systems. Springer, Berlin.
20. Ausloos, M. The Logistic Map and the Route to Chaos: From the Beginnings to Modern Applications XVI. 411, 2006.
21. Pastijn, H. (2006) Chaotic Growth with the Logistic Model of P.-F. Verhulst, Understanding Complex Systems. The Logistic Map and the Route to Chaos. 3-11.
22. Suansook, Y. and Paithoonwattanakij, K. (2009) Dynamic of Logistic Model at Fractional Order. IEEE International Symposium on Industrial Electronics.
23. Fridman, E., Fridman, L. and Shustin, E. (2000) Steady Modes in Relay Control Systems with Time Delay and Periodic Disturbances. Journal of Dynamic Systems, Measurement, and Control, 122, 732-737.http://dx.doi.org/10.1115/1.1320443
24. Davis, C.L. (2002) Modification of the Optimal Velocity Traffic Model to Include Delay due to Driver Reaction Time. Physica A, 319, 557-567. http://dx.doi.org/10.1016/S0378-4371(02)01457-7
25. Kuang, Y. (1993) Delay Differential Equations with Applications in Population Biology. Academic Press, Boston, San Diego, New York.
26. Epstein, I. and Luo, Y. (1991) Differential Delay Equations in Chemical Kinetics: Nonlinear Models: The Cross-Shaped Phase Diagram and the Originator. Journal of Chemical Physics, 95, 244-254.http://dx.doi.org/10.1063/1.461481
27. Mai-Duy, N., See, H. and Tran-Cong, T. (2009) A Spectral Collocation Technique Based on Integrated Chebyshev Polynomials for Biharmonic Problems in Irregular Domains. Applied Mathematical Modelling, 33, 284-299.http://dx.doi.org/10.1016/j.apm.2007.11.002
28. Bhrawy, A.H. and Alofi, A.S. (2013) The Operational Matrix of Fractional Integration for Shifted Chebyshev Polynomials. Applied Mathematics Letters, 26, 26-31. http://dx.doi.org/10.1016/j.aml.2012.01.027
29. Khader, M.M., Sweilam, N.H. and Mahdy, A.M.S. (2011) An Efficient Numerical Method for Solving the Fractional Diffusion Equation. Journal of Applied Mathematics and Bioinformatics, 1, 1-12.
30. Bell, W.W. (1968) Special Functions for Scientists and Engineers. Great Britain, Butler and Tanner Ltd, Frome and London.
31. El-Sayed, A.M.A., El-Mesiry, A.E.M. and El-Saka, H.A.A. (2007) On the Fractional-Order Logistic Equation. Applied Mathematics Letters, 20, 817-823. http://dx.doi.org/10.1016/j.aml.2006.08.013