Next Article in Journal
An Efficient Analytical Technique, for The Solution of Fractional-Order Telegraph Equations
Previous Article in Journal
Strong Convergence Theorems for Variational Inequalities and Common Fixed-Point Problems Using Relaxed Mann Implicit Iteration Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solution of Ambartsumian Delay Differential Equation with Conformable Derivative

by
Sayed M. Khaled
1,2,
Essam R. El-Zahar
3,4,* and
Abdelhalim Ebaid
1
1
Department of Studies and Basic Sciences, Faculty of Community, University of Tabuk, Tabuk 71491, Saudi Arabia
2
Department of Mathematics, Faculty of Sciences, Helwan University, Cairo 11795, Egypt
3
Department of Mathematics, College of Sciences and Humanities in Al-Kharj, Prince Sattam bin Abdulaziz University, Al-kharj 11942, Saudi Arabia
4
Department of Basic Engineering Science, Faculty of Engineering, Menofia University, Shebin El-Kom 32511, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(5), 425; https://doi.org/10.3390/math7050425
Submission received: 30 March 2019 / Revised: 5 May 2019 / Accepted: 6 May 2019 / Published: 13 May 2019

Abstract

:
This paper addresses the modelling of Ambartsumian equation using the conformable derivative as an application of the theory of surface brightness in astronomy. The homotopy perturbationmethod is applied to solve this model, where the approximate solution is given in terms of the conformable derivative order and the exponential functions. The present solution reduces to the corresponding one in the relevant literature as a special case. Moreover, a rapid rate of convergence has been achieved for the obtained approximate solutions. Furthermore, the accuracy of the obtained numerical results is validated via calculating the residual against the impeded parameters. It is shown graphically that the obtained residual approaches zero in various cases, which proves the efficiency of the current analysis.

1. Introduction

In this paper, we consider a generalized model of Ambartsumian equation (AE), a delay differential equation, in the form:
D t α z ( t ) = z ( t ) + 1 q z t q , 0 < α 1 , q > 1 ,
where q is a constant for the given model and α is the order of conformable derivative (CD). Equation (1) is subjectedto the initial condition:
z ( 0 ) = λ ,
where λ is also a constant. When α 1 , Equations (1) and (2) describe the surface brightness in the Milky Way, given by [1]:
z ( t ) = z ( t ) + 1 q z t q .
Moreover, Equation (3) was derived more than 25 years earlier by Ambartsumian [1] to describe the absorption of light by the interstellar matter. The existence and uniqueness of the standard model (Equations (2) and (3)) was analyzed by Kato and McLeod [2]. In the literature [3,4,5,6], the standard model Equations (2) and (3)) was solved via several analytical methods. Patade and Bhalekar [3] applied the Daftardar–Gejji and Jafari Method to obtain a power series solution and provided a theoretical analysis for convergence. Later, Bakodah and Ebaid [4] obtained the exact solution of the standard model (Equations (2) and (3)), i.e., as α 1 , by means of Laplace transform. In addition, Ebaid et al. [5] implemented the Adomian decomposition method (ADM) to derive an accurate approximate solution. Very recently, Alharbi and Ebaid [6] introduced a closed-form analytical solution in terms of exponential functions for the standard model (Equations (2) and (3)). They demonstrated that their new closed-form solution has many advantages over the published one by Patade and Bhalekar [3].
However, the present model (Equations (1) and (2)) is a generalized form of the standard AE (Equations (2) and (3)). In fractional calculus, Kumar et al. [7] solved the fractional model (Equations (1) and (2)) by means of the homotopy transform analysis method considering the Caputo’s fractional derivative definition.
However, the conformable derivative (CD) is one of the most prominent operators in this context. To solve the generalized model (Equations (1) and (2)) using the CD, several analytical approaches can be implemented such as the Adomian decomposition method (ADM) [8,9,10,11,12,13,14,15,16,17,18,19,20], the homotopy perturbation method (HPM) [21,22,23], the differential transform method (DTM)/Taylor expansion [24,25], and the the homotopy analysis method (HAM) [6]. In addition, many applications of the CD have been recently discussed by several authors [26,27,28,29]. Recent applications and advances of fractional operators in various fields can be also found in Refs. [30,31,32,33]. Therefore, the present paper extends the application of the CD to analyze the fractional model (Equations (1) and (2)). The objective of this paper is to investigate the system (Equations (1) and (2)) in view of the CD using the HPM. It will be declared that the HPM is an effective tool to deal with the current model. Moreover, it is shown that a few terms of the series is sufficient to obtain accurate numerical results, where a very small residual error is obtained.

