#### Abstract

In this research, an approximation symbolic algorithm is suggested to obtain an approximate solution of multipantograph system of type delay differential equations (DDEs) using a combination of Laplace transform and variational iteration algorithm (VIA). The corresponding convergence results are acquired, and an efficient algorithm for choosing a feasible Lagrange multiplier is designed in the solving process. The application of the Laplace variational iteration algorithm (LVIA) for the problems is clarified. With graphics and tables, LVIA approximates to a high degree of accuracy with a few numbers of iterates. Also, computational results of the considered examples imply that LVIA is accurate, simple, and appropriate for solving a system of multipantograph delay differential equations (SMPDDEs).

#### 1. Introduction

Many physical phenomena are formulated by delay differential equations which are similar to ordinary differential equations, but their development at assured time instant depends on past values. Pantograph equations are one of the most prominent kinds of functional differential equations with proportional delay and the pantograph kind equations have been studied extensively because of the various implementations in which these equations arise. The name of the pantograph is begun from being crafted by Ockendon and Tayler on the gathering of current by the pantograph head of an electric loco; this equation has appeared in many scientific models which is exceedingly applicable in physics, mathematics, engineering, and biology as in astrophysics, population models, probability theory, quantum mechanics, number theory, nonlinear dynamic system, electronic system, cell growth, and so forth (for further, see [1–3] and the references therein). Numerous examinations were done on the estimated arrangement of the referenced condition in the one-dimensional case. As of late, a collocation method based on the Genocchi delay operational matrix and the operational matrix of fractional derivative for solving generalized fractional pantograph equations is given in [4]. Some papers like [5, 6] have presented the solution for a system of multipantograph delay differential equations (SMPDDEs) with higher order by two different algorithms. A collocation method based on the Genocchi operational matrix for solving generalized pantograph equations is given in [7]. Then, again, the mathematical resolvability of another form of differential issues can be found in [8–18] and the references therein. The pantograph equation,is one of the important types of DDE that emerge in numerous scientific models which appear in dynamical systems, population studies, electrodynamics, and number theory. In particular, it was used by [3] to research how the electric current is composed through the pantograph of an electric locomotive, from where it gets its name. In this research, we are building on our work in [5] to develop the Laplace variational iteration algorithm (LVIA) for SMPDDEs:subject to the initial conditions:where , , are finite constants, are analytical functions such that , which verify all needed requirements for finding a unique solution, and , are functions that need to be found on the given interval. We were motivated to apply LVIA to find approximate solutions for SMPDDEs. By choosing a suitable value for the initial approximation, LVIA can easily be applied to the given problems. Moreover, the solution and its derivative are usable for each arbitrary point in the interval. LVIA provides a direct scheme for solving the problem with no physically unsuitable assumptions, discretization, linearization, transformation, or perturbation. Generally, in applications of LVIA to IVP (initial value problem) of SMPDDEs, one generally follows the following three proceedings:(a)Establishment of the correction functional(b)Identifying the Lagrange multipliers(c)Defining the initial iteration

LVIA is worth mentioning that the method is capable of decreasing the size of the computational work as compared to the classic methods while yet maintaining the elevation accuracy of the numerical results. The main feature of the variational iteration algorithm over decomposition restraint of Adomian is that the previous algorithm gives the solution of the problem without calculating Adomian’s polynomials which requires complex calculations. Accordingly, it is not affected by calculation round-off errors and one is not reputed with the requirement of large PC time and memory. Generally, by using LVIA, one iteration results in an accurate solution if the initial solution is closely chosen. The convergence of the method is a systematical debate in [8]. The remainder of this study is ordered as follows: in Section 2, we presented the basic idea of LVIA together with the analysis of the method. In Section 3, the symbolic approximate solutions for (2)-(3) can be provided using the extended LVIA. Depending on the above-mentioned, some of the numerical applications are given to explain the adequacy of LVIA in Section 3. The computations show that the approximate solutions can be achieved accurately and efficiently with a few iterations. Finally, Section 4 gives a brief conclusion.

#### 2. Methodology Basic Ideas of Laplace Variational Iteration Algorithm (LVIA)

