Second Kind Shifted Chebyshev Polynomials for Solving the Model Nonlinear ODEs

In this paper, we build the integral collocation method by using the second shifted Chebyshev polynomials. The numerical method solving the model non-linear such as Riccati differential equation, Logistic differential equation and Multi-order ODEs. The properties of shifted Chebyshev polynomials of the second kind are presented. The finite difference method is used to solve this system of equations. Several numerical examples are provided to confirm the reliability and effectiveness of the proposed method.


Introduction
In recent years, Chebyshev polynomials (family of orthogonal polynomials on the interval [−1, 1]) have become increasingly important in numerical analysis, from both theoretical and practical points of view.They have strong links with Fourier and Laurent series, with minimality properties in approximation theory and with discrete and continuous orthogonality in function spaces [1].These links have led to important applications, especially in spectral methods for ordinary and partial differential equations.There are four kinds of Chebyshev polynomials as in [2].The majority of books dealing with Chebyshev polynomials, contain mainly results of Chebyshev polynomials of all kinds ( ) ( ) ( ) , , and ( ) n W x and their numerous uses in different applications and research papers dealing with some types of these polynomials ( [3]- [8]) and other publications as ( [9] [10] [11] [12] [13]).However, there are only a limited researches of literature on shifted Chebyshev polynomials of the second kind   ( ) n U x may be generated by using the recurrence relations The analytical form of the Chebyshev polynomials of the second kind of degree n is given by: where π 2

Shifted Chebyshev Polynomials of the Second Kind ([2] [14])
In order to use these polynomials in Section 2.1 on the interval ( ) ( ) also there are important relation between the shifted and second kind Chebyshev polynomials as follows: ( ) these polynomials are orthogonal on the support interval [ ] 0,1 as the following inner product: with start values ( ) The analytical form of the shifted Chebyshev polynomials of the second kind The function which may be appear in solution of the model problem can be written as series of ( ) g x be a square integrable in [ ] 0,1 it can be expressed in terms of the shifted Chebyshev polynomials of the second kind as follows: where the coefficients , 0,1, i a i =  are given by: ( ) In practice, only the first ( ) terms of shifted Chebyshev polynomials of the second kind are considered in the approximate case.Then we have: Using the parctice shifted Chebyshev polynomials of the second kind to constraction the integral collocation method to give the N-th derivative of the unknown function ( ) u x as the following [15]: Using the integration we can obtain the lower-order derivatives and the function itself as follows ( ) ( ) from ( 4) and ( 9) we have We now collocate Equatuions ( 10)-( 14) at ( ) where Ω Ω Ω  are integrated matrices.

Integral Collocation Method for Solving Riccati, Logistic and Multi-Order Nonlinear ODEs [15]
In this section, we introduce the integral collocation method using shifted Chebyshev polynomials of the second kind for solving the Riccati, Logistic and multi-order nonlinear ODEs.

Model 1: Riccati Differential Equation [15]
( ) ( ) we also assume an initial condition The exact solution to this problem at 0 0 u = is ( ) The procedure of the implementation is given by the following steps: 1) Approximate the function ( ) u x using formula ( 9)-( 14) with 5 m = , as follows ( ) ( ) ( ) w x is defined in (14) as Then the Riccati differential Equation ( 16) is transformed to the following approximated form ( ) ( ) ( ) We now collocate Equation ( 19) at ( ) For suitable collocation points we use the roots of shifted Chebyshev polynomial 2) Also, by substituting from the initial condition (17) in (18) we can obtain ( ) an equation which gives the value of the constant 1 c as follows Equations ( 20) and (21) represent a system of non-linear algebraic equations which contains seven equations for the unknowns , 0,1, 2,3, 4,5 n a n = and 1 c . 3) Solve the resulting system using the Newton iteration method to obtain the unknowns , 0,1, 2,3, 4,5 Therefore, from Formula (19) we can obtain the approximate solution in the form ( ) ( ) ( ) The numerical results of the proposed problem ( 16) are given in Figure 1 with 5 m = in the interval [ ] 0,1 at 0 0 u = .From this Figure 1, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.We also assume an initial condition

Model 2: Logistic Differential Equation
The exact solution to this problem is given by ( ) ( ) The procedure of implementation is given by the following steps: 1) Approximate the function ( ) u x using formula ( 9)-( 14) with 5 m = Then the Logistic differential Equation ( 22) is transformed to the following approximated form ( ) ( ) ( ) ( ) ( ) We now collocate Equation (24) at ( ) x p = as ( ) ( ) ( ) ( ) ( ) For suitable collocation points we use roots of shifted Chebyshev polynomial ( ) 2) Also, by substituting from the initial condition (23) in (18) with 0 0.85 u = we can obtain ( ) an equation which gives the value of the constant 1 0.85 c = .Equation (25) represents a system of non-linear algebraic equations which contains six equations for the unknowns , 0, n a n = .3) Solve the resulting system using the Newton iteration method to obtain the unknowns , 0, Therefore, from Formula (18) we can obtain the approximate solution in the form ( ) ( ) ( ) The numerical results of the proposed problem ( 22) is given in Figure 2 with From this Figure 2, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.

Model Multi-Order Nonlinear ODEs ([16] [17])
Consider the following initial value problem ( [16] [17]) ( ) ( ) , n n w x w x and ( ) ( ) w x are defined as follows Then the multi-order ODE (27) can be written in the following approximated form ( ) ( ) ( ) We now collocate Equation (30) at ( ) For suitable collocation points we use roots of shifted Chebyshev polynomial ( ) 2) Also, by substituting from the initial conditions (28) in (29) we can obtain 3 k = of equations which give the values of the constants 1 2 , c c and 3 c .
The Equations ( 31) and (32) construct system of non-linear algebraic equations which contains seven equations for the unknowns 3) Solve the resulting system using the Newton iteration method to obtain the unknowns , 0,1, 2,3 n b n = .
Therefore, using the formula (29) we can find the required approximate solution in the following form:

Conclusion
In this paper, the Chebyshev polynomials of kind has been successfully applied to study the model nonlinear ODEs.The results show that Chebyshev polynomials of the second kind is an efficient and easy-to-use technique for finding exact and approximate solutions for nonlinear ordinary differential equations.The obtained approximate solutions using the suggested method is in excellent agreement with the exact solution and show that these approaches can be solved the problem effectively and illustrates the validity and the great potential of the proposed technique.
the second kind[2] are orthogonal polynomials of degree n in x defined on the [−1, 1] called shifted Chebyshev polynomials of the second kind means that the shifted Chebyshev polynomials of the second kind defined as ([2][14]): may be generated by using the recurrence relations

Figure 1 .
Figure 1.The behavior of the approximate solution and exact solution with m = 5.

Figure 2 .
Figure 2. The behavior of the approximate solution and exact solution with m = 5.
exact solution of the proposed problem (27).The numerical results of the proposed problem (27) are given in Figure3with 5 m = in the interval [ ] 0,10 .From this Figure3, since the obtained numerical solutions are in excellent agreement with the exact solution, so, we can conclude that the proposed technique is well for solving such class of ODEs.

Figure 3 .
Figure 3.The behavior of the approximate solution with m = 5.