Partial Variable Stability for a Class of Nonlinear Systems with Time Delay

Ran Huo¹, Xiaoli Wang^2*

¹College of Science, Inner Mongolia Agricultural University, Huhhot, China.
²College of Statistics and Mathematics, Inner Mongolia University of Finance and Economics, Huhhot, China.

DOI: 10.4236/am.2020.115027                

Abstract

This article first gives a new class of integral inequalities. Then, as an application, the nonlinear neutral differential system with multiple delays is considered, and the trivial solution of the nonlinear neutral system with multiple delays is obtained. Uniform asymptotic Lipschitz stability. Obviously, the above system is a generalization of the traditional differential system. The purpose of this paper is to study the dual stability of neutral differential equations with delays, including equal asymptotically Lipschitz stability and uniformly asymptotic Lipschitz stability. The author uses the method of integral inequality to establish a double stability criterion. As a result, the local stability of differential equations is widely used in theory and practice, such as dynamic systems and control systems.

Keywords

Nonlinear Neutral Systems, Double Stability, Lipschitz Asymptotic Stability, Integral Inequality

1. Introduction

In 1892, Lyapunov, a Russian mathematician, mechanician and physicist, proposed the notion of the stability of motion. He gave the general research methods in his doctoral dissertation “The general problem of the stability of motion” [1], in which he established the foundation of the stability theory. When studying nonlinear systems, especially studying dynamic systems or control systems, we cannot study the stability of all variables because of the technology difficulties, the limitation of practical conditions, or it is not necessary to study all variables considering the actual need. As a result, studying the partial stability of differential equations becomes more important. In addition, the partial stability is widely used in science and technology. For instance the absolute stability of famous Lurie adjusting systems can be changed into a problem of partial stability. In a word, it is of practical significance to study the partial stability of differential equations.

Since Bellman created a class of integral inequalities in 1958, integral inequalities have been greatly developed. The main results are:

In 1960, Li Yuesheng gave the following inequality in [2]:

$u\left(t\right)\le u_0+{\displaystyle {\int }_0^{t}g\left(s\right)u\left(s\right)\text{d}s}+{\displaystyle {\int }_0^{t}f\left(s\right)u^{\alpha }\left(s\right)\text{d}s}$

In 2005, Sligeng discussed the following inequality in [3]:

$\begin{array}{c}{u}_1\left(t\right)\le k_1+{\displaystyle {\int }_0^t{h}_1\left(s\right)u_1\left(s\right)\text{d}s}+{\displaystyle {\int }_0^t{h}_2\left(s\right)u_2\left(s\right){\text{e}}^{\mu s}\text{d}s}\\ +{\displaystyle {\int }_0^t{\bar{h}}_1\left(s\right)u_1^{\alpha }\left(s\right){\text{e}}^{-\left(\alpha -1\right)\mu s}\text{d}s}+{\displaystyle {\int }_0^t{\bar{h}}_2\left(s\right)u_2^{\alpha }\left(s\right){\text{e}}^{\mu s}\text{d}s}\end{array}$

$\begin{array}{c}{u}_2\left(t\right)\le k_2+{\displaystyle {\int }_0^t{h}_3\left(s\right)u_1\left(s\right){\text{e}}^{-\mu s}}\text{d}s+{\displaystyle {\int }_0^t{h}_4\left(s\right)u_2\left(s\right)\text{d}s}\\ +{\displaystyle {\int }_0^t{\bar{h}}_3\left(s\right)u_1^{\alpha }\left(s\right){\text{e}}^{-\alpha \mu s}}\text{d}s+{\displaystyle {\int }_0^t{\bar{h}}_4\left(s\right)u_2^{\alpha }\left(s\right)\text{d}s}\end{array}$

In 2009, the author discussed a new class of inequalities in [4].

Vorotnikov, V.I. [5] [6] considered the following system

$\{\begin{array}{l}\frac{\text{d}y}{\text{d}t}=A\left(t\right)y+B\left(t\right)z+Y\left(t,y,z\right)\\ \frac{\text{d}z}{\text{d}t}=C\left(t\right)y+D\left(t\right)z+Z\left(t,y,z\right)\end{array}$

and studied the double stability as $\Vert y\Vert +\Vert z\Vert \to 0$ and $\frac{\Vert Y\left(t,y,z\right)\Vert +\Vert Z\left(t,y,z\right)\Vert }{\Vert y\Vert +\Vert z\Vert }\to 0$.

In 2002, Wang Feng used the differential inequality of delay in article [7] to study the following delay system:

$\{\begin{array}{l}\frac{\text{d}y}{\text{d}t}=A\left(t\right)y+B\left(t\right)z+Y\left(t,y,z,y\left(t-\tau \right),z\left(t-\tau \right)\right)\\ \frac{\text{d}z}{\text{d}t}=C\left(t\right)y+D\left(t\right)z+Z\left(t,y,z,y\left(t-\tau \right),z\left(t-\tau \right)\right)\end{array}$

In 2006, Siligeng used the integral inequality extended in [3] in [8] to discuss the double stability of the following system to some variables:

$\{\begin{array}{l}\frac{\text{d}y}{\text{d}t}=A\left(t\right)y+f_1\left(t,y,z\right)\\ \frac{\text{d}z}{\text{d}t}=B\left(t\right)z+f_2\left(t,y,z\right)\end{array}$

In this paper the author consider a new class of the nonlinearly perturbed differential systems with time-delay

$\{\begin{array}{l}\frac{\text{d}y}{\text{d}t}=B\left(t\right)y+C\left(t\right)z+Y\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_1}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\hfill \\ \frac{\text{d}z}{\text{d}t}=D\left(t\right)y+E\left(t\right)z+Z\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_2}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\hfill \end{array}$

It is obvious that the above system is a generalization of the systems in [5] [6] [7] [8].

The aim of this paper is to investigate the double stability of neutural differential equations, including Uniform stability and Uniform Lipschitz stability. The author uses the method of differential inequalities with time-delay and integral inequalities to establish double stability criteria.

2. Preliminaries

Consider the following system:

$\frac{\text{d}x}{\text{d}t}=f\left(t,x\left(t\right),x\left(t-\tau \right),\stackrel{\dot{}}{x}\left(t-\tau \right)\right)$ (1)

where $x\in R^n$, $y=col\left(x_1,x_2,\cdots ,x_m\right)$, $z=col\left(x_{m+1},x_{m+2},\cdots ,x_n\right)$, $x=col\left(y,z\right)$, $f\left(t,0,0\right)\equiv 0$, $\tau$ is a nonnegative constant. Let $\varphi \left(t\right)$ be a continuous function, for $\forall t\in E_{t_0}=\left[t_0-\tau ,t_0\right]$.

Definition 1 [9] [10] [11] The trivial solution of system (1) has uniform stability and exponential asymptotic stability with respect to y if, for $\forall \epsilon >0$, $\forall t_0\in I$, $\exists \delta \left(\epsilon \right)>0$, and $\lambda >0$, when $\Vert \varphi \Vert <\delta$ (for $\forall t\in E_{t_0}$ ), such that

$\Vert y\left(t;t_0,\varphi \right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t;t_0,\varphi \right)\Vert <\epsilon \exp \left(-\lambda \left(t-t_0\right)\right),\forall t\ge t_0.$

Definition 2 [9] [10] [11] The trivial solution of system (1) has Lipschitz stability with respect to y if, there exists constants $M\left(t_0\right)>0$ and $\delta \left(t_0\right)>0$, when $\Vert \phi \Vert +\Vert \stackrel{\dot{}}{\phi }\Vert <\delta$ (for $\forall t\in E_{t_0}$ ), such that

$\Vert y\left(t;t_0,\varphi \right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t;t_0,\varphi \right)\Vert \le M\left(t_0\right)\left(\Vert \varphi \Vert +\Vert \stackrel{\dot{}}{\varphi }\Vert \right),\forall t\ge t_0\ge 0.$

Definition 3 [9] [10] [11] The trivial solution of system (1) has equi-exponential Lipschitz asymptotic stability with respect to y if, there exists $\lambda >0$, $K\left(t_0\right)>0$ and $\delta \left(t_0\right)>0$, when $\Vert \phi \Vert +\Vert \stackrel{\dot{}}{\phi }\Vert <\delta$ (for $\forall t\in E_{t_0}$ ), such that

$\Vert y\left(t;t_0,\varphi \right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t;t_0,\varphi \right)\Vert \le K\left(t_0\right)\left(\Vert \varphi \Vert +\Vert \stackrel{\dot{}}{\varphi }\Vert \right)\exp \left(-\lambda \left(t-t_0\right)\right),\forall t\ge t_0\ge 0.$

Definition 4 [9] [10] [11] The trivial solution of system (1) has uniform exponential Lipschitz asymptotic stability with respect to y if, K and $\delta >0$ in definition 3 are ndependent of $t_0$.

Lemma 1 [4] The following conditions are established on $t\ge t_0$ :

i) $k_1,k_2,\mu$ are non-negative constants;

ii) 1)

