# Solving Systems of Volterra Integral and Integrodifferential Equations with Proportional Delays by Differential Transformation Method

Authors: Şuayip Yüzbaşı; Nurbol Ismailov
Journal: Journal of Mathematics
 (2014)

---

## Abstract

In this paper, the differential transformation method is applied to the system of Volterra integral and integrodifferential equations with proportional delays. The method is useful for both linear and nonlinear equations. By using this method, the solutions are obtained in series forms. If the solutions of the problem can be expanded to Taylor series, then the method gives opportunity to determine the coefficients of Taylor series. Hence, the exact solution can be obtained in Taylor series form. In illustrative examples, the method is applied to a few types of systems.


---

## Body

## 1. Introduction

Integral and integrodifferential equations have found applications in engineering, physics, chemistry, and insurance mathematics [1–3]. In particular, functional-differential equations with proportional delays have described some models such as motion of particle in liquid and polymer crystallization which can be found in [4].There are a lot of methods of approach for solutions of systems of integral and integrodifferential equations. For example, the linear and nonlinear systems of integrodifferential equations have been solved by Haar functions [5]; Maleknejad and Tavassoli Kajani [6] used the hybrid Legendre functions, the Chebyshev polynomial method [7], the Bessel collocation method [8, 9], the Taylor collocation method [10], the homotopy perturbation method [11, 12], the variational iteration method [13], the differential transformation method [14], and the Taylor series method [15]. Biazar et al. [16] have obtained the solutions of systems of Volterra integral equations of the first kind by the Adomian method. In addition, the homotopy perturbation method has been used for systems of Abel’s integral equations [17]. On the other hand, the special systems of integral equations have been solved by the differential transformation method [18]. Katani and Shahmorad [19] have presented Romberg quadrature for the systems of Urysohn type Volterra integral equations. The nonlinear systems of Volterra integrodifferential equations with delay arguments have been studied by Yalçınbaş and Erdem [20].In this paper, we consider the system of Volterra integral and integrodifferential equations with proportional delays:(1)Fk(x,yk(l)(pkx),∫0rx‍Gk(x,t,yk(qkt))dt)=0,
where Fk, Gk are given functions, pk,qk,r∈(0,1], k=1,2,…,n, and l=0,1,…,m.


## 2. Differential Transformation Method

In 1987, the differential transformation method is introduced by Zhou [21] in the study of electric circuits. The method based on Taylor series and yields of differential transformation are difference equations which solutions give the exact values of derivatives of origin function at the given point. The method has been used for a wide class of problems [22–25]. The main advantage of differential transformation from Laplace and Fourier transformations is that it can be applied easily to linear equations with constant and variable coefficients and some nonlinear equations.The differential transformation of thekth derivative of function y(x) is defined by
(2)Y(k)=1k![dky(x)dxk]x=x0
and the inverse transformation is defined as follows:
(3)y(x)=∑k=0∞‍Y(k)(x-x0)k.The following theorems can be obtained from definitions (2) and (3).Theorem 1.
Assume thatF(k), G(k), H(k) are the differential transformations, at the x0=0, of the functions f(x), g(x), h(x), respectively; then one has the following.Iff(x)=dng(x)/dxn, then F(k)=((k+n)!/k!)G(k+n).Iff(x)=g(x)h(x), then F(k)=∑l=0k‍G(l)H(k-l).Iff(x)=xn, then F(k)=δ(k-n),δ is the Kronecker delta symbol.Iff(x)=eλx, then F(k)=λk/k!.Iff(x)=sin(ωx+λ), then F(k)=(ωk/k!)sin(kπ/2+λ).Iff(x)=cos⁡(ωx+λ), then F(k)=(ωk/k!)cos⁡(kπ/2+λ).Iff(x)=g(x)∫0x‍h(t)dt, then F(k)=∑l=0k-1‍(1/(k-l))G(l)H(k-l-1).Theorem 2.
Assume thatW(k),Y(k), and Yi(k), i=1,…,n, are the differential transformations, at the x0=0, of the functions w(x),y(x), and yi(x), respectively, and q,qi,r∈(0,1], i=1,2. Then, one has the following.Ifw(x)=y(qx), then W(k)=qkY(k).Ifw(x)=y1(q1x)y2(q2x), then W(k)=∑l=0k‍q1lq2k-lY1(l)Y2(k-l).Ifw(x)=dmy(qx)/dxm, then W(k)=((k+m)!/k!)qk+mY(k+m).Ifw(x)=∫0rx‍y(qt)dt, then W(k)=(rk/k)qk-1Y(k-1).Ifw(x)=∫0x‍y1(q1t)y2(q2t)dt, then W(k)=(1/k)∑l=0k-1‍q1lq2k-l-1Y1(l)Y2(k-l-1).Ifw(x)=y(qx)∫0x‍y1(q1t)y2(q2t)dt, then
(4)W(k)=∑l=0k-1‍∑s=0k-l-1‍1k-lqlq1sq2k-l-s-1Y(l)×Y1(s)Y2(k-l-s-1),where k∈N.The proofs of Theorems1 and 2 are given in [22, 25].


## 3. Illustrate Examples

