﻿Sinc-Collocation Method for Solving Linear and Nonlinear System of Second-Order Boundary Value Problems

Applied Mathematics
Vol. 3  No. 11 (2012) , Article ID: 24504 , 7 pages DOI:10.4236/am.2012.311225

Sinc-Collocation Method for Solving Linear and Nonlinear System of Second-Order Boundary Value Problems

Mohamed El-Gamel

Department of Mathematical Sciences, Faculty of Engineering, Mansoura University, Mansoura, Egypt

Email: gamel_eg@yahoo.com

Received July 23, 2012; revised October 8, 2012; accepted October 15, 2012

Keywords: Sinc Function; Collocation Method; System; Numerical Solution

ABSTRACT

Sinc methods are now recognized as an efficient numerical method for problems whose solutions may have singularities, or infinite domains, or boundary layers. This work deals with the sinc-collocation method for solving linear and nonlinear system of second order differential equation. The method is then tested on linear and nonlinear examples and a comparison with B-spline method is made. It is shown that the sinc-collocation method yields better results.

1. Introduction

Numerous problems in physics, chemistry, biology and engineering science are modelled mathematically by systems of ordinary differential equations, e.g. series circuits, mechanical systems with several springs attached in series lead to a system of differential equations (for example see [1,2]). However, many classical numerical methods used with second-order initial value problems cannot be applied to second-order boundary value problems (BVPs).

Most realistic systems of ordinary differential equations do not have exact analytic solutions, so approximation and numerical techniques must be used. There are many publications dealing with the linear system of second-order boundary value problems. They introduced various numerical methods. For instance, a finite difference method has been proposed in recent works [3-8]. For a nonlinear system of secondorder BVPs, there are few valid methods to obtain numerical solutions. Geng et al. have studied the numerical solution of a nonlinear system of second-order boundary value problems in the reproducing kernel space [9]. Lu considered the variational iteration method to solve a nonlinear system of second-order boundary value problems [10]. Recently, Bataineh et al. [11] represented modified homotopy method for solving systems of second-order boundary value problems. Sinc-collocation method was applied to solve nonlinear systems of second order boundary value problems in [12].

In this paper, we discuss the use of sinc-collocation method for solving a class of linear and non-linear system of differential equations

(1)

subject the boundary conditions

(2)

where, , , , and

, , , and, for, are analytic functions. It will always be assumed that (1) possesses a unique solution.

Numerical examples including regular, singular as well as singularly perturbed problems are considered. On the basis of these examples, the results reveal that the method is very effective and convenient.

The paper is organized into five sections. Section 2 contains notation, definitions and some results of sinc function theory. In Section 3, the sinc-collocation method is developed for linear second-order system of differential equation with homogeneous boundary conditions. The method is developed for nonlinear second-order system of differential equation in Section 4. Some numerical examples are presented in Section 5. Finally, Section 6 provides conclusions of the study.

2. Sinc Function

In recent years, a lot of attention has been devoted to the study of the sinc method to investigate various scientific models. The efficiency of the method has been formally proved by many researchers [13-22].

A general review of sinc function approximation is given in [23,24]. Hence, only properties of the sinc function that are used in the sequel.

If is defined on the real line, then for the Whittaker cardinal expansion of f is given by:

(3)

where, , and the mesh size is given by

(4)

where is suitably chosen and depends on the asymptotic behavior of. The n-th derivative of the function f at the sampling points can be approximated using a finite number of terms as:

where

In particular,

(5)

(6)

and

(7)

We note that

The interpolation formula for over takes the form

(8)

where the basis functions on are then given by

and the transformation function

(9)

transforms to the infinite range The interpolation points are then given by:

The n-th derivative of the function at points can be approximated using a finite number of terms as

(10)

Setting

and noting that

and

which will be used later.

3. System of Linear Second Order Equations

Consider a linear, system of linear second order equations of the form

(11)

We assume that and the solutions of (11) and (2), is approximated by the finite expansion of Sinc basis functions

(12)

and

(13)

where is the function for some fixed step size h. If we replace each term of (11) with its corresponding approximation given by the right-hand side of (10) and (8) we have

(14)

Substituting in (14) and applying the collocation to it, we eventually obtain the following theorem.

Theorem: If the assumed approximate solution of problem (11) and (2) is (12) and (13) , then the discrete sinc-collocation system for the determination of the unknown coefficients is given by

(15)

where

and

Now, since, and, we write the above equation in the form

(16)

We now rewrite these equations in matrix form. The system in (16) takes the matrix form

(17)

where

Now we have a linear system of equations of the unknown coefficients. We can obtain the coefficients of the approximate solution by solving this linear system. The system (17) may be easily solved by a variety of methods. In this paper we used the Q-R method. The solution c gives the coefficients in the approximate sinc-collocation solutions and of and.

4. System of Non-Linear Second Order Equations

Consider a nonlinear, system second-order equations of the form

(18)

where and are nonlinear functions of and is an analytic function and may be a polynomial or a rational function, or exponential. Due to the large number of different possibilities, our work will be focused mainly on the following forms

1)

2)

3)

or any analytic function of which has a power series expansion. We limit our study to the case

where is an integer, or a fraction.

We consider next applying of the sinc-collocation method to solve problem (18) and (2).

Lemma: The following relation holds

(19)

where N and h are now dependent on both and.

Replacing the terms of (18) with the appropriate representation defined in (8), (10) and (19) and applying the collocation to it, we eventually obtain the following theorem.

Theorem: If the assumed approximate solution of problem (18) and (2) is (12) and (13), then the discrete sinc-collocation system for the determination of the unknown coefficients is given by

(20)

Let be the -vector with j-th component given by. In this notation the system in (20) takes the matrix form

(21)

where

Now we have a nonlinear system of equations in the unknown coefficients. We can obtain the coefficients in the approximate solution by solving this nonlinear system by Newton’s method.

Starting from an initial estimate, the corrections are made using

Here, is the current iterate, and is the new iterate. A common numerical practice is to stop the Newton iteration whenever the distance between two iterates is less than a given tolerance, i.e. when

where the Euclidean norm is used. The solution gives the coefficients in the approximate sinc-collocation solution of.

5. Numerical Examples

In this section, some numerical examples are studied to demonstrate the accuracy of the present method. The results obtained by the method are compared with the exact solution of each example and are found to be in good agreement with each other.Comparison between sinc-collocation and other method shall be presented.

All computations were carried out using Matlab on a personal computer with a machine precision of. In all cases, is taken to be. The selection of a larger yields more accuracy, but at the expense of a lengthier computation. We report absolute error which is defined as

Example 1: [3,11] consider the linear system of second order boundary value problems

where

subject to the boundary conditions

whose exact solutions are

Maximum absolute errors for and are tabulated in Table 1 for the sinc-collocation method.

Maximum absolute error are tabulated in Table 2 for sinc-collocation together with the analogous results of N. Caglar and H. Caglar [3].

Example 2: [11] Now we turn to a nonlinear problem

where

subject to the boundary conditions

whose exact solutions are

The computational results are summarized in Table 3.

Example 3: Now we turn to a singular problem,

where

subject to the boundary conditions

whose exact solutions are

Table 1. Maximum absolute error for example 1.

Table 2. Comparison between maximum absolute error for example 1.

Table 3. Maximum absolute error for example 2.

The computational results are summarized in Table 4.

Example 4: Another example is also a singular problem

where

subject to the boundary conditions

whose exact solutions are

Maximum absolute errors for and are tabulated in Table 5 for the sinc-collocation method.

Example 5: Our final example is the singularly perturbed problem

Table 4. Maximum absolute error for example 3.

Table 5. Maximum absolute error for example 4.

where

subject to the boundary conditions

whose exact solutions are

Table 6 the maximum absolute errors obtained by using the sinc-collocation method for N = 40 and different.

6. Conclusions

This paper described an efficient method for solving the system of second-order boundary value problems. Our approach was based on the sinc-collocation method. Properties of the sinc-collocation method are utilized to reduce the computation of this problem to some linear or nonlinear algebraic equations. The method is computationally attractive and applications are demonstrated through illustrative examples. Numerical examples including regular, singular as well as singularly perturbed problems are presented. As expected, the accuracy increases as the number of terms N in the sinc expansion increases. The obtained results showed that this approach can solve the problem effectively.

The sinc-collocation method is a simple method with high accuracy for solving a large variety of linear and nonlinear system of differential equations. So it may be

Table 6. Maximum absolute error for example 5.

easily applied by researchers and engineers familiar with the sinc function. Extension of the method for solving systems of partial differential equations offers an excellent opportunity for future research.

REFERENCES

1. K.W. Tomantschger, “Series Solutions of Coupled Differential Equations with One Regular Singular Point,” Journal of Computational and Applied Mathematics, Vol. 140, No. 1-2, 2002, pp. 773-783. doi:10.1016/S0377-0427(01)00598-2
2. C. Wafo Soh and F. M. Mahomed, “Linearization Criteria for a System of Second-Order Ordinary Differential Equations,” International Journal of Non-Linear Mechanics, Vol. 36, No. 4, 2001, pp. 671-677. doi:10.1016/S0020-7462(00)00032-9
3. N. Caglar and H. Caglar, “B-Spline Method for Solving Linear System of Second-Order Boundary Value Problems,” Computers & Mathematics with Applications, Vol. 57, No. 5, 2009, pp. 757-762. doi:10.1016/j.camwa.2008.09.033
4. S. H. Chen, J. Hu, L. Chen and C. P. Wang, “Existence Results for n-Point Boundary Value Problem of Second Order Ordinary Differential Equations,” Journal of Computational and Applied Mathematics, Vol. 180, No. 2, 2005, pp. 425-432. doi:10.1016/j.cam.2004.11.010
5. X. Y. Cheng and C. K. Zhong, “Existence of Positive Solutions for a Second-Order Ordinary Differential System,” Journal of Mathematical Analysis and Applications, Vol. 312, No. 1, 2005, pp. 14-23. doi:10.1016/j.jmaa.2005.03.016
6. A. Lomtatidze and L. Malaguti, “On a Two-Point Boundary Value Problem for the Second-Order Ordinary Differential Equations with Singularities,” Nonlinear Analysis: Theory, Methods & Applications, Vol. 52, No. 6, 2003, pp. 1553-1567. doi:10.1016/S0362-546X(01)00148-1
7. H. Thompson and C. Tisdell, “Boundary Value Problems for Systems of Difference Equations Associated with Systems of Second-Order Ordinary Differential Equations,” Applied Mathematics Letters, Vol. 15, No. 6, 2002, pp. 761-766. doi:10.1016/S0893-9659(02)00039-3
8. H. Thompson and C. Tisdell, “The Nonexistence of Spurious Solutions to Discrete, Two-Point Boundary Value Problems,” Applied Mathematics Letters, Vol. 16, No. 1, 2003, pp. 79-84. doi:10.1016/S0893-9659(02)00147-7
9. F. Z. Geng and M. G. Cui, “Solving a Nonlinear System of Second-Order Boundary Value Problems,” Journal of Mathematical Analysis and Applications, Vol. 327, No. 2, 2007, pp. 1167-1181. doi:10.1016/j.jmaa.2006.05.011
10. J. F. Lu, “Variational Iteration Method for Solving a Nonlinear System of Second-Order Boundary Value Problems,” Computers & Mathematics with Applications, Vol. 54, No. 7-8, 2007, pp. 1133-1138. doi:10.1016/j.camwa.2006.12.060
11. A. Bataineh, M. S. M. Noorani and I. Hashim, “Modified Homotopy Analysis Method for Solving Systems of Second-Order BVPs,” Communications in Nonlinear Science and Numerical Simulation, Vol. 14, No. 2, 2009, pp. 430- 442. doi:10.1016/j.cnsns.2007.09.012
12. M. Dehghan and A. Saadatmandi, “The Numerical Solution of a Nonlinear System of Second-Order Boundary Value Problems Using the Sinc-Collocation Method,” Mathematical and Computer Modelling, Vol. 46, No. 11- 12, 2007, pp. 1434-1441. doi:10.1016/j.mcm.2007.02.002
13. B. Bialecki, “Sinc-Collocation Methods for Two-Point Boundary Value Problems,” IMA Journal of Numerical Analysis, Vol. 11, No. 3, 1991, pp. 357-375. doi:10.1093/imanum/11.3.357
14. M. El-Gamel and A. I. Zayed, “Sinc-Galerkin Method for Solving Nonlinear Boundary-Value Problems,” Computers & Mathematics with Applications, Vol. 48, No. 9, 2004, pp. 1285-1298. doi:10.1016/j.camwa.2004.10.021
15. M. El-Gamel, J. Cannon and A. Zayed, “Sinc-Galerkin Method for Solving Linear Sixth Order Boundary-Value Problems,” Mathematics of Computation, Vol. 73, No. 247, 2004, pp. 1325-1343.
16. M. El-Gamel, S. H. Behiry and H. Hashish, “Numerical Method for the Solution of Special Nonlinear FourthOrder Boundary Value Problems,” Applied Mathematics and Computation, Vol. 145, No. 2-3, 2003, pp. 717-734. doi:10.1016/S0096-3003(03)00269-8
17. M. El-Gamel and J. Cannon, “On the Solution of Second Order Singularly-Perturbed Boundary Value Problem by the Sinc-Galerkin Method,” Zeitschrift für Angewandte Mathematik und Physik (ZAMP), Vol. 56, No. 1, 2005, pp. 45-58. doi:10.1007/s00033-004-3002-6
18. A. Mohsen and M. El-Gamel, “A Sinc-Collocation Method for the Linear Fredholm Integro-Differential Equations,” Zeitschrift für Angewandte Mathematik und Physik (ZAMP), Vol. 58, No. 3, 2007, pp. 380-390. doi:10.1007/s00033-006-5124-5
19. A. Mohsen and M. El-Gamel, “On the Galerkin and Collocation Methods for Two-Point Boundary Value Problems Using Sinc Bases,” Computers & Mathematics with Applications, Vol. 56, No. 4, 2008, pp. 930-941. doi:10.1016/j.camwa.2008.01.023
20. A. Mohsen and M. El-Gamel, “On the Numerical Solution of Linear and Nonlinear Volterra Integral and IntegroDifferential Equations,” Applied Mathematics and Computation, Vol. 217, No. 7, 2010, pp. 3330-3337. doi:10.1016/j.amc.2010.08.065
21. R. Smith, G. Bogar, K. Bowers and J. Lund, “The SincGalerkin Method for Fourth-Order Differential Equations,” SIAM Journal on Numerical Analysis, Vol. 28, No. 3, 1991, pp. 760-788. doi:10.1137/0728041
22. G. Y. Yin, “Sinc-Collocation Method with Orthogonalization for Singular Poisson-Like Problem,” Mathematics of Computation, Vol. 62, No. 205, 1994, pp. 21-40. doi:10.1090/S0025-5718-1994-1203738-7
23. J. Lund and K. L. Bowers, “Sinc Methods for Quadrature and Differential Equations,” Society for Industry and Applied Mathematics (SIAM), Philadelphia, 1992. doi:10.1137/1.9781611971637
24. F. Stenger, “Numerical Methods Based on Sinc and Analytic Functions,” Springer, New York, 1993. doi:10.1007/978-1-4612-2706-9