$\begin{array}{c}{u}_1\left(t\right)\le k_1+{\displaystyle {\int }_{t_0}^{t}a\left(s\right)u_1\left(s\right)\text{d}s}+{\displaystyle {\int }_{t_0}^{t}b\left(s\right)u_2\left(s\right){\text{e}}^{\mu \left(s-t_0\right)}}\text{d}s\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right)u_1^{{\alpha }_i+1}\left(s\right){\text{e}}^{-{\alpha }_i\mu \left(s-t_0\right)}}\text{d}s}+{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right)u_2^{{\alpha }_i+1}\left(s\right){\text{e}}^{\mu \left(s-t_0\right)}}\text{d}s}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{e}_j\left(s\right){\displaystyle {\int }_{t_0}^s{f}_j\left(\tau \right)u_1\left(\tau \right)\text{d}\tau \text{d}s}}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{g}_j\left(s\right)\left[{\displaystyle {\int }_{t_0}^s{h}_j\left(\tau \right)u_2\left(\tau \right){\text{e}}^{\mu \left(\tau -t_0\right)}\text{d}\tau }\right]\text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{o}_k\left(s\right){\displaystyle {\int }_{t_0}^s{p}_k\left(\tau \right)u_1^{{\beta }_k+1}\left(\tau \right)}{\text{e}}^{-{\beta }_k\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}s}\end{array}$

$\begin{array}{l}+{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{q}_k\left(s\right){\displaystyle {\int }_{t_0}^s{r}_k\left(\tau \right)u_2^{{\beta }_k+1}\left(\tau \right)}{\text{e}}^{\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}s}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right)u_1^{{\alpha }_i}\left(s\right){\text{e}}^{-{\alpha }_i\mu \left(s-t_0\right)}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{v}_l\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{w}_l\left(\tau \right)u_1\left(\tau \right)\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right)u_2^{{\alpha }_i}\left(s\right)}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{x}_l\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{y}_l\left(\tau \right)u_2\left(\tau \right){\text{e}}^{\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right)u_1^{{\alpha }_i}\left(s\right){\text{e}}^{-{\alpha }_i\mu \left(s-t_0\right)}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{A}_m\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{B}_m\left(\tau \right)u_1^{{\gamma }_m+1}\left(\tau \right){\text{e}}^{-{\gamma }_m\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right)u_2^{{\alpha }_i}\left(s\right)}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{D}_m\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{E}_m\left(\tau \right)u_2^{{\gamma }_m+1}\left(\tau \right){\text{e}}^{\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}\theta \text{d}s}}\end{array}$

2)

$\begin{array}{c}{u}_2\left(t\right)\le k_2+{\displaystyle {\int }_{t_0}^t\bar{a}\left(s\right)u_1\left(s\right){\text{e}}^{-\mu \left(s-t_0\right)}\text{d}s}+{\displaystyle {\int }_{t_0}^t\bar{b}\left(s\right)u_2\left(s\right)\text{d}s}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(s\right)u_1^{{\alpha }_i+1}\left(s\right){\text{e}}^{-\left({\alpha }_i+1\right)\mu \left(s-t_0\right)}}\text{d}s}+{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(s\right)u_2^{{\alpha }_i+1}\left(s\right)\text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{\bar{e}}_j\left(s\right){\displaystyle {\int }_{t_0}^s{\bar{f}}_j\left(\tau \right)u_1\left(\tau \right){\text{e}}^{-\mu \left(\tau -t_0\right)}\text{d}\tau \text{d}s}}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{\bar{g}}_j\left(s\right)\left[{\displaystyle {\int }_{t_0}^s{\bar{h}}_j\left(\tau \right)u_2\left(\tau \right)\text{d}\tau }\right]\text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{\bar{o}}_k\left(s\right){\displaystyle {\int }_{t_0}^s{\bar{p}}_k\left(\tau \right)u_1^{{\beta }_k+1}\left(\tau \right)}{\text{e}}^{-\left({\beta }_k+1\right)\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}s}\end{array}$

