Abstract and Applied Analysis

Volume 2014 (2014), Article ID 796286, 7 pages

http://dx.doi.org/10.1155/2014/796286

Research Article

On Bernstein Polynomials Method to the System of Abel Integral Equations

A. Jafarian,¹ S. Measoomy Nia,¹ Alireza K. Golmankhaneh,² and D. Baleanu3,4,5

¹Department of Mathematics, Islamic Azad University, Urmia Branch, P.O. BOX 969, Urmia, Iran

²Department of Physics, Islamic Azad University, Urmia Branch, P.O. BOX 969, Urmia, Iran

³Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, P.O. Box 80204, Jeddah 21589, Saudi Arabia

⁴Department of Mathematics and Computer Science, Çankaya University, 06530 Ankara, Turkey

⁵Institute of Space Sciences, P.O. Box MG-23, 76900 Magurele-Bucharest, Romania

Received 20 November 2013; Revised 3 March 2014; Accepted 4 March 2014; Published 7 May 2014

Academic Editor: Hossein Jafari

Copyright © 2014 A. Jafarian et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper deals with a new implementation of the Bernstein polynomials method to the numerical solution of a special kind of singular system. For this aim, first the truncated Bernstein series polynomials of the solution functions are substituted in the given problem. Using some properties of these polynomials, the solution of the problem is reduced to solve a linear system of algebraic equations. In order to confirm the reliability and accuracy of the proposed method, some weakly Abel integral equations systems with comparisons are solved in detail as numerical examples.

1. Introduction

Many physical, biological, or engineering processes involve relating rates of change of various quantities according to physical or other principles. In some occasions, mathematical representation of these processes appears in a natural way in terms of singular integral equations. It should be noted that these types of equations contain much more mathematical information about the applied scientific problem aroused from natural phenomena. On the other hand, much singular integral equations are not easy to solve analytically. Hence, numerical methods based on solving these equations can be unbelievably powerful. Among the numerous numerical approaches that have been offered, let us mention some well-known methods. For example, Khan and Gondal [1] merged the Laplace transform and decomposition method and then presented a new mechanism for the solution of Abel type singular integral equations, that is to say, the two-step Laplace decomposition algorithm. Bougoffa et al. [2] used the Adomian decomposition method for solving linear and nonlinear Abel integral equations. Abdulkawi et al. [3] used the Chebyshev polynomials of the second kind for the numerical solution of the Cauchy type singular integral equations of the first kind, over a finite segment. The solution of a singular integral equation with logarithmic kernel in two disjoint intervals was obtained using function theoretic method by Banerjea and Dutta in [4]. The two interval conditions yielded to the two simple simultaneous algebraic equations in two unknowns, which were solved to recover the reflection and transmission coefficients in the scattering problem. Helsing in [5] presented a fast and stable solver based on the reduction, Nystrom discretization, composite quadrature, recursive compressed inverse preconditioning, and multipole acceleration techniques for the numerical solution of these integral equations on piecewise smooth curves. Maleknejad and Salimi Shamloo [6] utilized Laplace transform of the problem and operational matrices of piecewise constant orthogonal functions to solve singular Volterra integral equations. Moreover, Huang et al. [7] offered a stable approximate inversion of Abel integral equation by using the Taylor expansion method. This work transformed the equation approximately to a system of linear algebraic equations for the unknown function together with its derivatives. Akel and Hussein [8] proposed a convenient and efficient numerical method for the treatment of singular integral equations of the first and second kind. The method yielded a simple and closed form of the approximate solution by means of Sinc functions with smoothing transformations. Recently, we employed a structure of artificial neural networks for approximating solution of a fuzzy type Abel integral equation. Also, the application of Legendre wavelets as a basis function is offered for solving a special kind of singular system in [9]. Recently, [10–17], wavelets method has been employed for solving a class of system of nonlinear singular fractional Volterra integrodifferential equations. In this work, shifted Chebyshev polynomials and their properties were employed to derive a general procedure for forming the operational matrix of fractional derivative for Chebyshev wavelets.