To clarify the basic connotation of LVIA, firstly, we rewrite (2)-(3) as the next form:with the initial conditions:where , , , where , and . The idea of the VIA for (4) is to construct the next correction functional:where is the Lagrange multiplier which has a critical part in this study; it is recognized ideally by variety hypothesis and is the *n*-th approximate solution order for the exact solution , which will be gotten by utilizing starting approximation , which is acquired from (5). So, when tends to infinity, the approximate solution converges to . The numerical procedure of LVIA delineates how LVIA is utilized to approximate the solution of SMPDDEs. The basic steps included are given as follows:(i)The correction functional is acquired by taking the Laplace transform of (4):(ii) where , is used to indicate Laplace transform.(iii)To find the optimal value of , take the variation with respect to . Let (7) be stationary with respect to and ; then, we obtain(iv)Considering the terms as constrained variations, then we get(v)From (9), we derive the Lagrange multiplier as(vi)By using the inverse of Laplace transform for (7), we can obtain the succeeding approximations which giveand initial approximation can be determined by

Beginning with an initial approximation , we obtain the sequential approximations, and the exact solution can be acquired by using

Now, we display that the sequence given by (9) with converges to the exact solution of (4)-(5). To do this, we declare the following theorem.

Theorem 1 (see [9]). *Let , be the exact solution of (4) and be the solution of the sequencewith the initial approximation .*

If and , then the functional sequence defined by the above sequence converges to .

*Proof. *To prove this theorem, we follow the same proof of Theorem 1 in [9].

The accompanying considerable definitions given underneath are needed for this work and furthermore are needed for the convergence analysis in the next corollary.

Corollary 1. *Let and , be series for any iteration ; then the LVIA will attain approximately to the exact solution of (2)-(3). Now, we report some types of errors which are as follows:*(a)*Residual error defined by *(b)*Exact error (absolute error) which is defined by *(c)*Relative error , which is defined by *(d)*Consecutive error , which is defined by *

#### 3. Applications and Numerical Discussion

In this section, six numerical examples are given to illustrate the accuracy and the convergence of LVIA which is described in Section 2. In contrast, numerical results show that this method gives a good approximation to the exact solution for all possible values of , while the accuracy is in continuous increase by using only a few approximations. For comparison purposes, the solution intervals of problems are chosen generally the same as those in the references. Throughout this research, all the symbolical and computations results used Maple 18.0 software program.

*Example 1. *Regarding the SMPDDEs [19],subject to the initial conditions:where , with the exact solutions and .

We apply LVIA for solving system (15) and (16); we begin with selecting the initial conditions of the approximations such as and . Then, let us seek the approximate solutions and where is a positive integer number greater than or equal to zero. By taking for system (15)-(16), we obtainThe iteration formulas thus arewith the Lagrange multiplier , and by taking , we obtainTherefore,with the initial iteration . Now, in order to obtain the first approximation and of the LVIA solution for system (15) and (16), we put through (20) to getby using the first approximation and ; then, the second approximation and of the LVIA solution for system (15)-(16) can be written as follows:and by proceeding with the comparable style, the third and the fourth iterations for the LVIA solution of (15)-(16) lead to the following outcomes:To clarify the convergence of the solution of the proposed method to the exact solution , w.r.t. the number of iterations of the solutions, we set numerical values graphically. Figure 1(a) shows the exact solution , and , respectively. These plots detect that LVIA is an effective and appropriate method to solve SMPDDEs with fewer calculations and numbers of iterations. Figure 1(b) illustrates the efficiency of the method whether it has an exact solution or not by using the residual error values which shows that the residual error decreases as increases. Table 1 shows the CPU times of the present method for Example 1 to get 10^{th}-order approximate symbolic solutions with variable . The accuracy of the method can be spotted with the numerical results and by comparing the absolute errors of LVIA for and the Laplace decomposition algorithm (LDA) [19] is given in Table 2. It is obvious that, from this table, the approximate solutions are proven to be similar to the exact solution for all likely values of in [0,1].

**(a)**

**(b)**

*Example 2. *Consider SMPDDEs [20]:subject to the initial conditions:where with the exact solutions and . By the same procedures of Example 1, we apply the LVIA approach for solving (24)-(25), according to (11). Table 3 shows the CPU times of the LVIA for Example 2 to get 10^{th}-order approximate symbolic solutions with variable . To clarify the convergence of the solution of the proposed method to the exact solutions , with respect to the iteration number of the solution, we set numerical values graphically. Figure 2 shows the exact solution , and some iterated approximations , respectively. This graph and the results obtained in Table 4 detected that the LVIA is a very accurate and effective method to solve SMPDDEs with fewer computational and iteration steps.

