A rational Chebyshev (RC) spectral collocation technique is considered in this paper to solve high-order linear ordinary differential equations (ODEs) defined on a semi-infinite domain. Two definitions of the derivative of the RC functions are introduced as operational matrices. Also, a theoretical study carried on the RC functions shows that the RC approximation has an exponential convergence. Due to the two definitions, two schemes are presented for solving the proposed linear ODEs on the semi-infinite interval with the collocation approach. According to the convergence of the RC functions at the infinity, the proposed technique deals with the boundary value problem which is defined on semi-infinite domains easily. The main goal of this paper is to present a comparison study for differential equations defined on semi-infinite intervals using the proposed two schemes. To demonstrate the validity of the comparisons, three numerical examples are provided. The obtained numerical results are compared with the exact solutions of the proposed problems.

1. Introduction

In the spectral methods, the most common basis functions are Chebyshev polynomials (CPs), which play an important role in the interpolation problems. Many researchers considered CPs to solve differential equations in the finite domain (see [18] and [9]), but they often fail in the larger domain, also if the exact solution of the problem was in a rational form. For this reason, the rational Chebyshev (RC) functions are applied in the large domain where , which provide a major success in dealing with differential equations (DEs) defined in the open domain . Many researchers studied RC functions for treating man different problems of differential, integrodifferential equations (IDEs), partial, and some other physical-engineering problems as in [10, 11] and [12]. Abbasbandy et al. [13] applied the RC collocation method to get numerical solution of the magnetohydrodynamic flow (MHF) of an incompressible viscous fluid (VF) over a stretching sheet problem. Ramadan et al. in [14, 15] and [16], Yalçınbas et al. in [17], and Parand and Razzaghi in [18] are scrupulous in the use of RC functions to express the approximate solution of high-order ordinary differential equations (ODEs) by different spectral approaches. In [19], Parand and Razzaghi introduced RC functions for solving a population growth of a species within a closed system, named as a Volterra model, where the authors converted the Volterra population model first to an equivalent nonlinear ODE; the solution is approximated by the RC functions with the unknown coefficients. The authors of [20, 21] introduced the RC function approximation with the collocation technique for solving the natural convection of the Darcian fluid (DF) about a vertical full cone embedded in porous media (PM) with a prescribed wall temperature. Ramadan et al. [22, 23] obtained an approximate solution of the applied collocation method based on RC functions to treat high-order linear IDEs with variable coefficients. In [24], the authors applied the RC collocation approach for approximating nonlinear biomathematical problems, namely, the systematic logistic growth, the Lotka-Volterra system (prey-predator model), the simple two-species Lotka-Volterra model, and the prey-predator model with limit cycle (periodic behavior).

All the aforementioned work either relied on the RC functions as a basis defined of an open interval or used it to treat a specific application or used truncated matrices. The truncation in matrices was handled for the first time by us in [15], and it was also an application on open period equations. In this work, a comparison study for solving linear ODEs defined on semi-infinite domains using a spectral collocation method is presented. The equation under investigation is a high-order nonhomogeneous linear variable coefficient ordinary differential equation defined on semi-infinite domain . Two algorithms are considered by the collocation method with the RC functions as basis functions. The matrices of derivative for RC functions introduced with regular definition with truncation and improved one, which lead us to two different schemes. The collocation points that are used here are suggested by us, and they are defined on the interval which deal with the boundary value problems defend on or if one of the mixed conditions tends to infinity. So, we suggested new collocation points that are valid for dealing with a problem of this type.

2. Preliminaries

In this section, the definition and some properties of the RC functions are listed; also, the convergence for RC functions will be discussed.

The Chebyshev polynomials are an orthonormal system in the closed interval [-1, 1], where the weight function for is , and they may be generated using the recurrence formulae: with the initials

For more details about , see Ref. [8].

2.1. The RC Functions

The RC functions are orthonormal on the open interval defined as and they form an orthonormal set of functions with respect to the weight function , and they can be generated with the aid of the following recurrence formulae: with the initials and the property of the orthogonality is where for all and is the Kronecker delta function.

Let , and we note that is the eigen function of the singular Sturm-Liouville problem of the following form:

2.2. Function Spaces

In this subsection, the order of convergence for RC functions will be discussed; let us begin with assuming that represent the space functions, and the inner product is denoted here as such that

Subsequently, from the property of orthogonality (relation (6)), we get the fact that RC functions form a set of orthonormal basis of . Also, let us define the normed spaces and as where and is a positive integer constant, and we let be the Sturm-Liouville operator in (7), namely, it may be written as

Let be a positive integer such that , and .

Theorem 1. For any , and is a generic positive constant independent of any variable, and , then The proof of Theorem 1 is given in [25].

Since the set of RC functions is orthonormal and a complete set, we assume that is defined over the space ; then, it be may expanded in terms of RC functions as where