$\begin{array}{l}+{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{\bar{q}}_k\left(s\right){\displaystyle {\int }_{t_0}^s{\bar{r}}_k\left(\tau \right)u_2^{{\beta }_k+1}\left(\tau \right)}\text{d}\tau }\text{d}s}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(s\right)u_1^{{\alpha }_i}\left(s\right){\text{e}}^{-{\alpha }_i\mu \left(s-t_0\right)}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{\bar{v}}_l\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{\bar{w}}_l\left(\tau \right)u_1\left(\tau \right){\text{e}}^{-\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(s\right)u_2^{{\alpha }_i}\left(s\right)}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{\bar{x}}_l\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{\bar{y}}_l\left(\tau \right)u_2\left(\tau \right)\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(s\right)u_1^{{\alpha }_i}\left(s\right){\text{e}}^{-{\alpha }_i\mu \left(s-t_0\right)}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{\bar{A}}_m\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{\bar{B}}_m\left(\tau \right)u_1^{{\gamma }_m+1}\left(\tau \right){\text{e}}^{-\left({\gamma }_m+1\right)\mu \left(\tau -t_0\right)}\text{d}\tau }\text{d}\theta \text{d}s}}\\ +{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(s\right)u_2^{{\alpha }_i}\left(s\right)}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{\bar{D}}_m\left(\theta \right)}{\displaystyle {\int }_{t_0}^{\theta }{\bar{E}}_m\left(\tau \right)u_2^{{\gamma }_m+1}\left(\tau \right)\text{d}\tau }\text{d}\theta \text{d}s}}\end{array}$

where: $u_1\left(t\right),u_2\left(t\right),a\left(t\right),\bar{a}\left(t\right),b\left(t\right),\bar{b}\left(t\right),c_i\left(t\right),{\bar{c}}_i\left(t\right),d_i\left(t\right),{\bar{d}}_i\left(t\right)\left(i=1,2,\cdots ,I\right)$, $e_j\left(t\right),{\bar{e}}_j\left(t\right),f_j\left(t\right),{\bar{f}}_j\left(t\right),g_j\left(t\right),{\bar{g}}_j\left(t\right),h_j\left(t\right),{\bar{h}}_j\left(t\right)\left(j=1,2,\cdots ,J\right)$, $o_k\left(t\right),{\bar{o}}_k\left(t\right),p_k\left(t\right),{\bar{p}}_k\left(t\right),q_k\left(t\right),{\bar{q}}_k\left(t\right),r_k\left(t\right),{\bar{r}}_k\left(t\right)\left(k=1,2,\cdots ,K\right)$, $v_l\left(t\right),{\bar{v}}_l\left(t\right),w_l\left(t\right),{\bar{w}}_l\left(t\right),x_l\left(t\right),{\bar{x}}_l\left(t\right),y_l\left(t\right),{\bar{y}}_l\left(t\right)\left(l=1,2,\cdots ,L\right)$, $A_m\left(t\right),{\bar{A}}_m\left(t\right),B_m\left(t\right),{\bar{B}}_m\left(t\right),D_m\left(t\right),{\bar{D}}_m\left(t\right),E_m\left(t\right),{\bar{E}}_m\left(t\right)\left(m=1,2,\cdots ,M\right)$ are non-negative continuous function on $R_{+}$, and: ${\alpha }_i\left(i=1,2,\cdots ,I\right)$, ${\beta }_k\left(k=1,2,\cdots ,K\right)$, ${\gamma }_m\left(m=1,2,\cdots ,M\right)$ are all constants greater than 1, and $1\le {\alpha }_1\le {\alpha }_2\le \cdots \le {\alpha }_I$, $1\le {\beta }_1\le {\beta }_2\le \cdots \le {\beta }_K$, $1\le {\gamma }_1\le {\gamma }_2\le \cdots \le {\gamma }_M$ set: $\bar{\alpha }=\max \left({\alpha }_I,{\beta }_K,{\gamma }_M\right)$, $\underset{\_}{\alpha }=\min \left({\alpha }_1,{\beta }_1,{\gamma }_1\right)$

iii) let: $k=k_1+k_2$

$F\left(t\right)=\max \left\{a\left(t\right)+\bar{a}\left(t\right),b\left(t\right)+\bar{b}\left(t\right)\right\}$, $G_i\left(t\right)=\max \left\{c_i\left(t\right)+{\bar{c}}_i\left(t\right),d_i\left(t\right)+{\bar{d}}_i\left(t\right)\right\}$,

${\overline{\stackrel{¯}{e}}}_j\left(t\right)=\max \left\{e_j\left(t\right),{\bar{e}}_j\left(t\right)\right\}$, ${\overline{\stackrel{¯}{g}}}_j\left(t\right)=\max \left\{g_j\left(t\right),{\bar{g}}_j\left(t\right)\right\}$,

$H_j\left(t\right)=\max \left\{{\overline{\stackrel{¯}{e}}}_j\left(t\right),{\overline{\stackrel{¯}{g}}}_j\left(t\right)\right\}$, $N_j\left(t\right)=\max \left\{f_j\left(t\right)+{\bar{f}}_j\left(t\right),h_j\left(t\right)+{\bar{h}}_j\left(t\right)\right\}$,

${\overline{\stackrel{¯}{o}}}_k\left(t\right)=\max \left\{o_k\left(t\right),{\bar{o}}_k\left(t\right)\right\}$, ${\overline{\stackrel{¯}{q}}}_k\left(t\right)=\max \left\{q_k\left(t\right),{\bar{q}}_k\left(t\right)\right\}$,

$Q_k\left(t\right)=\max \left\{{\overline{\stackrel{¯}{o}}}_k\left(t\right),{\overline{\stackrel{¯}{q}}}_k\left(t\right)\right\}$, $R_k\left(t\right)=\max \left\{p_k\left(t\right)+{\bar{p}}_k\left(t\right),r_k\left(t\right)+{\bar{r}}_k\left(t\right)\right\}$,

${\overline{\stackrel{¯}{v}}}_l\left(t\right)=\max \left\{v_l\left(t\right),{\bar{v}}_l\left(t\right)\right\}$, ${\overline{\stackrel{¯}{x}}}_l\left(t\right)=\max \left\{x_l\left(t\right),{\bar{x}}_l\left(t\right)\right\}$,

$T_l\left(t\right)=\max \left\{{\overline{\stackrel{¯}{v}}}_l\left(t\right),{\overline{\stackrel{¯}{x}}}_l\left(t\right)\right\}$, $W_l\left(t\right)=\max \left\{w_l\left(t\right)+{\bar{w}}_l\left(t\right),y_l\left(t\right)+{\bar{y}}_l\left(t\right)\right\}$,

${\overline{\stackrel{¯}{A}}}_m\left(t\right)=\max \left\{A_m\left(t\right),{\bar{A}}_m\left(t\right)\right\}$, ${\overline{\stackrel{¯}{D}}}_m\left(t\right)=\max \left\{D_m\left(t\right),{\bar{D}}_m\left(t\right)\right\}$,

