Discrete Dynamics in Nature and Society
Volume 2015 (2015), Article ID 325364, 9 pages
http://dx.doi.org/10.1155/2015/325364

Research Article

Stability Analysis of One-Leg Methods for Nonlinear Neutral Delay Integrodifferential Equations

Yuexin Yu and Liping Wen

Department of Mathematics, Xiangtan University, Xiangtan 411105, China

Received 8 April 2015; Revised 16 June 2015; Accepted 24 June 2015

Academic Editor: Antonia Vecchio

Copyright © 2015 Yuexin Yu and Liping Wen. 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 paper is concerned with the numerical solution of nonlinear neutral delay integrodifferential equations (NDIDEs). The adaptation of one-leg methods is considered. It is proved that an A-stable one-leg method can preserve the global stability and a strongly A-stable one-leg method can preserve the asymptotic stability of the analytical solution of nonlinear NDIDEs. Numerical tests are given to confirm the theoretical results.

1. Introduction

In this paper, we consider the initial value problem (IVP) of nonlinear neutral delay integrodifferential equations:

where is a constant delay, is a given continuously differential function, and and are given continuous mappings.

Neutral delay integrodifferential equations (NDIDEs) arise widely in scientific and engineering fields such as physics, biology, medicine, economics, and control system (see [1–3] and the references therein). Generally speaking, it is difficult to obtain the analytical solutions of such equations. In view of this, people began to study the numerical solutions of the equations. For the special cases of NDIDEs, such as delay differential equations, delay integrodifferential equations, and neutral delay differential equations, the theory of computational methods has been studied by many authors and a great deal of interesting results have been found in recent 30 years. But, for NDIDEs, only a few results have been presented in the literature. In 2005, Zhao et al. [4] discussed the asymptotic stability of analytical solution and numerical solution (obtained by linear -methods and BDF methods) of linear neutral Volterra delay integrodifferential system:

where , , , , , and the matrix may be singular. Later, Xu and Zhao [5] further considered the asymptotic stability of Runge-Kutta methods for system (2). In 2008, Zhang and Vandewalle [6] dealt with the asymptotic stability of exact and discrete solutions of neutral multidelay integrodifferential equations. Sufficient conditions for the asymptotic stability of the analytical solution have been derived, and the asymptotic stability criteria of Runge-Kutta methods and linear multistep methods were constructed. Wu and Gan [7] investigated a test equation for one-dimension linear NDIDEs and got some delay-dependent stability results.

For the nonlinear NDIDEs (1), Yu and Li [8] discussed the global stability and asymptotic stability of algebraically stable Runge-Kutta methods. Recently, Hu and Huang [9] considered the analytical and numerical stability of nonlinear NDIDEs. Sufficient conditions for the analytical stability of nonlinear NDIDEs are derived, and they proved that any -stable linear multistep method can preserve the asymptotic stability of the analytical solution of nonlinear NDIDEs (1). For another case of nonlinear NDIDEs, namely, the NDIDEs of “Hale’s form”

Yu et al. [10] and Zhang et al. [11, 12] investigated the stability of Runge-Kutta methods and one-leg methods, respectively.

In this paper we are interested in the stability of one-leg methods for nonlinear NDIDEs (1). It is proved that an -stable one-leg method can preserve the global stability and a strongly -stable one-leg method can preserve the asymptotic stability of the analytical solution of nonlinear NDIDEs. Numerical tests are given to confirm the theoretical results in the end.

2. Problem Class and Its Stability

Let denote the inner product and the corresponding norm in space . Assume that the continuous mappings and in problem (1) satisfy the following conditions:

where , , , , and are real constants and . Furthermore, we also consider the function and assume that it is continuous and satisfies

where is a real constant.

Throughout this paper, we assume that problem (1) has a unique exact solution , and we use the symbol to denote the problem class consisting of all of problem (1) satisfying conditions (4)-(7).

Remark 1. When the right-hand side function of problem (1) does not possess the term , problem (1) degenerates into an IVP of delay integrodifferential equations (DIDEs):

