Journal of Applied Mathematics

Volume 2014 (2014), Article ID 134272, 8 pages

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

Research Article

Bernstein Collocation Method for Solving Nonlinear Fredholm-Volterra Integrodifferential Equations in the Most General Form

Ayşegül Akyüz-Daşcıoğlu,¹ Neşe İşler Acar,² and Coşkun Güler³

¹Department of Mathematics, Faculty of Arts & Sciences, Pamukkale University, 20070 Denizli, Turkey

²Department of Mathematics, Faculty of Arts & Sciences, Mehmet Akif Ersoy University, 15030 Burdur, Turkey

³Department of Mathematical Engineering, Faculty of Chemistry and Metallurgy, Yildiz Technical University, 34210 Istanbul, Turkey

Received 14 May 2014; Accepted 17 July 2014; Published 10 August 2014

Academic Editor: Naseer Shahzad

Copyright © 2014 Ayşegül Akyüz-Daşcıoğlu 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.

A collocation method based on the Bernstein polynomials defined on the interval is developed for approximate solutions of the Fredholm-Volterra integrodifferential equation (FVIDE) in the most general form. This method is reduced to linear FVIDE via the collocation points and quasilinearization technique. Some numerical examples are also given to demonstrate the applicability, accuracy, and efficiency of the proposed method.

1. Introduction

The quasilinearization method was introduced by Bellman and Kalaba [1] to solve nonlinear ordinary or partial differential equations as a generalization of the Newton-Raphson method. The origin of this method lies in the theory of dynamic programming. In this method, the nonlinear equations are expressed as a sequence of linear equations and these equations are solved recursively. The main advantage of this method is that it converges monotonically and quadratically to the exact solution of the original equations [2]. Therefore, the quasilinearization method is an effective approach for obtaining approximate solutions of nonlinear equations such as differential equations [3–7], functional equations [8, 9], integral equations [10–12], and integrodifferential equations [13–15].

In this paper, we consider the nonlinear FVIDE in the general form

under the initial

or boundary conditions

Here , , and are known functions, , , , , , , and are known constants, and is an unknown function.

Besides, we approximate to the nonlinear FVIDE (1) by the generalized Bernstein polynomials defined on the interval as

where denotes the generalized Bernstein basis polynomials of the form

For convenience, we set , if or .

Bernstein polynomials have many useful properties such as the positivity, continuity, differentiability, integrability, recursion’s relation, symmetry, and unity partition of the basis set over the interval . For more information about the Bernstein polynomials, see [16, 17]. Recently, these polynomials have been used for the numerical solutions of differential equations [4, 18, 19], integral equations [20–24], and integrodifferential equations [25–27].

Now, we give two main theorems for the generalized Berntein polynomials and their basis forms that were proved by Akyuz Dascioglu and Isler [4] as follows.

Theorem 1. If , for some integer , then

converges uniformly.

Proof. The above theorem can be easily proved by applying transformation to the theorem given on the interval [] by Phillips [28].

Theorem 2 (see [4]). There is a relation between generalized Bernstein basis polynomials matrix and their derivatives in the form

such that

Here the elements of matrix , , are defined by

In highlight of these theorems, a collocation method based on the generalized Bernstein polynomials, given in Section 2, is developed for the approximate solutions of the nonlinear FVIDE in the most general form (1) via the quasilinearization technique iteratively. In Section 3, some numerical examples are presented for exhibiting the accuracy and applicability of the proposed method. Finally, the paper ends with the conclusions in Section 4.

2. Method of the Solution

Our aim is to obtain a numerical solution of the nonlinear FVIDE in the general form (1) under conditions (2) or (3) in terms of the generalized Bernstein polynomials. For this, we firstly express this nonlinear equation as a sequence of linear FVIDEs via the quasilinearization technique iteratively. After that, using the collocation points yields the system of linear algebraic equations. This system represents a matrix equation given by the following theorem. Finally, solving this system with the conditions we get the desired approximate solution.

Theorem 3. Let be collocation points defined on the interval , and let the functions , , and be able to expand by Taylor series with respect to . Suppose that nonlinear FVIDE (1) has the generalized Bernstein polynomial solution. Then, the following matrix relation holds:

where is defined in Theorem 2, , , , and are    matrices, and and are matrices for , such that

Here is iteration index, is generalized Bernstein basis polynomials, and is given in the following proof.