$Y_m\left(t\right)=\max \left\{{\overline{\stackrel{¯}{A}}}_m\left(t\right),{\overline{\stackrel{¯}{D}}}_m\left(t\right)\right\}$, $Z_m\left(t\right)=\max \left\{B_m\left(t\right)+{\bar{B}}_m\left(t\right),E_m\left(t\right)+{\bar{E}}_m\left(t\right)\right\}$,

iv) Let:

$\begin{array}{c}\Delta \left(t\right)=F\left(t\right)+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{H}_j\left(t\right)}{\displaystyle {\int }_{t_0}^t{N}_j\left(s\right)\text{d}s}+2{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right)}{\displaystyle {\int }_{t_0}^t{R}_k\left(s\right)\text{d}s}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }\text{d}s}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}\end{array}$

$\begin{array}{c}\Phi \left(t\right)={\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}+{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right)}{\displaystyle {\int }_{t_0}^t{R}_k\left(s\right)\text{d}s}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }\text{d}s}\\ +2{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}\end{array}$

$\Gamma \left(t\right)={\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}$

${\Lambda }_1\left(t\right)=1-\bar{\alpha }{c}^{\bar{\alpha }}{\displaystyle {\int }_{t_0}^t\left[\Phi \left(\tau \right)+\Gamma \left(\tau \right)\right]\exp \left(\bar{\alpha }{\displaystyle {\int }_{t_0}^{\tau }\Delta \left(\sigma \right)\text{d}\sigma }\right)\text{d}\tau }$

$\Pi \left(t\right)=F\left(t\right)+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{H}_j\left(t\right){\displaystyle {\int }_{t_0}^t{N}_j\left(t\right)\text{d}t}}$

$\begin{array}{c}\Theta \left(t\right)={\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right){Щ}^{{\alpha }_i-\underset{\_}{\alpha }}}+{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right){\displaystyle {\int }_{t_0}^t{R}_k\left(s\right){Щ}^{{\beta }_k-\underset{\_}{\alpha }}\text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right){Щ}^{{\alpha }_i-\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l\left(s\right){\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }}}}\end{array}$

$\Sigma \left(t\right)={\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(t\right){\displaystyle {\int }_{t_0}^t{Z}_m\left(s\right){Щ}^{{\gamma }_m-\underset{\_}{\alpha }}\text{d}s}}$

${\Lambda }_2\left(t\right)=1-\underset{\_}{\alpha }{c}^{\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t\left[\Theta \left(\tau \right)+\Sigma \left(\tau \right)\right]\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_{t_0}^{\tau }\Pi \left(\sigma \right)\text{d}\sigma }\right)\text{d}\tau }$

And ${\displaystyle {\int }_{t_0}^{+\infty }\Delta \left(s\right)\text{d}s}<+\infty$,

${\Lambda }_1\left(t\right)>0$, $\left[1-\left(\bar{\alpha }-1\right)c^{\bar{\alpha }}{\displaystyle {\int }_{t_0}^t\Phi \left(s\right){\Lambda }_1^{-\frac{1}{\bar{\alpha }}}\left(s\right)\exp \left(\bar{\alpha }{\displaystyle {\int }_{t_0}^s\Delta \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]>0$

${\Lambda }_2\left(t\right)>0$, $\left[1-\left(\underset{\_}{\alpha }-1\right)c^{\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t\Theta \left(s\right){\Lambda }_2^{-\frac{1}{\underset{\_}{\alpha }}}\left(s\right)\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_{t_0}^s\Sigma \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]>0$

v) Assume:

$\Omega \left(t\right)\le k\exp \left({\displaystyle {\int }_{t_0}^t\Pi \left(s\right)\text{d}s}\right)\cdot {\left[1-\left(\underset{\_}{\alpha }-1\right)k^{\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t\Theta \left(s\right){\Lambda }_2^{-\frac{1}{\underset{\_}{\alpha }}}\left(s\right)\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_{t_0}^s\Pi \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]}^{-\frac{1}{\underset{\_}{\alpha }-1}}$ then: $u_1\left(t\right)\le \Omega \left(t\right){\text{e}}^{\mu \left(t-t_0\right)}$, $u_2\left(t\right)\le \Omega \left(t\right)$.

3. Main Results

Consider the following system

$\{\begin{array}{l}\begin{array}{c}\frac{\text{d}y}{\text{d}t}=B\left(t\right)y+C\left(t\right)z\\ +Y\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_1}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\end{array}\hfill \\ \begin{array}{c}\frac{\text{d}z}{\text{d}t}=D\left(t\right)y+E\left(t\right)z\\ +Z\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_2}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\end{array}\hfill \end{array}$ (2)

where $\tau \ge 0$ is a constant, initial condition is:

$x\left(t\right)=\varphi \left(t\right),\stackrel{\dot{}}{x}\left(t\right)=\stackrel{\dot{}}{\varphi }\left(t\right),t_0-\tau \le t\le t_0,$

$B\left(t\right)$ is an $m\times m$ matrix, $Y\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_1}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)$ is an $m\times 1$ matrix, $Z\left(s,y\left(s\right),z\left(s\right),{\displaystyle {\int }_0^t{h}_2}\left(s,y\left(s\right),z\left(s\right),y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)$ is an $\left(n-m\right)\times 1$ matrix, they are all continuous for $t\in I$ and satisfy the condition of existence and uniqueness theorem.

Set $Y\left(t,s\right)$ and $Z\left(t,s\right)$ satisfied:

$\{\begin{array}{l}\frac{\partial Y\left(t,s\right)}{\partial t}=B\left(t\right)Y\left(t,s\right)\\ Y\left(s,s\right)=I\end{array}$,

$\{\begin{array}{l}\frac{\partial Z\left(t,s\right)}{\partial t}=E\left(t\right)Z\left(t,s\right)\\ Z\left(s,s\right)=I\end{array}$

Theorem If (2) satisfies the following conditions:

i) $\Vert Y\left(t,s\right)\Vert \le m_1{\text{e}}^{-\lambda \left(t-s\right)}$, $\Vert Z\left(t,s\right)\Vert \le m_2$ ;

ii) $\begin{array}{l}\Vert Y\left(t,y,z,y\left(t-\tau \right),z\left(t-\tau \right),{\displaystyle {\int }_0^t{h}_1}\left(s,y,z,y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\Vert \\ \le a\left(t\right)\left(\Vert y\left(t-\delta \left(t\right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert \right)\\ +b\left(t\right)\left(\Vert z\left(t-\delta \left(t\right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert \right){\text{e}}^{-\epsilon \left(t-t_0\right)}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}\right){\text{e}}^{{\alpha }_i\epsilon \left(t-t_0\right)}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}\right){\text{e}}^{-\epsilon \left(t-t_0\right)}}\end{array}$

$\begin{array}{l}+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{e}_j\left(t\right){\displaystyle {\int }_{t_0}^t{f}_j\left(s\right)\left(\Vert y\left(s-\delta \left(s\right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(s-\delta \left(s\right)\right)\Vert \right)\text{d}s}}\\ +{\displaystyle \underset{j=1}{\overset{J}{\sum }}{g}_j\left(t\right){\displaystyle {\int }_{t_0}^t{h}_j\left(s\right)\left(\Vert z\left(s-\delta \left(s\right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(s-\delta \left(s\right)\right)\Vert \right){\text{e}}^{-\epsilon \left(s-t_0\right)}\text{d}s}}\\ +{\displaystyle \underset{k=1}{\overset{K}{\sum }}{o}_k\left(t\right){\displaystyle {\int }_{t_0}^t{p}_k\left(s\right)\left({\Vert y\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}+{\Vert \stackrel{\dot{}}{y}\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}\right){\text{e}}^{{\beta }_k\epsilon \left(s-t_0\right)}\text{d}s}}\\ +{\displaystyle \underset{k=1}{\overset{K}{\sum }}{q}_k\left(t\right){\displaystyle {\int }_{t_0}^t{r}_k\left(s\right)\left({\Vert z\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}+{\Vert \stackrel{\dot{}}{z}\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}\right){\text{e}}^{-\epsilon \left(s-t_1\right)}\text{d}s}}\end{array}$

$\begin{array}{l}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right){\text{e}}^{{\alpha }_i\epsilon \left(t-t_0\right)}}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{v}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{w}_l\left(\theta \right)\left(\Vert y\left(\theta -\delta \left(\theta \right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(\theta -\delta \left(\theta \right)\right)\Vert \right)\text{d}\theta \text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right)}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{x}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{y}_l\left(\theta \right)\left(\Vert z\left(\theta -\delta \left(\theta \right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(\theta -\delta \left(\theta \right)\right)\Vert \right){\text{e}}^{-\epsilon \left(\theta -t_0\right)}}\text{d}\theta \text{d}s}\end{array}$

$\begin{array}{l}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right){\text{e}}^{{\alpha }_i\epsilon \left(t-t_0\right)}}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{A}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{B}_m\left(\theta \right)\left({\Vert y\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}+{\Vert \stackrel{\dot{}}{y}\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}\right){\text{e}}^{{\gamma }_m\epsilon \left(\theta -t_0\right)}\text{d}\theta \text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right)}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{D}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{E}_m\left(\theta \right)\left({\Vert z\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}+{\Vert \stackrel{\dot{}}{z}\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}\right){\text{e}}^{-\epsilon \left(\theta -t_0\right)}}\text{d}\theta \text{d}s}\end{array}$

