Abstract and Applied Analysis

Volume 2013 (2013), Article ID 929478, 10 pages

http://dx.doi.org/10.1155/2013/929478

Research Article

Numerical Study of Two-Dimensional Volterra Integral Equations by RDTM and Comparison with DTM

Reza Abazari¹ and Adem Kılıçman²

¹Young Researchers and Elite Club, Ardabil Branch, Islamic Azad University, Ardabil, Iran

²Department of Mathematics and Institute of Mathematical Research, University Putra Malaysia (UPM), 43400 Serdang, Malaysia

Received 17 April 2013; Accepted 10 June 2013

Academic Editor: Santanu Saha Ray

Copyright © 2013 Reza Abazari and Adem Kılıçman. 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.

The two-dimensional Volterra integral equations are solved using more recent semianalytic method, the reduced differential transform method (the so-called RDTM), and compared with the differential transform method (DTM). The concepts of DTM and RDTM are briefly explained, and their application to the two-dimensional Volterra integral equations is studied. The results obtained by DTM and RDTM together are compared with exact solution. As an important result, it is depicted that the RDTM results are more accurate in comparison with those obtained by DTM applied to the same Volterra integral equations. The numerical results reveal that the RDTM is very effective, convenient, and quite accurate compared to the other kind of nonlinear integral equations. It is predicted that the RDTM can be found widely applicable in engineering sciences.

1. Introduction

Mathematical modeling of many problems in science, engineering, physics, and other disciplines leads to linear and nonlinear integrodifferential equations (IDE). The great use of mathematical models including integrodifferential equations is one of the main reasons obtaining the solutions of this kind of problems (see, e.g., [1–3] and the references therein). So, it is very important to get some information about the analytical solutions of these problems because these solutions give significant information about the character of the modeled event. But, in some cases, it is more difficult to obtain analytical solutions of these models. These are usually difficult to solve analytically, and in many cases the solution must be approximated. To approximate the solutions of these models, in recent years several numerical approaches have been proposed.

In this paper, we consider the following Volterra type of integral equation [4, 5]:

where and are continuous functions and has the following form:

The one-dimensional Volterra type of integral equation has been solved by many numerical methods, such as collocation methods [1], Taylor-series expansion methods [2], Gauss-type quadratures method [3], spectral methods [6], Chebyshev polynomial method [7], Tau method [8], sine-cosine wavelets method [9], Monte Carlo method [10], and Haar functions method [11].

But in two-dimensional cases, a small amount of work has been done (see, e.g., [12–14]). Very recently, Tari et al. in [4] employed the classic differential transform method for solving two-dimensional Volterra type of integral equations (1), and Jang in [5] improved the proofs of the presented theorems by Tari et al. in [4]. They derived fundamental properties of the differential transforms of some kernel functions in Volterra integral equations.

However, the classic differential transform method, introduced by Zhou [15], is based on the definition of the differential transform, which is a Taylor series. Thus, it requires a cumbersome calculation to obtain the basic properties of the differential transforms. Some of DTM applications are mentioned in [16–21].

Recently, Keskin and Oturanç introduced a reduced form of DTM as reduced DTM (RDTM) and applied it to approximate some PDE [22] and factional PDEs [23]. More recently, Abazari and Ganji [24] extended RDTM to study the partial differential equation with proportional delay in and shrinking in and showed that, as a special advantage of RDTM rather than DTM, the reduced differential transform recursive equations produce exactly all the Poisson series coefficients of solutions based on the initial condition as weighted function, whereas the differential transform recursive equations produce exactly all the Taylor series coefficients of solutions.

Here, we suggest the RDTM, for the approximating of the solutions of the two-dimensional Volterra integral equations (1) with the same kernel functions in [4, 5]. In order to demonstrate the effectiveness of the RDTM, the illustrative examples for the same kernel function of references [5] are presented. These examples show that the RDTM produces exactly all the POISSON series coefficients (see Remark 5) of the exact solutions, whereas, the classic DTM produces exactly all the TAYLOR series coefficients of the exact solutions. As an important result, notwithstanding the simplicity and robustness of RDTM, it is depicted that the RDTM results are more accurate in comparison with those obtained by classic DTM.

2. Basic Definitions