*Example 3. *Consider the two-dimensional nonlinear pantograph equations [19]:the initial conditions in which the system subjected toThe exact solutions are and . By the same procedures of Example 1, we apply LVIA for solving (26)-(27), according to (11). To clarify the convergence of to , with respect to the iteration number of the solution, we present numerical results of Example 3 graphically where Figures 3(a) and 3(b) show the exact solution and the approximate solution of LVIA for , respectively. These plots clarify that LVIA is a good accurate method for solving such SMPDDEs with fewer calculations and numbers of iterations. Also, Figure 3(c) displays the absolute errors obtained by LVIA for of (26)-(27). From Figure 4(c), the approximate solutions at are identical to the exact solution for all values of in [0,1]. These plots clarify that LVIA is an efficient method for solving such a system with a few iterations. We note that the error decreases while the value increases.

**(a)**

**(b)**

**(c)**

*Example 4. *Consider the 2-dimensional nonlinear SMPDDEs [21]:subject to the initial functions:The exact solutions are and . By using the given initial functions and also using the same procedures of Example 1, we apply the LVIA approach for solving (28)-(29), according to (11). Without loss of generality, we will test the accuracy of LVIA for Example 4 using error analysis of , , for (28)-(29) with a step size of 0.2 as well as comparison among the absolute errors, relative errors, consecutive errors, and residual errors of 10^{th}-order approximate LVIA solutions shown in Table 5. Also, to illustrate the convergence of to , with respect to the order of for the solutions, we ready numerical results of Example 4 graphically where Figure 4 shows the exact solutions , and approximate solutions , and , respectively. These plots clarify that LVIA is an efficient method for solving such a system with a few iterations. We note that the error decreases while the value increases. The results show that the LVIA provides us with the precise approximate solutions of (28)-(29). Moreover, we can control the error also by evaluating more iterations.

*Example 5. *Consider 3-dimensional SMPDDEs [22]:subject to the initial conditions:where and with the exact solutions , , and . Using the LVIA method for solving (30)-(31), at first we select the initial values of the approximations such as and ; then, let us seek the approximate solutions , , and . Taking for (30)-(31), we obtainThe iteration formulas thus arewith the Lagrange multiplier , and by taking , we obtainwith the initial iterations above. Now, in order to obtain the first approximation , , and of LVIA solution for (30)-(31), we put through (34) to getby using the first approximations , and and then the second approximations , of LVIA solution for (30)-(31) which can be written in the formUsing a similar approach, we get the 4^{th}-order approximations , of the LVIA solution for (30) and (31) leading to the following results:To clarify the convergence of to , with respect to the number of iterations of the solution, we present numerical solutions of Example 5 graphically, where Figure 5 shows the exact solutions , and the approximate solution of LVIA for and , respectively. We will test the accuracy of LVIA for Example 5 by using error analysis of , , for (30)-(31) with a step size of 0.2 as well as comparison among the absolute errors, relative errors, consecutive errors residual errors, and computing times (CPU) of 4^{th}-order approximate LVIA solutions shown in Table 6.

*Example 6. *Consider 3-dimensional nonlinear SMPDDEs [23]:subject to the initial conditions:The exact solutions are , , and . By comparing between the absolute errors of LVIA for four iterations approximations and Adomian decomposition method (ADM) [23] which are shown in Table 7, it is obvious that, from this table, the approximate solutions are identical to the exact solutions for all possible values of in . Figure 6 shows that the residual error is the difference between both sides of equation (38), indicating that the two sides are identical with the increase of , which shows the efficiency of the proposed method.

#### 4. Conclusion

The main objective of this paper is to adapt Laplace variational iteration algorithm to investigate SMPDDEs. This work aims to show the power of LVIA by reducing the numerical calculation which does not need any perturbations, discretization, and other restrictive hypotheses which may change the framework for the problem being solved. LVIA gives fast convergent successive approximations. We approve that the accuracy for LVIA gives it numerous wider applicability. Additionally, the results obtained by LVIA have the ability to be successful and furthermore helpful in linear and nonlinear cases since they do not require much computational work. It is obvious from the tables and figures that the absolute errors relief as rises. In any case, the numerical results show that LVIA is a reliable technique for SMPDDEs. LVIA is a relatively new approach to provide a symbolic approximation, in a rapidly convergent series by using the initial conditions only and without tedious calculations which show the solution in a closed form, for linear and nonlinear problems and it changes SMPDDEs to the recurrence sequence of functions. Though the LVIA leads to convergent solutions for a wider region as in Examples 1 and 2, but in certain situations, the series converges very quickly in a very small zone or neighborhood of the boundary points, and it has a very slow convergence rate in the wider and/or outer zone.

#### Abbreviations

VIA: | Variational iteration algorithm |

LVIA: | Laplace variational iteration algorithm |

DDEs: | Delay different |