The stability of numerical methods for DIDEs has been investigated in [13–17].

Remark 2. When the right-hand side function of problem (1) does not possess the integral term, problem (1) degenerates into an IVP of neutral delay differential equations (NDDEs):

The stability of numerical methods for NDDEs has been researched in [18–22].

Remark 3. When the right-hand side function of problem (1) does not possess the term and the integral term, problem (1) degenerates into an IVP of delay differential equations (DDEs):

The stability results of numerical methods for DDEs can be found in [23–28] and so forth.

For problems of the class , Hu and Huang derived the following stability results (see [9]).

Theorem 4. Suppose problem (1) belongs to the class satisfying . Then one has

where

and denotes the solution of any given perturbed problem of (1):

where is a given continuously differential function.

Theorem 5. Suppose problem (1) belongs to the class satisfying . Then one has

Inequality (11) characterizes the stability property and relation (14) characterizes the asymptotic stability property of problem (1), respectively.

3. Stability Analysis of One-Leg Methods for NDIDEs

Consider using a one-leg -step method (for ordinary differential equations)

to solve problem (1); we have

where , is an arbitrarily given positive integer, , is the translation operator, , is an approximation to , for , and is an approximation to , which can be computed by the repeated trapezoidal rule:

is an approximation to , which is obtained by using the following formula:

where for and and are generating polynomials which are assumed to have real coefficients and no common divisor. We also assume , , and .

Similarly, applying the same method to perturbed problem (13), we have

where and are approximations to and , respectively, for , and is an approximation to , which can be computed by

where for .

For a real symmetric positive matrix , the norm is defined by

Theorem 6. Assume that one-leg method (15) is -stable. Then the numerical solutions and , obtained by using corresponding method (16) to problems (1) and (13) which belong to the class with , respectively, satisfy the global stability inequality

for all , where depends only on the method, , , , , , and , and

Proof. Let

where denotes the integer part; then .
Since -stability is equivalent to -stability (cf. [29]), there is a real symmetric positive definite matrix such that, for any real sequence , the following inequality holds:

where . Therefore, we can easily obtain (cf. [29])

Using condition (4), we have

When , that is, , (27) leads to

When , that is, , using conditions (4)–(7), (27) leads to

where, here and below, we define equal to 0 for . Combining (28) and (29) yields

Substituting (30) into (26) and using condition , we obtain

By induction, (31) gives

Let and denote the maximum and minimum eigenvalues of the matrix , respectively. Then, we have

Hence

Therefore, there is a real constant depending only on the method, , , , , , and such that inequality (22) holds and this completes the proof of Theorem 6.

Remark 7. It is well known that many one-leg methods, such as implicit Euler method, the second-order BDF formula method, and one-leg -methods (), are all -stable. Therefore, in terms of Theorem 6, the corresponding methods are globally stable for solving the nonlinear NDIDEs of the class which satisfies the condition .

Next, we further discuss the asymptotic stability of the one-leg method. One-leg method (15) is called strongly -stable if it is -stable and the modulus of any root of is strictly less than .

Theorem 8. Assume that one-leg method (15) is strongly -stable. Then the numerical solutions and , obtained by using corresponding method (16) to problems (1) and (13) which belong to the class with , respectively, satisfy

for all . Relation (35) characterizes the asymptotic stability property of method (16).

Proof. In terms of the proof of Theorem 6, we have

Since , it is easily obtained from (36) that

By analogy with the proof of Theorem   in [27], we have

and this completes the proof of Theorem 8.

Remark 9. It is well known that many one-leg methods, such as implicit Euler method, the second-order BDF formula method, and one-leg -methods (), are all strongly -stable. Therefore, in terms of Theorem 8, the corresponding methods are asymptotically stable for solving the nonlinear NDIDEs of the class which satisfies the condition .

4. Numerical Experiments

Example 1. Consider the one-dimensional parabolic problem with neutral type

After application of the numerical method of lines, we obtain the following NDIDEs:

where is the spatial step, is a natural number such that , , , and . Then, problem (40) belongs to the class with

