A Legendre-Gauss collocation method for neutral functional-differential equations with proportional delays

In this paper, we present a unified framework for analyzing the spectral collocation method for neutral functional-differential equations with proportional delays using shifted Legendre polynomials. The proposed collocation technique is based on shifted Legendre-Gauss quadrature nodes as collocation knots. Error analysis and stability of the proposed algorithm are theoretically investigated under several mild conditions. The accuracy of the proposed method has been compared with a variational iteration method, a one-leg θ-method, a particular Runge-Kutta method, and a reproducing kernel Hilbert space method. Numerical results show that the proposed methods are of high accuracy and are efficient for solving such an equation. Also, the results demonstrate that the proposed method is a powerful algorithm for solving other delay differential equations.


Introduction
One of the fundamental classes of delay differential equations (DDEs) is that of neutral functional-differential equations (NFDEs) with proportional delays. Such equations arise in many areas of science and engineering and play an important role in the modeling of real-life phenomena in other fields of science (cf. [-]). For these reasons, NFDEs have received much attention in the last decades. The principal difficulty in studying DDEs lies in their special transcendental nature. Obviously, most of NFDEs cannot be solved by the well-known exact methods. Therefore, it is highly desirable to design accurate numerical approaches to approximate the solutions of NFDEs.
The theory of DDEs with multiple delays has been analyzed by many authors, and we briefly review some of them. In [] Jackiewicz and Lo proposed and developed the Adams predictor-corrector method for the numerical solution of NFDEs. In [], the authors investigated the Adomian decomposition method for solving a special class of DDEs and established the convergence of this approach.
Some analytical and numerical solutions of a family of DDEs were presented in []. The authors of [] investigated the Runge-Kutta (RK) method for NFDEs with different proportional delays and proved the stability of RK for this equation. Wang and Li [] proposed the one-leg-methods for solving nonlinear neutral functional differential equations. http://www.advancesindifferenceequations.com/content/2013/1/63 We also refer to the articles [-] and the references therein for applying the waveform relaxation method to solve neutral type functional-differential systems.
Over the years, it was found that the spectral methods were a valid method to obtain approximations for differential equations. The solution of a DDE globally depends on its history due to the delay variable, so a global spectral method could be a good candidate for numerical DDEs. In this direction, Ishiwata and Muroya [] proposed the rational approximation algorithm based on Legendre polynomials and Ishiwata et al. [] developed the collocation approximation for solving DDEs. Meanwhile, in [] the authors discussed the tau approach for solving a class of NFDEs, and the solution was expanded by a shifted Legendre polynomial, and the unknown expansion coefficients were obtained by a segmented Lanczos-tau formulation. Our fundamental goal of this paper is to develop a suitable way to approximate the neutral functional-differential equations with proportional delays on the interval [, T] by using the shifted Legendre polynomials. In other words, we propose the spectral shifted Legendre-Gauss collocation (SLC) method to find the solution u N (x). For suitable collocation points, we use the (Nm + ) nodes of the shifted Legendre-Gauss interpolation on the interval [, T]. These equations together with m initial conditions generate (N + ) algebraic equations which can be solved. It should be noted that the basic requirement for using any spectral base (e.g., Legendre polynomials) is the smoothness of the solution of the considered problem. This may be guaranteed by the smoothness of known functions in the neutral equation. By these assumptions (i.e., being smooth), exponential convergence behavior of spectral approximations is exhibited in any test problem.
The paper is organized as follows. Section  is devoted to preliminaries needed hereafter. In Section , we design the shifted Legendre-Gauss collocation technique for NFDEs with proportional delays. The error analysis and stability of the solution are provided in Section  under several mild conditions. In Section , we present some numerical results demonstrating the efficiency of the suggested numerical algorithm. Concluding remarks are given in Section .

