Collocation Method via Jacobi Polynomials for Solving Nonlinear Ordinary Differential Equations

We extend a collocation method for solving a nonlinear ordinary differential equation ODE via Jacobi polynomials. To date, researchers usually use Chebyshev or Legendre collocation method for solving problems in chemistry, physics, and so forth, see the works of Doha and Bhrawy 2006, Guo 2000, and Guo et al. 2002 . Choosing the optimal polynomial for solving every ODEs problem depends on many factors, for example, smoothing continuously and other properties of the solutions. In this paper, we show intuitionally that in some problems choosing other members of Jacobi polynomials gives better result compared to Chebyshev or Legendre polynomials.


Introduction
The Jacobi polynomials with respect to parameters α > −1, β > −1 see, e.g., 1, 2 are sequences of polynomials P α,β n x ; n 0, 1, 2, . . .satisfying the following relation  These polynomials are eigenfunctions of the following singular Sturm-Liouville equation: A consequence of this is that spectral accuracy can be achieved for expansions in Jacobi polynomials so that 1.4 The Jacobi polynomials can be obtained from Rodrigue's formula as Furthermore, we have that The Jacobi polynomials are normalized such that An important consequence of the symmetry of weight function w x and the orthogonality of Jacobi polynomial is the symmetric relation     x .
x .that is, the Jacobi polynomials are even or odd depending on the order of the polynomial.In this form the polynomials may be generated using the starting form which is obtained from Rodrigue's formula as follows: Following the two seminal papers of Doha 3, 4 let f x be an infinitely differentiable function defined on −1, 1 ; then we can write 1.12 and, for the qth derivative of f x , Then, a q n 2 −q ∞ i 0 n i q λ − 1 q C n i,n α q, β q, α, β a n i q , ∀n ≥ 0, q ≥ 1, 1.14 where For the proof of the above, see 3 .The formula for the expansion coefficients of a generalorder derivative of an infinitely differentiable function in terms of those of the function is available for expansions in ultraspherical and Jacobi polynomials in Doha 5 .Another interesting formula is where

1.18
For the proof see 4 .Chebyshev, Legendre, and ultraspherical polynomials are particular cases of the Jacobi polynomials.These polynomials have been used both in the solution International Journal of Mathematics and Mathematical Sciences of boundary value problems 6 and in computational fluid dynamics 7, 8 .For the ultraspherical coefficients of the moments of a general-order derivative of an infinitely differentiable function, see 5 .Collocation method is a kind of spectral method that uses the delta function as a test function.Test functions have an important role because these functions are applied to obtain the minimum value for residual by using inner product.Attention to this point is very important since Tau and collocation methods do not give the best approximation solution of ODEs.Due to adding the boundary conditions, we are forced to eliminate several conditions of orthogonality properties, and this decreases the accuracy of the ODEs solution.
Previous explanations do not mean that the accuracy of the collocation method is less than that of the Galerkin method.Many studies on the collocation method have recently appeared in the literature on the numerical solution of boundary value problems 9-12 .The collocation solution is a piecewise continuous polynomial function.The error has been analyzed for different conditions by different authors.Frank and Ueberhuber 13 showed equivalence between solutions of collocation method and fixed points of Iterated Defect Correction IDeC method.They proved that IDeC method can be regarded as an efficient scheme for solving collocation equations.C ¸elik 14, 15 investigated the corrected collocation method for the approximate computation of eigenvalues of Sturm-Liouville and periodic Sturm-Liouville problems by a Truncated Chebyshev Series TCS .On the other hand, some results on Jacobi approximation were used for analyzing the p-version of finite element method, boundary element method, spectral method for axisymmetrical domain, and various rational spectral methods 16-21 .The rest of this paper is organized as follows.In Section 2, the method of solution is developed.In Section 3, we report our numerical findings and demonstrate the accuracy of the proposed scheme by considering a numerical example.Finally, a brief conclusion is presented in Section 4.

Jacobi Collocation Method
We consider the class of ODEs with the following form where s, l, m ≥ 0, P i x , Q i x , R i x , and F x are defined on the interval a ≤ x ≤ b.The most general boundary conditions are where the real coefficients γ i , α i , δ, and μ are appropriate constants and the solution is expressed in the form

2.4
Hence, there is no loss of generality in taking a −1 and b 1.To obtain such a solution, the Jacobi collocation method can be used.Collocation points can be taken as the roots of P α,β n x 22 .Until now we have not the explicit formulas that determine these roots.Therefore, we use the numerical method to determine these roots, for example, the Newton method.Substituting the Jacobi collocation points in 2.1 , the following expression can be obtained: x j i F x j , j 0, 1, 2, 3, . .., n.

2.5
The above system of equations constructs n 1 nonlinear equations.We impose the boundary conditions to these equations as follows:

2.6
Equations 2.5 and 2.6 construct the n 3 nonlinear equations and so we eliminate two equations from 2.1 .Thus we have n 1 equations and solve this system with Newton iteration method.

Illustrative Numerical Example
Consider the following nonlinear boundary value problem:

International Journal of Mathematics and Mathematical Sciences
The exact solution is y sin x .We solve it with the above-mentioned method for different values of α and β , and we show the results.We compute the numerical solution for 4 terms in Figure 1.
At first we obtain the roots of P α,β n x .For this goal we must determine the α and β parameters, and Table 1 shows the roots.
To apply the method we assume that the solution is in the following form: 3.3 By substituting y in 3.1 , we obtain a system of equations, and, with attention to Gaussian nodes, we must solve the set of nonlinear equations.The boundary conditions are imposed on the system of equations with elimination of two equations, and we add these boundary conditions to the system as follows: x j 0, j 0, 1, 2, 3.

3.4
The above system can be solved with Newton iteration method.We show our results for every parameter as in Table 2. Representation of the error by increasing the number of terms in our polynomial solution is shown in Table 3.It is observable that as the number of terms increases, the absolute error decreases.

Conclusion
In this paper, we show that this method is very accurate, even for small value of n i.e., n 4 .We see that Chebyshev polynomial the cases c 1 and c 2 and Legendre polynomial the cases a 1 and a 2 have not obtained the best approximation see the case α −0.4,β 0.6 .We observe that the error is symmetric when α β see the cases a 2 , b 2 , c 2 , e 2 , and f 2 .When we increase the number of terms in our computation, we monitor that the error decreases and our computed results are the best.
the certain conditions r 0, 1, 2, . . . .Any finite range a ≤ w ≤ b can be transformed to the basic range −1 ≤ x ≤ 1 with the change of variable w b − a x b a 2 .

Table 2 :
The solution of our example for different values of α and β.

Table 3 :
Representation of the error by increasing the number of terms in our polynomial solution.