Abstract
The application of the sincGalerkin method to an approximate solution of secondorder singular Dirichlettype boundary value problems were discussed in this study. The method is based on approximating functions and their derivatives by using the Whittaker cardinal function. The differential equation is reduced to a system of algebraic equations via new accurate explicit approximations of the inner products without any numerical integration which is needed to solve matrix system. This study shows that the sincGalerkin method is a very effective and powerful tool in solving such problems numerically. At the end of the paper, the method was tested on several examples with secondorder Dirichlettype boundary value problems.
Keywords:
sincGalerkin method; sinc basis functions; Dirichlettype boundary value problems; LU decomposition method1 Introduction
Sinc methods were introduced by Frank Stenger in [1] and expanded upon by him in [2]. Sinc functions were first analyzed in [3] and [4]. An extensive research of sinc methods for twopoint boundary value problems can be found in [5,6]. In [7,8], parabolic and hyperbolic problems were discussed in detail. Some kind of singular elliptic problems were solved in [9], and the symmetric sincGalerkin method was introduced in [10]. Sinc domain decomposition was presented in [1113] and [14]. Iterative methods for symmetric sincGalerkin systems were discussed in [15,16] and [17]. Sinc methods were discussed thoroughly in [18]. Applications of sinc methods can also be found in [19,20] and [21]. The article [22] summarizes the results obtained to date on sinc numerical methods of computation. In [14], a numerical solution of a Volterra integrodifferential equation by means of the sinc collocation method was considered. The paper [2] illustrates the application of a sincGalerkin method to an approximate solution of linear and nonlinear secondorder ordinary differential equations, and to an approximate solution of some linear elliptic and parabolic partial differential equations in the plane. The fully sincGalerkin method was developed for a family of complexvalued partial differential equations with timedependent boundary conditions [19]. Some novel procedures of using sinc methods to compute solutions to three types of medical problems were illustrated in [23], and sincbased algorithm was used to solve a nonlinear set of partial differential equations in [24]. A new sincGalerkin method was developed for approximating the solution of convection diffusion equations with mixed boundary conditions on halfinfinite intervals in [25]. The work which was presented in [26] deals with the sincGalerkin method for solving nonlinear fourthorder differential equations with homogeneous and nonhomogeneous boundary conditions. In [27], sinc methods were used to solve secondorder ordinary differential equations with homogeneous Dirichlettype boundary conditions.
2 Sinc functions preliminaries
Let C denote the set of all complex numbers, and for all
For
For various values of k, the sinc basis function
Figure 1. The basis functions
Figure 2. Central sinc basis function
If a function
is called the Whittaker cardinal expansion of f whenever this series converges. The infinite strip
In general, approximations can be constructed for infinite, semiinfinite and finite intervals. Define the function
which is a conformal mapping from
This is shown in Figure 3.
Figure 3. The relationship between the eyeshaped domain
For the sincGalerkin method, the basis functions are derived from the composite translated sinc functions
for
the evenly spaced nodes
A list of conformal mappings may be found in Table 1[6].
Figure 4. Three adjacent members
Table 1. Conformal mappings and nodes for some subintervals ofR
Definition 2.1 Let
and
Definition 2.2 Let
where
and those on the boundary of
The proof of following theorems can be found in [2].
Theorem 2.1Let Γ be
where
For the sincGalerkin method, the infinite quadrature rule must be truncated to a finite sum. The following theorem indicates the conditions under which an exponential convergence results.
Theorem 2.2If there exist positive constantsα, βandCsuch that
then the error bound for the quadrature rule (2.14) is
The infinite sum in (2.14) is truncated with the use of (2.16) to arrive at the inequality (2.17). Making the selections
where
We used Theorems 2.1 and 2.2 to approximate the integrals that arise in the formulation of the discrete systems corresponding to a secondorder boundary value problem.
Theorem 2.3Letϕbe a conformal onetoone map of the simply connected domain
3 The sincGalerkin method for singular Dirichlettype boundary value problems
Consider the following problem:
with Dirichlettype boundary condition
where P, Q and F are analytic on D. We consider sinc approximation by the formula
The unknown coefficients
Let
where w is the weight function. For the secondorder problems, it is convenient to take [2].
For Eq. (3.1), we use the notations (2.21)(2.23) together with the inner product that, given (3.5) [2], showed to get the following approximation formulas:
where
Using (3.5), (3.8)(3.11), we obtain a linear system of equations for
The
Let
With these notations, the discrete system of equations in (3.5) takes the form:
Theorem 3.1Letcbe anmvector whosejth component is
Now we have a linear system of
4 Examples
Three examples were given in order to illustrate the performance of the sincGalerkin
method to solve a singular Dirichlettype boundary value problem in this section.
The discrete sinc system defined by (3.18) was used to compute the coefficients
Figure 5. Approximation to the exact solution: the red colored curve displays the exact solution and the green one is the approximate solution of Eq. (4.1).
Figure 6. Approximation to the exact solution: the red colored curve displays the exact solution and the green one is the approximate solution of Eq. (4.2).
Figure 7. Approximation to the exact solution: the red colored curve displays the exact solution and the green one is the approximate solution of Eq. (4.3).
Table 2. The numerical results for the approximate solutions obtained by sincGalerkin in comparison
with the exact solutions of Eq. (4.1) for
Table 3. The numerical results for the approximate solutions obtained by sincGalerkin in comparison
with the exact solutions of Eq. (4.2) for
Table 4. The numerical results for the approximate solutions obtained by sincGalerkin in comparison
with the exact solutions of Eq. (4.3) for
Example 4.1 Consider the following singular Dirichlettype boundary value problem on the interval
The exact solution of (4.1) is
We choose the weight function according to [2],
Example 4.2 Let us have the following form of a singular Dirichlettype boundary value problem
on the interval
The problem has an exact solution like
where
Example 4.3 The following problem is given on the interval
where the exact solution of (4.3) is
In this case,
5 Conclusion
The sincGalerkin method was employed to find the solutions of secondorder Dirichlettype boundary value problems on some closed real interval. The main purpose was to find the solution of boundary value problems which arise from the singular problems. The examples show that the accuracy improves with increasing number of sinc grid points N. We have also developed a very efficient and rapid algorithm to solve secondorder Dirichlettype BVPs with the sincGalerkin method on the Maple computer algebra system. All of the above computations and graphical representations were prepared by using Maple.
We give the Maple code in the Appendix section.
Appendix: Maple code which we developed for the sincGalerkin approximation
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
AS proposed main idea of the solution schema by using Sinc Method for linear BVPs. He developed computer algorithm and worked on theoretical aspect of problem. MK searched the materials about study and compared with other techniques, contributed with his experience on Nonlinear Approximation methods.
References

Stenger, F: Approximations via Whittaker’s cardinal function. J. Approx. Theory. 17, 222–240 (1976). Publisher Full Text

Stenger, F: A sincGalerkin method of solution of boundary value problems. Math. Comput.. 33, 85–109 (1979)

Whittaker, ET: On the functions which are represented by the expansions of the interpolation theory. Proc. R. Soc. Edinb.. 35, 181–194 (1915)

Whittaker, JM: Interpolation Function Theory, Cambridge University Press, London (1935)

Lund, J: Symmetrization of the sincGalerkin method for boundary value problems. Math. Comput.. 47, 571–588 (1986). Publisher Full Text

Lund, J, Bowers, KL: Sinc Methods for Quadrature and Differential Equations, SIAM, Philadelphia (1992)

Lewis, DL, Lund, J, Bowers, KL: The spacetime sincGalerkin method for parabolic problems. Int. J. Numer. Methods Eng.. 24, 1629–1644 (1987). Publisher Full Text

McArthur, KM, Bowers, KL, Lund, J: Numerical implementation of the sincGalerkin method for secondorder hyperbolic equations. Numer. Methods Partial Differ. Equ.. 3, 169–185 (1987). Publisher Full Text

Bowers, KL, Lund, J: Numerical solution of singular Poisson problems via the sincGalerkin method. SIAM J. Numer. Anal.. 24(1), 36–51 (1987). Publisher Full Text

Lund, J, Bowers, KL, McArthur, KM: Symmetrization of the sincGalerkin method with block techniques for elliptic equations. IMA J. Numer. Anal.. 9, 29–46 (1989). Publisher Full Text

Lybeck, NJ: Sinc domain decomposition methods for elliptic problems. PhD thesis, Montana State University, Bozeman, Montana (1994)

Lybeck, NJ, Bowers, KL: Domain decomposition in conjunction with sinc methods for Poisson’s equation. Numer. Methods Partial Differ. Equ.. 12, 461–487 (1996). Publisher Full Text

Morlet, AC, Lybeck, NJ, Bowers, KL: The Schwarz alternating sinc domain decomposition method. Appl. Numer. Math.. 25, 461–483 (1997). Publisher Full Text

Morlet, AC, Lybeck, NJ, Bowers, KL: Convergence of the sinc overlapping domain decomposition method. Appl. Math. Comput.. 98, 209–227 (1999). Publisher Full Text

Alonso, N, Bowers, KL: An alternatingdirection sincGalerkin method for elliptic problems. J. Complex.. 25, 237–252 (2009). Publisher Full Text

Ng, M: Fast iterative methods for symmetric sincGalerkin systems. IMA J. Numer. Anal.. 19, 357–373 (1999). Publisher Full Text

Ng, M, Bai, Z: A hybrid preconditioner of banded matrix approximation and alternatingdirection implicit iteration for symmetric sincGalerkin linear systems. Linear Algebra Appl.. 366, 317–335 (2003)

Stenger, F: Numerical Methods Based on Sinc and Analytic Functions, Springer, New York (1993)

Koonprasert, S: The sincGalerkin method for problems in oceanography. PhD thesis, Montana State University, Bozeman, Montana (2003)

McArthur, KM, Bowers, KL, Lund, J: The sinc method in multiple space dimensions: model problems. Numer. Math.. 56, 789–816 (1990)

Stenger, F: Numerical methods based on Whittaker cardinal, or sinc functions. SIAM Rev.. 23, 165–224 (1981). Publisher Full Text

Stenger, F: Summary of sinc numerical methods. J. Comput. Appl. Math.. 121, 379–420 (2000). Publisher Full Text

Stenger, F, O’Reilly, MJ: Computing solutions to medical problems via sinc convolution. IEEE Trans. Autom. Control. 43, 843 (1998). Publisher Full Text

Narasimhan, S, Majdalani, J, Stenger, F: A first step in applying the sinc collocation method to the nonlinear Navier Stokes equations. Numer. Heat Transf., Part B. 41, 447–462 (2002). Publisher Full Text

Mueller, JL, Shores, TS: A new sincGalerkin method for convectiondiffusion equations with mixed boundary conditions. Comput. Math. Appl.. 47, 803–822 (2004). Publisher Full Text

ElGamel, M, Behiry, SH, Hashish, H: Numerical method for the solution of special nonlinear fourthorder boundary value problems. Appl. Math. Comput.. 145, 717–734 (2003). Publisher Full Text

Lybeck, NJ, Bowers, KL: Sinc methods for domain decomposition. Appl. Math. Comput.. 75, 4–13 (1996)