With reference to the articles [16–21], the basic definitions of two-dimensional differential transform method (DTM) and their reduced form (RDTM) are introduced in the following two subsections, respectively.

2.1. Two-Dimensional DTM

Consider a function of two variables , and suppose that it can be represented as a product of two single-variable functions, that is, . On the basis of the properties of the one-dimensional differential transform, the function can be represented as

where is called the spectrum of .

The basic definitions and operations for two-dimensional differential transform are introduced as follows.

Definition 1. If is analytic and continuously differentiable with respect to time in the domain of interest, then

where the spectrum function is the transformed function, which is also called -function in brief.

The differential inverse transform of   is defined as

Combining (4) and (5), it can be obtained that

When are taken as , then (5) can be expressed as

In real applications, the function is represented by a finite series of (7) that can be written as

and (7) implies that is negligibly small. Usually, the values of and are decided by convergency of the series coefficients.

From the above definitions, it can be found that the concept of the two-dimensional differential transform is derived from the two-dimensional Taylor series expansion. With (4) and (5), the fundamental mathematical operations performed using the two-dimensional differential transform may be readily obtained, and these are listed in Table 1. (See [4, 5, 15, 16].)

Table 1: The fundamental operations of two-dimensional differential transform method.

Original function
Transformed function

Recently, Jang [5] extended the two-dimensional DTM on (1) as follows.

Theorem 2. Assume that , , , and are the differential transforms of the functions , , , and , respectively; then we have the following:

(a)if , then

(b)if , then

(c)if , then

Proof. See [5].

Theorem 3. Assume that , , and are the differential transforms of the functions , , and , respectively; then we have the following:

(a)if , then

(b)if , then

Proof. See [5].

2.2. Two-Dimensional Reduced DTM (RDTM)

Consider a function of two variables , and suppose that it can be represented as a product of two single-variable functions, that is, . Based on the properties of one-dimensional differential transform, the function can be represented as

where is called the spectrum of .

Remark 4. The poisson function series generates a multivariate Taylor series expansion of the input expression , with respect to the variables , to order , using the variable weights .

Remark 5. The relationship introduced in (14) is the poisson series form of the input expression , with respect to the variables and , to order , using the variable weights .

Similarly on previous section, the basic definitions of two-differential reduced differential transformation are introduced as follows.

Definition 6. If is analytical function in the domain of interest, then the spectrum function

is the reduced transformed function of .

Similarly on previous sections, the lowercase respects the original function while the uppercase stands for the reduced transformed function. The differential inverse transform of is defined as

Combining (15) and (16), it can be obtained that

In real applications, the function is represented by a finite series of (16), around , and can be written as

and (18) implies that is negligibly small. Usually, the values of and are decided by convergency of the series coefficients. From the above proposition, it can be found that the concept of the reduced two-dimensional differential transform is derived from the two-dimensional differential transform method. With (15) and (16), the fundamental mathematical operations performed by reduced two-dimensional differential transform can readily be obtained and listed in Table 2.

Table 2: The fundamental operations of two-dimensional RDTM.

Original function
Reduced transformed function

Similarly on previous subsection, we can extend the RDTM on Volterra integral equations (1) as follow.

Theorem 7. Assume that , , , and are the reduced differential transforms of the functions , , , and , respectively; then we have the following:

(a)if , then

(b)if , then

Proof. (a) According to the fundamental operations of two-dimensional RDTM listed in Table 2 and from Leibnitz formula, we get

therefore

and then, from using (15), for , we get

(b) Analogous to part (a), we get

therefore

and then from using (15), for , we get

Theorem 8. Assume that , , and are the reduced differential transforms of the functions , , and , respectively; then we have the following:

(a)if , then

(b)if , then

Proof. (a) By following the same manner as in the Theorem 7, we get

then

therefore

and then from using (15), for , we get

(b) Analogous to part (a), we get

then

therefore

and then from using (15), for , we get

and therefore

3. Numerical Results of DTM and RDTM

In this section, the reduced differential transform technique is described to solve a class of Volterra integral equations (1) with kernel functions of (2). In order to demonstrate the effectiveness of the RDTM, the illustrative examples for the same kernel function of [5] are presented. In each example, the numerical results of DTM, RDTM, and their comparisons with exact solution are given in separate tables. The results of the test examples show that the RDTM results are more powerful than DTM results.

