báo cáo hóa học:" Research Article Solvability of Nonautonomous Fractional Integrodifferential Equations with Infinite Delay" ppt

18 258 0
báo cáo hóa học:" Research Article Solvability of Nonautonomous Fractional Integrodifferential Equations with Infinite Delay" ppt

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

Hindawi Publishing Corporation Advances in Difference Equations Volume 2011, Article ID 806729, 18 pages doi:10.1155/2011/806729 Research Article Solvability of Nonautonomous Fractional Integrodifferential Equations with Infinite Delay Fang Li School of Mathematics, Yunnan Normal University, Kunming 650092, China Correspondence should be addressed to Fang Li, fangli860@gmail.com Received 4 September 2010; Revised 19 October 2010; Accepted 29 October 2010 Academic Editor: Toka Diagana Copyright q 2011 Fang Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. We study the existence and uniqueness of mild solution of a class of nonlinear nonautonomous fractional integrodifferential equations with infinite delay in a Banach space X. The existence of mild solution is obtained by using the theory of the measure of noncompactness and Sadovskii’s fixed point theorem. An application of the abstract results is also given. 1. Introduction The Cauchy problem for various delay equations in Banach spaces has been receiving more and more attention during the past decades cf., e.g., 1–15. This paper is concerned with existence results for nonautonomous fractional integrodifferential equations with infinite delay in a Banach space X d q u  t  dt q  −A  t  u  t   f  t,  t 0 K  t, s  u  s  ds, u t  ,t∈  0,T  , u  t   φ  t  ,t∈  −∞, 0  , 1.1 where T>0, 0 <q<1, {At} t∈0,T is a family of linear operators in X, K ∈ CD, R   with D  {t, s ∈ R 2 :0≤ s ≤ t ≤ T} and sup t∈0,T  t 0 K  t, s  ds < ∞, 1.2 2 Advances in Difference Equations f : 0,T × X ×P → X, u t : −∞, 0 → X defined by u t θut  θ for θ ∈ −∞, 0, φ belongs to the phase space P,andφ00. The fractional derivative is understood here in the Riemann-Liouville sense. In recent years, the fractional differential equations have been proved to be good tools in the investigation of many phenomena in engineering, physics, economy, chemistry, aerodynamics, electrodynamics of complex medium and they have been studied by many researchers cf., e.g., 13, 14, 16, 17 and references therein. Moreover, many phenomena cannot be described through classical differential equations but the integral and integrodifferential equations in abstract spaces in fields like electronic, fluid dynamics, biological models, and chemical kinetics. So many significant works on this topic have been appeared cf., e.g., 10, 15, 18–25 and references therein. In this paper, we study the existence of mild solution of 1.1 and obtain the existence theorem based on the measures of noncompactness without the assumptions that the nonlinearity f satisfies a Lipschitz type condition and the semigroup {exp−tAs} generated by −Ass ∈ 0,T is compact see Theorem 3.1 . An example is given to show an application of the abstract results. 2. Preliminaries Throughout this paper, we set J 0,T, a compact interval in R. We denote by X a Banach space, LX the Banach space of all linear and bounded operators on X,andCJ, X the space of all X-valued continuous functions on J.Weset Gu  t  :  t 0 K  t, s  u  s  ds, G ∗ : sup t∈J  t 0 K  t, s  ds < ∞. 2.1 Next, we recall the definition of the Riemann-Liouville integral. Definition 2.1 see 26. The fractional arbitrary order integral of the function g ∈ L 1 R  , R of order ν>0 is defined by I ν g  t   1 Γ  ν   t 0  t − s  ν−1 g  s  ds, 2.2 where Γ is the Gamma function. Moreover, I ν 1 I ν 2  I ν 1 ν 2 , for all ν 1 ,ν 2 > 0. Remark 2.2. 1 I ν : L 1 0,T → L 1 0,Tsee 26, 2 obviously, for g ∈ L 1 J, R, it follows from Definition 2.1 that  t 0  η 0  t − η  q−1  η − s  γ−1 g  s  ds dη  B  q, γ   t 0  t − s  qγ−1 g  s  ds, 2.3 where Bq, γ is a beta function. See the following definition about phase space according to Hale and Kato 27. Advances in Difference Equations 3 Definition 2.3. A linear space P consisting of functions from R − into X, with seminorm · P , is called an admissible phase space if P has the following properties. 1 If x : −∞,T → X is continuous on J and x 0 ∈P, then x t ∈Pand x t is continuous in t ∈ J,and  x  t   ≤ L  x t  P , 2.4 where L ≥ 0 is a constant. 2 There exist a continuous function C 1 t > 0 and a locally bounded function C 2 t ≥ 0 in t ≥ 0, such that  x t  P ≤ C 1  t  sup s∈  0,t   x  s    C 2  t   x 0  P , 2.5 for t ∈ 0,T and x as in 1. 3 The space P is complete. Remark 2.4. Equation 2.4 in 1 is equivalent to φ0≤Lφ P , for all φ ∈P. Next, we consider the properties of Kuratowski’s measure of noncompactness. Definition 2.5. Let B be a bounded subset of a seminormed linear space Y . The Kuratowski’s measure of noncompactnessfor brevity, α-measure of B is defined as α  B   inf  d>0:B has a finite cover by sets of diameter ≤ d  . 2.6 From the definition we can get some properties of α-measure immediately, see 28. Lemma 2.6 see 28. Let A and B be bounded subsets of X,then 1 αA ≤ αB,ifA ⊆ B; 2 αAα A,whereA denotes the closure of A; 3 αA0 if and only if A is precompact; 4 αλA|λ|αA,λ ∈ R; 5 αA ∪ Bmax{αA,αB}; 6 αA  B ≤ αAαB,whereA  B  {x  y : x ∈ A, y ∈ B}; 7 αA  xα A, for any x ∈ X. For H ⊂ CJ, X, we define  t 0 H  s  ds    t 0 u  s  ds : u ∈ H  , for t ∈ J, 2.7 where Hs{us ∈ X : u ∈ H}. 4 Advances in Difference Equations The following lemma will be needed. Lemma 2.7. If H ⊂ CJ, X is a bounded, equicontinuous set, then i αHsup t∈J αHt, ii α  t 0 Hsds ≤  t 0 αHsds,fort ∈ J. For a proof refer to 28. Lemma 2.8 see 29. If {u n } ∞ n1 ⊂ L 1 J, X and there exists an m ∈ L 1 J, R   such that u n t≤ mt,a.e.t ∈ J,thenα{u n t} ∞ n1  is integrable and α   t 0 u n sds  ∞ n1  ≤ 2  t 0 α  { u n  s  } ∞ n1  ds. 2.8 We need to use the following Sadovskii’s fixed point theorem here, see 30. Definition 2.9. Let P be an operator in Banach space X.IfP is continuous and takes bounded sets into bounded sets, and αPB <αB for every bounded set B of X with αB > 0, then P is said to be a condensing operator on X. Lemma 2.10 Sadovskii’s fixed point theorem 30. Let P be a condensing operator on Banach space X.IfP H ⊆ H for a convex, closed, and bounded set H of X,thenP has a fixed point in H. In this paper, we denote that C is a positive constant, and assume that a family of closed linear operators {At} t∈J satisfying the following. A1 The domain DA of {At} t∈J is dense in the Banach space X and independent of t. A2 The operator Atλ −1 exists in LX for any λ with Re λ ≤ 0and     A  t   λ  −1    ≤ C | λ |  1 ,t∈ J. 2.9 A3 There exist constants γ ∈ 0, 1 and C such that     A  t 1  − A  t 2  A −1  s     ≤ C | t 1 − t 2 | γ ,t 1 ,t 2 ,s∈ J. 2.10 Under condition A2, each operator −As, s ∈ J, generates an analytic semigroup exp−tAs, t>0, and there exists a constant C such that   A n  s  exp  −tA  s    ≤ C t n , 2.11 where n  0, 1, t>0, s ∈ J see 31. Advances in Difference Equations 5 Let Ω be set defined by Ω  u :  −∞,T  −→ X such that u | −∞,0 ∈Pand u | J ∈ C  J, X   . 2.12 According to 16, a mild solution of 1.1 can be defined as follows. Definition 2.11. A function u ∈ Ω satisfying the equation u  t   ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ φ  t  ,t∈  −∞, 0  ,  t 0 ψ  t − η, η  f  η, Gu  η  ,u η  dη   t 0  η 0 ψ  t − η, η  ϕ  η, s  f  s, Gu  s  ,u s  ds dη, t ∈ J 2.13 is called a mild solution of 1.1, where ψ  t, s   q  ∞ 0 θt q−1 ξ q  θ  exp  −t q θA  s  dθ , 2.14 and ξ q is a probability density function defined on 0, ∞ such that its Laplace transform is given by  ∞ 0 e −σx ξ q  σ  dσ  ∞  j0  −x  j Γ  1  qj  , 0 <q≤ 1,x>0, ϕ  t, τ   ∞  k1 ϕ k  t, τ  , 2.15 where ϕ 1  t, τ    A  t  − A  τ  ψ  t − τ, τ  , ϕ k1  t, τ    t τ ϕ k  t, s  ϕ 1  s, τ  ds, k  1, 2, 2.16 To our purpose the following conclusions will be needed. For the proofs refer to 16. Lemma 2.12 see 16. The operator-valued functions ψt−η, η and Atψt−η,η are continuous in uniform topology in the variables t, η,where0 ≤ η ≤ t − ε, 0 ≤ t ≤ T, for any ε>0. Clearly,   ψ  t − η, η    ≤ C  t − η  q−1 . 2.17 Moreover, we have   ϕ  t, η    ≤ C  t − η  γ−1 . 2.18 6 Advances in Difference Equations 3. Main Results We need the hypotheses as follows: H1 f : J × X ×P → X satisfies f·,v,w : J → X is measurable for all v, w ∈ X ×P and ft, ·, · : X ×P → X is continuous for a.e. t ∈ J, and there exist a positive function μ· ∈ L p J, R  p>1/q > 1 and a continuous nondecreasing function W : 0, ∞ → 0, ∞, such that   f  t, v, w    ≤ μ  t  W   v    w  P  ,  t, v, w  ∈ J × X ×P, 3.1 and set T p,q : T q−1/p , H2 for any bounded sets D 1 ⊂ X, D 2 ⊂P,and0≤ τ ≤ s ≤ t ≤ T, α  ψ  t − s, s  f  s, D 1 ,D 2   ≤ β 1  t, s  α  D 1   β 2  t, s  sup −∞<θ≤0 α  D 2  θ  , α  ψ  t − s, s  ϕ  s, τ  f  τ,D 1 ,D 2   ≤ β 3  t, s, τ  α  D 1   β 4  t, s, τ  sup −∞<θ≤0 α  D 2  θ  , 3.2 where sup t∈J  t 0 β i t, sds : β i < ∞ i  1, 2,sup t∈J  t 0  s 0 β i t, s, τdτ ds : β i < ∞ i  3, 4 and D 2 θ{wθ : w ∈ D 2 }, H3 there exists M,with0< M<1 such that C  1  CB  q, γ  T p,q,γ M p,q  G ∗  C ∗ 1    μ   L p J,R   lim inf τ →∞ W  τ  τ < M, 3.3 where M p,q :p − 1/pq − 1 p−1/p , C ∗ 1  sup 0≤η≤T C 1 η and T p,q,γ  max{T p,q , T p,qγ }. Theorem 3.1. Suppose that H1–H3 are satisfied, and if 4G ∗ β 1  2β 3 β 2  2β 4  < 1,then for 1.1 there exists a mild solution on −∞,T. Proof. Consider the operator Φ : Ω → Ω defined by  Φu  t   ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ φ  t  ,t∈  −∞, 0  ,  t 0 ψ  t − η, η  f  η, Gu  η  ,u η  dη   t 0  η 0 ψ  t − η, η  ϕ  η, s  f  s, Gu  s  ,u s  ds dη, t ∈ J. 3.4 It is easy to see that Φ is well-defined. Let x· : −∞,T → X be the function defined by x  t   ⎧ ⎨ ⎩ φ  t  ,t∈  −∞, 0  , 0,t∈ J. 3.5 Advances in Difference Equations 7 Let ut xtyt, t ∈ −∞,T. It is easy to see that y satisfies y 0  0and y  t    t 0 ψ  t − η, η  f  η, G  x  η   y  η  , x η  y η  dη   t 0  η 0 ψ  t − η, η  ϕ  η, s  f  s, G  x  s   y  s   , x s  y s  ds dη, t ∈ J 3.6 if and only if u satisfies u  t    t 0 ψ  t − η, η  f  η, Gu  η  ,u η  dη   t 0  η 0 ψ  t − η, η  ϕ  η, s  f  s, Gu  s  ,u s  ds dη, t ∈ J 3.7 and utφt, t ∈ −∞, 0. Let Y 0  {y ∈ Ω : y 0  0}. For any y ∈ Y 0 ,   y   Y 0  sup t∈J   y  t       y 0   P  sup t∈J   y  t    . 3.8 Thus Y 0 , · Y 0  is a Banach space. We define the operator  Φ : Y 0 → Y 0 by   Φyt0, t ∈ −∞, 0 and   Φy   t    t 0 ψ  t − η, η  f  η, G  x  η   y  η  , x η  y η  dη   t 0  η 0 ψ  t − η, η  ϕ  η, s  f  s, G  x  s   y  s   , x s  y s  ds dη, t ∈ J. 3.9 Obviously, the operator Φ has a fixed point if and only if  Φ has a fixed point. So it turns out to prove that  Φ has a fixed point. Let {y k } k∈N be a sequence such that y k → y in Y 0 as k →∞. Since f satisfies H1, for almost every t ∈ J,weget f  t, G  x  t   y k  t   , x t  y k t  −→ f  t, G  x  t   y  t   , x t  y t  , as k −→ ∞ . 3.10 8 Advances in Difference Equations For t ∈ −∞,T, we can prove that  Φ is continuous. In fact,      Φy k   t  −   Φy   t     ≤  t 0    ψ  t − η, η   f  η, G  x  η   y k  η   , x η  y k η  − f  η, G  x  η   y  η  , x η  y η      dη   t 0  η 0    ψ  t − η, η  ϕ  η, s   f  s, G  x  s   y k  s   , x s  y k s  − f  s, G  x  s   y  s   , x s  y s      ds dη. 3.11 Let C ∗ 2  sup 0≤η≤T C 2 η, and noting 2.4, 2.5, we have   x t  y t   P ≤  x t  P    y t   P ≤ C 1  t  sup 0≤τ≤t  x  τ    C 2  t   x 0  P  C 1  t  sup 0≤τ≤t   y  τ     C 2  t    y 0   P  C 2  t    φ   P  C 1  t  sup 0≤τ≤t   y  τ    ≤ C ∗ 2   φ   P  C ∗ 1 sup 0≤τ≤t   y  τ    . 3.12 Moreover,   G  x  t   y  t     ≤  t 0 K  t, τ    x  τ   y  τ    dτ   t 0 K  t, τ  ·   y  τ    dτ. 3.13 Noting that y k → y in Y 0 , we can see that there exists ε>0 such that y k − y≤ε, for k sufficiently large. Therefore, we have    f  t, G  x  t   y k  t   , x t  y k t  − f  t, G  x  t   y  t   , x t  y t     ≤ μ  t   W     G  x  t   y k  t          x t  y k t    P   W    G  x  t   y  t        x t  y t   P   ≤ μ  t   ω 1 k  t   ω 2 k  t   , 3.14 Advances in Difference Equations 9 where ω 1 k  t   W  G ∗ ε  G ∗   y   Y 0  C ∗ 2   φ   P  C ∗ 1 ε  C ∗ 1   y   Y 0  , ω 2 k  t   W  G ∗   y   Y 0  C ∗ 2   φ   P  C ∗ 1   y   Y 0  . 3.15 In view of 2.17 and the Lebesgue Dominated Convergence Theorem ensure that  t 0    ψ  t − η, η   f  η, G  x  η   y k  η   , x η  y k η  − f  η, G  x  η   y  η  , x η  y η      dη ≤ C  t 0  t − η  q−1    f  η, G  x  η   y k  η   , x η  y k η  − f  η, G  x  η   y  η  , x η  y η    dη −→ 0, as k −→ ∞ . 3.16 Similarly,by 2.17 and 2.18 , we have  t 0  η 0    ψ  t − η, η  ϕ  η, s   f  s, G  x  s   y k  s   , x s  y k s  − f  s, G  x  s   y  s   , x s  y s      ds dη ≤ C 2  t 0  η 0  t − η  q−1  η − s  γ−1    f  s, G  x  s   y k  s   , x s  y k s  − f  s, G  x  s   y  s   , x s  y s    ds dη −→ 0, as k −→ ∞ . 3.17 Therefore, we deduce that lim k →∞     Φy k −  Φy    Y 0  0. 3.18 This means that  Φ is continuous. We show that  Φ maps bounded sets of Y 0 into bounded sets in Y 0 . For any r>0, we set B r  {y ∈ Y 0 : y Y 0 ≤ r}.Now,fory ∈ B r ,by3.12, 3.13,andH1, we can see   f  t, G  x  t   y  t   , x t  y t    ≤ μ  t  W  M 1  , 3.19 where M 1 : G ∗ r  C ∗ 2 φ P  C ∗ 1 r. 10 Advances in Difference Equations Then for any y ∈ B r ,by2.17, 2.18, 3.19,andRemark 2.2, we have      Φy   t     ≤  t 0   ψ  t − η, η  f  η, G  x  η   y  η  , x η  y η    dη   t 0  η 0   ψ  t − η, η  ϕ  η, s  f  s, G  x  s   y  s   , x s  y s    ds dη ≤ C  t 0  t − η  q−1 μ  η  W  M 1  dη  C 2  t 0  η 0  t − η  q−1  η − s  γ−1 μ  s  W  M 1  ds dη  M 2  C  t 0  t − η  q−1 μ  η  dη  C 2 B  q, γ   t 0  t − η  qγ−1 μ  η  dη  , 3.20 where M 2  WM 1 . Noting that the H ¨ older inequality, we have  t 0  t − η  q−1 μ  η  dη  t pq−1/p M p,q   μ   L p J,R   ≤ T p,q M p,q   μ   L p J,R   ,  t 0  t − η  γq−1 μ  η  dη ≤ T p,qγ M p,qγ   μ   L p J,R   . 3.21 Thus      Φy   t     ≤ M 2 M p,q T p,q,γ  C  C 2 B  q, γ     μ   L p J,R   : r. 3.22 This means  ΦB r  ⊂ B r . Next, we show that there exists k ∈ N such that  ΦB k  ⊂ B k . Suppose contrary that for every k ∈ N there exist y k ∈ B k and t k ∈ J such that   Φy k t k  >k. However, on the other hand, similar to the deduction of 3.20 and noting    f  t, G  x  t   y k  t   , x t  y k t     ≤ μ  t  W  G ∗ k  C ∗ 2   φ   P  C ∗ 1 k  , 3.23 we have k<      Φy k   t k     ≤  M 2  C  t k 0  t k − η  q−1 μ  η  dη  C 2 B  q, γ   t k 0  t k − η  qγ−1 μ  η  dη  ≤  M 2 M p,q T p,q,γ  C  C 2 B  q, γ     μ   L p J,R   , 3.24 where  M 2  WG ∗ k  C ∗ 2 φ P  C ∗ 1 k. [...]... integrodifferential equations with infinite delay,” Applied Mathematics Letters, vol 17, no 4, pp 473– 477, 2004 11 J Liang and T.-J Xiao, Solvability of the Cauchy problem for infinite delay equations, ” Nonlinear Analysis: Theory, Methods & Applications, vol 58, no 3-4, pp 271–297, 2004 12 J Liang and T.-J Xiao, “Solutions to nonautonomous abstract functional equations with infinite delay,” Taiwanese Journal of Mathematics,... M Mophou and G M N’Gu´ r´ kata, “A note on a semilinear fractional differential equation of ee neutral type with infinite delay,” Advances in Difference Equations, vol 2010, Article ID 674630, 8 pages, 2010 14 G M Mophou and G M N’Gu´ r´ kata, “Existence of mild solutions of some semilinear neutral ee fractional functional evolution equations with infinite delay,” Applied Mathematics and Computation, vol... solutions of infinite delay evolution equations, ” Journal of Mathematical Analysis and Applications, vol 286, no 2, pp 705–712, 2003 6 J Liang and T J Xiao, “Functional-differential equations with infinite delay in Fr´ chet space,” Sichuan e Daxue Xuebao, vol 26, no 4, pp 382–390, 1989 7 J Liang and T J Xiao, “Functional-differential equations with infinite delay in Banach spaces,” International Journal of Mathematics... existence of solutions to integral equations with infinite delay in Banach spaces,” Nonlinear Analysis: Theory, Methods & Applications, vol 71, no 12, pp e1442– e1447, 2009 16 M M El-Borai, “The fundamental solutions for fractional evolution equations of parabolic type,” Journal of Applied Mathematics and Stochastic Analysis, no 3, pp 197–211, 2004 17 G M Mophou and G M N’Gu´ r´ kata, “Existence of the... mild solution for some fractional differential ee equations with nonlocal conditions,” Semigroup Forum, vol 79, no 2, pp 315–322, 2009 18 J.-H Chen, T.-J Xiao, and J Liang, “Uniform exponential stability of solutions to abstract Volterra equations, ” Journal of Evolution Equations, vol 9, no 4, pp 661–674, 2009 19 J Liang and T.-J Xiao, “Semilinear integrodifferential equations with nonlocal initial conditions,”... vol 192, no 1, pp 93–110, 2003 Advances in Difference Equations 17 3 Faming Guo, B Tang, and F Huang, “Robustness with respect to small delays for exponential stability of abstract differential equations in Banach spaces,” The ANZIAM Journal, vol 47, no 4, pp 555–568, 2006 4 J H Liu, “Periodic solutions of infinite delay evolution equations, ” Journal of Mathematical Analysis and Applications, vol 247,... supported by the NSF of Yunnan Province no 2009ZC054M References 1 Z Fan, “Existence and continuous dependence results for nonlinear differential inclusions with infinite delay,” Nonlinear Analysis: Theory, Methods & Applications, vol 69, no 8, pp 2379–2392, 2008 2 A Favini and L Vlasenko, “Degenerate non-stationary differential equations with delay in Banach spaces,” Journal of Differential Equations, vol... 497–508, 1991 8 J Liang and T J Xiao, “Solutions to abstract functional-differential equations with infinite delay,” Acta Mathematica Sinica, vol 34, no 5, pp 631–644, 1991 9 J Liang and T.-J Xiao, “The Cauchy problem for nonlinear abstract functional differential equations with infinite delay,” Computers & Mathematics with Applications, vol 40, no 6-7, pp 693–703, 2000 10 J Liang, T.-J Xiao, and J van... dependent Desch-Schappacher type perturbations of Volterra integral equations, ” Integral Equations and Operator Theory, vol 44, no 4, pp 494–506, 2002 26 S G Samko, A A Kilbas, and O I Marichev, Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach, Yverdon, Switzerland, 1993 27 J K Hale and J Kato, “Phase space for retarded equations with infinite delay,” Funkcialaj Ekvacioj, vol... vol 21, no 1, pp 11–41, 1978 28 J Bana´ and K Goebel, Measures of Noncompactness in Banach Spaces, vol 60 of Lecture Notes in Pure s and Applied Mathematics, Marcel Dekker, New York, NY, USA, 1980 18 Advances in Difference Equations 29 H.-P Heinz, “On the behaviour of measures of noncompactness with respect to differentiation and integration of vector-valued functions,” Nonlinear Analysis: Theory, Methods . in Difference Equations Volume 2011, Article ID 806729, 18 pages doi:10.1155/2011/806729 Research Article Solvability of Nonautonomous Fractional Integrodifferential Equations with Infinite Delay Fang. existence and uniqueness of mild solution of a class of nonlinear nonautonomous fractional integrodifferential equations with infinite delay in a Banach space X. The existence of mild solution is obtained. and integrodifferential equations with infinite delay,” Applied Mathematics Letters, vol. 17, no. 4, pp. 473– 477, 2004. 11 J. Liang and T J. Xiao, Solvability of the Cauchy problem for infinite delay equations, ”

Ngày đăng: 21/06/2014, 11:20

Từ khóa liên quan

Tài liệu cùng người dùng

Tài liệu liên quan