Example 1.
Let us consider the following linear system of Volterra integrodifferential equations with proportional delays and separable kernels:(5)y1(x)-y1′(x3)+9y2′(x3)=4cos⁡x+8∫0x‍cos⁡(x-t)y2(t3)dt+19∫0x‍y1(t3)dty2(x)+y1′(x)-4y2′(x)=-2cos⁡x2-2+3∫0x‍sin(x-t)y1(t2)dt+∫0x‍y2(t2)dt,
with the initial conditions y1(0)=1 and y2(0)=0.The differential transformation of the last system is(6)Y1(k)-k+13k+1Y1(k+1)+9k+13k+1Y2(k+1)=4k!cos⁡πk2+8∑l=0k-1∑s=0k-l-1Y2(k-l-s-1)3k-l-s-1(k-l)l!s!cos⁡(π2[l-s])+Y1(k-1)3k+1k,Y2(k)+(k+1)Y1(k+1)-4(k+1)Y2(k+1)=-12k-1k!cos⁡πk2-2δ(k)+3∑l=0k-1∑s=0k-l-1Y2(k-l-s-1)2k-l-s-1(k-l)l!s!sin(π2[l-s])+Y2(k-1)2k-1k.Substitutingx=0 in Example 1, we obtain values of first derivatives of unknown functions; that is, Y1(1)=0 and Y2(1)=1.Fork=1 in (6), we have the following system:
(7)Y1(1)-29Y1(2)+2Y2(2)=8Y2(0)+19Y1(0),Y2(1)+2Y1(2)-8Y2(2)=Y2(0).Solving the last system, we getY1(2)=-1/2 and Y2(2)=0.Substitutingk=2 in (6), we obtain the following system of equations:
(8)19Y1(3)-Y2(3)-23=-12,-3Y1(3)+12Y2(3)+2=0.Solving the last system with two unknown, we haveY1(3)=0 and Y2(3)=-1/6.Continue this process and use inverse transformation; we gety1(x)=cos⁡x and y2(x)=sinx which are the exact solutions of Example 1.Example 2.
Consider the following system of nonlinear Volterra integrodifferential equations with proportional delays:(9)y1(x)+y2′(x)=x-14x2+∫0x‍y1(t2)dt,y2(x)+y1′(x)=ex-3(x-3)ex/3+∫0x‍y1(t)y2(t3)dt,
with the initial conditions y1(0)=0 and y2(0)=1.Analogously, forx=0 in (9), we get y1′(0)=1 and y2′(0)=1.Applying the differential transformation to (9), we have the following system of difference equations:
(10)Y1(k)+(k+1)Y2(k+1)=1k!+δ(k-1)-14δ(k-2)+1k2k-1Y1(k-1),Y2(k)+(k+1)Y1(k+1)=1k!-3∑l=0k‍δ(l-1)13k-l(k-l)!+13k-2k!+1k∑l=0k-1‍Y1(l)Y2(k-l-1)3k-l-1.Fork=1 in (10), we get the following system:
(11)Y1(1)+2Y2(2)=2+Y1(0),Y2(1)+2Y1(2)=1.Solving the last system, we obtainY1(2)=0 and Y2(2)=1/2=1/2!.Substitutingk=2 in (10) and solving corresponding system, we have Y1(3)=0 and Y2(3)=1/6=1/3!, and for k≥3, Y1(k)=0 and Y2(k)=1/k!. Then using (3), we gain y1(x)=x and y2(x)=ex which are the exact solutions of system (9).Example 3.
Consider the following system of nonlinear Volterra integral equations:(12)y1(x)=x3324+89x290-∫0x/3‍[y1(t2)-y2(t5)]dt,y2(x)=x+8x327+∫0x‍[y1(t3)-y22(t)]dt,
with the initial condition y1(0)=y2(0)=0.Now, applying the differential transformation to (12), we get the following system of difference equations:
(13)Y1(k)=δ(k-3)324+89δ(k-2)90-Y1(k-1)3k2k-1k+Y2(k-1)3k5k-1k,Y2(k)=δ(k-1)+8δ(k-3)27+Y1(k-1)3k-1k-1k∑l=0k-1‍Y2(l)Y2(k-l-1).From initial conditions, we haveY1(0)=Y2(0)=0.Using (13), we get the following values:
(14)Y1(1)=0,Y1(2)=1,Y2(1)=1,Y2(2)=0,
and for k≥3, Y1(k)=Y2(k)=0.Using (3), we have y1(x)=x and y2(x)=x2 which are exact solutions of system (12).Example 4.
In last we consider the linear system with variable coefficients of Volterra integral equations with proportional delays:(15)x2y1(x4)+2sinxy2(x3)=sin4x3+sin2x3+2∫0x/2‍xy1(t)dt,sinx3y1(x)+xy2(x2)=xcos⁡x2+∫0x/3‍xy2(t)dt,
with exact solutions y1(x)=x and y2(x)=cos⁡x.Applying DTM, we have(16)∑l=0k‍δ(l-2)Y1(k-l)14k-l+2∑l=0k‍1l!3k-lsinlπ2Y2(k-l)=4k3kk!sinkπ2+2k3kk!sinkπ2+2∑l=0k‍δ(l-1)Y1(k-l-1)1(k-l)2k-l,∑l=0k‍1l!3lsinlπ2Y1(k-l)+∑l=0k‍δ(l-1)Y2(k-l)12k-l=∑l=0k‍δ(l-1)1(k-l)!2k-lcos⁡(k-l2π)+∑l=0k-1‍δ(l-1)Y2(k-l-1)1(k-l)3k-l.Solving the last system we havey1(x)=x and y2(x)=1-(x2/2)+(x4/4)-⋯ which are Taylor series of exact solutions of Example 4.


## 4. Conclusions

In this study, the differential transformation method has been presented for solving system of integral and integrodifferential equations with proportional delays. The major benefits of method from integral transformations are that the method can be applied for linear equations with variable coefficients and nonlinear equations and the method gives the exact solutions in series forms.



---

*Source: 725648-2014-08-12.xml*