where the inner product is standard inner product. We take for the numerical method of lines; thus the condition (<0) is satisfied, which means the analytical solution of problem (40) is stable and asymptotically stable.
We use the 2-step one-leg methods of order 2:

which is -stable and strongly -stable, for solving problem (40) and its perturbed problem, where the initial function of the perturbed problem is
As a comparison, we also use the 2-step one-leg method of order 3:

which is not -stable, for solving problem (40) and its perturbed problem. We denote the numerical solutions of problem (40) and its perturbed problems and , where and are approximations to and , respectively. The values obtained by different methods are listed in Figure 1.

(a)

(b)

Figure 1: Values obtained by method (42) (a) or by method (43) (b) with .

Example 2. Consider the following neutral delay integrodifferential equations, which are the models for population dynamics (see [1]):

We choose , , , and . Then, problem (44) belongs to the class with

Thus the condition (<0) is satisfied, which means the analytical solution of problem (44) is stable and asymptotically stable. We use the 2-step one-leg method (42) for solving problem (44) and its perturbed problem, where the initial functions of problem (44) and its perturbed problem are , and , for , respectively. We denote the numerical solutions of problem (44) and its perturbed problems and , where and are approximations to and , respectively. The values obtained by method (42) with different stepsize are listed in Figure 2.
From Figures 1-2, one can see that the values (or ) obtained by method (42) are bounded and tend to zero. These coincide with Theorems 6 and 8. However, for method (43), the situation is inverse as one can see that the values are divergent as .

(a)

(b)

Figure 2: Values obtained by method (42) with (a) and (b).

5. Conclusion

In the present paper, the adaptation of one-leg methods is applied for solving nonlinear NDIDEs. It is proved that an -stable one-leg method can preserve the global stability and a strongly -stable one-leg method can preserve the asymptotic stability of the analytical solution of nonlinear NDIDEs. Investigating the stability of other numerical methods, for example, multistep Runge-Kutta methods, for solving nonlinear NDIDEs will be our future work.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank an anonymous referee for his many valuable suggestions. This work was supported by the NSF of China (11171282, 11371302) and Scientific Research Fund of Hunan Provincial Education Department.

References

G. A. Bocharov and F. A. Rihan, “Numerical modelling in biosciences using delay differential equations,” Journal of Computational and Applied Mathematics, vol. 125, no. 1-2, pp. 183–199, 2000.

J. K. Hale and S. M. V. Lunel, Introduction to Functional Differential Equations, Springe, Berlin, Germany, 1993.

V. Kolmanovskii and A. Myshkis, Introduction to the Theory and Applications of Functional Differential Equations, Kluwer Academic, Dordrecht, The Netherlands, 1999.

J. J. Zhao, Y. Xu, and M. Z. Liu, “Stability analysis of numerical methods for linear neutral Volterra delay-integro-differential system,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1062–1079, 2005.

Y. Xu and J. J. Zhao, “Stability of Runge-Kutta methods for neutral delay-integro-differential-algebraic system,” Mathematics and Computers in Simulation, vol. 79, no. 3, pp. 571–583, 2008.

C. Zhang and S. Vandewalle, “Stability criteria for exact and discrete solutions of neutral multidelay-integro-differential equations,” Advances in Computational Mathematics, vol. 28, no. 4, pp. 383–399, 2008.

S. F. Wu and S. Q. Gan, “Analytical and numerical stability of neutral delay integro-differential equations and neutral delay partial differential equations,” Computers & Mathematics with Applications, vol. 55, no. 11, pp. 2426–2443, 2008.

Y.-X. Yu and S.-F. Li, “Stability analysis of Runge-Kutta methods for nonlinear neutral delay integro-differential equations,” Science in China—Series A: Mathematics, vol. 50, no. 4, pp. 464–474, 2007.

P. Hu and C. M. Huang, “Analytical and numerical stability of nonlinear neutral delay integro-differential equations,” Journal of the Franklin Institute, vol. 348, no. 6, pp. 1082–1100, 2011.

