In this research, a new approach is presented for solving delay differential equations (DDEs) which is a blend of Sumudu transform and variational iteration method (VIM). A general Lagrange multiplier is used to construct a correction functional. This is done with an uncommon Sumudu transform alongside variational theory. A few numerical cases were solved to demonstrate methodology of this new approach. Objective of this research is to reduce the complexity of computational work compared to the conventional approaches. It can be concluded that the amount of evaluation is reduced but at the same time the results are comparable as in the previous works.
1. Introduction
In recent years, especially, for several decades, applications of powerful tools for solving numerical and analytical cases have attracted attentions of scientists from all over the world [1]. Delay differential equations (DDEs) are a kind of functional differential equation having a widespread range of uses in the arena of science, technology, and engineering which acquires numerical/analytical solutions. It had already been applied in control theory for many years and recently it is being used extensively in many biological models. DDEs are form of differential equation in which the derivative of the unknown function at a particular period is provided in terms of values of the function at earlier periods. Introducing delays in models enhances the vitality of these models and allows an accurate depiction of actual occurrences. DDEs arise commonly in various physical occurrences. To be specific, they are essential once ordinary differential equation (ODE) based models are unsuccessful. In disparate ODEs, where initial conditions are stated at preliminary point, DDEs need history of the system over the delayed interval which are then provided as initial conditions. In line with this, delay systems turn out to be intricate and multifaceted in nature which in turn complicates the process of analyzing DDEs analytically and therefore requires a numerical approach. The delay differential equation in its simplest form is(1)v′it=ft,vit,viqt,in which i=1,2,3,…,n and q is a constant delay.
Application of a few numerical methods that were introduced turns out to be very useful to solve these types of equations as in numerical simulations such as the variational iteration method and the optimal perturbation method. DDE problems have constantly led to an infinite spectrum of frequencies. Hence, few approaches have been implemented in solving them such as approximation, asymptotic solutions, numerical method, and graphical approaches. In recent times, Anakira et al. employed the Optimal Homotopy Asymptotic Method (OHAM) in solving linear and nonlinear DDE [2] while Alomari et al. used the Homotopy Analysis Method (HAM) to find the solution for DDE [3]. Among all these methods applied to solve differential equations, one of them is the variational iteration method (VIM) which was firstly initiated by Ji-Huan He which could be seen throughout [4–8]. Subsequent works [9–17] reflect the flexibility, consistency, and effectiveness of the procedure in VIM. The application of VIM to differential equations generally involves discovering a correction functional, finding the Lagrange multiplier and deciding on a good initial approximation. Application of VIM has remained successful on initial as well as boundary value problems, Schrodinger equations, integrodifferential equation, fractional differential equation, chaotic Chen system, coupled sine-Gordon equation, general Riccati differential equations, fractional heat- and wave-like equations, and many more.
Various authors have made efforts in the developments of the VIM with the wide application of the method. Particularly, Wu [13, 16] used the Laplace transform to compute the Lagrange multipliers, which overcomes main shortcomings in implementation of the VIM to fractional equations. Sumudu transform is a simple variant of the Laplace transform and is essentially identical with the Laplace. Sumudu transform has many interesting properties that make it easy to visualize making it an ideal transform for control engineering and applied mathematicians. Recently, Sumudu transform is implemented in some well-known analytical approaches [19], where the coupling of homotopy perturbation method (HPM) and Sumudu transform is used to make the process of the solution simpler and improves the solution’s accurateness.
A new modified variational iteration method was found, inspired and driven by Wu’s thoughts and combining with the Sumudu transform [20]. The new approach is based on variational iteration theory and Sumudu transform. In this paper, the basic motivation is extending this new reliable approach for the solution of linear and nonlinear DDEs which are normally challenging to analyze due to their multifaceted nature and boundless dimensionality.
2. Sumudu Variational Iteration Method (SVIM)
First, let us take the general nonlinear differential equation [21] to illustrate main idea of VIM for DDEs,(2)dmvtdtm+Rvt+Nvt=ft,with the following initial conditions:(3)vk0=v0k,k=0,1,2,…,m-1,in which v=v(t), R is a linear operator, N is a nonlinear operator, f(t) is a knowncontinuous function, and dmv(t)/dtm is the term of the maximum order derivative.
Basic idea of the VIM is building a correction functional for (2) of formula (4)vn+1t=vnt+∫λt,τdmvndtm+Rv~n+Nv~n-fτdτ.The successive approximation vn,n≥1, could be attained by finding λτ, a general Lagrange multiplier that could be known optimally with variational theory. The function v~n is considered as a restricted variation indicating δv~n=0. At first, integration by parts is done to determine the Lagrange multiplier which enables the consecutive approximations, vnt of the exact solution v(t) to be attained by means of a good initial approximation v0(t). Initial conditions in (3) typically give the initial approximation. Finally, as n→∞,vnt converges to the exact solution vt.
2.1. Combination of VIM and Sumudu Transform (SVIM)
The entire procedure of Lagrange multipliers is expressed as a case of algebraic equation where its solution fx=0 could be found by (5)xn+1=xn+λfxn.Optimality condition for the extreme δxn+1/δxn=0 advances to(6)λ=-1f′xn,in which δ represents the traditional variational operator. Implementing the initial point x0 provided, the approximate solution xn+1 could be determined via next iterative scheme, with (5) and (6)(7)xn+1=xn-fxnf′xn,f′x0≠0,n=0,1,2,….The formula above (7) is the famous Newton-Raphson formula that possesses a quadratic convergence.
In this article, we outspread the idea in finding the unknown Lagrange multiplier. Key step is the application of Sumudu transform into (2) with its fundamental properties in [22]. Then, the linear equation will be converted into an algebraic equation below: (8)u-mv-s-u-mv0-…-u-1vm-10+SRv+Nv=Sft.Hence, the algorithm of the Sumudu Variational Iteration Method (SVIM) is given below:
Applying Sumudu transform to (2) gives the correction functional as (9)v-n+1u=v-nu+λuu-mv-u-u-mv0-…-u-1vm-10+SRv+Nv-Sft,with the notation used S to indicate Sumudu transform.
Considering the terms SRv~n+Nv~n-S(ft) as restricted variations, we let (9) be stationary with respect to v-n(10)δv-n+1u=δv-nu+λu1umδv-nu.From (9), we define Lagrange multiplier as (11)λu=-um.
Succeeding approximations can then be attained with the application of inverse Sumudu transform S-1 into (8) which gives(12)vn+1t=vnt-S-1umu-mv-u-u-mv0-…-u-1vm-10+SRvn+Nvn-Sft,=S-1v0+…+um-1vm-10+S-1umSRvn+Nvn-ft,
with initial approximation(13)v0t=S-1v0+…+um-1vm-10=v0+v′0t+…+tm-1vm-10m-1!.Equation (12) shows that the first iteration in traditional VIM is made up by the Taylor series.
3. Numerical Applications
Over this segment, application of the above procedure is done to some nonlinear DDEs.
Example 1.
Consider a first order nonlinear DDE (14)v′t=1-2v2t2,v0=0.The exact solution is known as vt=sint.
Taking the Sumudu transform, we obtain (15)Vuu-v0u=S1-2v2t2.The iteration formula thus is(16)Vn+1u=Vnu+λuVuu-v0u-S1-2v2t2,and its Lagrange multiplier(17)λu=-um.Applying inverse Sumudu transform gives(18)vn+1t=vnt-S-1uVuu-v0u-S1-2vn2t2,=vnt+S-1uS1-2vn2t2.Therefore,(19)vn+1t=vnt+S-1uS1-2vn2t2;with initial approximation v0t=0 and applying the iteration formula (19) above, we attain(20)v1t=t,v2t=t-t36,v3t=t-t36+t5120-t78064,v4t=t-t36+t5120-t75040+61t923224320-67t113406233600+t1312881756160-0.125146226×10-12t15.From these approximations, it can be seen that the solution tends to form the Taylor series expansion of sin(t). In order to attest numerically whether or not the suggested approach maintains the accurateness, numerical solutions of the approximation up to v4 were evaluated. The absolute values of LVIM and SVIM are compared in Table 1 while Figure 1 displays behavior of the error between these two methods in Example 1. The outcome is in good agreement with each other. Terms of sequences obtained from SVIM are computed using the Maple package.
Absolute error comparison representing Example 1.
Absolute Error
T
Laplace VIM [18]
Sumudu VIM
0.0
0
0
0.1
0.0980877
0.0980877
0.2
0.1951794
0.1951794
0.3
0.290284
0.290284
0.4
0.3824387
0.3824387
0.5
0.4707035
0.4707035
0.6
0.554168
0.554168
0.7
0.632003
0.632003
0.8
0.703398
0.703398
0.9
0.767623
0.767623
1.0
0.824018
0.824018
Graphs of the absolute errors for some values of t with LVIM and SVIM for Example 1.
Example 2.
Let us look at second order linear DDE for a second example.(21)v′′t=34vt+vt2-t2+2,v0=0,v′0=0.The exact solution is given by vt=t2.
Taking the Sumudu transform, we obtain (22)Vn+1u=Vnu+λuVuu2-v0u2-v′0u-S34vt+vt2-t2+2,and its Lagrange multiplier(23)λu=-u2.Applying inverse Sumudu transform gives(24)vn+1t=vnt+S-1u2S34vt+vt2-t2+2,with initial approximation v0t=0. Thus, from (24), we attain(25)v1t=t2-t412,v2t=t2-13t65760,v3t=t2-4.09323×10-5t8,v4t=t2-3.428788×10-7t10.The convergence of the approximations to the exact solution t2 could be observed from the solution. The absolute values of LVIM and SVIM are compared in Table 2 while Figure 2 displays behavior of error between these two methods in Example 2. The values demonstrate that the outcomes are in excellent agreement with those of the other approaches.
Absolute error comparison representing Example 2.
Absolute Error
t
Laplace VIM [18]
Sumudu VIM
0.0
0
0
0.1
0
0
0.2
0
0
0.3
0
0
0.4
0
0
0.5
3E-10
3E-10
0.6
1.6E-09
2.1E-09
0.7
7.3E-09
9.7E-09
0.8
2.78E-08
3.68E-08
0.9
9.01E-08
1.196E-07
1.0
2.585E-07
3.429E-07
Graphs of the absolute errors for some values of t with LVIM and SVIM for Example 2.
Example 3.
Lastly, a third order nonlinear DDE is considered(26)v′′′t=-1+2v2t2,v0=0,v′0=1,v′′0=0.Exact solution is known as vt=sint.
Succeeding the processes of earlier examples, we attain the following consecutive approximations:(27)v1t=t-t36+t5120,v2t=t-t36+t5120-t75040+t9362880-t1145619200+7.904061×10-11t13,v3t=t-t36+t5120-t75040+t9362880-t1139916800+1.605904×10-10t13-7.635961×10-13t15,v4t=t-t36+t5120-t75040+t9362880-t1139916800+1.605904×10-10t13-7.647106×10-13t15+2.811514×10-15t17.That converges to exact solution vt=sin(t) when n→∞. Approximate solutions obtained with SVIM could be seen to have a good agreement with the other method. The absolute values of LVIM and SVIM are compared in Table 3 while Figure 3 displays the behavior of the error among these two methods in Example 3.
Absolute error comparison representing Example 3.
Absolute Error
t
Laplace VIM [18]
Sumudu VIM
0.0
0
0
0.1
0.098088117
0.098088117
0.2
0.195178731
0.195178731
0.3
0.290284207
0.290284207
0.4
0.382437042
0.382437042
0.5
0.470699039
0.470699039
0.6
0.554170473
0.554170473
0.7
0.632000687
0.632000687
0.8
0.703394091
0.703394091
0.9
0.76761991
0.76761991
1.0
0.824018986
0.824018985
Graphs of the absolute errors for some values of t with LVIM and SVIM for Example 3.
4. Conclusion
The Sumudu transform is a simple variant of the Laplace transform and is essentially identical with the Laplace. It has many interesting properties that make it easy to visualize and hence making the process of the solution simpler. Sumudu-Lagrange multiplier is derived from Sumudu transform plus integrating with procedures of VIM to obtain the approximate solutions to DDEs. Lagrange multipliers, λ(u) which are defined in (11), could be known optimally with this different approach of variational theory. A new modification of the VIM was attained. This recommended procedure was obtained by not having applied any linearization, discretization, or impractical rules. This new technique provides more convincing or accurate sequence of results that converges quickly in physical problems. In this article, the SVIM was successfully applied in solving linear and nonlinear DDEs. In order to attest numerically whether or not the suggested approach maintains the accurateness, numerical solutions of the approximation up to v4 were evaluated. The absolute values of LVIM and SVIM are compared in the tables while Figures 1, 2, and 3 display behavior of the error between these two methods in Example 1 to Example 3. The outcomes are in good agreement with each other. It is worth stating that this method is efficacious in minimizing number of computations contrasting with traditional methods in spite of sustaining great accurateness of the numerical outcome.
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
AnakiraN. R.JameelA.AlomariA. K.SaabanA.AlmahameedM.HashimI.Approximate Solutions of Multi-Pantograph Type Delay Differential Equations Using Multistage Optimal Homotopy Asymptotic Method201850322123210.5614/j.math.fund.sci.2018.50.3.1Ratib AnakiraN.AlomariA. K.HashimI.Optimal homotopy asymptotic method for solving delay differential equations201320131149890210.1155/2013/4989022-s2.0-84889031592AlomariA. K.NooraniM. S.NazarR.Solution of delay differential equation by means of homotopy analysis method2009108239541210.1007/s10440-008-9318-zMR25514812-s2.0-70349884138HeJ.Variational iteration method—a kind of non-linear analytical technique: Some examples199934469970810.1016/S0020-7462(98)00048-1Zbl051378912-s2.0-0000092673HeJ.Variational iteration method—some recent results and new interpretations2007207131710.1016/j.cam.2006.07.009MR2332941Zbl1119.650492-s2.0-34250668369HeJ. H.WuX. H.Variational iteration method: new development and applications2007547-888189410.1016/j.camwa.2006.12.083MR2395625Zbl1141.653722-s2.0-34748870677HeJ. H.WuG. C.AustinF.The variational iteration method which should be followed201011130HerisanuN.MarincaV.A modified variational iteration method for strongly nonlinear problems201012183192OdibatZ. M.MomaniS.Application of variational iteration method to nonlinear differential equations of fractional order200671273410.1515/IJNSNS.2006.7.1.27MR31103302-s2.0-30344464250Zbl1401.65087GohS. M.NooraniM. S.HashimI.On solving the chaotic Chen system: a new time marching design for the variational iteration method using Adomian's polynomial201054224526010.1007/s11075-009-9333-9MR2646079Zbl1190.651892-s2.0-77952008913Molliq RY.NooraniM. S. M.HashimI.Variational iteration method for fractional heat- and wave-like equations2009103185418692-s2.0-6054911241410.1016/j.nonrwa.2008.02.026Zbl1172.35302NoorM. A.Mohyu-DinS. T.Variational iteration method for solving higher-order nonlinear boundary value problems using Hes polynomials200892141156WuG.Challenge in the variational iteration method—a new approach to identification of the Lagrange multipliers201325217517810.1016/j.jksus.2012.12.0022-s2.0-84875236585BatihaB.NooraniM. S. M.HashimI.Approximate analytical solution of the coupled sine-Gordon equation using the variational iteration method20077654454482-s2.0-3624893555710.1088/0031-8949/76/5/007Zbl1131.35068BatihaB.NooraniM. S. M.HashimI.Application of variational iteration method to a general Riccati equation20072562759277010.12988/imf.2007.07248Zbl1137.65372WuG.BaleanuD.Variational iteration method for fractional calculus—a universal approach by Laplace transform201311910.1186/1687-1847-2013-18MR3020186WazwazA.-M.The variational iteration method for analytic treatment of linear and nonlinear ODEs2009212112013410.1016/j.amc.2009.02.003MR2519268BialaT. A.AsimO. O.AfolabiY. O.A combination of the laplace transform and the variational iteration method for the analytical treatment of delay differential equations2014133ElbelezeA. A.KılıçmanA.TaibB. M.Homotopy perturbation method for fractional black-scholes european option pricing equations using sumudu transform20132013752485210.1155/2013/524852MR3049813LiuY.ChenW.A new iterational method for ordinary equations using sumudu transform20161210.22606/aan.2016.12004GoswamiP.AlqahtaniR. T.Solutions of fractional differential equations by Sumudu transform and variational iteration method20169419441951MR345269110.22436/jnsa.009.04.48BelgacemF. B. M.KaraballiA.Sumudu transform fundamental properties investigations and applications200620062310.1155/JAMSA/2006/9108391083MR2221002Zbl1115.440012-s2.0-33749525486