Rational Chebyshev bases and Galerkin method are used to obtain the approximate solution of a system of high-order integro-differential equations on the interval [0,∞). This method is based on replacement of the unknown functions by their truncated series of rational Chebyshev expansion. Test examples are considered to show the high accuracy, simplicity, and efficiency of this method.
1. Introduction
In recent years, there has been a growing interest in the system of integrodifferential equations (IDE), which arise frequently in many applied areas which include engineering, mechanics, physics, chemistry, astronomy, biology, economics, potential theory, electrostatics, and so forth [1–8]. The systems of integrodifferential equations are generally difficult to solve analytically, thus finding efficient computational algorithms for obtaining numerical solution is required.
There are various techniques for solving systems of IDE, for example, operational Tau method [9, 10], Adomian decomposition method [11], Galerkin method [12], rationalized Haar functions method [13], He’s homotopy perturbation method (HPM) [14, 15], and Ghebyshev polynomial [16].
A number of problems arising in science and engineering are set in semi-infinite domains. One can apply different spectral methods that are used to solve problems in semi-infinite domains. The first approach is using Laguerre polynomials [17–20]. The second approach is replacing semi-infinite domain with [0,L] interval by choosing L, sufficiently large. This method is named domain truncation [21]. The third approach is reformulating original problem in semi-infinite domain to singular problem in bounded domain by variable transformation and then using the Jacobi polynomials to approximate the resulting singular problem [22]. The fourth approach of spectral method is based on rational orthogonal functions.
Boyd [23] defined a new spectral basis, named rational Chebyshev functions on the semi-infinite interval, by mapping to the Chebyshev polynomials. Guo et al. [24] introduced a new set of rational Legendre functions which are mutually orthogonal in L2(0,+∞). They applied a spectral scheme using the rational Legendre functions for solving the Korteweg-de Vries equation on the half line. Boyd et al. [25] applied pseudospectral methods on a semi-infinite interval and compared rational Chebyshev, Laguerre, and mapped Fourier sine.
The authors of [26–29] applied spectral method to solve nonlinear ordinary differential equations on semi-infinite intervals. Their approach was based on a rational Tau method. They obtained the operational matrices of derivative and product of rational Chebyshev and Legendre functions and then they applied these matrices together with the Tau method to reduce the solution of these problems to the solution of system of algebraic equations.
Zarebnia and Ali Abadi [30] used Sinc-Collocation method for solving system of nonlinear second-order integrodifferential equations of the Fredholm type. Rational second (third) kind Chebyshev (RSC) functions, for the first time, were proposed by Tavassoli Kajani and Ghasemi Tabatabaei [31] to find the numerical solution of Lane-Emden equation.
This paper outlines the application of rational second kind Chebyshev functions and Galerkin method to the following system of linear high-order integrodifferential equations on the interval [0,∞). Two problems of such equations are solved to make clear the application of the proposed method. One has
(1.1)∑i=1l(∑j=0mνpij(x)yi(j)(x)+λip∫ab(kip(x,t)yi(t))dt)=fp(x),p=1,2,…,l,yi(j)(0)=yij,i=1,2,…,l,j=0,1,…,m-1,x∈[0,∞),
where 0≤a<b<∞.
2. Properties of RSC Functions
In this section, we present some properties of RSC functions.
2.1. RSC Functions
The second kind Chebyshev polynomials Un(x),n≥0, are orthogonal in the interval [-1,1] with respect to the weight function 1-x2 and we find that Un(x) satisfies the recurrence relation [32]
(2.1)U0(x)=1,U1(x)=2x,Un(x)=2xUn-1(x)-Un-2(x),n≥2.
The RSC functions are defined by [31, 33]
(2.2)Rn(x)=Un(x-1x+1),
thus RSC functions satisfy
(2.3)R0(x)=1,R1(x)=2(x-1x+1),Rn(x)=2(x-1x+1)Rn-1(x)-Rn-2(x),n≥2.
2.2. Function Approximation
Let w(x)=4x/(x+1)3 denotes a nonnegative, integrable, real-valued function over the interval I=[0,+∞). We define
(2.4)Lw2(I)={y:I→ℝ∣yismeasurableand∥y∥w<∞},
where
(2.5)∥y∥w=(∫0∞|y(x)|2w(x)dx)1/2
is the norm induced by the scalar product
(2.6)〈y,z〉w=∫0∞y(x)z(x)w(x)dx.
Thus {Rn(x)}n≥0 denote a system which is mutually orthogonal under (2.6), that is,
(2.7)∫0∞Rn(x)Rm(x)w(x)dx=π2δnm,
where δnm is the Kronecker delta function. This system is complete in Lw2(I); as a result, any function y∈Lw2(I) can be expanded as follows:
(2.8)y(x)=∑k=0∞ykRk(x),
with
(2.9)yk=2π〈y,Rk〉w.
The yk’s are the expansion coefficients associated with the family {Rk(x)}. If the infinite series in (2.8) is truncated, then it can be written as
(2.10)y(x)≃yN(x)=∑k=0NykRk(x)=YTR(x),
where Y=[y0,y1,…,yN]T and R(x)=[R0(x),R1(x),…,RN(x)]T.
We can also approximate the function k(x,t) in Lw2(I×I) as follows:
(2.11)k(x,t)≃kN(x,t)=RT(x)KR(t),
where K is an (N+1)×(N+1) matrix that
(2.12)Kij=2π2〈Ri(x),〈k(x,t),Rj(t)〉〉,i,j=0,1,…,N.
Moreover, from recurrence relation in (2.3) we have
(2.13)R(0)=[1,-2,3,-4,…,(-1)N(N+1)]T=e1.
2.3. Product Integration of the RSC Functions
We also use the matrix Pab as follows:
(2.14)Pab=∫abR(t)RT(t)dt.
To illustrate the calculation Pab we choose a=0 and b=1, then we obtain
(2.15)P01=[12-4ln211-16ln228-40ln21673-80ln2⋯2-4ln212-16ln230-44ln22003-96ln23743-180ln2⋯11-16ln230-44ln22033-96ln23803-184ln2332915-320ln2⋯28-40ln22003-96ln23803-184ln2334415-320ln2539615-520ln2⋯1673-80ln23743-180ln2332915-320ln2539615-520ln21945135-800ln2⋯⋮⋮⋮⋮⋮⋱].
2.4. Operational Matrix of Derivative
The derivative of the vector R(x) defined in (2.10) can be approximated by
(2.16)R′(x)≃DR(x),
where D is named the n×n operational matrix of derivative. Differentiating (2.3) we get
(2.17)R0′(x)=0,R1′(x)=54R0(x)-R1(x)+14R2(x),Rn′(x)=(R1(x)Rn-1(x))′-Rn-2′(x),n≥2.
By using (2.17) the matrix D can be calculated. The matrix D is a lower Hessenberg matrix and can be expressed as D=D1+D2, where D1 is a tridiagonal matrix which is obtained from
(2.18)D1=diag(-2+7i4,-i,i4),i=0,…,n-1,
and the dij elements of matrix D2 are obtained from
(2.19)dij={0,i≤j+1,(-1)i+j+1(2j),i>j+1.
2.5. The Product Operational Matrix
The following property of the product of two rational Chebyshev vectors will also be used:
(2.20)R(x)RT(x)Y=Y~R(x),
where Y~ is called (N+1)×(N+1) product operational matrix for the vector Y. Using (2.20) and the orthogonal property, the elements Y~ij,i=0,…,N,j=0,…,N of the matrix Y~ can be calculated from
(2.21)Y~ij=2π∑k=0Nckgijk,
where gijk is given by
(2.22)gijk=∫0∞Ri(x)Rj(x)Rk(x)w(x)dx.
3. Solving System of Integrodifferential Equations over Semi-Infinite Interval
Consider the following system of integrodifferential equations:
(3.1)∑i=1l(∑j=0mνpij(x)yi(j)(x)+λip∫ab(kip(x,t)yi(t))dt)=fp(x),p=1,2,…,l,yi(j)(0)=yij,i=1,2,…,l,j=0,1,…,m-1,x∈[0,∞).
Using (2.10) and (2.11) to approximate yi, fp, kip, and νpij when i,p=1,2,…,l and j=0,…,m, we have
(3.2)yi(x)≃YiTR(x),fp(x)≃FpTR(x),kip(x,t)≃RT(t)KipR(x),νpij(x)≃VpijTR(x).
According to the operational matrix of derivative we can approximate yi(j) as
(3.3)yi(j)(x)≃YiTR(j)(x)≃YiTDjR(x),yi(j)(0)≃YiTDjR(0)=YiTDje1.
With substituting these approximations in (3.1) we have
(3.4)∑i=1l(∑j=0mYiTDjR(x)RT(x)Vpij+λip∫abYiTR(t)RT(t)KipR(x)dt)=FpTR(x),YiTDje1=yij,i=1,2,…,l,j=0,1,…,m-1,p=1,2,…,l.
Then using (2.20) we obtain
(3.5)∑i=1l(∑j=0mYiTDjV~pijR(x)+λipYiT(∫abR(t)RT(t)dt)KipR(x))=FpTR(x),YiTDje1=yij,i=1,2,…,l,j=0,1,…,m-1,p=1,2,…,l,
which can be simplified using (2.14)
(3.6)∑i=1l(∑j=0mYiTDjV~pij+λipYiTPabKip)=FpT,p=1,2,…,l,YiTDje1=yij,i=1,2,…,l,j=0,1,…,n-1.
By solving this linear system of algebraic equations we can find vectors Yi, i=1,2,…,n, and then approximate the solutions yi(x) as
(3.7)yi(x)≃YiTR(x).
4. Numerical ExamplesExample 4.1.
Consider the following system of linear integrodifferential equations:
(4.1)1x+1y1′′(x)+y1′(x)+y2′′(x)+∫01y1(t)+y2(t)(x+1)2(t+1)2dt=x-154(x+1)3,y1′(x)-2y2′(x)+y2′′(x)+24∫01(y1(t)(x+1)2(t+1)2+y2(t)(x+1)3(t+1)3)dt=8x-1(x+1)3,y1(0)=1,y1′(0)=0,y2(0)=-1,y2′(0)=2.The exact solution of this example is y1(x)=1 and y2(x)=(x-1)/(x+1).
We solved Example 4.1 using the present method with N=3, and we obtained Y1T=[1,0,0,0] and Y2T=[0,0.5,0,0], which imply
(4.2)y1(x)=Y1TR(x)=1,y2(x)=Y2TR(x)=x-1x+1
that are the exact solutions.
Example 4.2.
Next, consider the following system of integrodifferential equations with the exact solution y1(x)=e-x and y2(x)=e-2x:
(4.3)y1′′(x)+y1′(x)+y2′′(x)+2y2′(x)+∫016e-t-x(y1(t)+y2(t))dt=(5-3e-2-2e-3)e-x,4e-xy1′′(x)-y2′′(x)+∫0112e-2(t+x)(y1(t)+y2(t))dt=(5-4e-3-3e-4)e-2x,y1(0)=1,y1′(0)=-1,y2(0)=1,y2′(0)=-2.
We solved this example by using the method described in Section 3 for N=11 and N=13. Results are shown in Table 1 and Figures 1 and 2. The errors for large values of x are shown in Table 2. It is seen that the proposed method provides accurate results even for large values of x.
Numerical results of Example 4.2.
x
y1(x)
y2(x)
N=11
N=13
Exact
N=11
N=13
Exact
0.0
1.00000
1.00000
1.00000
1.00000
1.00000
1.00000
0.2
0.82435
0.81643
0.81873
0.67203
0.66964
0.67032
0.4
0.67955
0.66654
0.67032
0.45203
0.44825
0.44933
0.6
0.56024
0.54411
0.54881
0.30436
0.29994
0.30119
0.8
0.46199
0.44410
0.44933
0.20516
0.20062
0.20190
1.0
0.38119
0.36238
0.36788
0.13844
0.13413
0.13534
1.2
0.31473
0.29559
0.30119
0.09351
0.08965
0.09072
1.4
0.26007
0.24100
0.24660
0.06321
0.05991
0.06081
1.6
0.21512
0.19637
0.20190
0.04274
0.04006
0.04076
1.8
0.17819
0.15987
0.16530
0.02886
0.02681
0.02732
2.0
0.14787
0.13003
0.13534
0.01942
0.01798
0.01832
2.2
0.12300
0.10562
0.11080
0.01297
0.01211
0.01228
2.4
0.10261
0.08565
0.09072
0.00855
0.00822
0.00823
2.6
0.08589
0.06932
0.07427
0.00550
0.00565
0.00552
2.8
0.07218
0.05596
0.06081
0.00339
0.00395
0.00370
3.0
0.06094
0.04503
0.04970
0.00193
0.00285
0.00248
Absolute error for large values of x.
N=11
x=50
x=100
x=200
x=400
e(y1(x)-Y1TR(x))
3.1×10-4
7.9×10-5
1.1×10-6
5.3×10-7
e(y2(x)-Y2TR(x))
4.5×10-3
2.4×10-3
1.2×10-4
6.3×10-5
Graph of the exact and numerical solutions of y1(x) for N=13; symbols correspond to the numerical solution.
Graph of the exact and numerical solutions of y2(x) for N=13; symbols correspond to the numerical solution.
5. Conclusion
The fundamental goal of this paper has been to construct an approximation to the solution of the integrodifferential equations system in a semi-infinite interval. In the above discussion, the Galerkin method with RSC functions, which have the property of orthogonality, is employed to achieve this goal. Advantages of this method is that we do not reform the problem to a finite domain, and with a small value of N accurate results are obtained. There is a good agreement between obtained results and exact values that demonstrates the validity of the present method for this type of problems and gives the method a wider applicability.
Abdul JerriJ.1999New York, NY, USAJohn Wiley & SonsMR1800272BoydJ. P.Evaluating of Dawson's integral by solving its differential equation using orthogonal rational Chebyshev functions2008204291491910.1016/j.amc.2008.07.039MR2460332ZBL1157.65331LinzP.19857Philadelphia, Pa, USASociety for Industrial and Applied Mathematics10.1137/1.9781611970852MR796318MaleknejadK.MahmoudiY.Numerical solution of linear Fredholm integral equation by using hybrid Taylor and block-pulse functions2004149379980610.1016/S0096-3003(03)00180-2MR2033162ZBL1038.65147RenY.ZhangB.QiaoH.A simple Taylor-series expansion method for a class of second kind integral equations19991101152410.1016/S0377-0427(99)00192-2MR1715520ZBL0936.65146RashedM. T.Numerical solution of functional differential, integral and integro-differential equations2004156248549210.1016/j.amc.2003.08.021MR2087525ZBL1061.65146WangW.An algorithm for solving the high-order nonlinear Volterra-Fredholm integro-differential equation with mechanization2006172112310.1016/j.amc.2005.01.116MR2197883ZBL1088.65118YalçınbaşS.SezerM.The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials20001122-329130810.1016/S0096-3003(99)00059-4MR1758715ZBL1023.65147AbbasbandyS.TaatiA.Numerical solution of the system of nonlinear Volterra integro-differential equations with nonlinear differential part by the operational Tau method and error estimation2009231110611310.1016/j.cam.2009.02.014MR2532653ZBL1170.65101EbadiG.RahimiM. Y.ShahmoradS.Numerical solution of the system of nonlinear Fredholm integro-differential equations by the operational tau method with an error estimation2007146546554MR2374110ZBL1178.65146KhaniA.Mohseni MoghadamM.ShahmoradS.Numerical solution of special class of systems of non-linear Volterra integro-differential equations by a simple high accuracy method2008342141152MR2477999MaleknejadK.Tavassoli KajaniM.Solving linear integro-differential equation system by Galerkin methods with hydrid functions2004159360361210.1016/j.amc.2003.10.046MR2097737MaleknejadK.MirzaeeF.AbbasbandyS.Solving linear integro-differential equations system by using rationalized Haar functions method2004155231732810.1016/S0096-3003(03)00778-1MR2077051ZBL1056.65144BiazarJ.GhazviniH.EslamiM.He's homotopy perturbation method for systems of integro-differential equations20093931253125810.1016/j.chaos.2007.06.001MR2512929ZBL1197.65106YusufoğluE.An efficient algorithm for solving integro-differential equations system20071921515510.1016/j.amc.2007.02.134MR2385568ZBL1193.65234Akyüz-DaşcıoğluA.SezerM.Chebyshev polynomial solutions of systems of higher-order linear Fredholm-Volterra integro-differential equations2005342668870110.1016/j.jfranklin.2005.04.001MR2166749ZBL1086.65121GuoB.-Y.ShenJ.Laguerre-Galerkin method for nonlinear partial differential equations on a semi-infinite interval200086463565410.1007/PL00005413MR1794346ZBL0969.65094MadayY.Pernaud-ThomasB.VandevenH.Shock-Fitting techniques for solving hyperbolic problems with spectral methods19856192-s2.0-0022954035ZBL0608.76053SiyyamH. I.Laguerre tau methods for solving higher-order ordinary differential equations20013217318210.1023/A:1010141309991ZBL1024.65049ShenJ.Stable and efficient spectral methods in unbounded domains using Laguerre functions20003841113113310.1137/S0036142999362936MR1786133ZBL0979.65105BoydJ. P.20002ndBerlin, GermanySpringerMR1874071GuoB.-Y.Jacobi spectral approximations to differential equations on the half line200018195112MR1741177ZBL0948.65071BoydJ. P.Orthogonal rational functions on a semi-infinite interval1987701638810.1016/0021-9991(87)90002-7MR888932ZBL0614.42013GuoB.-Y.ShenJ.WangZ.-Q.A rational approximation and its applications to differential equations on the half line200015211714710.1023/A:1007698525506MR1827574ZBL0984.65104BoydJ. P.RanganC.BucksbaumP. H.Pseudospectral methods on a semi-infinite interval with application to the hydrogen atom: a comparison of the mapped Fourier-sine method with Laguerre series and rational Chebyshev expansions2003188156742-s2.0-003834644510.1016/S0021-9991(03)00127-XZBL1028.65086ParandK.RazzaghiM.Rational Chebyshev tau method for solving higher-order ordinary differential equations2004811738010.1080/00207160310001606061aMR2032164ZBL1047.65052SiyyamH. I.Laguerre Tau methods for solving higher-order ordinary differential equations20013217318210.1023/A:1010141309991ParandK.RazzaghiM.Rational Chebyshev tau method for solving Volterra's population model2004149389390010.1016/j.amc.2003.09.006MR2033170ZBL1038.65149ParandK.RazzaghiM.Rational legendre approximation for solving some physical problems on semi-infinite intervals20046953533572-s2.0-244263937310.1238/Physica.Regular.069a00353ZBL1063.65146ZarebniaM.Ali AbadiM. G.Numerical solution of system of nonlinear second-order integro-differential equations201060359160110.1016/j.camwa.2010.05.005MR2665659ZBL1201.65138Tavassoli KajaniM.Ghasemi TabatabaeiF.Rational Chebyshev approximations for solving Lane-Emde equation of index mProceeding of the International Conference on Computational and Applied MathematicsMarch 2011Bangkok, Thailand840844AbramowitzM.StegunI. A.197210thNew York, NY, USADoverDadkhah TiraniM.Ghasemi TabatabaeiF.Tavassoli KajaniM.Rational second (third) kind Chebyshev approximations for solving Volterras population modelProceeding of the International Conference on Computational and Applied MathematicsMarch 2011Bangkok, Thailand835839