Proof. Firstly, by applying the quasilinearization method to the nonlinear FVIDE (1), we obtain a sequence of linear FVIDEs:

where the expressions , , and represent partial differentiation of the functions , , and with respect to , and these are defined, respectively, as

Here is a reasonable initial approximation of the function , and is always considered known and is obtained from previous iteration. The recurrence relation (12) can now be written compactly in the form

denoting :

Notice that (14) is a linear FVIDE with variable coefficients, since is known function of . is an unknown function that has the Bernstein polynomial solution; also this function and its derivatives can be expressed by

By utilizing Theorem 2 and collocation points, above relation becomes

Substituting the collocation points and the relation (17) into (14), we obtain a linear algebraic system:

Here and are denoted by

and are defined, respectively, in Theorems 2 and 3.

For , the system (18) can be written compactly in the matrix form

so that

where matrices are clearly

Hence, the proof is completed.

Corollary 4. For , the nonlinear FVIDE (1) is reduced to the th order nonlinear differential equation

and by utilizing Theorem 3, this equation can be written as

Here the matrices , , and are defined as above, and elements of the matrix are denoted by

Corollary 5. For , the nonlinear FVIDE (1) is reduced to the nonlinear Fredholm-Volterra integral equation in the form

such that . From Theorem 3, this equation has the iteration matrix form as

Here , , , , and elements of these matrices are denoted as follows:

Now we can solve the nonlinear FVDIE (1) under the initial (2) or boundary (3) conditions as follows.

Step 1. Firstly, we use Theorem 3 for the nonlinear FVDIE (1) and determine the matrices in (10). This matrix equation is a system of linear algebraic equations with -unknown coefficients . Let the augmented matrix corresponding equation (10) be denoted by .

Step 2. We need to choose the first iteration function for calculating the and . Notice that this function can be obtained in a variety of ways. For instance, it can be obtained from the physical situation for engineering problems. However, a very rough choice for the first iteration function such as initial value is enough for the procedure to converge. We can also consider that the first iteration function as the highest degree polynomial satisfied the given conditions (2) or (3).

Step 3. From expression (17), initial (2) and boundary (3) conditions can be written in the matrix forms, respectively,

where the matrices are

Besides, (29) can be denoted by the augmented matrices and .

Step 4. To obtain the solution of nonlinear FVIDE (1) under the given conditions, we insert the elements of the row matrices or to the end of the augmented matrix . In this way, we have the new augmented matrix , that is, rectangular matrix.

Step 5. If , then unknown coefficients are uniquely determined for each iteration . This kind of systems can be solved by the Gauss Elimination, Generalized Inverse, and QR factorization methods.

3. Numerical Results

Four numerical examples are given to illustrate the applicability, accuracy, and efficiency of the proposed method. All results are computed by using an algorithm written in Matlab 7.1. Besides, in the tables, the absolute and -norm errors are computed numerically on the collocation points , by the folowing formulas:

such that is a Bernstein approximation for the th iteration function and is an exact solution.

Example 1. Consider the nonlinear Volterra integrodifferential equation

under the initial condition . Exact solution of this equation is . Here and .

Let be the first iteration function. From Theorem 3, matrix relation of the above problem can be written as

where , , , and ( is an identity matrix with dimensional ), because of the . Elements of the matrices become

A numerical comparison of the proposed method with the Hybrid method [29] is given in Table 1. The obtained numerical results for two different initial approximations are also listed in Table 1. It can be seen that the computational results of the proposed method are better and more effective for smaller values and iterations than the other method, and the choice of the higher degree polynomial for the initial approximation leads to better results.

Table 1: Comparison of the -norm errors for Example 1.

Presented method
Hybrid method [29]

Example 2. Consider the nonlinear Fredholm-Volterra integrodifferential equation

with the initial condition that exact solution is . Here the required functions and constants are denoted by , , , , and .

We have two choices satisfying the initial condition for the first iteration functions such that and . Let the first iteration function be , because of the higher degree. From Theorem 3, matrix relation of the above problem is

where , , , , and because of the . Elements of these matrices are as follows:

Table 2 contains a numerical comparison of the proposed method between the numerical method based on fixed point theorem [30] and direct method by using triangular functions [31]. The table reveals that convergence of the presented method is faster and more accurate than the others.

Table 2: Comparison of the absolute errors for Example 2.

Presented method
Direct method [31]
Fixed point method [30]

,
,
,
,

Example 3. Consider the third-order nonlinear Fredholm integrodifferential equation

with the boundary conditions

The exact solution of the above equation is . Here and .

Let the first iteration function be . From Theorem 3, matrix relation of the above problem can be denoted by

such that , , , and . Here elements of matrices are, respectively,

In Table 3, absolute errors of the proposed method are given for different values and iterations . The table shows that the presented method converges quite rapidly for increasing values and iterations . Besides, the absolute error of the homotopy analysis method [32] given with figure is approximately for iteration . Therefore, we can say that the proposed method has more effective numerical results than the other methods.

Table 3: Absolute errors for Example 3.

Example 4. Consider the fourth-order Volterra integrodifferential equation

with the boundary conditions

that exact solution is . Here the functions are and .

Let the first iteration function be . From Theorem 3, matrix form of the above problem can be written as

where , , , and . Here elements of matrices are, respectively,

In Table 4, -norm errors of the proposed method are compared with the -norm errors computed via the Hybrid of Block-Pulse functions and Lagrange interpolation polynomials [29] for different values . While the results of other methods are effective, the computational results of the proposed method are more rapid, effective, for smaller values .

Table 4: Comparison of the -norm errors for Example 4.

Presented method
Hybrid method [29]

6

8

10

12

15

4. Conclusions

In general, nonlinear integrodifferential equations can not be solved analytically. For this reason, numerical solutions of nonlinear equations are needed. With the presented method, we have reduced the nonlinear FVIDE (1) to a sequence of linear equations depending on the collocation points and the iteration function, and then, combining the conditions with obtained linear matrix equation, we have the generalized Bernstein polynomials solution. Besides, the proposed method is valid for both nonlinear differential and integral equations; this is explained in Section 2. If unknown function is continuous on the interval , then the Bernstein collocation method can be used for solving these kinds of equations via the quasilinearization technique. For convenience, the first iteration function can be choosen as initial value. Besides, this function can be select as the higher degree polynomials satisfied the given conditions for providing the better approximation. These polynomials can be obtained from interpolation or the least square approximation methods. The applicability and accuracy of the proposed method have been tested with some numerical examples. These results have shown that the presented method converges rapidly for all problems. Consequently, all these positive implications lead to applicability of the proposed method for numerical solutions of any other kinds of problems including nonlinear equations.

Conflict of Interests

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

Acknowledgment

This work is supported by Scientific Research Project Coordination Unit of Pamukkale University with no. 2012FBE036.

References

R. E. Bellman and R. E. Kalaba, Quasilinearization and Nonlinear Boundary Value Problems, Elsevier, New York, NY. USA, 1965.

E. L. Stanley, Quasilinearization and Invariant Imbedding, Academic Press, New York, NY, USA, 1968.

R. P. Agarwal and Y. M. Chow, “Iterative methods for a fourth order boundary value problem,” Journal of Computational and Applied Mathematics, vol. 10, no. 2, pp. 203–217, 1984.

A. Akyuz Dascıoglu and N. Isler, “Bernstein collocation method for solving nonlinear differential equations,” Mathematical & Computational Applications, vol. 18, no. 3, pp. 293–300, 2013.

A. Charles and J. Baird, “Modified quasilinearization technique for the solution of boundary-value problems for ordinary differential equations,” Journal of Optimization Theory and Applications, vol. 3, no. 4, pp. 227–242, 1969.

V. B. Mandelzweig and F. Tabakin, “Quasilinearization approach to nonlinear problems in physics with application to nonlinear ODEs,” Computer Physics Communications, vol. 141, no. 2, pp. 268–281, 2001.

J. I. Ramos, “Piecewise quasilinearization techniques for singular boundary-value problems,” Computer Physics Communications, vol. 158, no. 1, pp. 12–25, 2004.

B. Ahmad, R. Ali Khan, and S. Sivasundaram, “Generalized quasilinearization method for nonlinear functional differential equations,” Journal of Applied Mathematics and Stochastic Analysis, vol. 16, no. 1, pp. 33–43, 2003.

Z. Drici, F. A. McRae, and J. V. Devi, “Quasilinearization for functional differential equations with retardation and anticipation,” Nonlinear Analysis: Theory, Methods & Applications, vol. 70, no. 4, pp. 1763–1775, 2009.