Infinite series expression (14) represents as a spectral truncated approximation as follows

The order of convergence for the RC function approximation will be obtained using several orthonormal projections. From (16), it is clear that is the orthogonal projection of onto with respect to the inner product (10). For all of the above, especially Theorem 1, the following theorem is presented and contains the order of convergence of RC functions.

Theorem 2. For any function such that , where , there exists a positive constant independent of such that The complete proof of Theorem 2 is found in [13] (or see Ref. [25] for more details); this theorem shows that the RC approximation has exponential convergence.

2.3. Operational Matrix

This subsection introduces the form of operational matrix for the RC functions; the derivative of the vector can be expressed by where is operational differentiation matrix. The elements of are found by differentiating (4) and using , then also using the multiplication relation:

The approximation sign in (18) made by a truncation to the last column of (by consideration that , for ) to get an invertible square operational matrix (see Ref. [14]). The structure of the matrix is obtained as a lower-Heisenberg matrix. The matrix can be expressed as , where is a tridiagonal matrix which is obtained from and the entire elements of matrix are , obtained from

In addition to and for .

Consequently, the -order derivative of the row matrix , which is given in (18), is obtained as

And we note here that , , for . Definitions (18) and (23) were introduced for the first time in [18, 19], and many works have used them; see, for example, Refs. [2022], [23, 24], and [17].

3. The Improved Differentiation of the RC Functions

In the present section, an improved definition of differentiation for the RC functions is introduced. There was a need to find an improvement to definition (23), because in the higher derivatives, when using this definition, a week approximation is obtained. The truncated definition (18) and the -order derivative (23) give us a regular truncated differentiation of the RC functions (RRC). Generally, the derivative of the rational or fractional functions increases the order of the denominator (in contrast to polynomials that reduce the order at differentiation), so the truncation increases as the order of the derivative increases (more than one column in Equation (18)). For example, the fourth-order derivative of the vector at , the truncated terms in the row vector are (the last four terms) . This will lead to unsatisfied approximating in high-order DEs using the presented RRC definition (18). Therefore, an improved definition of the derivative of RC functions will be proposed next.

First, a vector will be inserted into (18) to treat the truncated terms, which will improve the regular definition. This manner will be called an improved derivative of the RC functions and will be indicated by IRC.

Additionally, the first-order derivative of the row matrix is where

Theorem 3. The derivative of the matrix vector , in terms of itself, is defined as where

Proof. Using the assumption (24) as the first derivative, and differentiating (24), then we get or and by induction, we get the derivative as relation (26).☐

As a special case, if , it leads us to the regular RC definition in relations (18) and (23). The introduced definition (24) and the -order derivative (26) give us an improved differentiation of the RC functions.

4. Problem Statement

In this study, the form of high-order ODEs which represents a linear nonhomogeneous with variable coefficients defined on a semi-infinite domain is which forms -order ODEs; the previous forms of DEs are subjected to the following conditions: where the and are well-defined functions on , and and are constants (initial value problem), where may tend to (boundary value problems).

Now, we consider that the approximate solution according to (16) for the exact solution of equation (30) in the vector form as where

5. Fundamental Relation-Based Matrix Forms

In the beginning, we provide the fundamental matrix relation of the solution of (30) by two schemes using the RC collocation approach.

Assuming that the solutionof (30) can be expressed as relation (32), which is a truncated RC series, thenand itsderivativeare written in the matrix forms (32) and (33) such thatand, where,are the RC coefficients to be determined later.

Now, let the collocation points as and at the end points , namely,

Permanently, the RC functions are convergent at both end points 0 and ; in addition, the presence of infinity in the collocation points does not cause a failure in the substitution.

Hence, upon substituting these points (35) into (30), one obtains

The matrix form of the obtained system (37) is written farther as where

5.1. The RRC Scheme

From (23), we know that the first scheme RRC gives us a derivative of RC functions from the order; thus, equation (38) takes the form where

Hence, from (38) and (40), one obtains the fundamental matrix equation for (30) as

Also, we obtain the matrix forms corresponding to condition (31) as follows: setting in (33), we get the fundamental matrix form corresponding to the condition (31):

5.2. The IRC Scheme

We studied the improved and regular differentiating RC functions in the preceding section. Now, we deduce the fundamental matrix relation by the IRC scheme.

Substituting relation (26) into (40), we get where

Hence, from (38) and (44), the fundamental matrix equation for (30) is obtained as

Similarly, the matrix form corresponding to the condition (31) using (26) is obtained as for so that .

6. Method of Solution

Due to the collocation method, the regular (42) and the improved (46) fundamental matrix equations for the proposed problem (30) correspond to a system of algebraic equations with equations for the unknown RC coefficients .

One writes matrix equations (42) and (46) compactly as or in the augmented form as

Equations (43) and (47) obtain the matrix form for the condition (31); also, they are written compactly as so that and for RRC are defined by while and for IRC are defined by