iv) $\begin{array}{l}\Vert Z\left(t,y,z,y\left(t-\tau \right),z\left(t-\tau \right),{\displaystyle {\int }_0^t{h}_2}\left(s,y,z,y\left(s-\tau \right),\stackrel{\dot{}}{z}\left(s-\tau \right)\right)\text{d}s\right)\Vert \\ \le \bar{a}\left(t\right)\left(\Vert y\left(t-\delta \left(t\right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert \right){\text{e}}^{\epsilon \left(t-t_0\right)}\\ +\bar{b}\left(t\right)\left(\Vert z\left(t-\delta \left(t\right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert \right)\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}\right){\text{e}}^{\left({\alpha }_i+1\right)\epsilon \left(t-t_0\right)}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i+1}\right)}\end{array}$

$\begin{array}{l}+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{\bar{e}}_j\left(t\right){\displaystyle {\int }_{t_0}^t{\bar{f}}_j\left(s\right)\left(\Vert y\left(s-\delta \left(s\right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(s-\delta \left(s\right)\right)\Vert \right){\text{e}}^{\epsilon \left(s-t_0\right)}\text{d}s}}\\ +{\displaystyle \underset{j=1}{\overset{J}{\sum }}{\bar{g}}_j\left(t\right){\displaystyle {\int }_{t_0}^t{\bar{h}}_j\left(s\right)\left(\Vert z\left(s-\delta \left(s\right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(s-\delta \left(s\right)\right)\Vert \right)\text{d}s}}\\ +{\displaystyle \underset{k=1}{\overset{K}{\sum }}{\bar{o}}_k\left(t\right){\displaystyle {\int }_{t_0}^t{\bar{p}}_k\left(s\right)\left({\Vert y\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}+{\Vert \stackrel{\dot{}}{y}\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}\right){\text{e}}^{\left({\beta }_k+1\right)\epsilon \left(s-t_0\right)}\text{d}s}}\\ +{\displaystyle \underset{k=1}{\overset{K}{\sum }}{\bar{q}}_k\left(t\right){\displaystyle {\int }_{t_0}^t{\bar{r}}_k\left(s\right)\left({\Vert z\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}+{\Vert \stackrel{\dot{}}{z}\left(s-\delta \left(s\right)\right)\Vert }^{{\beta }_k+1}\right)\text{d}s}}\end{array}$

$\begin{array}{l}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right){\text{e}}^{{\alpha }_i\epsilon \left(t-t_0\right)}}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{v}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{w}_l\left(\theta \right)\left(\Vert y\left(\theta -\delta \left(\theta \right)\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(\theta -\delta \left(\theta \right)\right)\Vert \right){\text{e}}^{\epsilon \left(\theta -t_0\right)}\text{d}\theta \text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right)}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{\bar{x}}_l\left(s\right){\displaystyle {\int }_{t_0}^s{\bar{y}}_l\left(\theta \right)}}}\left(\Vert z\left(\theta -\delta \left(\theta \right)\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(\theta -\delta \left(\theta \right)\right)\Vert \right)\text{d}\theta \text{d}s\end{array}$

$\begin{array}{l}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{c}}_i\left(t\right)\left({\Vert y\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{y}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right){\text{e}}^{{\alpha }_i\epsilon \left(t-t_0\right)}}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{\bar{A}}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{\bar{B}}_m\left(\theta \right)\left({\Vert y\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}+{\Vert \stackrel{\dot{}}{y}\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}\right){\text{e}}^{\left({\gamma }_m+1\right)\epsilon \left(\theta -t_0\right)}\text{d}\theta \text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{\bar{d}}_i\left(t\right)\left({\Vert z\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}+{\Vert \stackrel{\dot{}}{z}\left(t-\delta \left(t\right)\right)\Vert }^{{\alpha }_i}\right)}\\ \cdot {\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{\bar{D}}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{\bar{E}}_m\left(\theta \right)\left({\Vert z\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}+{\Vert \stackrel{\dot{}}{z}\left(\theta -\delta \left(\theta \right)\right)\Vert }^{{\gamma }_m+1}\right)\text{d}\theta \text{d}s}}\end{array}$