2. Analysis of the HPM

The HPM assumes the solution as an infinite series using an artificial parameter p. To implement this method to deal with the system (Equations (1) and (2)), we rewrite Equation (1) in the form:
D t α z ( t ) = z ( t ) + p 1 q z t q ,
where p is an embedding parameter, which is used to construct the homotopy series solution as
z ( t ) = n = 0 p n z n ( t ) .
Inserting Equation (5) into Equation (4), we have
D t α z 0 ( t ) + z 0 ( t ) + n = 0 p n + 1 D t α z n + 1 ( t ) + z n + 1 ( t ) 1 q z n t q = 0 ,
which leads to the following systems of initial value problems:
D t α z 0 ( t ) + z 0 ( t ) = 0 , z 0 ( 0 ) = λ ,
D t α z n + 1 ( t ) + z n + 1 ( t ) = 1 q z n t q , z n + 1 ( 0 ) = 0 , n 0 .
The conformable derivative of arbitrary order α , 0 < α 1 , of a function z ( t ) : [ 0 , ) R is defined by [26,27,28,29]
D t α z ( t ) = t 1 α d z d t .
Accordingly, the systems in Equations (7) and (8) become
t 1 α z 0 ( t ) + z 0 ( t ) = 0 , z 0 ( 0 ) = λ ,
t 1 α z n + 1 ( t ) + z n + 1 ( t ) = 1 q z n t q , z n + 1 ( 0 ) = 0 , n 0 .
Equation (11) is defined as nth-order systems. The solution of the zeroth-order system (i.e., the system of order zero) in Equation (10) is obtained as
z 0 ( t ) = λ Exp t α α .
From Equations (11) and (12), the first-order system becomes
t 1 α z 1 ( t ) + z 1 ( t ) = λ q Exp q α t α α , z 1 ( 0 ) = 0 .
To solve the first-order ordinary differential equation in Equation (13), we rewrite such equation in the form:
z 1 ( t ) + t α 1 z 1 ( t ) = λ q t α 1 Exp q α t α α .
The solution of Equation (14) subject to the initial condition z 1 ( 0 ) = 0 is given in terms of the integrating factor Exp t α α as
Exp τ α α z 1 ( τ ) τ = 0 t = λ q 0 t τ α 1 Exp τ α α Exp q α τ α α d τ ,
or
Exp t α α z 1 ( t ) = λ q 1 q α Exp t α α 1 q α 1 .
Therefore, z 1 ( t ) is finally given by
z 1 ( t ) = λ q 1 q α Exp q α t α α Exp t α α .
The second-order system takes the from of Equation (11):
t 1 α z 2 ( t ) + z 2 ( t )    = 1 q z 1 t q ,             = λ q 2 1 q α Exp q 2 α t α α Exp q α t α α , z 2 ( 0 ) = 0 .
Proceeding as above, the solution of the system in Equation (18) can be solved as follows:
Exp τ α α z 2 ( τ ) τ = 0 t    = λ q 2 1 q α 0 t τ α 1 [ Exp q α τ α α Exp τ α α Exp τ α α ] d τ ,             = λ q 2 1 q α 0 t τ α 1 Exp τ α α ( 1 q 2 α ) Exp τ α α ( 1 q α ) d τ ,             = λ q 2 1 q α I 1 I 2 ,
where I 1 and I 2 are defined, respectively, by
I 1 = 0 t τ α 1 Exp τ α α ( 1 q 2 α ) d τ ,
and
I 2 = 0 t τ α 1 Exp τ α α ( 1 q α ) d τ .
The two integrals above can be explicitly calculated as
I 1 = 1 1 q 2 α Exp t α α ( 1 q 2 α ) 1 ,
and
I 2 = 1 1 q α Exp t α α ( 1 q α ) 1 .
From Equation (19), we obtain z 2 ( t ) in the form:
z 2 ( t ) = λ q 2 1 q α Exp t α α I 1 I 2 .
Implementing Equations (22) and (23), the expression Exp t α α I 1 I 2 in Equation (24) reads
Exp t α α I 1 I 2 = 1 1 q 2 α Exp q 2 α t α α + 1 1 q α 1 1 q 2 α ×           Exp t α α 1 1 q α Exp q α t α α ,           = 1 1 q 2 α Exp q 2 α t α α + q α 1 q 2 α Exp t α α           1 1 q α Exp q α t α α ,
where the following identity:
1 1 q α 1 1 q 2 α = q α 1 q 2 α ,
is used to obtain Equation (25). From Equations (24) and (25), we finally get z 2 ( t ) as
z 2 ( t ) = λ q 2 1 q α 1 q 2 α Exp q 2 α t α α + q α Exp t α α 1 + q α Exp q α t α α .
In the same manner, the solutions of higher-order systems can be provided. Now, the HPM gives the n-term approximate analytic solution ρ n ( t ) when p 1 as
ρ n ( t ) = i = 0 n 1 z i ( t ) .
Accordingly, we have the following three-term approximate analytic solution:
ρ 3 ( t )    = λ Exp t α α + λ q 1 q α Exp q α t α α Exp t α α +       Exp q 2 α t α α + q α Exp t α α 1 + q α Exp q α t α α ×       λ q 2 1 q α 1 q 2 α .
At this point, the approximate solution in Equation (29) can be validated as α 1 . In that case, the approximate solution in Equation (29) reduces to
ρ 3 ( t )    = λ Exp t + λ q 1 Exp t / q Exp t + λ q 1 q 2 1 ×       q Exp t / q 2 + Exp t q + 1 Exp t / q ,
which agrees with the third-order approximate solution in [5] for the model in Equations (1)–(3). Hence, the present results are applicable for 0 < α 1 and more general than those in relevant literature [3,4,5,6] when α 1 . The applicability and validity of the obtained results are illustrated in the next section, where the convergence of the approximate solutions in Equation (28) is valid in the whole domain t 0 for all values of q > 1 . Moreover, the present approximate numerical results are validated by calculating the residual | R E n ( t ) | defined by
R E n ( t ) = t 1 α ρ n ( t ) + ρ n ( t ) 1 q ρ n t q , n 1 ,
by using the n-term approximate solution to estimate z ( t ) .