In many problems in science and engineering, we have some unknown functions which are too complicated to be determined. The Bernstein polynomials method is one of the earliest analytic-numeric algorithms for approximating the unknowns in different kinds of mathematical problems. This is an extremely useful way of expressing a complicated function in terms of simple polynomial. The only requirement is that the given function should be smooth. In other words, at a point of interest it must be possible to differentiate the function as often as we please. In our earlier works, the linear integral equations system was handled by using these kinds of polynomials [18]. Moreover, in [19–22] Bernstein polynomials and Bernstein operational matrices were employed for solving some types of singular integral and integrodifferential equations. In this paper, we are going to apply a development of the Bernstein series method to approximate solution of a system of weakly Abel integral equations which has been considered in the standard form:

where is a real constant, and are predetermined real-valued functions, and is the vector solution that will be determined. For more details on the integral equations, the reader is referred to the book [23]. Supposedly the present problem has a unique solution [23], a comfort computational form of the above system is obtained when the Bernstein polynomials are used as basis functions. Now, putting (for ) transforms the resulting system to a system of linear algebraic equations for the unknown functions. Then, a numerical solution can be obtained by solving the resulting system using a standard rule. If more and more terms are used from the Bernstein series, then the polynomial representations better and better approximate the unknowns. It is clear that the -th order series solution converges to the exact solution if the unknown functions are polynomials of degree up to . Here is an outline of the paper. Section 2 intends to describe how to find approximate solution of a weakly Abel integral equations system by using the present approach. In Section 3, the error analysis of approximate solution corresponding to the singular integral equations system is given. To demonstrate the efficiency and reliability of the method, some numerical test examples with comparisons are given in Section 4. Section 5 concludes the paper.

2. Description of the Method

The Bernstein polynomials representation is extremely important from many points of view. This approach can be used to approximate more complicated functions in terms of simpler polynomials. In this section, our aim is to show how the Bernstein polynomials method can be applied to the numerical solution of singular system (1). First, assume that the functions and are continually differentiable in the interval of interest. We know that the Bernstein approximation of function is the polynomial:

where is a constant coefficient and is the Bernstein basis polynomial of degree for each positive integer , which is defined by Mandal and Bhattacharya [24]:

It is clear that the value of approaches as gets larger. Suppose we carry this argument and try to determine polynomial approximations of degree for the solutions in the given problem. To do this, first is substituted for in system (1). Now, we get

or equivalently

where

After some simplifications, (5) is reduced to the following form:

where

For suitable choice of the distinct values (for ) in (e.g., ), we achieve the following linear algebraic system:

where

Therefore, system (9) with linear equations of unknowns (for ) can be rewritten in the matrix form as follows:

or equivalently

with the block matrices:

Consequently, system (11) can be solved using a standard rule to obtain the unknown constants (for ). Ultimately, the unknown function can be easily approximated using the aforementioned assumptions.

3. Convergence Analysis

In this section, we intend to prove that the presented numerical method converges to the exact solution of system (1).

Theorem 1. Let (for ) be polynomials of degree that their numerical coefficients are produced by solving the linear system (11). There exists an integer such that, for all , these polynomials converge to the exact solution of the integral equations system (1).

Proof. We know that the Bernstein polynomials for the function that is continuous on interval converge uniformly to this function [24]. Consider again system (1). In other words, can be expanded as a uniformly convergent Bernstein series in :

Based on the suggested process, Abel integral equations system (1) can be transformed to the following equivalent infinitely systems of linear equations for unknown :

with

where , , and were defined in the previous section. For the above system, the unique solution can be expressed as

where . Alternatively, the above system can be rewritten as

where

Consequently, one finds that the vector composed of the first elements of the exact solution vector must satisfy the following relation:

In addition, based on the analysis of the foregoing section, the unique solution of (11) is denoted as

Subtracting (21) from (20) yields

where . Expanding the right-hand side of (22), the left-hand side of this is expressed by

where and are the elements of and , respectively. Thus,

follows from the well-known Cauchy-Schwarz inequality. Since and , so we can conclude that and proof is completed.

4. Numerical Examples

In this section, the above-mentioned method is illustrated with the help of three illustrative examples which include second kind Volterra integral equations system with Abel kernels. Note that the empirical results obtained here will be compared with the ones achieved from the Taylor expansion method (TEM) of degree [7].

Example 1. Consider the following singular integral equations system:

where

with the exact vector solution . In this example, we aim to apply the present technique to approximate solution of the problem by taking . The solutions in series forms are given by