where： $a\left(t\right),\bar{a}\left(t\right),b\left(t\right),\bar{b}\left(t\right),c_i\left(t\right),{\bar{c}}_i\left(t\right),d_i\left(t\right),{\bar{d}}_i\left(t\right)\left(i=1,2,\cdots ,I\right)$,

$e_j\left(t\right),{\bar{e}}_j\left(t\right),f_j\left(t\right),{\bar{f}}_j\left(t\right),g_j\left(t\right),{\bar{g}}_j\left(t\right),h_j\left(t\right),{\bar{h}}_j\left(t\right)\left(j=1,2,\cdots ,J\right)$,

$o_k\left(t\right),{\bar{o}}_k\left(t\right),p_k\left(t\right),{\bar{p}}_k\left(t\right),q_k\left(t\right),{\bar{q}}_k\left(t\right),r_k\left(t\right),{\bar{r}}_k\left(t\right)\left(k=1,2,\cdots ,K\right)$,

$v_l\left(t\right),{\bar{v}}_l\left(t\right),w_l\left(t\right),{\bar{w}}_l\left(t\right),x_l\left(t\right),{\bar{x}}_l\left(t\right),y_l\left(t\right),{\bar{y}}_l\left(t\right)\left(l=1,2,\cdots ,L\right)$,

$A_m\left(t\right),{\bar{A}}_m\left(t\right),B_m\left(t\right),{\bar{B}}_m\left(t\right),D_m\left(t\right),{\bar{D}}_m\left(t\right),E_m\left(t\right),{\bar{E}}_m\left(t\right)\left(m=1,2,\cdots ,M\right)$ are non-negative continuous monotonic non-increasing functions on $R_{+}$, and: ${\alpha }_i\left(i=1,2,\cdots ,I\right)$, ${\beta }_k\left(k=1,2,\cdots ,K\right)$, ${\gamma }_m\left(m=1,2,\cdots ,M\right)$ are all constants greater than 1, $1\le {\alpha }_1\le {\alpha }_2\le \cdots \le {\alpha }_I$, $1\le {\beta }_1\le {\beta }_2\le \cdots \le {\beta }_K$, $1\le {\gamma }_1\le {\gamma }_2\le \cdots \le {\gamma }_M$ let: $\bar{\alpha }=\max \left({\alpha }_I,{\beta }_K,{\gamma }_M\right)$, $\underset{\_}{\alpha }=\min \left({\alpha }_1,{\beta }_1,{\gamma }_1\right)$

iv) Set:

$\begin{array}{c}\Delta \left(t\right)=F\left(t\right)+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{H}_j\left(t\right)}{\displaystyle {\int }_{t_0}^t{N}_j\left(s\right)\text{d}s}+2{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right)}{\displaystyle {\int }_{t_0}^t{R}_k\left(s\right)\text{d}s}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l\left(s\right)}{\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }\text{d}s}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}\end{array}$

$\begin{array}{c}\Phi \left(t\right)={\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}+{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right)}{\displaystyle {\int }_{t_0}^t{R}_k\left(s\right)\text{d}s}+{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l(s)}{\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }\text{d}s}\\ +2{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right)}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}\end{array}$

$\Gamma \left(t\right)={\displaystyle {\int }_{t_0}^t{\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(s\right)}{\displaystyle {\int }_{t_0}^s{Z}_m\left(\sigma \right)\text{d}\sigma }\text{d}s}$

${\Lambda }_1\left(t\right)=1-\bar{\alpha }{c}^{\bar{\alpha }}{\displaystyle {\int }_{t_0}^t\left[\Phi \left(\tau \right)+\Gamma \left(\tau \right)\right]\exp \left(\bar{\alpha }{\displaystyle {\int }_{t_0}^{\tau }\Delta \left(\sigma \right)\text{d}\sigma }\right)\text{d}\tau }$

$\Pi \left(t\right)=F\left(t\right)+{\displaystyle \underset{j=1}{\overset{J}{\sum }}{H}_j\left(t\right){\displaystyle {\int }_{t_0}^t{N}_j\left(t\right)\text{d}t}}$

$\begin{array}{c}\Theta \left(t\right)={\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right){Щ}^{{\alpha }_i-\underset{\_}{\alpha }}}+{\displaystyle \underset{k=1}{\overset{K}{\sum }}{Q}_k\left(t\right){\displaystyle {\int }_{t_0}^t{R}_k\left(s\right){Щ}^{{\beta }_k-\underset{\_}{\alpha }}\text{d}s}}\\ +{\displaystyle \underset{i=1}{\overset{I}{\sum }}{G}_i\left(t\right){Щ}^{{\alpha }_i-\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{l=1}{\overset{L}{\sum }}{T}_l\left(s\right){\displaystyle {\int }_{t_0}^s{W}_l\left(\sigma \right)\text{d}\sigma }}}}\end{array}$

$\Sigma \left(t\right)={\displaystyle \underset{m=1}{\overset{M}{\sum }}{Y}_m\left(t\right){\displaystyle {\int }_{t_0}^t{Z}_m\left(s\right){Щ}^{{\gamma }_m-\underset{\_}{\alpha }}\text{d}s}}$

${\Lambda }_2\left(t\right)=1-\underset{\_}{\alpha }{c}^{\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t\left[\Theta \left(\tau \right)+\Sigma \left(\tau \right)\right]\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_{t_0}^{\tau }\Pi \left(\sigma \right)\text{d}\sigma }\right)\text{d}\tau }$

and ${\displaystyle {\int }_{t_0}^{+\infty }\Delta \left(s\right)\text{d}s<+\infty }$,