Example 9. In the first example, consider the following two-dimensional Volterra integral equation [5]:

(a)DTM: Jang [5] solved this equation by using DTM and obtained the following five-term DTM solution:

(b)RDTM: from Volterra integral equation (38), it is easy to see that the , and therefore RDTM version is

By applying the RDTM properties listed in Theorem 8, on Volterra integral equation (38), for , we get

where is the reduced differential transform of . After expanding the RDTM recurrence equations (41), with initial value of (40), for , the first five terms of   are obtained as follows:

In the same manner, the rest of the components can be obtained by using the recursive equations (41). Substituting the quantities (41) in (18), the approximation solution of Volterra integral equation (38) in the Poisson series form is

which is the same as the first five terms of the Poisson series of the exact solution . The numerical results obtained with RDTM are presented in Table 3, in comparison with the classic DTM solution of [5] and the exact solution , for some points of the intervals and .

Table 3: Comparisons of the exact solution , with obtained by classic DTM [5] and obtained by RDTM at some test points in Example 9.

Classic DTM [5]
Reduced DTM

0.2
0.1
+0.4591040413
+0.4591115320

+0.4591040577

0.4
−1.0258569894
−1.0257575253

−1.0257906614

0.7
−2.7833269096
−2.7813944067

−2.7814532650

1
−4.6912625375
−4.6755840000

−4.6756686352

0.5
0.1
+2.0517139787
+2.0516661667

+2.0517139939

0.4
+0.5566653819
+0.5573213333

+0.5567258924

0.7
−1.1728156344
−1.1698785000

−1.1711316158

1
−2.9974949866
−2.9817708333

−2.9836944470

0.8
0.1
+3.8966576735
+3.8929643413

+3.8966576865

0.4
+2.3456312688
+2.3439851520

+2.3456822669

0.7
+0.5994989973
+0.6014245867

+0.6008906123

1
−1.1895390524
−1.1754026667

−1.1783733000

1
0.1
+5.2041732481
+5.1881855000

+5.2041732592

0.4
+3.5825396760
+3.5680853333

+3.5825824990

0.7
+1.7949988863
+1.7840151667

+1.7961453619

1
+0.0000000000
−0.0000000000
−0.0000000000
+0.0090050384

Example 10. In the second example, consider the following two-dimensional Volterra integral equation [5]:

(a)DTM: the approximation solution of this equation is also obtained by DTM in [5] as follows:

(b)RDTM: it is easy to see that the , and therefore RDTM version is

By applying the RDTM on nonlinear Volterra integral equation (44), for , we get

where is the reduced differential transform of . After expanding the RDTM recurrence equations (47), with initial value of (46), for , the first five terms of are obtained as follows:

In the same manner, the rest of the components were obtained by using the recursive equations (47). Substituting the quantities (48) in (18), the approximation solution of Volterra integral equation (44) in the Poisson series form is

which is the same as the first five terms of the Poisson series of the exact solution . The numerical results obtained with reduced DTM are presented in Table 4, in comparison with the classic DTM solution of [5] and the exact solution , for some points of the intervals and .

Table 4: Comparisons of the exact solution , with obtained by classic DTM [5] and obtained by reduced DTM at some test points in Example 10.

Classic DTM [5]
Reduced DTM

0.2
0.1
0.2955202067
0.2955202184

0.2955202070

0.4
0.5646424734
0.5646439552

0.5646439183

0.7
0.7833269096
0.7833750550

0.7833749959

1
0.9320390860
0.9325020000

0.9325019239

0.5
0.1
0.5646424734
0.5646461680

0.5646424741

0.4
0.7833269096
0.7833397500

0.7833299139

0.7
0.9320390860
0.9321460859

0.9321309857

1
0.9974949866
0.9983398437

0.9983208230

0.8
0.1
0.7833269096
0.7834038828

0.7833269106

0.4
0.9320390860
0.9322215424

0.9320433813

0.7
0.9974949866
0.9978859380

0.9976224907

1
0.9738476309
0.9752880000

0.9749626963

1
0.1
0.8912073601
0.8915382743

0.8912073612

0.4
0.9854497300
0.9861662222