3. Results

This section is devoted to the domain of applicability and validity of the obtained approximate solutions in the previous section for the conformable form of the AE. The convergence of the Homotopy perturbation method (HPM) was discussed in detail for ODEs (ordinary differential equations) by Ayati and Biazar [23]. For FPDEs (fractional partial differential equations), the convergence of the HPM was recently analyzed by Touchent et al. [34] and Sene and Fall [33]. Ayati and Biazar [23] proved that the series in Equation (28) is convergent in the limit if ( 0 ϵ i < 1 ) such that ϵ i = z i + 1 z i , i N , where f ( t ) = max 0 t 1 f ( t ) . For illustration of convergence, we consider α = 0.8 , λ = 1 , and q = 1.4 . Accordingly, we have ϵ 1 = 0.2988 , ϵ 2 = 0.3853 , ϵ 3 = 0.1932 , ϵ 4 = 0.1089 , ϵ 5 = 0.0657 , ϵ 6 = 0.0414 , ϵ 7 = 0.0268 , ϵ 8 = 0.0175 , and ϵ 9 = 0.0120 , which proves the convergence of the present approximate solutions, even if the domain 0 t 1 is extended. If the domain is enlarged to be 0 t τ , τ > 1 , we still have 0 ϵ i < 1 , i N .
The convergence of the sequence ρ n ( t ) in Equation (28) is demonstrated in Figure 1, where the approximate solutions ρ 6 ( t ) , ρ 8 ( t ) and ρ 10 ( t ) are depicted at α = 0.8 , λ = 1 , and q = 1.6 . This figure shows that the approximate solutions using a few terms approach to a certain curve as the number of terms increases, even for large values of the independent variable. This may explain how the convergence of approximate solution is estimated, in addition to the above analysis of convergence. Furthermore, it is shown below that the approximate solution ρ 10 ( t ) leads to highly accurate numerical results, especially, when q 2 .
The variation of ρ 10 ( t ) versus t is displayed in Figure 2 at various values of the arbitrary order α of the CD. To validate the current results, the residual | R E 10 | is plotted versus t in Figure 3 and Figure 4 at two particular values of the delay parameter q, q = 1.6 and q = 2 , respectively, for various values of α . The obtained results in these two figures reveal that | R E 10 | is small when q = 1.6 , while it rapidly decreases as q increases as observed from Figure 4 when q = 2 . In addition, the profile of | R E 10 | against two sub-domains of the parameter q ( > 1 ) is introduced through Figure 5 ( 1.1 q 1.6 ) and Figure 6 ( 5 q 10 ) in the domain t [ 0 , 50 ] at λ = 1 and α = 0.8 . It is noticed from these figures that the maximum values of | R E 10 | in these sub-domains are, respectively, 6 × 10 2 and 2 × 10 16 . In addition, it can be concluded that the residual becomes very small and rapidly decreases (approaches zero) for relatively higher values of q, i.e., when q 2 .
At α = 0.8 and q = 2 , the impact of the initial condition λ ( 0 λ 5 ) on the variation of | R E 10 | is depicted in Figure 7. Figure 8 shows the profile of | R E 10 | verses the arbitrary order α ( 0.5 α 1 ). Finally, at a particular value of the independent variable t = 100 , the behavior of the residual | R E 10 | is presented in Figure 9 versus both λ ( 0 λ 5 ) and q ( 2 q 10 ). The obtained results in Figure 7, Figure 8 and Figure 9 confirm that the present 10-term approximate solution ρ 10 ( t ) is highly accurate. Therefore, the preceding discussion shows the effectiveness of the HPM in accurately and quickly solving the generalized form of AE in view of the CD.