As shown, the method is quite accurate when the solutions are polynomials of degree up to .

Example 2. Consider the system of Abel integral equations:

which has the exact solutions and . For suitable distinct points (for ) in in which is taken near 0 and , the numerical results obtained here and the Taylor expansion quadrature method are presented in Table 1. Also, the exact and approximate solutions are compared in Figure 1.

Table 1: Numerical results for Example 2.

Exact solution
Proposed method
TEM

0.2153
1.1072
0.215290
1.107242
0.215169
1.107682

0.3259
0.8541
0.325884
0.854115
0.325780
0.854328

0.4276
0.7555
0.427565
0.755497
0.427509
0.755881

0.5293
0.7046
0.529335
0.704608
0.529611
0.704209

0.6350
0.6752
0.635031
0.675237
0.635295
0.675831

0.7470
0.6575
0.747040
0.657498
0.747548
0.657637

0.8672
0.6468
0.867187
0.646770
0.867735
0.646340

0.9971
0.6406
0.997089
0.640586
0.997818
0.640385

1.1383
0.6375
1.138298
0.637503
1.138803
0.637804

1.2924
0.6366
1.292388
0.636619
1.292765
0.636391

(a) and

(b) and

Figure 1: Comparison of the exact and approximated solutions for Example 2.

Example 3. Let us consider the Abel integral equations system:

where

with the exact solutions , , and . In Table 2, the numerical results obtained here for are compared with the exact solutions. It is clear that errors between the exact and approximate values of and are zero. Similarly, the series solutions of order 3 with Taylor approximation converge to the exact solutions and . Also, the exact and approximate solutions of are compared in Figure 2 for and .

Table 2: Numerical results for Example 3.

Exact solution
Approximate solution

3.3000
0.1500
1.0488
3.3000
0.1500
1.0488

3.6000
0.6000
1.0954
3.6000
0.6000
1.0951

3.9000
1.3500
1.1402
3.9000
1.3500
1.1389

4.2000
2.4000
1.1832
4.2000
2.4000
1.1810

4.5000
3.7500
1.2247
4.5000
3.7500
1.2198

4.8000
5.4000
1.2649
4.8000
5.4000
1.2590

5.1000
7.3500
1.3038
5.1000
7.3500
1.2898

5.4000
9.6000
1.3416
5.4000
9.6000
1.3310

5.7000
12.150
1.3784
5.7000
12.150
1.3687

6.0000
15.000
1.4142
6.0000
15.000
1.4050

Figure 2: , , and for Example 3.

5. Conclusions

In this paper, we presented a useful numerical method that originated mainly from the Bernstein polynomials for solving Abel integral equations systems. This method converted the present problem to a system of linear algebraic equations for unknown functions. Having determined the unknown Bernstein coefficients of the solution functions, the series solutions were produced immediately. An interesting feature of this method is that the exact solution is derived if it is a polynomial of degree or less than . On the other hand, if an exact solution is not obtainable, then the obtained series can be used for numerical purposes. In this case, more terms must be evaluated to the higher accuracy level. Additionally, the proposed technique has been compared with the Taylor expansion method. The obtained numerical results from analyzed examples illustrated that, in applications involving computations with polynomials, the Bernstein form offers efficient algorithm versus the traditional Taylor rule, for many basic functions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

M. Khan and M. A. Gondal, “A reliable treatment of Abel's second kind singular integral equations,” Applied Mathematics Letters, vol. 25, no. 11, pp. 1666–1670, 2012.

L. Bougoffa, R. C. Rach, and A. Mennouni, “A convenient technique for solving linear and nonlinear Abel integral equations by the Adomian decomposition method,” Applied Mathematics and Computation, vol. 218, no. 5, pp. 1785–1793, 2011.

M. Abdulkawi, Z. K. Eshkuvatov, and N. M. A. Nik Long, “A note on the numerical solution of singular integral equations of cauchy type,” World Academy of Science, Engineering and Technology, vol. 3, no. 9, pp. 893–896, 2009.

S. Banerjea and B. Dutta, “Solution of a singular integral equation and its application to water wave problems,” Journal of the Indian Institute of Science, vol. 86, no. 3, pp. 265–278, 2006.