${\Lambda }_1\left(t\right)>0$, $\left[1-\left(\bar{\alpha }-1\right)c^{\bar{\alpha }}{\displaystyle {\int }_{t_0}^t\Phi \left(s\right){\Lambda }_1^{-\frac{1}{\bar{\alpha }}}\left(s\right)\exp \left(\bar{\alpha }{\displaystyle {\int }_{t_0}^s\Delta \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]>0$

${\Lambda }_2\left(t\right)>0$, $\left[1-\left(\underset{\_}{\alpha }-1\right)c^{\underset{\_}{\alpha }}{\displaystyle {\int }_{t_0}^t\Theta \left(s\right){\Lambda }_2^{-\frac{1}{\underset{\_}{\alpha }}}\left(s\right)\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_{t_0}^s\Sigma \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]>0$

v) Set:

$\Omega \left(t\right)=k\exp \left({\displaystyle {\int }_0^t\Pi \left(s\right)\text{d}s}\right)\cdot {\left[1-\left(\underset{\_}{\alpha }-1\right)k^{\underset{\_}{\alpha }}{\displaystyle {\int }_0^t\Theta \left(s\right){\Lambda }_2^{-\frac{1}{\underset{\_}{\alpha }}}\left(s\right)\exp \left(\underset{\_}{\alpha }{\displaystyle {\int }_0^s\Pi \left(\tau \right)\text{d}\tau }\right)\text{d}s}\right]}^{-\frac{1}{\underset{\_}{\alpha }-1}}$ then:

1) when $\lambda >\epsilon$, the trivial solution of (2) is $LS$, $G{E}_{q}ELAS$ with respect to y;

2) when $\lambda =\epsilon$, the trivial solution of (2) is, $GUELAS$ with respect to.

Proof Apply constant variation method to system (2), it can be deduced that:

$y\left(t\right)=Y\left(t,t_0\right)y_0+{\displaystyle {\int }_{t_0}^{t}Y\left(t,s\right)F_1\left(s,x\left(s-\delta \left(s\right)\right),{\displaystyle {\int }_{t_0}^s{h}_1\left(\tau ,x\left(\tau -\delta \left(\tau \right)\right)\right)\text{d}\tau }\right)\text{d}s}$ (3)

$z\left(t\right)=Z\left(t,t_0\right)z_0+{\displaystyle {\int }_{t_0}^{t}Z\left(t,s\right)F_2\left(s,x\left(s-\delta \left(s\right)\right),{\displaystyle {\int }_{t_0}^s{h}_2\left(\tau ,x\left(\tau -\delta \left(\tau \right)\right)\right)\text{d}\tau }\right)\text{d}s}$ (4)

By the condition of the theory, available from (3):

(5)

however, set

$u_1\left(t\right)=\left(\Vert y\left(t\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t\right)\Vert \right){\text{e}}^{\lambda \left(t-t_0\right)}$, $u_2\left(t\right)=\Vert z\left(t\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(t\right)\Vert$

${\phi }_1^{\left(1\right)}=\underset{-\delta \le t\le 0}{\sup }\Vert {\phi }_1\left(t\right)\Vert$, ${\phi }_2^{\left(1\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_1\left(t\right)\Vert }^{{\alpha }_i+1}$

${\phi }_3^{\left(1\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_1\left(t\right)\Vert }^{{\beta }_k+1}$, ${\phi }_4^{\left(1\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_1\left(t\right)\Vert }^{{\gamma }_m+1}$

${\phi }^{\left(1\right)}=\max \left\{{\phi }_1^{\left(1\right)},{\phi }_2^{\left(1\right)},{\phi }_3^{\left(1\right)},{\phi }_4^{\left(1\right)},{\phi }_2^{\left(1\right)}{\phi }_4^{\left(1\right)}\right\}$

${\phi }_1^{\left(2\right)}=\underset{-\delta \le t\le 0}{\sup }\Vert {\phi }_2\left(t\right)\Vert$, ${\phi }_2^{\left(2\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_2\left(t\right)\Vert }^{{\alpha }_i+1}$

${\phi }_3^{\left(2\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_2\left(t\right)\Vert }^{{\beta }_k+1}$, ${\phi }_4^{\left(2\right)}=\underset{-\delta \le t\le 0}{\sup }{\Vert {\phi }_2\left(t\right)\Vert }^{{\gamma }_m+1}$

${\phi }^{\left(2\right)}=\max \left\{{\phi }_1^{\left(2\right)},{\phi }_2^{\left(2\right)},{\phi }_3^{\left(2\right)},{\phi }_4^{\left(2\right)},{\phi }_2^{\left(2\right)}{\phi }_4^{\left(2\right)}\right\}$

Set:

$\Vert u_1\left(t\right)\Vert =\{\begin{array}{l}\max \left\{{\phi }_1^{\left(1\right)},\max \left(u\left(\xi \right)\right)\right\},0\le \xi \le t\\ {\phi }_1^{\left(1\right)},-\delta \le t\le 0\end{array}$

$\Vert u_2\left(t\right)\Vert =\{\begin{array}{l}\max \left\{{\phi }_1^{\left(2\right)},\max \left(u\left(\xi \right)\right)\right\},0\le \xi \le t\\ {\phi }_1^{\left(2\right)},-\delta \le t\le 0\end{array}$

Obviously, $\Vert u_1\left(t\right)\Vert$, $\Vert u_2\left(t\right)\Vert$ are monotonous, and defined by, $u_1\left(t\right)$, $u_2\left(t\right)$, we have

$u_1\left(t-\delta \left(t\right)\right)\le {\phi }_1^{\left(1\right)}$,

$u_2\left(t-\delta \left(t\right)\right)\le {\phi }_1^{(\; 2\; )}$

So substituting (5) into (2) gives:

$\begin{array}{c}{u}_1\left(t\right)\\ \le {\phi }^{\left(1\right)}+M_2{\displaystyle {\int }_{t_0}^{t}a\left(s\right)\Vert u_1\left(s\right)\Vert \text{d}s}+M_2{\displaystyle {\int }_{t_0}^{t}b\left(s\right)\Vert u_2\left(s\right)\Vert {\text{e}}^{\left(\lambda -\epsilon \right)\left(s-t_0\right)}\text{d}s}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right){\Vert u_1\left(s\right)\Vert }^{{\alpha }_i+1}{\text{e}}^{-{\alpha }_i\left(\lambda -\epsilon \right)\left(s-t_0\right)}\text{d}s}}+M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right){\Vert u_2\left(s\right)\Vert }^{{\alpha }_i+1}{\text{e}}^{\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}\text{d}s}}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{e}_j\left(s\right){\displaystyle {\int }_{t_0}^s{f}_j\left(\theta \right)\Vert u_1\left(\theta \right)\Vert \text{d}\theta \text{d}s}}}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{j=1}{\overset{J}{\sum }}{g}_j\left(s\right){\displaystyle {\int }_{t_0}^s{h}_j\left(\theta \right)\Vert u_2\left(\theta \right)\Vert {\text{e}}^{\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}\text{d}\theta \text{d}s}}}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{o}_k\left(s\right){\displaystyle {\int }_{t_0}^s{p}_k\left(\theta \right){\Vert u_1\left(\theta \right)\Vert }^{{\beta }_k+1}}{\text{e}}^{-{\beta }_k\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}\text{d}\theta \text{d}s}}\end{array}$