4. Conclusions

The homotopy perturbation method is applied in this paper to solve the fractional form of the Ambartsumian equation in view of the conformable derivative. The obtained approximate solution is expressed in terms of the exponential functions and the arbitrary order of the conformable derivative. The present approximate solution reduces to the corresponding one in the literature when the order of the conformable derivative tends to unity. Moreover, the impacts of the initial condition, the delay parameter, and the arbitrary order of the conformable derivative on the residual error are discussed in detail. It is also declared that the residual achieved by using only ten terms is very small for moderate values of the delay parameter while such residual approaches zero as the the delay parameter increases.The obtained results are applicable for 0 < α 1 and more general than those in relevant literature [3,4,5,6] in which α 1 . The obtained residual errors in Figure 7, Figure 8 and Figure 9 are less than 3 × 10 6 , which confirms the accuracy of the approximate solution. In addition, the obtained residual error in Figure 6 is less than 2 × 10 16 which is highly accurate.

Author Contributions

The authors participated equally in this work. All the authors read and approved the final manuscript.

Funding

The authors would like to acknowledge financial support for this work from the Deanship of Scientific Research (DSR), University of Tabuk, Tabuk, Saudi Arabia, under Grant no. S-015-1439.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ambartsumian, V.A. On the fluctuation of the brightness of the milky way. Doklady Akad Nauk USSR 1994, 44, 223–226. [Google Scholar]
  2. Kato, T.; McLeod, J.B. The functional-differential equation y′(x) = ay(λx) + by(x). Bull. Am. Math. Soc. 1971, 77, 891–935. [Google Scholar]
  3. Patade, J.; Bhalekar, S. On Analytical Solution of Ambartsumian Equation. Natl. Acad. Sci. Lett. 2017, 40, 291–293. [Google Scholar] [CrossRef]
  4. Bakodah, H.O.; Ebaid, A. Exact Solution of Ambartsumian Delay Differential Equation and Comparison with Daftardar-Gejji and Jafari Approximate Method. Mathematics 2018, 6, 331. [Google Scholar] [CrossRef]
  5. Ebaid, A.; Al-Enazi, A.; Albalawi, B.Z.; Aljoufi, M.D. Accurate Approximate Solution of Ambartsumian Delay Differential Equation via Decomposition Method. Math. Comput. Appl. 2019, 24, 7. [Google Scholar] [CrossRef]
  6. Alharbi, F.M.; Ebaid, A. New analytic solution for Ambartsumian equation. J. Math. Syst. Sci. 2018, 8, 182–186. [Google Scholar]
  7. Kumar, D.; Singh, J.; Baleanu, D.; Rathore, S. Analysis of a fractional model of the Ambartsumian equation. Eur. Phys. J. Plus 2018, 133, 259. [Google Scholar] [CrossRef]
  8. Wazwaz, A.M. Adomian decomposition method for a reliable treatment of the Bratu-type equations. Appl. Math. Comput. 2005, 166, 652–663. [Google Scholar] [CrossRef]
  9. Ebaid, A. Approximate analytical solution of a nonlinear boundary value problem and its application in fluid mechanics. Z. Naturforschung A 2011, 6, 423–426. [Google Scholar] [CrossRef]
  10. Ebaid, A. A new analytical and numerical treatment for singular two-point boundary value problems via the Adomian decomposition method. J. Comput. Appl. Math. 2011, 235, 1914–1924. [Google Scholar] [CrossRef] [Green Version]
  11. Aly, E.H.; Ebaid, A.; Rach, R. Advances in the Adomian decomposition method for solving two-point nonlinear boundary value problems with Neumann boundary conditions. Comput. Math. Appl. 2012, 63, 1056–1065. [Google Scholar] [CrossRef] [Green Version]
  12. Chun, C.; Ebaid, A.; Lee, M.; Aly, E. An approach for solving singular two point boundary value problems: Analytical and numerical treatment. ANZIAM J. 2012, 53, 21–43. [Google Scholar] [CrossRef]
  13. Wazwaz, A.M.; Rach, R.; Duan, J.S. Adomian decomposition method for solving the Volterra integral form of the Lane-Emden equations with initial values and boundary conditions. Appl. Math. Comput. 2013, 219, 5004–5019. [Google Scholar] [CrossRef]
  14. Ebaid, A.; Aljoufi, M.D.; Wazwaz, A.M. An advanced study on the solution of nanofluid flow problems via Adomian’s method. Appl. Math. Lett. 2015, 46, 117–122. [Google Scholar] [CrossRef]
  15. Alshaery, A.; Ebaid, A. Accurate analytical periodic solution of the elliptical Kepler equation using the Adomian decomposition method. Acta Astronaut. 2017, 140, 27–33. [Google Scholar] [CrossRef]
  16. Wazwaz, A.M.; Raja, M.A.Z.; Syam, M.I. Reliable treatment for solving boundary value problems of pantograph delay differential equation. Rom. Rep. Phys. 2017, 69, 102. [Google Scholar]
  17. Aljohani, A.F.; Rach, R.; El-Zahar, E.; Wazwaz, A.M.; Ebaid, A. Solution of the hyperbolic kepler equation by Adomian’s asymptotic decomposition method. Rom. Rep. Phys. 2018, 70, 14. [Google Scholar]
  18. Bakodah, H.O.; Ebaid, A.; Wazwaz, A.M. Analytical and numerical treatment of Falkner-Skan equation via a transformation and Adomian’s method. Rom. Rep. Phys. 2018, 70, 17. [Google Scholar]
  19. Gaber, A.A.; Ebaid, A. Analytical study on the slip flow and heat transfer of nanofluids over a stretching sheet using Adomian’s method. Rom. Rep. Phys. 2018, 70, 15. [Google Scholar]
  20. Bakodah, H.O.; Ebaid, A. The Adomian decomposition method for the slip flow and heat transfer of nanofluids over a stretching/shrinking sheet. Rom. Rep. Phys. 2018, 70, 115. [Google Scholar]
  21. Golmankhaneh, A.K.; Golmankhaneh, A.K.; Baleanu, D. Homotopy perturbation method for solving a system of Schrodinger-Korteweg-De Vries equations. Rom. Rep. Phys. 2011, 63, 609–623. [Google Scholar]
  22. Patra, A.; Ray, S.S. Homotopy perturbation sumudu transform method for solving convective radial fins with temperature-dependent thermal conductivity of fractional order energy balance equation. Int. J. Heat Mass Transf. 2014, 76, 162–170. [Google Scholar] [CrossRef]
  23. Ayati, Z.; Biazar, J. On the convergence of Homotopy perturbation method. J. Egypt. Math. Soc. 2014, 23, 424–428. [Google Scholar] [CrossRef]
  24. Ebaid, A. A reliable aftertreatment for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 528–536. [Google Scholar] [CrossRef]
  25. Ebaid, A.; Momani, S.; Chang, S.-H. On the periodic solutions of the nonlinear oscillators. J. Vibroeng. 2014, 16, 1–14. [Google Scholar]
  26. Feng, Q. A new approach for seeking coefficient function solutions of conformable fractional partial differential equations based on the Jacobi elliptic equation. Chin. J. Phys. 2018, 56, 2817. [Google Scholar] [CrossRef]
  27. Ebaid, A.; Masaedeh, B.; El-Zahar, E. A New Fractional Model for the Falling Body Problem. Chin. Phys. Lett. 2017, 34, 020201. [Google Scholar] [CrossRef]
  28. Yaslan, H.C. Numerical solution of the conformable space-time fractional wave equation. Chin. J. Phys. 2018, 56, 2916. [Google Scholar] [CrossRef]
  29. Rezazadeh, H.; Tariq, H.; Eslami, M.; Mirzazadeh, M.; Zhou, Q. New exact solutions of nonlinear conformable time-fractional Phi-4 equation. Chin. J. Phys. 2018, 56, 2805. [Google Scholar] [CrossRef]
  30. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998; Volume 198. [Google Scholar]
  31. Hristov, J. Approximate solutions to fractional sub-diffusion equations. Eur. Phys. J. Spec. Top. 2011, 193, 229–243. [Google Scholar] [CrossRef]
  32. Dos Santos, M. Non-Gaussian distributions to random walk in the context of memory kernels. Fractal Fract. 2018, 2, 20. [Google Scholar] [CrossRef]
  33. Sene, N.; Fall, A.N. Homotopy Perturbation-Laplace Transform Method and Its Application to the Fractional Diffusion Equation and the Fractional Diffusion Reaction Equation. Fractal Fract. 2019, 3, 14. [Google Scholar] [CrossRef]
  34. Ait Touchent, K.; Hammouch, Z.; Mekkaoui, T.; Belgacem, F. Implementation and Convergence Analysis of Homotopy Perturbation Coupled With Sumudu Transform to Construct Solutions of Local-Fractional PDEs. Fractal Fract. 2018, 2, 22. [Google Scholar] [CrossRef]