Preliminaries
Let P i (x) be the standard Legendre polynomial of degree i, then we have that If α > -, β > - and P (α,β) k (x) is the kth degree of the Jacobi polynomial, then for an integer m, the mth-order derivative of Legendre polynomials is Let w (α,β) (x) = (x) α ( + x) β , then we define the weighted space Ł  w (α,β) (-, ) as usual, equipped with the following inner product and norm: The set of Legendre polynomials forms a complete Ł  (-, )-orthogonal system, and Let T > , then the shifted Legendre polynomial of degree k on the interval (, T) is defined by L T,k (x) = P k ( x T -). Thus L T,k (x) can be obtained by With the aid of the standard Legendre polynomials, we get Next, we define the space Ł  (, T) in the usual way, with the following inner product and norm: The set of shifted Legendre polynomials is a complete Ł  (, T)-orthogonal system. Moreover, due to (), we have

Shifted Legendre-Gauss collocation method
In this section, we develop a spectral Legendre-Gauss collocation approach to analyze the following NFDEs with proportional delays: with the initial conditions Here, a(x) and b n (x) (n = , , . . . , m -) are given analytical functions, and γ n , β, λ n are constants with  < γ n <  (n = , , . . . , m). Now we introduce the Legendre-Gauss quadratures in two different intervals (-, ) and (, T) that will be used in the sequel. Denote by x N,j (x T,N,j ),  ≤ j ≤ N , and N,j ( T,N,j ), http://www.advancesindifferenceequations.com/content/2013/1/63  ≤ i ≤ N , the nodes and Christoffel numbers of the standard (shifted) Legendre-Gauss-Lobatto quadratures on the intervals (-, ), (, T), respectively. Then one can clearly deduce that [] and if S N (, T) denotes the set of all polynomials of degree at most N , then it follows that We define the discrete inner product and norm as follows: Obviously, Thus, for any u ∈ S N (, T), the norms u N and u coincide. We set The following relation for the qth derivative of shifted Legendre polynomials L T,k (x) will be needed for our main results []: where .

(   )
In virtue of (), the high-order derivative of shifted Legendre polynomials with proportional delay can be written as where C q (k, i) is defined in (). http://www.advancesindifferenceequations.com/content/2013/1/63 The shifted Legendre-Gauss collocation method for solving () and () is to seek Now, we derive an efficient algorithm for tackling () and (). Let us expand the numerical approximation in terms of shifted Legendre polynomials We first approximate u(x) and u n (x) as Eq. (). By substituting these approximations in Then, by virtue of (), we deduce that To find the solution u N (x), we collocate Eq. (). Using () and (), at the (Nm + ) shifted Legendre roots, yields () http://www.advancesindifferenceequations.com/content/2013/1/63 Next, Eq. (), after using (), can be written as Thus Eq. () with relation () generate (N + ) of a set of algebraic equations which can be solved for the unknown coefficients a j , j = , , , . . . , N , by using any standard solver technique.

Error analysis and stability of the solution
This section is divided into two subsections. The first subsection is related to presenting a bound for the error of the proposed method; meanwhile, in the second subsection, the stability of the numerical solution is investigated briefly.

Error bound of the method
In this part, the error bound of the method will be provided under several mild conditions such as solution boundedness of the main neutral differential equation. However, some definitions and lemmas should be provided for clarifying the main theorem of this subsection.
where ξ L p denotes the usual Lebesgue norm, and ξ (x) stands for any finite dimensional norm in n .

Lemma . Given a function
, there exists a polynomial u N (x) of degree less than or equal to N such that where C is a constant independent of N , m is the order of smoothness of ξ , and C  = ξ W m,∞ . Here, u N (x) with the smallest norm ξ (x)u N (x) L ∞ is called the N th order best polynomial approximation of ξ (x) in the norm of L ∞ .
Note that if ξ ∈ C ∞ , then m = ∞. This implies that u N (x) converges to ξ at a spectral rate, i.e., it is faster than any given polynomial rate. Moreover, we denote the set of continuous Again, we consider Eq. () with the initial conditions (). For clarity of presentation, we assume that m = . A similar procedure can be applied for higher values of m. Therefore, ()-() can be written as follows: Integrating the aforementioned equation in the interval [, x] yields Accordingly, we can rewrite the above-obtained equation in the following form: . In the following theorem, we show that the approximate solution which was expressed in terms of Legendre polynomials converges to the exact solution under several mild conditions. Proof Since the known functions a(x), g(x), andb  (x) are approximated in terms of shifted Legendre polynomials, then the obtained solution is an approximated polynomial in the form of u N (x). Our aim is to find an upper bound for the associated error between u(x) and u N (x) for (). We rewrite () in the following form:

Theorem . Consider Eq. () again. Assume u(x) and u N (x) are the exact and approximate solutions of (). Also, let the approximations of a(x), g(x), andb
According to the assumptions, one reaches the following equation: According to the trigonometric inequality, we have

Subtraction of the above equations yields
Since  ≤ γ  ≤  and x ∈ [, T], the above inequality reduces to a simple form as follows: In other words, For clarity of presentation, we rewrite the above inequality in the following form: ). This conclusion is made because of the smoothness of a(x), g(x), andb  (x). These assumptions . This completes the proof. where p * is the best approximation of N th order of u and

Theorem . Let us consider u(x) ∈ C[, T] together with its Lagrange interpolation polynomial I N u(x) which is based on the shifted Gaussian points; then
are the Lebesgue constant and the Lebesgue function, respectively. Also, F j (x) stands for the jth Lagrange polynomial which is constructed by using shifted Gaussian points.
Proof Since u(x) ∈ C[, T], the best approximating polynomial for u, say p * , exists and we immediately have However, both p * and I N u are N th order polynomials The Lebesgue function and the Lebesgue constant both depend only on the choice of interpolation nodes. The above theorem indicates that a good set of grid points is one which minimizes the Lebesgue constant for interpolation. We know that on equally spaced interpolation points, the Lebesgue constant blows up exponentially with the increasing of N (as made apparent by the classical Runge phenomenon). On the contrary, the best possible Lebesgue constant among all distributions of N prohibits only a logarithmic growth with N . Fortunately, the main families of Gaussian points (Gauss, Gauss-Radau, Gauss-Lobatto) for the Legendre or Chebyshev weights have Lebesgue constants that grow logarithmically or sub-linearly with N .

Numerical results
To illustrate the effectiveness of the proposed method in the present paper, several test examples are carried out in this section. Comparisons of the results obtained by the present method with those obtained by other methods reveal that the present method is very effective and convenient. We consider the following examples.
Example  Consider the first-order neutral functional-differential equation with proportional delay used in [] which has the exact solution xe -x .
We apply the proposed shifted Legendre-Gauss collocation method at t =  to (). In Table  The logarithmic graph of absolute coefficients of shifted Legendre polynomials of a neutral functional-differential equation with proportional delay is shown in Figure , which shows that the method has exponential convergence rate. Absolute errors obtained by the SLC method with N =  are plotted in Figure . Example  Consider the following first-order neutral functional-differential equation with proportional delay [, , , , ]:  In Table  Example  Consider the following second-order neutral functional-differential equation with proportional delays [, , , , ]:   In Table      Example  Consider the following third-order neutral functional-differential equation with proportional delays: and the function f (x) is chosen such that the analytical solution u(x) = log( + sin( π x  )).
In Table 

Conclusions and future works
In this paper, we have proposed a new collocation method based on the shifted Legendre polynomials to numerically solve the neutral functional-differential equations with studying the error analysis of the proposed method. The comparison of the obtained results with those based on other methods shows that the present method is a powerful mathematical tool for finding the numerical solutions of such equations. High accuracy in long computational intervals and the stability of the approximated solutions encourage us to apply a similar method for solving other applied mathematics problems (see, for instance, []) in the future.