Abstract

The modified generalized Laguerre-Gauss collocation (MGLC) method is applied to obtain an approximate solution of fractional neutral functional-differential equations with proportional delays on the half-line. The proposed technique is based on modified generalized Laguerre polynomials and Gauss quadrature integration of such polynomials. The main advantage of the present method is to reduce the solution of fractional neutral functional-differential equations into a system of algebraic equations. Reasonable numerical results are achieved by choosing few modified generalized Laguerre-Gauss collocation points. Numerical results demonstrate the accuracy, efficiency, and versatility of the proposed method on the half-line.

1. Introduction

The fractional calculus which means the calculus of derivatives and integral of any order started to be used as a powerful tool in various disciplines from science and engineering [1, 2]. The models involving fractional derivatives and integrals have memory; therefore it has proven to be very suitable for the description of memory and hereditary properties of various processes [311]. Interested readers can also see [1221].

Neutral functional-differential equations play an important role in the mathematical modeling of several phenomena. It is well known that most of delay differential equations cannot be solved exactly. Therefore, numerical methods would be presented and developed to get approximate solutions of these equations. The rational approximation scheme has been proposed by Ishiwata and Muroya [22], for approximating the solution of delay differential equations. Moreover, the authors of [23] investigated the variational iteration method to find a closed form analytical solution of a class of the neutral functional-differential equations. Recently, the authors of [24] proposed the Chebyshev cardinal functions in combination of operational matrices for numerical solution of delay differential equations which arise in electrodynamics. Very recently, Tohidi et al. [25] proposed and developed an efficient collocation method which depends on Bernoulli operational matrix for numerical solution of a class of delay differential equations. Meanwhile, Doha et al. [26] proposed an accurate approximate solution of generalized pantograph equations on the half-line using Jacobi rational-Gauss collocation method. In this paper we propose an approximate solution of a class of delay differential equations, namely, fractional neutral functional-differential equations (FNFDEs) with proportional delay.

In the last decade or so, comprehensive research has been accomplished on the development of numerical algorithms which are numerically stable for both linear and nonlinear FDEs. Tripathi et al. [27] presented a new operational matrix of hat functions to solve linear FDEs. The spectral tau method was proposed in [28] to achieve an accurate solution of linear and nonlinear FDEs subject to multipoint conditions. In [29], Ma and Huang developed spectral collocation method for solving linear fractional integrodifferential equations. Yang and Huang [30] analyzed and developed the Jacobi collocation scheme for pantograph integrodifferential equations with fractional orders in finite interval. In [31], Yin et al. proposed a new fractional-order Legendre function with spectral method to solve partial FDEs; based on the operational matrix of these functions, the same authors developed their approach in combination with variational iteration formula to solve a class of FDEs; see [32]. More recently, the Jacobi Galerkin method was extended in [33] to solve stochastic FDEs, for fractional differential equations in unbounded domains. The operation matrices of fractional derivatives and fractional integrals of generalized Laguerre polynomials were investigated for solving multiterm FDEs on a semi-infinite interval; see [11, 12]. The generalized Laguerre spectral tau and collocation techniques were given in [12] to solve linear and nonlinear FDEs on the half- line. These spectral techniques were developed and generalized by using the modified generalized Laguerre polynomials in [3436].

Our fundamental goal of this paper is to develop a suitable way to approximate the neutral fractional functional-differential equations with proportional delays on the interval using the modified generalized Laguerre polynomials; we propose the spectral modified generalized Laguerre-Gauss collocation (MGLC) method to find the solution . The modified generalized Laguerre spectral collocation (MGLC) approximation, which is more reliable, is employed to obtain approximate solution of neutral fractional functional-differential equations with proportional delays of order and initial conditions. For suitable collocation points we use the nodes of the modified generalized Laguerre-Gauss interpolation on . These equations together with initial conditions generate algebraic equations which can be solved. Finally, the accuracy of the proposed methods is demonstrated by test problems; numerical results are presented in which the usual exponential convergence behavior of spectral approximations is exhibited.

This paper is organized as follows. In Section 2 we give an overview of modified generalized Laguerre polynomials and their relevant properties needed hereafter, and in Section 3, the way of constructing the collocation technique for neutral fractional functional-differential equations with proportional delays is described using the modified generalized Laguerre polynomials. In Section 4, we present some numerical results exhibiting the accuracy and efficiency of our numerical algorithms. The last section offers some obtained conclusions.

2. Some Basic Preliminaries

We give some definitions and properties of fractional derivatives and modified generalized Laguerre polynomials.

2.1. The Fractional Derivative in the Caputo Sense

The two most commonly used definitions are the Riemann-Liouville operator and the Caputo operator.