0.9854546786

0.7
0.9916648105
0.9928385451

0.9918098802

1
0.9092974268
0.9118055556

0.9105512242

Example 11. In the third example, consider the following two-dimensional Volterra integral equation [5]:

(a)DTM: the approximation solution of this equation is also obtained by DTM in [5] as follows

(b)RDTM: it is easy to see that the , and therefore RDTM version is

By applying the RDTM on nonlinear Volterra integral equation (50), for , we get

where is the reduced differential transform of . After expanding the RDTM recurrence equations (53), with initial value of (52), for , the first five terms of are obtain as follows:

In the same manner, the rest of the components were obtained by using the recursive equations (47). Substituting the quantities (48) in (18), the approximation solution of Volterra integral equation (44) in the Poisson series form is

which is same as the first five terms of the Poisson series of the exact solution . The numerical results obtained with reduced DTM are presented in Table 5, in comparison with the classic DTM solution of [5] and the exact solution , for some points of the intervals and .

Table 5: Comparisons of the exact solution , with obtained by classic DTM [5] and obtained by reduced DTM at some test points in Example 11.

Classic DTM [5]
Reduced DTM

0.2
0.1
−0.1001667500
−0.1001667561

−0.1001667498

0.4
+0.2013360025
+0.2013367851

+0.2013368189

0.7
+0.5210953055
+0.5211116472

+0.5211117115

1
+0.8881059822
+0.8881853333

+0.8881854339

0.5
0.1
−0.4107523258
−0.4107529453

0.4
−0.1001667500
−0.1001714167

0.7
+0.2013360025
+0.2013887643

+0.2014033478

1
+0.5210953055
+0.5215820313

+0.5216052465

0.8
0.1
−0.7585837018
−0.7585783977

−0.7585837006

0.4
−0.4107523258
−0.4108535808

−0.4107476947

0.7
−0.1001667500
−0.1002690359

−0.1000423588

1
+0.2013360025
+0.2019546667

+0.2023226727

1
0.1
−1.0265167257
−1.0264561563

−1.0265167241

0.4
−0.6366535821
−0.6370106667

−0.6366473802

0.7
−0.3045202934
−0.3051720521

−0.3043519610

1
+0.0000000000
−0.0000000000
−0.0000000000
+0.0013512390

4. Conclusions

In this study, we presented the definition and operation of both two-dimensional differential transformation method (DTM) and their reduced form, the so-called reduced-DTM (RDTM) for finding the solutions of a class of Volterra integral equations. For illustration purposes, we consider three different examples. It is worth pointing out that both DTM and RDTM have convergence for the solutions; actually, the accuracy of the series solution increases when the number of terms in the series solution is increased. From the computational process of DTM and RDTM, we find that the RDTM is easier to apply. In other words, it is obvious that DTM has very complicated computational process rather than RDTM. The RDTM reduces the computational difficulties of the DTM and all the calculations can be made with simple manipulations MATLAB. Actually, as a special advantage of RDTM rather than DTM, the reduced differential transform recursive equations produce exactly all the Poisson series coefficients of solutions, whereas the differential transform recursive equations produce exactly all the Taylor series coefficients of solutions. The reliability of the RDTM and the reduction in the size of computational domain give this method a wider applicability. For small value of , , in Tables 3, 4, and 5, we find that the RDTM has a smaller error than DTM. Also, for large values of , , we may increase the accuracy of the series solution by computing more terms, which is quite easy using MATLAB.

References

H. Brunner, Collocation Methods for Volterra Integral and Related Functional Differential Equations, vol. 15 of Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, Cambridge, UK, 2004.

P. Huabsomboon, B. Novaprateep, and H. Kaneko, “On Taylor-series expansion methods for the second kind integral equations,” Journal of Computational and Applied Mathematics, vol. 234, no. 5, pp. 1466–1472, 2010.

H. Kaneko and Y. Xu, “Gauss-type quadratures for weakly singular integrals and their application to Fredholm integral equations of the second kind,” Mathematics of Computation, vol. 62, no. 206, pp. 739–753, 1994.

A. Tari, M. Y. Rahimi, S. Shahmorad, and F. Talati, “Solving a class of two-dimensional linear and nonlinear Volterra integral equations by the differential transform method,” Journal of Computational and Applied Mathematics, vol. 228, no. 1, pp. 70–76, 2009.

