Solving a Quadratic Riccati Differential Equation, Multi-Pantograph Delay Differential Equations, and Optimal Control Systems with Pantograph Delays

: An effective algorithm for solving quadratic Riccati differential equation (QRDE), multi-pantograph delay differential equations (MPDDEs), and optimal control systems (OCSs) with pantograph delays is presented in this paper. This technique is based on Genocchi polynomials (GPs). The properties of Genocchi polynomials are stated, and operational matrices of derivative are constructed. A collocation method based on this operational matrix is used. The ﬁndings show that the technique is accurate and simple to use.


Introduction
Riccati differential equations (RDEs) play significant role in many fields of applied science [1]. For example, a one-dimensional static Schrödinger equation [2][3][4]. The applications of this equation found not only in random processes, optimal control, and diffusion problems [1] but also in stochastic realization theory, optimal control, network synthesis and financial mathematics. Now, RDEs attracted much attention. Recently, various iterative methods are employed for the numerical and analytical solution of functional equations such as Adomian's decomposition method (ADM) (see [5,6]), homotopy perturbation method (HPM) [7], variational iteration method (VIM) [8], and differential transform method (DTM) [9].
The GPs are non-orthogonal polynomials, which were first applied to solve fractional calculus problem (FCP) involving differential equation [10], this GPs were successfully applied to solve different kinds of problems in numerical analysis, system of Volterra integro-differential equation [11] and fractional Klein-Gordon equation [12], differential topology (differential structures on spheres), theory of modular forms (Eisenstein series).
In this paper, a new operational matrix of fractional order derivative based on Genocchi polynomials is introduced to provide approximate solutions of QRDE. Although the method is very easy to utilize and straightforward, the obtained results are satisfactory (see the numerical results).
The outline of this sequel is as follows: In Section 2, Some basic preliminaries are stated. Explanation of the problem is explained in Section 3. Some numerical results are provided in Section 4. A remark is provided about MPDDEs and OCSs. Numerical applications for solving MPDDEs are stated in Section 5. Finally, Section 6 will give a conclusion briefly.

Some Basic Preliminaries
Genocchi numbers (G n ) and Genocchi polynomials (G n (x)) have been extensively studied in various papers, (see [13]). The classical Genocchi polynomials G n (x) are usually defined by the following form where also, we have

Explanation of the Problem
Firstly, Riccati differential equation (RDE) is considered where p(x), q(x) and r(x) are continuous, x 0 , x f and α are arbitrary constants, and y(x) is unknown function. Now, the collocation method based on Genocchi operational matrix of derivatives to solve numerically RDEs is presented.
Our strategy is utilizing Genocchi polynomials (GPs) to approximate the solution y(x) by y N (x) is as given below.
then, the k-th derivative of y N (x) can be stated as by Equations (3) and (9), we have to obtain y N (x), one may use the collocation points x j = j−1 N , j = 1, 2, . . . , N − 1. These equations can be solved by Maple 15 software.
Proof. (The proof is coming in [14], but we state again here). One may set

Numerical Applications
In this section, some results are given to demonstrate the quality of the sated technique in approximating the solution of RDEs. Example 1. First, the following RDE is considered (see [15]) One may achieve y approx (x) = 0.4836486196 + 1.959259361x + 0.1873135074x 2 − 0.5349351716x 3 with this technique by n = 4. The approximate and exact solution for y(x) are shown in Figure 1. Table 1 demonstrates the absolute error of the this technique.
One may obtain y approx (x) = 1 − 4x − 2.864375404 × 10 −14 x 2 + 1.687538998 × 10 −14 x 3 with this method by n = 4. The approximate and exact solution for y(x) are shown in Figure 2. Table 2 demonstrates the absolute error of the this technique and stated technique in [15].   [15]) One may obtain y approx (x) = 0.9999999998 + 1.026769538x + 0.3911797260x 2 + 0.3003325646x 3 with this method by n = 4. The approximate and exact solutions for y(x) are shown in Figure 3. Table 3 demonstrates the absolute error of the this technique.   Remark 1. Delay differential equations (DDEs) are defined as distributed delay systems. DDEs are encountered in various practical systems such as engineering, and in the modeling of feeding system (see [16]). Many researchers used various polynomials for solving DDEs. Orthogonal functions were utilized for solving OCSs with time delay ( [17]). Also Chebyshev polynomials (ChPs) were used to solving time-varying systems with distributed time delay. The stated technique in [18] is based on expanding all time functions in terms of ChPs. The Bezier technique is utilized for solving DDEs and switched systems (see [15]). Using Bessel polynomials, pantograph equations were solved in [19].
MPDDEs are in various applications such as astrophysics, number theory, nonlinear dynamical systems (NDSs), quantum mechanics and cell growth, probability theory on algebraic structures, and etc. Properties of the analytic solution of MPDDEs as well as numerical techniques have been studied by several researchers. For example, there are treated in [20].
In this sequel, a new operational matrix of fractional order derivative based on GPs is introduced to provide approximate solutions of MPDDEs and optimal control systems with pantograph delays.
Our strategy is utilizing GPs to approximate the solution u i (x) by u N i (x) is as given below.
also G(x) satisfies in Equations (4) and (5), then, the k-th derivative of u i (x) can be stated as by Equations (8) and (9), we have to obtain u i (x), one may use the collocation points x j = j−1 N , j = 1, 2, . . . , N − 1.

Numerical Applications for Solving MPDDEs
In this section, some findings are given to demonstrate the quality of the sated technique in approximating the solution of MPDDEs and optimal control systems with pantograph delays.

Example 4.
The following time-varying system described by (see [21]) One may obtain with this method by n = 5. The approximate and exact solution for u(x) are shown in Figure 4. Table 4 demonstrates the absolute error of the this technique.
with this technique by n = 4. The approximate and exact u(x) are shown in Figure 5. Example 6. The following optimal control system with pantograph delay is considered (see [21]) One may obtain u 1,approx (x) = −1.11022302462516 × 10 −16 with this technique by n = 5. The approximate and exact u 2 (x) is shown in Figure 6. Example 7. First, the following two-dimensional pantograph equations is considered (see [22]) One may achieve with this technique by n = 3. The approximate and exact u 1,approx (x) and u 2,appox (x) are shown in Figures 7 and 8. Table 5 demonstrates the absolute error of the this technique.

Conclusions
In this paper, GPs stated for solving the RDEs, also GPs stated for solving the MPDDEs and optimal control systems with pantograph delays. The stated technique is computationally attractive. Some results are included to explain the validity of this technique. The presented approximate solutions are more accurate compared to the references as it is shown in the tables. By stated technique, the high orders of convergence obtained when it achieved accurate solutions even for small values of n.