Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 RESEARCH Open Access Finite-time stability analysis of fractional singular time-delay systems Denghao Pang* and Wei Jiang * Correspondence: [email protected] School of Mathematical Sciences, Anhui University, Hefei, 230039, China Abstract This paper studies the finite-time stability of fractional singular time-delay systems. First, by the method of the steps, we discuss the existence and uniqueness of the solutions for the equivalent systems to the fractional singular time-delay systems. Furthermore, we give the Mittag-Leffler estimation of the solutions for the equivalent systems and obtain the sufficient conditions of the finite-time stability for the original systems. MSC: 34K20 Keywords: finite-time stability; singular systems; time delay; fractional calculus; Mittag-Leffler estimation; generalized Gronwall inequality 1 Introduction In the past years or so, fractional calculus has attracted many physicists, mathematicians, and engineers, and notable contributions have been made to both the theory and the applications of fractional differential equations (see [–]). Moreover, the different techniques have been applied to investigate the stability of various fractional dynamical systems, such as the principle of contraction mappings [], the Lyapunov direct method [], linear matrix inequalities [], Gronwall inequalities [–] and fixed-point theorems []. At the same time, we notice that large numbers of practical systems, such as economic systems, power systems and so on, are singular differential systems which are also named differential-algebraic systems or descriptor systems. Such systems have some particular properties including regularity and impulse behavior which does not need to be considered in normal systems. In [–], the authors discuss singular systems with or without delay and obtain some important results. However, in the previous literature, there are few results on the stability of fractional singular systems, especially the fractional singular systems with time delay. In this regard, it is necessary and important to study the stability problems for fractional singular dynamical systems. Motivated by this consideration, in this paper, we investigate the stability of fractional singular dynamical systems with state delay via the generalized Gronwall approach. In this paper, we consider the following fractional singular time-delay system: E(c Dα x(t)) = Ax(t) + Bx(t – τ ), x(t) = ϕ(t), t ∈ [–τ , ], t ∈ [, T], (.) ©2014 Pang and Jiang; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 2 of 11 where c Dα denotes the Caputo fractional derivative of order < α ≤ ; the vector function x(t) ∈ Rn is a state vector; A, B, E ∈ Rn×n are constant matrices; E ∈ Rn×n is a singular matrix i.e. rank(E) = q < n; the constant parameter τ > represents the delay argument and ϕ(t) is a given sufficiently often differentiable function on [–τ , ]. The organization of this paper is as follows. In Section , we summarize some notations and give preliminary results which will be used in this paper. In Section , we present our main results. 2 Preliminaries and lemmas For completeness, in this section, we firstly demonstrate and study the definitions and some fundamental results of fractional calculus which can be found in [–]. Definition . (see []) The Euler gamma function is defined as ∞ (z) = z ∈ C, e–t t z– dt, (.) where C denotes the complex plane. Definition . (see []) The fractional integral of order α with the lower limit zero for any function f (t) ∈ C([, +∞), R), t ≥ is defined as t n α f (t – rh) = (t – θ )α– f (θ) dθ , I f (t) = lim h h→ (α) r r= α α α > , (.) nh=t where α r = α(α+)···(α+r–) , r! (·) is the gamma function. Definition . (see []) The Riemann-Liouville derivative of order α with the lower limit zero for any function f (t) ∈ C([, +∞), R), t ≥ is defined as l dn D f (t) = (n – α) dt n t α (t – θ )n–α– f (θ) dθ , n – < α < n. (.) Definition . (see []) The Caputo derivative of order α for any function f (t) ∈ C n ([, +∞), R), t ≥ , is defined as c Dα f (t) = (n – α) t (t – θ )n–α– f (n) (θ) dθ = I n–α f (n) (t), n – < α < n. (.) Remark . (see []) (i) The Laplace transform of the Caputo derivative is n– sα–k– f (k) () (n – < α ≤ n); L c Dα f (t); s = sα F(s) – (.) k= (ii) The Caputo fractional derivative is a linear operator satisfying the relation c Dα λf (t) + μg(t) = λc Dα f (t) + μc Dα g(t), where λ and μ are scalars. (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 3 of 11 Lemma . (see []) Let < α < , then we have I α c Dα x(t) = x(t) – x(). (.) Definition . (see []) The Mittag-Leffler function in two parameters is defined as Eα,β (z) = ∞ k= zk , (αk + β) z ∈ C, (.) where α > , β > , and z ∈ C. Remark . (see []) λk (zα )k z (i) For β = , Eα, (λzα ) = Eα (λzα ) = ∞ k= (αk+) , and E, (z) = e , z ∈ C; (ii) for β = , the matrix extension of the aforementioned Mittag-Leffler function has Ak (t α )k the following representation: Eα (At α ) = ∞ k= (αk+) , z ∈ C and c α D Eα (At α ) = AEα (At α ); (iii) we have the Laplace transform of the Mittag-Leffler function in two parameters (k) ±at α ; s = L t αk+β– Eα,β k!sα–β (sα ∓ a) k+ Re(s) > |a|/α , (.) where Re(s) denotes the real parts of s. Next, we introduce some fundamental definitions and lemmas about singular systems. Definition . (see []) For any given two matrices E, A ∈ Rn×n , the pencil (E, A) is called regular if there exists a constant scalar λ ∈ C such that |λE + A| = , or the polynomial |sE – A| ≡ . Lemma . (see []) The pencil (E, A) is regular if and only if two nonsingular matrices Q, P may be chosen such that QEP = diag(In , N), QAP = diag(A , In ), (.) where n + n = n; A ∈ Rn ×n ; N ∈ Rn ×n is nilpotent; In , In are identity matrices. Remark . (see []) (i) N ∈ Rn ×n is nilpotent (the nilpotent index is denoted by h), and we have N h = and N h– = , (.) where h is also called the index of the matrix pair (E, A); (ii) the system (.) will be termed regular if the pencil (E, A) is regular. In the following, we present the first equivalent form (FE) of system (.) by the coordinate transformation, which is also called the standard decomposition of a singular system. Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 4 of 11 For convenience, we denote c Dα x(t) by x(α) (t), from Lemma . and Remark ., we deduce the following statement. Assume that the system (.) is regular throughout this paper, there exist two nonsingular matrices Q and P such that the system (.) is a restricted system equivalent to ⎧ (α) ⎪ ⎪ x (t) = A x (t) + B x (t – τ ) + B x (t – τ ), ⎪ ⎪ ⎨ x (t) = ϕ (t), –τ ≤ t ≤ ; (α) ⎪ ⎪ ⎪ Nx (t) = x (t) + B x (t – τ ) + B x (t – τ ), ⎪ ⎩ x (t) = ϕ (t), –τ ≤ t ≤ , t ≥ , (a) (.) t ≥ , (b) with the coordinate transformation x¯ = [x /x ] = P– x, ϕ¯ = [ϕ /ϕ ] = P– ϕ; x , ϕ ∈ Rn ; x , ϕ ∈ Rn , (.) and QEP = diag(In , N), QAP = diag(A , In ), B QBP = B B , B (.) where ϕ ∈ Rn ; ϕ ∈ Rn ; n + n = n; A , B ∈ Rn ×n ; N, B ∈ Rn ×n ; B ∈ Rn ×n ; N is nilpotent. The following definitions and lemmas will play important roles in our next analysis. Definition . (see []) If X and Y are normed linear spaces, an operator T : X → Y is linear if T(αx + βx ) = αT(x ) + βT(x ), (.) for all x , x in X and scalars α and β. Remark . (see []) We say the linear operator T is a bounded linear operator from X to Y if there is a finite constant C such that Tx Y ≤ C x X for all x in X. Lemma . (see []) If T : X → Y is a linear operator from a normed linear space X to a normed linear space Y , the following are equivalent: (i) T is bounded; (ii) T is continuous; (iii) T is continuous at . Lemma . (see []; Generalized Gronwall Inequality) Suppose x(t), a(t) are nonnegative and local integrable on ≤ t < T; some T ≤ ∞, and g(t) is a nonnegative, nondecreasing continuous function defined on ≤ t < T; g(t) ≤ C , where C is a constant, α > with t x(t) ≤ a(t) + g(t) (t – s)α– x(s) ds (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 5 of 11 on this interval. Then t ∞ [g(t)(α)]n (t – s)nα– a(s) ds, x(t) ≤ a(t) + g(t) (nα) n= ≤ t < T. (.) Lemma . (see []) Under the hypothesis of Theorem ., let a(t) be a nondecreasing function on [, T). Then x(t) ≤ a(t)Eα g(t)(α)t α , (.) where Eα is the Mittag-Leffler function. 3 Main results In this section, we discuss some problems of the singular fractional time-delay system (.). Let Dα be the Caputo fractional differential operator of order < α ≤ , Dnα f (t) = α α α α – D D · · · D f (t) and T = (ND – I) . It is not difficult to verify the following: n T = – I + NDα + N Dα + · · · + N h– D(h–)α , (.) where I ∈ Rn ×n is an identity matrix. Theorem . The fractional differential operator T is bounded i.e. there exists a positive constant M such that for ∀x (t) we have Tx (t) ≤ Mx (t). (.) Proof Obviously, the fractional differential operator T is linear. According to Lemma ., we are only necessary to show that T is continuous at . Let any sequences xn (t) → , yn (t) → y (t), and yn (t) = Txn (t), all we finally need to do is to show that y (t) = . According to yn (t) = Txn (t) = (NDα – I)– xn (t), we have NDα – I yn (t) = xn (t), (.) and for n → ∞, NDα – I y (t) = ⇒ ⇒ NDα y (t) = y (t) y (t) = NDα y (t) = N Dα y (t) = · · · = N h– D(h–)α y (t) = N h Dhα y (t). (.) Combining Remark . and (.) yields y (t) = . Therefore, the operator T is bounded. To give the solution of systems (.), let us define a new function. (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 6 of 11 Definition . (see []) Let α obey ( ≤ α < ), the function δ α (t) = ( – α) t δ(θ) dθ (t – θ )α (.) is called an α – δ function, where δ(t) is the Dirac delta function. Remark . (see []) The Laplace transformation of the α –δ function is L{δ α (t); s} = sα– . Theorem . If the system (.) is regular, the solution for the system (.) exists uniquely. Proof From Remark ., we know that the pencil (E, A) is regular if the system (.) is regular. By the coordinate transformation, the system (.) is equivalent to the system (.). For t ∈ [, τ ], then t – τ ∈ [–τ , ], the system (.)(a) may be written as x(α) (t) = A x (t) + B ϕ (t – τ ) + B ϕ (t – τ ), t ∈ [, τ ]. (.) Let f (t) = B ϕ (t – τ ) + B ϕ (t – τ ). Obviously, if f (t) is the known function, then (.) may be written as x(α) (t) = A x (t) + f (t). (.) Applying the Laplace transformation on both sides of (.) and using (.) yield sα X (s) – sα– x () = A X (s) + F (s), – – X (s) = sα I – A sα– x () + sα I – A F (s). (.) Applying the Laplace inverse transformation on both sides of (.) and using (.) yield x (t) = Eα, A t α x () + t (t – θ )α– Eα,α A (t – θ )α f (θ) dθ , t ∈ [, τ ]. (.) As for the system (.)(b), it may be rewritten as Nx(α) (t) = x (t) + B ϕ (t – τ ) + B ϕ (t – τ ), t ∈ [, τ ]. (.) Similarly, let f (t) = B ϕ (t – τ ) + B ϕ (t – τ ), and f (t) is the known sufficiently often differentiable function, then (.) may be written as Nx(α) (t) = x (t) + f (t). Taking the Laplace transformation on both sides of (.), we have sα N – I X (s) = sα– Nx () + F (s), (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 7 of 11 – α– X (s) = sα N – I s Nx () + F (s) =– h– i N i sα sα– Nx () + F (s) i= =– h N i– iα– s x () – i= h– N i siα F (s). (.) i= According to Remark ., the inverse Laplace transformation of X (s) yields x (s) = – h N i– δ iα x () – i= h– N i Diα f (t), t ∈ [, τ ]. (.) i= Obviously, by the method of steps, once the solution x¯ (t) of the system (.) on [, τ ] is known, continuing the above process, we can easily obtain the solution x¯ (t) of the system (.) on [τ , τ ], [τ , τ ], . . . . Thus the solution x¯ (t) of the system (.) on [, T] exists uniquely. Furthermore, we give the following theorems as regards the Mittag-Leffler estimation of the solution and finite-time stability for this singular system. Let us denote by C([a, b]) the space of all continuous real functions defined on [a, b] and by C([a, b], Rn ) the Banach space of continuous functions mapping the interval [a, b] into Rn with the topology of uniform convergence. Let C = C([–τ , ], Rn ), [a, b] = [–τ , ], and designate the norm of an element ϕ in C by ϕ = sup ϕ(t). (.) –τ ≤t≤ Let X = C([–τ , T], Rn ) and x(t) = ϕ(t), t ∈ [–τ , ] be equipped with the norm x(t) := sup x(t), ≤t≤T xt := x(t + θ ) := sup x(t + θ ), –τ ≤θ ≤ ∀x ∈ X. (.) Definition . (see []) The system given by (.) satisfying the initial condition x(t) = ϕ(t), for t ∈ [–τ , ] is finite-time stable w.r.t. {t , δ, , J}, δ < , J = [t , t + T] if and only if ϕ < δ (.) implies x(t) < , ∀t ∈ J. (.) Theorem . If x¯ (t) = [x (t)/x (t)] is a solution of the system (.), then there exist positive constants a and b such that (i) a = + M B + M B ; (ii) b > A + B + B ; (iii) ¯x(t) ≤ a ϕ E ¯ α (bt α ), ∀t ∈ J = [, T]. Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 8 of 11 Proof According to Lemma ., the system (.)(a) may be rewritten in the form of the equivalent Volterra integral equation x (t) = x () + (α) t (t – s)α– A x (s) + B x (s – τ ) + B x (s – τ ) ds t (t – s)α– A x (s) ds (α) t (t – s)α– B x (s – τ ) + B x (s – τ ) ds, + (α) = ϕ () + t ≥ . (.) Using the appropriate property of the norm · on (.), it follows that x (t) ≤ ϕ () + + (α) (α) t t |t – s|α– A x (s) ds |t – s|α– B x (s – τ ) + B x (s – τ ) ds t |t – s|α– A x¯ (s) ds ≤ ϕ ¯ + (α) t |t – s|α– B x¯ (s – τ ) + B x¯ (s – τ ) ds, + (α) t ≥ . (.) As for the system (.)(b), we have NDα x (t) = x (t) + B x (t – τ ) + B x (t – τ ), – x (t) = NDα – I B x (t – τ ) + B x (t – τ ) = TB x (t – τ ) + TB x (t – τ ). (.) Applying the appropriate property of the norm · and Theorem ., we have x (t) ≤ M B x (t – τ ) + M B x (t – τ ) ≤ M B + M B x¯ (t – τ ), t ≥ . (.) Combining (.) and (.) yields x¯ (t) ≤ x (t) + x (t) ≤ ϕ ¯ + M B + M B x¯ (t – τ ) t + |t – s|α– A x¯ (s) ds (α) t |t – s|α– B + B x¯ (s – τ ) ds, + (α) t ≥ . (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 9 of 11 For ≤ t ≤ τ , ¯x(t – τ ) ≤ ϕ , ¯ (.) can be written as x¯ (t) ≤ + M B + M B ϕ ¯ t |t – s|α– A x¯ (s) ds + (α) t |t – s|α– B + B x¯ (s – τ ) ds, + (α) ≤ t ≤ τ. (.) From Definition ., we know that I α f (t) is an increasing function of t, if f (t) > . So t t α– |t – s|α– A ¯x(s) ds and (α) [ B + B ] ¯x(s – τ ) ds are both in |t – s| (α) creasing functions with regard to t. Taking into account (.) and (.) yields ¯xt ≤ + M B + M B ϕ ¯ t + |t – s|α– A ¯xs ds (α) t |t – s|α– B + B ¯xs ds + (α) ¯ ≤ + M B + M B ϕ t |t – s|α– A + B + B ¯xs ds, + (α) ≤ t ≤ τ. (.) Let a = + M B + M B and b = A + B + B , we can see the function a(t) in Lemma . to be a(t) = + M B + M B ϕ ¯ = a ϕ , ¯ (.) obviously, it is nondecreasing. An application of the corollary of the generalized Gronwall inequality (.) yields ¯ α b t α , ¯xt ≤ a ϕ E ≤ t ≤ τ. (.) Similarly, the same argument implies the following estimate: ¯xt ≤ a ¯xτ Eα b (t – τ )α , τ ≤ t ≤ τ + τ , τ ≥ . (.) From Definition ., we know that the Mittag-Leffler function Eα (t) is an increasing function with regard to t. Therefore, there exists b > b such that Eα (bτ α ) > Eα (b τ α ) and Eα (b(t–τ )α )Eα (b τ α ) < a . Eα (bt α ) Equations (.) and (.) suggest the following general expression: ¯xt ≤ a ϕ E ¯ α bt α , ≤ t ≤ nτ ≤ T. (.) To prove (.) by induction we have to show that it holds for n = because of (.) and if it holds for n = k, then it also holds for n = k + . Indeed, for t ∈ [τ , (k + )τ ], so that t – τ ∈ [, kτ ], on the one hand, using (.), we have ¯xt ≤ a ¯xt–τ Eα b τ α . (.) Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 Page 10 of 11 On the other hand, using (.) we obtain ¯xt–τ ≤ a ϕ E ¯ α b(t – τ )α . (.) Taking into account (.) and (.) we conclude that ¯xt ≤ a a ϕ E ¯ α b(t – τ )α Eα b τ α aEα (b(t – τ )α )Eα (b τ α ) = a ϕ E ¯ α bt α Eα (bt α ) ≤ a ϕ E ¯ α bt α . (.) That is, x¯ (t) ≤ ¯xt ≤ a ϕ E ¯ α bt α . (.) The proof is completed. Theorem . The fractional singular time-delay system given by (.) is finite-time stable w.r.t. {, δ, , J}, δ < , if the following condition is satisfied: a P Eα bt α ≤ , δ ∀t ∈ J = [, T]. (.) Proof From the coordinate transformation (.), we have x(t) = P¯x(t) = P x (t)/x (t) , ϕ(t) = Pϕ(t) ¯ = P ϕ (t)/ϕ (t) . (.) From Theorem . we obtain x(t) ≤ P x¯ (t) ≤ a P ϕ Eα bt α . (.) Hence, using Definition . and the basic condition of Theorem ., it follows that x(t) < , ∀t ∈ J = [, T]. The proof is completed. (.) Competing interests The authors declare that they have no competing interests. Authors’ contributions The authors contributed equally to the manuscript. Both of them read and approved the final manuscript. Acknowledgements The authors are sincerely thankful to the reviewers for their valuable suggestions and insightful comments. This research was jointly supported by National Natural Science Foundation of China (no. 11371027 and no. 11471015), Program of Natural Science of Colleges of Anhui Province (KJ2013A032, KJ2011A020) and the Special Research Fund for the Doctoral Program of the Ministry of Education of China (20093401110001). Received: 11 May 2014 Accepted: 26 September 2014 Published: 13 Oct 2014 Pang and Jiang Advances in Difference Equations 2014, 2014:259 http://www.advancesindifferenceequations.com/content/2014/1/259 References 1. Miller, KS, Boss, B: An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York (1993) 2. Podlubny, I: Fractional Differential Equations. Academic Press, San Diego (1999) 3. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006) 4. Kaczorek, T: Selected Problems of Fractional Systems Theory. Springer, Berlin (2011) 5. Weitzner, H, Zaslavsky, GM: Some applications of fractional equations. Commun. Nonlinear Sci. Numer. Simul. 8(3-4), 273-281 (2003) 6. Machado, JT, Kiryakova, V, Mainardi, F: Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 16(3), 1140-1153 (2011) 7. Machado, JAT, Costa, AC, Quelhas, MD: Fractional dynamics in DNA. Commun. Nonlinear Sci. Numer. Simul. 16(8), 2963-2969 (2011) 8. Zhou, Y, Jiao, F, Li, J: Existence and uniqueness for fractional neutral differential equations with infinite delay. Nonlinear Anal. 71(7-8), 3249-3256 (2009) 9. Zhou, X, Jiang, W, Hu, L: Controllability of a fractional linear time-invariant neutral dynamical system. Appl. Math. Lett. 26(4), 418-424 (2013) 10. Deng, W: Smoothness and stability of the solutions for nonlinear fractional differential equations. Nonlinear Anal. 72(3-4), 1768-1777 (2010) 11. Li, Y, Chen, Y, Podlubny, I: Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput. Math. Appl. 59(5), 1810-1821 (2010) 12. Sabatier, J, Moze, M, Farges, C: LMI stability conditions for fractional order systems. Comput. Math. Appl. 59(5), 1594-1609 (2010) 13. Lazarevi´c, MP: Finite time stability analysis of PDa fractional control of robotic time-delay systems. Mech. Res. Commun. 33(2), 269-279 (2006) 14. Zhang, X: Some results of linear fractional order time-delay system. Appl. Math. Comput. 197(1), 407-411 (2008) 15. Lazarevi´c, MP, Spasi´c, AM: Finite-time stability analysis of fractional order time-delay systems: Gronwall’s approach. Math. Comput. Model. 49(3-4), 475-481 (2009) 16. Pang, D, Jiang, W: Finite-time stability of neutral fractional time-delay systems via generalized Gronwall’s inequality. Abstr. Appl. Anal. 2014, Article ID 610547 (2014) 17. Wang, J, Lv, L, Zhou, Y: New concepts and results in stability of fractional differential equations. Commun. Nonlinear Sci. Numer. Simul. 17(6), 2530-2538 (2012) 18. Dai, L: Singular Control Systems. Springer, Berlin (1989) 19. Fridman, E: Stability of linear descriptor systems with delay: a Lyapunov-based approach. J. Math. Anal. Appl. 273(1), 24-44 (2002) 20. Jiang, W: Eigenvalue and stability of singular differential delay systems. J. Math. Anal. Appl. 297(1), 305-316 (2004) 21. Campbell, SL, Linh, VH Stability criteria for differential-algebraic equations with multiple delays and their numerical solutions. Appl. Math. Comput. 208(2), 397-415 (2009) 22. Liu, X, Zhong, S, Ding, X: A Razumikhin approach to exponential admissibility of switched descriptor delayed systems. Appl. Math. Model. 38(5-6), 1647-1659 (2014) 23. MacCluer, BD: Elementary Functional Analysis. Springer, New York (2009) 24. Ye, H, Gao, J, Ding, Y: A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 328(2), 1075-1081 (2007) 25. Jiang, W: The constant variation formulae for singular fractional differential systems with delay. Comput. Math. Appl. 59(3), 1184-1190 (2010) 10.1186/1687-1847-2014-259 Cite this article as: Pang and Jiang: Finite-time stability analysis of fractional singular time-delay systems. Advances in Difference Equations 2014, 2014:259 Page 11 of 11
© Copyright 2024