J. Helsing, “A fast and stable solver for singular integral equations on piecewise smooth curves,” SIAM Journal on Scientific Computing, vol. 33, no. 1, pp. 153–174, 2011.

K. Maleknejad and A. Salimi Shamloo, “Numerical solution of singular Volterra integral equations system of convolution type by using operational matrices,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 500–505, 2008.

L. Huang, Y. Huang, and X. Li, “Approximate solution of Abel integral equation,” Computers and Mathematics with Applications, vol. 56, no. 7, pp. 1748–1757, 2008.

M. S. Akel and H. S. Hussein, “Numerical treatment of solving singular integral equations by using Sinc approximations,” Applied Mathematics and Computation, vol. 218, no. 7, pp. 3565–3573, 2011.

J. Biazar and H. Ebrahimi, “A new technique for systems of Abel-Volterra integral equations,” International Journal of the Physical Sciences, vol. 7, no. 1, pp. 89–99, 2012.

M. H. Heydari, M. R. Hooshmandasl, F. Mohammadi, and C. Cattani, “Wavelets method for solving systems of nonlinear singular fractional Volterra integro-differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 1, pp. 37–48, 2014.

S. A. Yousefi and M. Behroozifar, “Operational matrices of Bernstein polynomials and their applications,” International Journal of Systems Science, vol. 41, no. 6, pp. 709–716, 2010.

M. Alipour and D. Rostamy, “Solving nonlinear fractional differential equations by Bernstein polynomials operational matrices,” The Journal of Mathematics and Computer Science, vol. 5, no. 3, pp. 185–196, 2012.

D. Rostamy, M. Alipour, H. Jafari, and D. Baleanu, “Solving multi-term orders fractional differential equations by operational matrices of BPs with convergence analysis,” Romanian Reports in Physics, vol. 65, no. 2, pp. 334–349, 2013.

D. Baleanu, M. Alipour, and H. Jafari, “The Bernstein operational matrices for solving the fractional quadratic riccati differential equations with the Riemann-Liouville derivative,” Abstract and Applied Analysis, vol. 2013, Article ID 461970, 7 pages, 2013.

M. Alipour and D. Baleanu, “Approximate analytical solution for nonlinear system of fractional differential equations by BPs operational matrices,” Advances in Mathematical Physics, vol. 2013, Article ID 954015, 9 pages, 2013.

M. Alipour, D. Rostamy, and D. Baleanu, “Solving multi-dimensional fractional optimal control problems with inequality constraint by Bernstein polynomials operational matrices,” Journal of Vibration and Control, vol. 19, no. 16, pp. 2523–2540, 2013.

M. Alipour and D. Rostamy, “BPs operational matrices for solving time varying fractional optimal control problems,” The Journal of Mathematics and Computer Science, vol. 6, pp. 292–304, 2013.

A. Jafarian, S. Measoomy Nia, A. K. Golmankhaneh, and D. Baleanu, “Numerical solution of linear integral equations system using the Bernstein collocation method,” Advances in Difference Equations, vol. 2013, article 123, 2013.

V. K. Singh, R. K. Pandey, and O. P. Singh, “New stable numerical solutions of singular integral equations of Abel type by using normalized Bernstein polynomials,” Applied Mathematical Sciences, vol. 3, no. 5, pp. 241–255, 2009.

O. P. Singh, V. K. Singh, and R. K. Pandey, “A stable numerical inversion of Abel's integral equation using almost Bernstein operational matrix,” Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 111, no. 1, pp. 245–252, 2010.

S. Dixit, O. P. Singh, and S. Kumar, “A stable numerical inversion of generalized Abels integral equation,” Applied Numerical Mathematics, vol. 62, no. 5, pp. 567–579, 2012.

O. R. Işik, M. Sezer, and Z. Güney, “Bernstein series solution of a class of linear integro-differential equations with weakly singular kernel,” Applied Mathematics and Computation, vol. 217, no. 16, pp. 7009–7020, 2011.

A. M. Wazwaz, Linear and Nonlinear Integral Equations: Methods and Applications, Higher Education Press, Beijing, China; Springer, Berlin, Germany, 2011.

B. N. Mandal and S. Bhattacharya, “Numerical solution of some classes of integral equations using Bernstein polynomials,” Applied Mathematics and Computation, vol. 190, no. 2, pp. 1707–1716, 2007.