F. Caliò, F. Munoz, and E. Marchetti, “Direct and iterative methods for the numerical solution of mixed integral equations,” Applied Mathematics and Computation, vol. 216, no. 12, pp. 3739–3746, 2010.

K. Maleknejad and E. Najafi, “Numerical solution of nonlinear Volterra integral equations using the idea of quasilinearization,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 93–100, 2011.

S. G. Pandit, “Quadratically converging iterative schemes for nonlinear volterra integral equations and an application,” Journal of Applied Mathematics and Stochastic Analysis, vol. 10, no. 2, pp. 169–178, 1997.

B. Ahmad, “A quasilinearization method for a class of integro-differential equations with mixed nonlinearities,” Nonlinear Analysis: Real World Applications, vol. 7, no. 5, pp. 997–1004, 2006.

B. Ahmad, R. A. Khan, and S. Sivasundaram, “Generalized quasilinearization method for integro-differential equations,” Nonlinear Studies, vol. 8, no. 3, pp. 331–341, 2001.

P. Wang, Y. Wu, and B. Wiwatanapaphee, “An extension of method of quasilinearization for integro-differential equations,” International Journal of Pure and Applied Mathematics, vol. 54, no. 1, pp. 27–37, 2009.

R. T. Farouki and V. T. Rajan, “Algorithms for polynomials in Bernstein form,” Computer Aided Geometric Design, vol. 5, no. 1, pp. 1–26, 1988.

G. G. Lorentz, Bernstein Polynomials, Chelsea, New York, NY, USA, 2nd edition, 1986.

E. H. Doha, A. H. Bhrawy, and M. A. Saker, “Integrals of Bernstein polynomials: an application for the solution of high even-order differential equations,” Applied Mathematics Letters, vol. 24, no. 4, pp. 559–565, 2011.

O. R. Işık, M. Sezer, and Z. Güney, “A rational approximation based on Bernstein polynomials for high order initial and boundary values problems,” Applied Mathematics and Computation, vol. 217, no. 22, pp. 9438–9450, 2011.

S. Bhattacharya and B. N. Mandal, “Use of Bernstein polynomials in numerical solutions of Volterra integral equations,” Applied Mathematical Sciences, vol. 2, no. 33–36, pp. 1773–1787, 2008.

K. Maleknejad, E. Hashemizadeh, and R. Ezzati, “A new approach to the numerical solution of Volterra integral equations by using Bernstein's approximation,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 2, pp. 647–655, 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.

A. Shirin and M. S. Islam, “Numerical solutions of Fredholm integral equations using Bernstein polynomials,” Journal of Scientific Research, vol. 2, no. 2, pp. 264–272, 2010.

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–8, pp. 241–255, 2009.

O. R. Işık, 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.

Y. Ordokhani and S. Davaei far, “Application of the Bernstein polynomials for solving the nonlinear Fredholm integro-differential equations,” Journal of Applied Mathematics & Bioinformatics, vol. 1, no. 2, pp. 13–31, 2011.

B. M. Pandya and D. C. Joshi, “Solution of a Volterra's population model in a Bernstein polynomial basis,” Applied Mathematical Sciences, vol. 5, no. 69, pp. 3403–3410, 2011.

G. M. Phillips, Interpolation and Approximation by Polynomials, Springer, New York, NY, USA, 2003.

H. R. Marzban and S. M. Hoseini, “Solution of nonlinear Volterra-Fredholm integrodifferential equations via hybrid of block-pulse functions and Lagrange interpolating polynomials,” Advances in Numerical Analysis, vol. 2012, Article ID 868279, 14 pages, 2012.

M. I. Berenguer, D. Gámez, and A. J. López Linares, “Fixed point techniques and Schauder bases to approximate the solution of the first order nonlinear mixed Fredholm—Volterra integro-differential equation,” Journal of Computational and Applied Mathematics, vol. 252, pp. 52–61, 2013.

E. Babolian, Z. Masouri, and S. Hatamzadeh-Varmazyar, “Numerical solution of nonlinear Volterra-Fredholm integro-differential equations via direct method using triangular functions,” Computers & Mathematics with Applications, vol. 58, no. 2, pp. 239–247, 2009.

A. Shidfar, A. Molabahrami, A. Babaei, and A. Yazdanian, “A series solution of the nonlinear Volterra and Fredholm integro-differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 2, pp. 205–215, 2010.