Definition 1. The Riemann-Liouville fractional integral operator of order is defined as

Definition 2. The Caputo fractional derivatives of order are defined as where is the classical differential operator of order .

For the Caputo derivative we have where and are the ceiling and floor functions, respectively, while and .

Caputo’s fractional differentiation is a linear operation, similar to the integer-order differentiation where and are constants.

2.2. Properties of Modified Generalized Laguerre Polynomials

Let and let be a weight function on in the usual sense. Define equipped with the following inner product and norm:

Next, let be the modified generalized Laguerre polynomial of degree for and is defined by Moreover, for and , we have where and .

The set of modified generalized Laguerre polynomials is the -orthogonal system; namely, where is the Kronecker function and .

The modified generalized Laguerre polynomials of degree on the interval are given by where .

The special value will be of important use later.

Since the analytic form of the modified generalized Laguerre polynomials of degree is given by (11), by using (4), (5), and (11) we get

Now, approximating by terms of modified generalized Laguerre series, we have where is given from (18) with , and

Employing (13)–(15) we get where , and

A function may be expressed in terms of modified generalized Laguerre polynomials as In practice, only the first terms of modified generalized Laguerre polynomials are considered. Then we have

2.3. Modified Generalized Laguerre-Gauss and Gauss-Radau Quadratures

Let and , , be the zeros of and , respectively. They are arranged in ascending order. Denote by ,  , , , the corresponding Christoffel numbers such that where For the Gauss-Radau weights, we have Note that the earlier two types of quadratures have close relations:

In the next section, we will extend the application of spectral methods (see [3742]) by using the modified generalized Laguerre polynomials to solve fractional neutral functional-differential equations on a semi-infinite domain.

3. Fractional Neutral Functional-Differential Equations

In this section, we use the modified generalized Laguerre pseudospectral method based on the modified generalized Laguerre-Gauss quadrature points as the collocation nodes, to solve numerically the following model problem: with the initial conditions Here, and are given analytical functions; , , and , , , denote given constants with . By using the modified generalized Laguerre-Gauss collocation method, we can approximate the fractional neutral functional-differential equations with proportional delays on a semi-infinite domain directly, without any artificial boundary and variable transformation. Let us first introduce some basic notation that will be used in the sequel. We set and we define the discrete inner product and norm as follows: where and are the nodes and the corresponding weights of the modified generalized Laguerre-Gauss quadrature formula on the interval , respectively. Obviously, Thus, for any , the norms and coincide.

Associating with this quadrature rule, we denote by the modified generalized Laguerre-Gauss interpolation:

The modified generalized Laguerre-Gauss collocation method for solving (24) and (25) is to seek , such that We now derive the algorithm for solving (24) and (25). To do this, let

We first approximate and , as (31). By substituting these approximations in (24), we get Making use of (16), we deduce that Also, by substituting (31) into (25) we obtain

Now, collocating (33) at the modified generalized Laguerre-Gauss interpolation points yields

Next (34), after using (12), can be written as

Finally, (35) with relation (36) generates set of algebraic equations which can be solved for the unknown coefficients , , by using any standard solver technique.

4. Numerical Results

In this section, we implement the proposed method to solve two numerical examples from the fractional neutral functional-differential equations with proportional delays in semi-infinite intervals to demonstrate the accuracy and capability of the proposed algorithm. The absolute errors in the given tables are the values of at selected points.

Example 1. Consider the following fractional neutral functional-differential equation with proportional delay: where
and the exact solution is given by .

Table 1 lists the results obtained by the modified generalized Laguerre collocation method in terms of maximum absolute errors at with , , , , and . In the case of ,  , the approximate solution by the presented method is shown in Figure 1, to make it easier to compare with the exact solution.

Example 2. Consider the following fractional neutral functional-differential equation with proportional delay: subject to with exact solution .

In Table 2, we list the absolute errors obtained by the modified generalized Laguerre collocation method, with several values of , and at . Meanwhile, Figure 2 presents the MGLC solution with , at and exact solution, which are found to be in excellent agreement.

5. Conclusions

We have presented the modified generalized Laguerre-Gauss collocation scheme for approximating the solution of fractional neutral functional-differential equations with proportional delays on the half-line. The technique was based on modified generalized Laguerre polynomials and Gauss quadrature integration of such polynomials, which is defined on the half-line. In the previous section, two numerical results were given to demonstrate the accuracy, efficiency, and versatility of the proposed method on the half-line. Indeed, reasonable numerical results were achieved by choosing few modified generalized Laguerre-Gauss collocation points.

Conflict of Interests

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

Acknowledgment

This paper was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under Grant no. (4-135-35-RG). The authors, therefore, acknowledge DSR technical and financial support.