B. Jang, “Comments on “Solving a class of two-dimensional linear and nonlinear Volterra integral equations by the differential transform method”,” Journal of Computational and Applied Mathematics, vol. 233, no. 2, pp. 224–230, 2009.

Y.-J. jiang, “On spectral methods for Volterra-type integro-differential equations,” Journal of Computational and Applied Mathematics, vol. 230, no. 2, pp. 333–340, 2009.

A. Akyüz-Daşcıoğlu and M. Sezer, “Chebyshev polynomial solutions of systems of higher-order linear Fredholm-Volterra integro-differential equations,” Journal of the Franklin Institute. Engineering and Applied Mathematics, vol. 342, no. 6, pp. 688–701, 2005.

S. M. Hosseini and S. Shahmorad, “Numerical solution of a class of integro-differential equations by the tau method with an error estimation,” Applied Mathematics and Computation, vol. 136, no. 2-3, pp. 559–570, 2003.

M. T. Kajani, M. Ghasemi, and E. Babolian, “Numerical solution of linear integro-differential equation by using sine-cosine wavelets,” Applied Mathematics and Computation, vol. 180, no. 2, pp. 569–574, 2006.

R. Farnoosh and M. Ebrahimi, “Monte Carlo method for solving Fredholm integral equations of the second kind,” Applied Mathematics and Computation, vol. 195, no. 1, pp. 309–315, 2008.

K. Maleknejad and F. Mirzaee, “Numerical solution of integro-differential equations by using rationalized Haar functions method,” Kybernetes, vol. 35, no. 10, pp. 1735–1744, 2006.

K. E. Atkinson, The Numerical Solution of Integral Equations of the Second Kind, vol. 4 of Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, Cambridge, UK, 1997.

H. Guoqiang and W. Jiong, “Extrapolation of Nystrom solution for two dimensional nonlinear Fredholm integral equations,” Journal of Computational and Applied Mathematics, vol. 134, no. 1-2, pp. 259–268, 2001.

G. Han and R. Wang, “Richardson extrapolation of iterated discrete Galerkin solution for two-dimensional Fredholm integral equations,” Journal of Computational and Applied Mathematics, vol. 139, no. 1, pp. 49–63, 2002.

J. K. Zhou, Differential Transformation and Its Application for Electrical CircuIts, Huazhong University Press, Wuhan, China, 1986.

R. Abazari and A. Borhanifar, “Numerical study of the solution of the Burgers and coupled Burgers equations by a differential transformation method,” Computers & Mathematics with Applications, vol. 59, no. 8, pp. 2711–2722, 2010.

A. Borhanifar and R. Abazari, “Exact solutions for non-linear Schrödinger equations by differential transformation method,” Journal of Applied Mathematics and Computing, vol. 35, no. 1-2, pp. 37–51, 2011.

A. Borhanifar and R. Abazari, “Numerical study of nonlinear Schrödinger and coupled Schrödinger equations by differential transformation method,” Optics Communications, vol. 283, no. 10, pp. 2026–2031, 2010.

R. Abazari and M. Abazari, “Numerical simulation of generalized Hirota-Satsuma coupled KdV equation by RDTM and comparison with DTM,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 2, pp. 619–629, 2012.

R. Abazari, “Solution of Riccati types matrix differential equations using matrix differential transform method,” Journal of Applied Mathematics & Informatics, vol. 27, pp. 1133–1143, 2009.

R. Abazari and R. Abazari, “Numerical study of some coupled PDEs by using differential transformation method,” World Academy of Science, Engineering and Technology, vol. 66, pp. 52–59, 2010.

Y. Keskin and G. Oturanç, “Reduced differential transform method for partial differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 6, pp. 741–749, 2009.

Y. Keskin and G. Oturanç, “The reduced differential transform method: a new approach to factional partial differential equations,” Nonlinear Science Letters A, vol. 1, pp. 207–217, 2010.

R. Abazari and M. Ganji, “Extended two-dimensional DTM and its application on nonlinear PDEs with proportional delay,” International Journal of Computer Mathematics, vol. 88, no. 8, pp. 1749–1762, 2011.