Hence, the approximate solution of (30) under the condition (31) may be obtained by exchanging the rows of matrices (50) by the last (or first) rows of the matrix (49), then getting the required augmented matrix as

If rank is equal to rank , then the algebraic system has a solution, and if the two ranks are equal to , then the solution is unique, the inverse matrix method is used here to solve the system, and one may write the matrix equation (49) as

Therefore, the RC coefficients , are uniquely determined.

7. Test Examples

In the present section, three numerical test examples are given to explain the applicability of the proposed two techniques. Using own codes written in MATHEMATICA 10.0. package, the numerical results and figures are presented, as shown in the illustrative comparison tables.

The absolute error is given to compare the efficiency of the proposed schemes, given by , and evaluated at selected points for some . The error norms and calculated in an interval are given by

for the step size along the interval . All numeric calculations are carried out on a regular machine Intel(R) Core(TM) i7 CPU, 3.2 GHz.

Example 1. Let us assume the following fourth-order boundary value problem with when .
We have

Thus, for , the numeric collocation points according to (35) are

The fundamental matrix equation of problem using RRC is while the fundamental matrix equation of problem using IRC is where , and are matrices with a given size of , for this example at ,

Then, we obtain the augmented matrix (49) with respect to RRC as

Then, we solve the equation (49) to find the RC coefficients in the matrix form: while the augmented matrix (49) with respect to IRC as

Similarly, we compute the RC coefficients in the matrix form:

Substituting by the RC coefficients of the RRC scheme in Equation (32), we get an approximate solution of this problem, while, if we use the RC coefficients of the IRC scheme in Equation (32), we get the exact solution . We note that the RRC scheme gives the exact solution at while the IRC scheme gives the exact solution at , which shows that the IRC scheme obtains better accuracy than the RRC due to the truncation found in RRC (18). In addition, the CPU time used for the RRC and IRC schemes at different is listed in Table 1. Table 1 shows that the time used for calculating the solution by RRC is less than the time for IRC, because of the added truncated terms in coding. Therefore, the time for getting the exact solution by IRC (better accuracy at ) is close to that given by RRC (at ). The following examples show the comparison of the two schemes strongly with numeric table’s results.

Example 2. Consider Laguerre-eigen problem [26] has the exact eigen solutions where are the well-known Laguerre polynomials of degree ; Equation (66) is a boundary value problem, so the conditions are “natural” at both endpoints where and when . By applying the proposed two schemes as in the previous example, we obtain the solution , which is the exact eigen solution when . For , the RRC scheme gives the exact eigen solution, while the IRC is satisfied with to find this solution. In Tables 2 and 3, the resulting values for and 30 using the present two schemes together with the exact values of the solution are tabulated with different values of . The error reduces when the series increased. The numeric values of the error norms and in interval is given in Tables 4 and 5 with different values of . Additionally, the absolute errors for the two schemes when are plotted in Figure 1. In Table 6, the CPU time used for RRC and IRC schemes at different shows that IRC needs time more than RRC because of the calculation of the added terms of .

Example 3. Consider Whittaker’s equation eigen problem [26] of the form where is the eigen value; it represents a special case of Whittaker’s equation. The exact given solution of (68) is where is an integer and is the well-known associated Laguerre polynomials of order and degree . If (68) gives a special case or linear Whittaker’s equation as The exact solution is with the boundary conditions and when .

By the same procedure, the RC collocation method using the proposed two schemes is applied to solve (69) with the subjected boundary conditions. In Table 7, the numerical result for and 30 using the proposed schemes is compared with the exact values of . The computation of and on interval is given in Table 8. In Table 9, the CPU time used for RRC and IRC schemes at different shows that IRC needs time more than RRC.

8. Conclusion

A rational Chebyshev (RC) spectral collocation technique is considered in this paper to solve high-order ordinary differential equations (ODEs) defined on a semi-infinite domain using the proposed two schemes. Two definitions of the derivative of the RC functions are introduced, namely, the regular and improved definitions. Due to the two definitions, two schemes are presented for solving the proposed ODEs with variable coefficients in the semi-infinite interval. According to the convergence of the RC functions at the infinity, the proposed technique deals with the boundary value problem which is defined on semi-infinite domains easily. Furthermore, an intriguing advantage of this approach is the ability to find the analytical exact solutions if the equation has a solution in a rational function form. To demonstrate the applicability of the proposed approach, three illustrative examples are given. The calculated numerical values and comparisons proved that the improved scheme is better with more calculation than the regular scheme which is based on the truncation in the definition. The method may extend to the case of nonlinear DEs with variable coefficients, which the authors are investigating.

Data Availability

The authors declare that there is no data associated with this research.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors carried out the proofs and conceived of the study. All authors read and approved of the final form of the manuscript.


The researchers would like to thank the Deanship of Scientific Research, Qassim University, for funding the publication of this project.