Y. X. Yu, L. P. Wen, and S. F. Li, “Nonlinear stability of Runge-Kutta methods for neutral delay integro-differential equations,” Applied Mathematics and Computation, vol. 191, no. 2, pp. 543–549, 2007.

C. J. Zhang, T. T. Qin, and J. Jin, “An improvement of the numerical stability results for nonlinear neutral delay-integro-differential equations,” Applied Mathematics and Computation, vol. 215, no. 2, pp. 548–556, 2009.

C. J. Zhang and Y. Y. He, “The extended one-leg methods for nonlinear neutral delay-integro-differential equations,” Applied Numerical Mathematics, vol. 59, no. 6, pp. 1409–1418, 2009.

T. Koto, “Stability of Runge-Kutta methods for delay integro-differential equations,” Journal of Computational and Applied Mathematics, vol. 145, no. 2, pp. 483–492, 2002.

H. Brunner and R. Vermiglio, “Stability of solutions of delay functional integro-differential equations and their discretizations,” Computing. Archives for Scientific Computing, vol. 71, no. 3, pp. 229–245, 2003.

C. J. Zhang and S. Vandewalle, “Stability analysis of Runge-Kutta methods for nonlinear Volterra delay-integro-differential equations,” IMA Journal of Numerical Analysis, vol. 24, no. 2, pp. 193–214, 2004.

C. J. Zhang and S. Vandewalle, “General linear methods for Volterra integro-differential equations with memory,” SIAM Journal on Scientific Computing, vol. 27, no. 6, pp. 2010–2031, 2006.

C. M. Huang, “Stability of linear multistep methods for delay integro-differential equations,” Computers & Mathematics with Applications, vol. 55, no. 12, pp. 2830–2838, 2008.

G.-D. Hu and T. Mitsui, “Stability analysis of numerical methods for systems of neutral delay-differential equations,” BIT Numerical Mathematics, vol. 35, no. 4, pp. 504–515, 1995.

Y. Liu, “Numerical solution of implicit neutral functional-differential equations,” SIAM Journal on Numerical Analysis, vol. 36, no. 2, pp. 516–528, 1999.

A. Bellen, N. Guglielmi, and M. Zennaro, “Numerical stability of nonlinear delay differential equations of neutral type,” Journal of Computational and Applied Mathematics, vol. 125, no. 1-2, pp. 251–263, 2000.

R. Vermiglio and L. Torelli, “A stable numerical approach for implicit non-linear neutral delay differential equations,” BIT Numerical Mathematics, vol. 43, no. 1, pp. 195–215, 2003.

W.-S. Wang, S.-F. Li, and K. Su, “Nonlinear stability of Runge-Kutta methods for neutral delay differential equations,” Journal of Computational and Applied Mathematics, vol. 214, no. 1, pp. 175–185, 2008.

L. Torelli, “Stability of numerical methods for delay differential equations,” Journal of Computational and Applied Mathematics, vol. 25, no. 1, pp. 15–26, 1989.

A. Bellen and M. Zennaro, “Strong contractivity properties of numerical methods for ordinary and delay differential equations,” Applied Numerical Mathematics, vol. 9, no. 3–5, pp. 321–346, 1992.

K. J. In 't Hout, “Stability analysis of Runge-Kutta methods for systems of delay differential equations,” IMA Journal of Numerical Analysis, vol. 17, no. 1, pp. 17–27, 1997.

M. Zennaro, “Asymptotic stability analysis of Runge-Kutta methods for nonlinear systems of delay differential equations,” Numerische Mathematik, vol. 77, no. 4, pp. 549–563, 1997.

C. M. Huang, S. F. Li, H. Y. Fu, and G. N. Chen, “Stability and error analysis of one-leg methods for nonlinear delay differential equations,” Journal of Computational and Applied Mathematics, vol. 103, no. 2, pp. 263–279, 1999.

A. Bellen and M. Zennaro, Numerical Methods for Delay Differential Equations, Clarendon Press, Oxford, UK, 2003.

G. Dahlquist, “-stability is equivalent to -stability,” BIT Numerical Mathematics, vol. 18, no. 4, pp. 384–401, 1978.