$\begin{array}{l}+M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{k=1}{\overset{K}{\sum }}{q}_k\left(s\right){\displaystyle {\int }_{t_0}^s{r}_k\left(\theta \right){\Vert u_2\left(\theta \right)\Vert }^{{\beta }_k+1}}{\text{e}}^{\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}\text{d}\theta \text{d}s}}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right){\Vert u_1\left(s\right)\Vert }^{{\alpha }_i}{\text{e}}^{-{\alpha }_i\left(\lambda -\epsilon \right)\left(s-t_0\right)}}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{v}_l\left(\sigma \right)}{\displaystyle {\int }_{t_0}^{\sigma }{w}_l\left(\theta \right)\Vert u_1\left(\theta \right)\Vert \text{d}\theta \text{d}\sigma \text{d}s}}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right){\Vert u_2\left(s\right)\Vert }^{{\alpha }_i}}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{l=1}{\overset{L}{\sum }}{x}_l\left(\sigma \right)}{\displaystyle {\int }_{t_0}^{\sigma }{y}_l\left(\theta \right)\Vert u_2\left(\theta \right)\Vert {\text{e}}^{\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}}\text{d}\theta \text{d}\sigma \text{d}s}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{c}_i\left(s\right){\Vert u_2\left(s\right)\Vert }^{{\alpha }_i}{\text{e}}^{-{\alpha }_i\left(\lambda -\epsilon \right)\left(s-t_0\right)}}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{A}_m\left(\sigma \right)}{\displaystyle {\int }_{t_0}^{\sigma }{B}_m\left(\theta \right){\Vert u_1\left(\theta \right)\Vert }^{{\gamma }_m+1}}{\text{e}}^{-{\gamma }_m\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}\text{d}\theta \text{d}\sigma \text{d}s}\\ +M_2{\displaystyle {\int }_{t_0}^t{\displaystyle \underset{i=1}{\overset{I}{\sum }}{d}_i\left(s\right){\Vert u_2\left(s\right)\Vert }^{{\alpha }_i}}}{\displaystyle {\int }_{t_0}^s{\displaystyle \underset{m=1}{\overset{M}{\sum }}{D}_m\left(\sigma \right)}{\displaystyle {\int }_{t_0}^{\sigma }{E}_m\left(\theta \right){\Vert u_2\left(\theta \right)\Vert }^{{\gamma }_m+1}{\text{e}}^{\left(\lambda -\epsilon \right)\left(\theta -t_0\right)}}\text{d}\theta \text{d}\sigma \text{d}s}\end{array}$

Similarly available:

where $M_2=\max \left\{M_1,M_1^{{\alpha }_i+1},M_1^{{\beta }_k+1},M_1^{{\gamma }_m+1},M_1^{{\alpha }_i+{\gamma }_m+1}\right\}$.

So it can be obtained from Lemma: $u_1\left(t\right)\le k{\text{e}}^{\left(\lambda -\epsilon \right)\left(t-t_0\right)}\Omega \left(t\right)$, $u_2\left(t\right)\le k\Omega \left(t\right)$ here $k={\text{e}}^{\lambda t_0}\varphi$, $\varphi =\max \left({\phi }^{\left(1\right)},{\phi }^{\left(2\right)}\right)$, $\Omega \left(t\right)$ As the lemma states, then:

$\Vert y\left(t\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t\right)\Vert \le M_3\phi {\text{e}}^{-\epsilon \left(t-t_0\right)}\left(\lambda >\epsilon \right)$, $\Vert y\left(t\right)\Vert +\Vert \stackrel{\dot{}}{y}\left(t\right)\Vert \le M_4\phi {\text{e}}^{-\epsilon \left(t-t_0\right)}\left(\lambda =\epsilon \right)$

$\Vert z\left(t\right)\Vert +\Vert \stackrel{\dot{}}{z}\left(t\right)\Vert \le M_5\phi$ (6)

here: $M_3={\text{e}}^{\left(\lambda -\epsilon \right)t_0}\Omega \left(t\right)$ ; $M_4=\Omega \left(t\right)$ ; $M_5={\text{e}}^{\lambda t_0}\Omega \left(t\right)$.

Notice the theorem conditions, we have $M_1$ is a constant that has nothing to do with $t_0$, $M_2$ and $M_3$ are constants that has nothing to do with $t_0$.

Therefore, when $\lambda >\epsilon$, (6) means the trivial solution of (2) is $LS$, $G{E}_{q}ELAS$ with respect to y; when $\lambda =\epsilon$, (6) means the trivial solution of (2) is $LS$, $GUELAS$ with respect to y.

Note: The differential system discussed in this paper is the time-differential form of the ordinary differential system in [5] [6]. The time-differential system in [7] is generalized to a neutral system, and the Lipschitz stability in [8] is further extended to equi-exponential Lipschitz asymptotic stability and uniform exponential Lipschitz asymptotic stability and added global results

4. Conclusion

In this paper, we use the method of integral inequalities to establish double stability criteria. As a result, studying the partial stability of differential equations becomes more important. In addition, the partial stability of differential equations is widely used in science and technology.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Lyapunov, A.M. (1892) The General Problem of the Stability of Motion. Doctor dissertation, University of Moscow, Moscow.
[2]	Li, Y.S. (1960) Basic Inequalities and Uniqueness of Solutions of Differential Equations. Journal of Jilin University, 1, 7-12.
[3]	Si, L.G. (2005) A Note on a Class of Integral Inequalities. Journal of Inner Mongolia Normal University, 34, 127-128.
[4]	Huo, R., Wang, X.L. and Si, L.G. (2009) A Generalization of a Class of Integral Inequality. Journal of Inner Mongolia Agricultral University (Natural Science Edition), 30, 267-273.
[5]	Vorotnikov, V.I. (1999) On the Theory of Lyapunov Stablity in Critical Cases. Dokl. Akad. Nank., 367, 481-484.
[6]	Vorotnikov, V.I. (1999) On Problems of Stablity with Respect to Some of the Variables. Prikl. Mat. Meth., 63, 736-745. https://doi.org/10.1016/S0021-8928(99)00088-X
[7]	Wang, F. (2002) Stability of Differential Equations on Some Variables. Inner Mongolia Normal University, Hohhot.
[8]	Si, L.G. (2006) The Stability of a Class of Nonlinear Differential Systems with Respect to Partial Variables. Journal of Inner Mongolia Normal University, 35, 1-6.
[9]	Si, L.G. (1974) The Boundness and Stablity of Solutions for Nonlinear Neutral Differential Systems with Variable Delay. Acta Mathematica Sinica, 17, 197-204.
[10]	Zhang, Y. (1994) K-Stablity of Nonlinear System with Time-Delay. Scientia Sinica Mathematica, 37, 247-255.
[11]	Si, L.G. (2001) Differential Inequality and Differential Quation with Time-Delay, Inner Mongolia People’s Press, Hohhot.