Figure 1. Logarithmic plot of the approximate solutions.
Figure 1. Logarithmic plot of the approximate solutions.
Mathematics 07 00425 g001
Figure 2. The variation of the approximate solution versus the conformable derivative order α .
Figure 2. The variation of the approximate solution versus the conformable derivative order α .
Mathematics 07 00425 g002
Figure 3. Effect of the conformable derivative order α on the residual at λ = 1 , q = 1.6 .
Figure 3. Effect of the conformable derivative order α on the residual at λ = 1 , q = 1.6 .
Mathematics 07 00425 g003
Figure 4. Effect of the conformable derivative order α on the residual at λ = 1 , q = 2 .
Figure 4. Effect of the conformable derivative order α on the residual at λ = 1 , q = 2 .
Mathematics 07 00425 g004
Figure 5. Residual at α = 0.8 , λ = 1 , 1.1 q 1.6 .
Figure 5. Residual at α = 0.8 , λ = 1 , 1.1 q 1.6 .
Mathematics 07 00425 g005
Figure 6. Residual at α = 0.8 , λ = 1 , 5 q 10 .
Figure 6. Residual at α = 0.8 , λ = 1 , 5 q 10 .
Mathematics 07 00425 g006
Figure 7. Residual at α = 0.8 , q = 2 , 0 λ 5 .
Figure 7. Residual at α = 0.8 , q = 2 , 0 λ 5 .
Mathematics 07 00425 g007
Figure 8. Residual at λ = 1 , q = 2 , 0.5 α 1 .
Figure 8. Residual at λ = 1 , q = 2 , 0.5 α 1 .
Mathematics 07 00425 g008
Figure 9. Residual at α = 0.8 , 2 q 10 , 0 λ 5 , t = 100.
Figure 9. Residual at α = 0.8 , 2 q 10 , 0 λ 5 , t = 100.
Mathematics 07 00425 g009

Share and Cite

MDPI and ACS Style

Khaled, S.M.; El-Zahar, E.R.; Ebaid, A. Solution of Ambartsumian Delay Differential Equation with Conformable Derivative. Mathematics 2019, 7, 425. https://doi.org/10.3390/math7050425

AMA Style

Khaled SM, El-Zahar ER, Ebaid A. Solution of Ambartsumian Delay Differential Equation with Conformable Derivative. Mathematics. 2019; 7(5):425. https://doi.org/10.3390/math7050425

Chicago/Turabian Style

Khaled, Sayed M., Essam R. El-Zahar, and Abdelhalim Ebaid. 2019. "Solution of Ambartsumian Delay Differential Equation with Conformable Derivative" Mathematics 7, no. 5: 425. https://doi.org/10.3390/math7050425

APA Style

Khaled, S. M., El-Zahar, E. R., & Ebaid, A. (2019). Solution of Ambartsumian Delay Differential Equation with Conformable Derivative. Mathematics, 7(5), 425. https://doi.org/10.3390/math7050425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop