**ABSTRACT**

Neutral Delay Differential Systems (NDDSs) arise in many areas of various mathematical modeling. Infectious diseases, population dynamics, physiological and pharmaceutical kinetics and chemical kinetics, the navigational control of ships and mechanical systems, chemical process simulation and optimal control are the main field concerning with NDDSs. The purpose of this study was to present an extension of the algebraic formulation of the Operational Tau Method (OTM) for the numerical solution of NDDSs. The proposed method converts the delay parts of the desired NDDS to some operational matrices. Then the NDDS reduces to a set of algebraic equations. Some orthogonal bases including shift Chebyshev and shifted Legendre polynomials are used to decrease the volume of computations. Two illustrative linear and nonlinear experiments are included to show the high accuracy and efficiency of the proposed method.

PDF Abstract XML References Citation

**Received:**March 18, 2011;

**Accepted:**May 04, 2011;

**Published:**June 07, 2011

####
**How to cite this article**

*Journal of Applied Sciences, 11: 2585-2591.*

**DOI:**10.3923/jas.2011.2585.2591

**URL:**https://scialert.net/abstract/?doi=jas.2011.2585.2591

**INTRODUCTION**

Spectral methods and spectral modelings provide a computational approach which achieved substantial popularity in the last three decades (Ghafarian *et al*., 2011; Taiwo and Abubakar, 2011; Tasci, 2003). Tau method is one of the most important spectral methods which is extensively applied for numerical solution of many problems. This method was invented by Lanczos (1938) for solving Ordinary Differential Equations (ODEs) and then the expansion of the method were done for many different problems such as Partial Differential Equations (PDEs) (Matos * et al*., 2004; Kong and Wu, 2009; Doha and Abd-Elhameed, 2005), Integral Equations (IEs) (El-Daou and Khajah, 1997), Integro-Differential Equations (IDEs) (Pour-Mahmoud *et al*., 2005) and etc. (Rahimi-Ardabili and Shahmorad, 2007; Garcia-Olivares, 2002; Parand and Razzaghi, 2004).

A time delay phenomenon is encountered in a wide variety of scientific and engineering applications, including circuit analysis, computer-aided design, real-time simulation of mechanical systems, chemical process simulation, optimal control, population dynamics and vibrating masses attached to an elastic bar (Hale and Verduyn Lunel, 1993; Taiwo and Odetunde, 2010; Rao *et al*., 2011; Vanani and Aminataei, 2009, 2010; Kolmanovskii and Myshkis, 1992; Salamon, 1984).

In this study, we are interested in solving NDDSs with an operational approach of the Tau method. Because in the Tau method, we are dealing with a system of equations wherein the matrix of unknown coefficients is sparse and can be easily invertible. Also, the delay parts appearing in the equation are replaced by their operational Tau representation. Then, we obtain a system of algebraic equations wherein its solution is easy.

**Operational tau method:** In this section, we state some preliminaries and notations using in this study.

For any integrable functions Ψ (x) and φ (x) on (a, b), we define the scalar product <, > by:

where, and ω (x) is a weight function. Let L^{2}_{ω} [a, b] be the space of all functions f: [a, b] →R with ||f||^{2}_{ω} < ∞.

The main idea of the method is to seek a polynomial to approximate u (x) ε L^{2}_{ω} [a,b]. Let ΦX be a set of arbitrary orthogonal polynomial bases defined by a lower triangular matrix Φ and X_{x} = [1, x, x^{2}, ...]^{T}.

**Lemma 1:** Suppose that u (x) is a polynomial as , then we have:

(1) |

(2) |

and:

(3) |

where, u = [u_{0}, u_{1}, ....,u_{n}, ....], x_{a} = [1, a, a^{2}, ...]^{T}, a ε and M, N and P are infinite matrices with only nonzero elements:

**Proof:** (Liu and Pan, 1999).

Let us consider:

(4) |

To be an orthogonal series expansion of the exact solution where, is a vector of unknown coefficients, Φ X_{x} is an orthogonal basis for polynomials in .

In the Tau method, the aim is to convert the linear and nonlinear terms to an algebraic system using some operational matrices. Therefore, we state the following lemma.

**Lemma 2:** Let X_{x} = [1, x, x^{2}, ...]^{T}, u = [u_{0}, u_{1}, u_{2}, ....] be infinite vectors and Φ = [φ_{0}|φ_{1}|φ_{2}|...], φ_{i} are infinite columns of matrix Φ. Then, we have:

(5) |

where, U is an upper triangular matrix as:

(6) |

In addition, if we suppose that u (x) = uΦ X_{x} represents a polynomial, then for any positive integer p, the relation:

(7) |

is valid.

**Proof:** We have:

Therefore:

If we call the last upper triangular coefficient matrix as U, then we have:

Now, in order to prove Eq. 7, we apply induction. For p = 1, it is obvious that u (x) = u Φ X_{x}. For p = 2 we rewrite u^{2} (x) = u Φ X_{x} u Φ X_{x}, = u Φ (X_{x} u Φ X_{x}) and using Eq. 5, we have:

Therefore, Eq. 7 is hold for p = 2. Now, suppose that Eq. 7 is hold for p = k, then we must prove that the relation is valid for s = k + 1 . Thus:

So, Eq. 7 is proved.

**Application on NDDSs:** Let us consider the following NDDS:

(8) |

Where:

(9) |

is the state vector and:

(10) |

Such that and are delay functions; A (x), B (x) and C (x) are m-dimensional matrices which their entries are complex functions of x. Also:

(11) |

(12) |

represent the initial vector function and known vector function, respectively.

Now the aim is to write u_{i} (α_{i} (x)) and (β_{i} (x)) i = 0, 1..., m; in operational forms. Using Eq. 4, we have:

(13) |

We know that X_{x} = [1, x, x^{2},...]^{T}, therefore:

By approximating each power of α_{i} (x) as we obtain:

If be the last coefficient matrix then we have:

(14) |

Substituting Eq. 14 in Eq. 13 we get:

(15) |

Also, from Eq. 1 and 4, it is obvious that:

In the same manner from Eq. 13 to 15, there exist the coefficient matrices Δ_{i} such that:

(16) |

In next step we desire to approximate each elements of matrices A (x), B (x) and C (x) in operational forms. Since, each elements of A (x), B (x) and C (x) are smooth functions therefore we can approximate them as follows:

(17) |

Substituting above equations in Eq. 8 and using Eq. 2, we obtain:

(18) |

Therefore:

(19) |

In the same way, we have:

(20) |

(21) |

The vectors U (x) and F (x) also can be considered as:

(22) |

(23) |

Thus Eq. 8 is replaced by the following algebraic system:

(24) |

So, the residual matrix R (x) of Eq. 8, can be written as:

(25) |

Where:

Now, we set the residual matrix or we use the following inner products:

(26) |

For supplementary conditions of Eq. 8 we have:

(27) |

Therefore, imposing supplementary conditions and setting , a system of algebraic equations is obtained. Since, somewhere we require finite terms of approximation, then we must truncate the series solution to finite number of terms. This is the so-called operational Tau method which is applicable for finite, infinite, regular and irregular domains.

**Some shifted orthogonal polynomials:** We have considered OTM based on arbitrary orthogonal polynomials. Orthogonal functions can be used to obtain a good approximation for transcendental functions. Since shifted Chebyshev and Legendre polynomials are more applicable orthogonal functions for a wide range of problems therefore, we consider them, briefly.

**Shifted chebyshev polynomials:** The Chebyshev polynomials are defined on [-1, 1] as:

(28) |

or:

and shifted Chebyshev polynomials are defined as:

(29) |

Now, we consider the following lemma.

**Lemma 3:**. Suppose that T and T* are coefficient matrices of Chebyshev polynomials {T_{i} (x)| x ε [-1, 1], i = 0, 1, 2,...} and shifted Chebyshev polynomials {T*_{i} (x)| x ε [a, b], i = 0, 1, 2,...} , respectively. Hence, we have:

T* = TQ |

Where:

with v = 2/b-a and w = a + b/a-b.

**Proof:** Definition of T states that:

We know that , thus:

If we let Q to be the last coefficient matrix, then:

so

T* = TQ |

Therefore, the lemma is valid.

**Shifted legendre polynomials:** The Legendre polynomials on [-1, 1] are defined as:

(30) |

and we define shifted Legendre polynomials as:

(31) |

In a similar manner with lemma 3 we can prove P* = PC, where, P and P* are coefficient matrices of Legendre and shifted Legendre polynomials, respectively.

**Illustrative numerical experiments:** In this section, two experiments of NDDSs are given to illustrate the efficiency of the method. In all experiments, we consider the shifted Chebyshev and Legendre polynomials as basis functions and have compared the obtained results with the exact solutions. The computations associated with the experiments discussed above were performed in Maple 14 on a PC with a CPU of 2.4 GHz.

**Experiment 1:** Consider the following NDDS (Vanani and Aminataei, 2010):

The exact solution in the interval [0, 1] is:

We have solved this experiment by OTM for n = 20 with shifted Chebyshev and Legendre bases and compared with the exact solution. Results are given in Table 1 and 2 for u_{1} (x), u_{2} (x) and u_{3} (x), respectively.

From the numerical results in Table 1 and 2, it is easy to conclude that obtained results by OTM are in good agreement with the exact solution. Also, during the running of programs we find out the run time of OTM is 0.952 sec. Therefore, the algorithm of OTM is fast.

**Experiment 2:** Consider the following nonlinear NDDS (Vanani and Aminataei, 2009):

The exact solution is:

We have solved this experiment by OTM for n = 20 with shifted Chebyshev and Legendre bases. Results are given in Table 3 and 4 for u_{1} (x), u_{2} (x), u_{3} (x) and u_{4} (x), respectively. Numerical results in Tables 3 and 4 illustrate a good agreement between OTM solutions and exact solutions. In this experiment, the run time of OTM is 1.607 sec. Again, we can conclude that OTM is a fast method.

Table 1: | Exact and approximate solutions of u_{1} (x) and u_{2} (x) of experiment 1 |

Table 2: | Exact and approximate solution of u_{3} (x) of experiment 1 |

Table 3: | Exact and approximate solution of u_{1} (x) and u_{2} (x) of experiment 2 |

Table 4: | Exact and approximate solution of u_{1} (x) and u_{2} (x) of experiment 2 |

**CONCLUSION **

In the present study, OTM is proposed for solving NDDSs. Reducing the NDDSs to algebraic equations is the first characteristic of the proposed method. The main idea of the proposed method is to convert the NDDS including linear and nonlinear terms to an algebraic system to simplify the computations. Arbitrary orthogonal polynomial bases were applied as basis functions to reduce the volume of computations. Furthermore, this method yields the desired accuracy only in a few terms in a series form of the exact solution. All of these advantages of the OTM to solve nonlinear problems assert the method as a convenient, reliable and powerful tool.

####
**REFERENCES**

- Ghafarian, P., M.T. Zamanian and T. Haghroosta, 2011. Design of barotropic spectral model for predicting of vorticity field. Res. J. Environ. Sci., 5: 386-393.

CrossRef - Taiwo, O.A. and A. Abubakar, 2011. An application method for the solution of second order non linear ordinary differential equations by chebyshev polynomials. Asian J. Applied Sci., 4: 255-262.

CrossRef - Tasci, F., 2003. On the approximative solution of boundary value problems by collocation. J. Applied Sci., 3: 210-215.

CrossRefDirect Link - Matos, J., M.J. Rodrigues and P.B. Vasconcelos, 2004. New implementation of the Tau method for PDEs. J. Comput. Applied Math., 164-165: 555-567.

CrossRef - Kong, W. and X. Wu, 2009. Chebyshev tau matrix method for Poisson-type equations in irregular domain. J. Comput. Applied Math., 228: 158-167.

CrossRef - Doha, E.H. and W.M. Abd-Elhameed, 2005. Accurate spectral solutions for the parabolic and elliptic partial differential equations by the ultraspherical tau method. J. Comput. Applied Math., 181: 24-45.

CrossRefDirect Link - El-Daou, M.K. and H.G. Khajah, 1997. Iterated solutions of linear operator equations with the Tau method. Math. Comput., 66: 207-213.

Direct Link - Pour-Mahmoud, J., M.Y. Rahimi-Ardabili and S. Shahmorad, 2005. Numerical solution of the system of Fredholm integro-differential equations by the Tau method. Applied Math. Comput., 168: 465-478.

CrossRefDirect Link - Rahimi-Ardabili, M.Y. and S. Shahmorad, 2007. Iterative numerical solution of non-linear integro-differential equations by the Tau method. Applied Math. Comput., 193: 514-522.

CrossRef - Garcia-Olivares, A., 2002. Analytical approximants of time-dependent partial differential equations with Tau methods. Math. Comput. Simul., 61: 35-45.

CrossRef - Parand, K. and M. Razzaghi, 2004. Rational chebyshev tau method for solving volterra`s population model. Applied Math. Comput., 149: 893-900.

CrossRef - Taiwo, O.A. and O.S. Odetunde, 2010. On the numerical approximation of delay differential equations by a decomposition method. Asian J. Math. Stat., 3: 237-243.

CrossRefDirect Link - Rao, T.P., K.S. Rao and C.L. Usha, 2011. Stochastic modeling of blood glucose levels in type-2 diabetes mellitus. Asian J. Math. Stat., 4: 56-65.

CrossRefDirect Link - Vanani, S.K. and A. Aminataei, 2009. Multiquadric approximation scheme on the numerical solution of delay differential systems of neutral type. Math. Comput. Modell., 49: 234-241.

CrossRef - Liu, K.M. and C.K. Pan, 1999. The automatic solution to systems of ordinary differential equations by the tau method. Comput. Math. Appl., 38: 197-210.

CrossRef