ON THE STABILITY OF THE QUASI-LINEAR IMPLICIT EQUATIONS IN HILBERT SPACES

Similar documents
ALMOST PERIODIC SOLUTIONS OF HIGHER ORDER DIFFERENTIAL EQUATIONS ON HILBERT SPACES

CESARO OPERATORS ON THE HARDY SPACES OF THE HALF-PLANE

MAIN ARTICLES. In present paper we consider the Neumann problem for the operator equation

DISCRETE METHODS AND EXPONENTIAL DICHOTOMY OF SEMIGROUPS. 1. Introduction

Conservative Control Systems Described by the Schrödinger Equation

BOUNDARY VALUE PROBLEMS IN KREĬN SPACES. Branko Ćurgus Western Washington University, USA

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR FOURTH-ORDER BOUNDARY-VALUE PROBLEMS IN BANACH SPACES

On (h, k) trichotomy for skew-evolution semiflows in Banach spaces

THE PERRON PROBLEM FOR C-SEMIGROUPS

1.17 : Consider a continuous-time system with input x(t) and output y(t) related by y(t) = x( sin(t)).

NONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality

ESTIMATES ON THE NUMBER OF LIMIT CYCLES OF A GENERALIZED ABEL EQUATION

EQUIVALENCE OF CONDITIONS FOR CONVERGENCE OF ITERATIVE METHODS FOR SINGULAR LINEAR SYSTEMS

arxiv: v1 [math.sp] 29 Jan 2018

c 1995 Society for Industrial and Applied Mathematics Vol. 37, No. 1, pp , March

Chapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors

Mathematical Journal of Okayama University

ON THE SOLUTION OF DIFFERENTIAL EQUATIONS WITH DELAYED AND ADVANCED ARGUMENTS

ASYMPTOTIC REPRESENTATIONS OF SOLUTIONS OF THE NONAUTONOMOUS ORDINARY DIFFERENTIAL n-th ORDER EQUATIONS

Sectorial Forms and m-sectorial Operators

APPROXIMATE CONTROLLABILITY OF DISTRIBUTED SYSTEMS BY DISTRIBUTED CONTROLLERS

LINEAR MAPS ON M n (C) PRESERVING INNER LOCAL SPECTRAL RADIUS ZERO

On the Well-Posedness of the Cauchy Problem for a Neutral Differential Equation with Distributed Prehistory

A quantitative Fattorini-Hautus test: the minimal null control time problem in the parabolic setting

A NOTE ON QUASI ISOMETRIES II. S.M. Patel Sardar Patel University, India

COMPLETENESS THEOREM FOR THE DISSIPATIVE STURM-LIOUVILLE OPERATOR ON BOUNDED TIME SCALES. Hüseyin Tuna

Analysis Preliminary Exam Workshop: Hilbert Spaces

Memoirs on Differential Equations and Mathematical Physics

arxiv: v1 [math.sp] 29 Jan 2018

Sung-Wook Park*, Hyuk Han**, and Se Won Park***

Stability of an abstract wave equation with delay and a Kelvin Voigt damping

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.

A Spectral Characterization of Closed Range Operators 1

2. Review of Linear Algebra

HYPO-EP OPERATORS 1. (Received 21 May 2013; after final revision 29 November 2014; accepted 7 October 2015)

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is

CONVERGENCE OF THE STEEPEST DESCENT METHOD FOR ACCRETIVE OPERATORS

Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control

A Concise Course on Stochastic Partial Differential Equations

Nonlinear Programming Algorithms Handout

ON OPERATORS WITH AN ABSOLUTE VALUE CONDITION. In Ho Jeon and B. P. Duggal. 1. Introduction

Spectral analysis for rank one perturbations of diagonal operators in non-archimedean Hilbert space

A G Ramm, Implicit Function Theorem via the DSM, Nonlinear Analysis: Theory, Methods and Appl., 72, N3-4, (2010),

Left invertible semigroups on Hilbert spaces.

ON THE NUMERICAL RANGE OF AN OPERATOR

ME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms

Math Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.

The i/s/o resolvent set and the i/s/o resolvent matrix of an i/s/o system in continuous time

Weaker hypotheses for the general projection algorithm with corrections

2 Statement of the problem and assumptions

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE

EXISTENCE RESULTS FOR OPERATOR EQUATIONS INVOLVING DUALITY MAPPINGS VIA THE MOUNTAIN PASS THEOREM

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

Stolz angle limit of a certain class of self-mappings of the unit disk

GORDON TYPE THEOREM FOR MEASURE PERTURBATION

EQUIVALENCE OF TOPOLOGIES AND BOREL FIELDS FOR COUNTABLY-HILBERT SPACES

T.8. Perron-Frobenius theory of positive matrices From: H.R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton 2003

On the simplest expression of the perturbed Moore Penrose metric generalized inverse

THE IMPLICIT FUNCTION THEOREM FOR CONTINUOUS FUNCTIONS. Carlos Biasi Carlos Gutierrez Edivaldo L. dos Santos. 1. Introduction

On m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry

ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.

Applied Mathematics Letters. Comparison theorems for a subclass of proper splittings of matrices

Quasireversibility Methods for Non-Well-Posed Problems

EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ

LOCAL FIXED POINT THEORY INVOLVING THREE OPERATORS IN BANACH ALGEBRAS. B. C. Dhage. 1. Introduction

AN EFFECTIVE METRIC ON C(H, K) WITH NORMAL STRUCTURE. Mona Nabiei (Received 23 June, 2015)

On the bang-bang property of time optimal controls for infinite dimensional linear systems

The Fine Spectra of the Difference Operator Δ over the Sequence Spaces l 1 and bv

On intermediate value theorem in ordered Banach spaces for noncompact and discontinuous mappings

GLOBAL ATTRACTOR FOR A SEMILINEAR PARABOLIC EQUATION INVOLVING GRUSHIN OPERATOR

Block companion matrices, discrete-time block diagonal stability and polynomial matrices

PERTURBATION OF A FREDHOLM COMPLEX BY INESSENTIAL OPERATORS

Two-Step Iteration Scheme for Nonexpansive Mappings in Banach Space

Finite groups determined by an inequality of the orders of their elements

arxiv: v2 [math.fa] 8 Jan 2014

Scalar Asymptotic Contractivity and Fixed Points for Nonexpansive Mappings on Unbounded Sets

On Solutions of Evolution Equations with Proportional Time Delay

Characterization of half-radial matrices

Applied Mathematics Letters

************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan

SOME INEQUALITIES FOR THE EUCLIDEAN OPERATOR RADIUS OF TWO OPERATORS IN HILBERT SPACES

ON HARMONIC FUNCTIONS ON SURFACES WITH POSITIVE GAUSS CURVATURE AND THE SCHWARZ LEMMA

1 Math 241A-B Homework Problem List for F2015 and W2016

On perturbations in the leading coefficient matrix of time-varying index-1 DAEs

arxiv: v1 [math.fa] 26 May 2012

EXISTENCE AND APPROXIMATION OF SOLUTIONS OF SECOND ORDER NONLINEAR NEUMANN PROBLEMS

NOTE ON HILBERT-SCHMIDT COMPOSITION OPERATORS ON WEIGHTED HARDY SPACES

THE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS

THE NEARLY ADDITIVE MAPS

Linear maps preserving AN -operators

Research Article The Solution by Iteration of a Composed K-Positive Definite Operator Equation in a Banach Space

Periodic solutions of differential equations with three variable in vector-valued space

Existence And Uniqueness Of Mild Solutions Of Second Order Volterra Integrodifferential Equations With Nonlocal Conditions

GENERALIZED POPOVICIU FUNCTIONAL EQUATIONS IN BANACH MODULES OVER A C ALGEBRA AND APPROXIMATE ALGEBRA HOMOMORPHISMS. Chun Gil Park

The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation

COMMON FIXED POINT THEOREMS FOR A PAIR OF COUNTABLY CONDENSING MAPPINGS IN ORDERED BANACH SPACES

16 1 Basic Facts from Functional Analysis and Banach Lattices

Weaker assumptions for convergence of extended block Kaczmarz and Jacobi projection algorithms

Transcription:

Khayyam J. Math. 5 (219), no. 1, 15-112 DOI: 1.2234/kjm.219.81222 ON THE STABILITY OF THE QUASI-LINEAR IMPLICIT EQUATIONS IN HILBERT SPACES MEHDI BENABDALLAH 1 AND MOHAMED HARIRI 2 Communicated by J. Brzdęk Abstract. We use the generalized theorem of Liapounov to obtain some necessary and sufficient conditions for the stability of the stationary implicit equation Ax (t) = Bx(t), t, where A and B are bounded operators in Hilbert spaces. The achieved results can be applied to the stability for the quasi-linear implicit equation Ax (t) = Bx(t) + θ(t, x(t)), t. 1. Introduction Consider the abstract implicit differential equation Ax (t) = Bx(t) + θ(t, x(t)), t, (1.1) where A and B are two linear bounded operators on a Hilbert space H and θ(, ) is a continuous function from [, ) H to H. The operator A is not necessarily invertible. The equation (1.1) has been considered in various forms by many authors as Favini and Yagi [5], Rutkas [6], Vlasenko [7] and others. In the present paper, we study the stationary implicit equation Ax (t) = Bx(t), t, (1.2) and also the quasi-linear implicit equation (1.1), with the initial condition x() = x. Date: Received: 1 May 218 ; Revised: 1 October 218; Accepted: 29 November 218. Corresponding author. 21 Mathematics Subject Classification. Primary 34D2; Secondary 47A1, 34K32. Key words and phrases. Exponential stability, operator theory, implicit equations. 15

16 M. BENABDALLAH, M. HARIRI In [2] the authors obtained results concerning the stability of the degenerate difference systems that is similar to (1.1). Some practical examples of (1.1) can be found in [5, 6, 7]. The organization of this paper is as follows: in Section 2, we introduce some preliminaries and expand the famous Liapounov general theorem [4], which has an important role in this paper. In section 3, we present our main results concerning the exponential stability of the solution for the quasi-linear implicit equation (1.1). We use the following definitions. Definition 1.1. The equation (1.2) is called exponentially stable, if there exist two constants M > and α < such that, for any solution x(t), we have x(t) Me αt x for any t. (1.3) Definition 1.2. The equation (1.2) is said to be well-posed, if it satisfies the following properties: (i) for any solution x(.) such that x() = x =, then x(t) = for all t ; (ii) it generates an evolution semigroup of bounded operators S(t) : x x(t) for all t. The operators S(t) are defined on the set D = {x } of the admissible initial vectors x. Definition 1.3 (see [6]). The complex number λ C is called a regular value of the pencil λa B, if the resolvent (λa B) 1 exists and is bounded. The set of all regular values is denoted by ρ(a, B) and its complement σ(a, B) = C\ρ(A, B) is called the spectrum of the pencil λa B. The set of all eigenvalues of the pencil λa B is denoted by σ p (A, B) = { λ C : v ; (λa B)v = }. 2. Stationary implicit equation For the stationary implicit equation (1.2), we can obtain the following criterion for the exponential stability. Theorem 2.1. The equation (1.2) is exponentially stable if and only if it is wellposed. Proof. Suppose that (1.2) is exponentially stable. Then, it has a unique solution x(t). In fact, if x =, then by (1.3), we obtain x(t), and consequently x(t) = for all t. On the other hand, we have x(t) = S(t)x Me αt x. It means that the operator S(t) is bounded and S(t) Me αt. So, (1.2) is well-posed. Conversely, if (1.2) is well-posed, then one obtains ln S(t) ω = lim < t t (see [4, p. 26]), where ω is the strict Liapounov exponent of Φ(t) = S(t). More precisely, ω is the greatest lower bound of the set of real numbers ϱ for which

EXPONENTIAL STABILITY 17 there exists a positive constant N ϱ (see [3, pp. 8 9]) such that Φ(t) = S(t) N ϱ e ϱt for all t, ω ϱ <. Hence which achieves the proof. x(t) = S(t)x N ϱ e ϱt x, Theorem 2.2. If (1.2) is exponentially stable, then all eigenvalues of the pencil λa B are in the half-plane Reλ α, where α is the constant defined in (1.3). Proof. Suppose that there exists an eigenvalue λ σ p (A, B) such that Reλ > α. Then (λ A B)v =, where v is the corresponding eigenvector. Consequently, y(t) = e λ t v is a solution of (1.2) verifying the condition y() = v, and we have y(t) = e λ t v = e (Reλ )t v > e αt y(). So, the solution y(t) does not satisfy (1.3) and consequently (1.2) is not exponentially stable. Remark 2.3. If (1.2) is exponentially stable, then all the eigenvalues of the pencil λa B are inside the left half-plane, that is since α <, where α is given by (1.3). σ p (A, B) {λ C : Reλ < }, We can now extend the generalized Liapounov theorem [4] for the spectrum of the bounded operator T, to the spectrum of the pencil λa B of the bounded operators A and B on a Hilbert space H, using the spectral theory of the pencil of operators and an appropriate conformal mapping as follows. Theorem 2.4. A necessary condition, for the spectrum σ(a, B) of the pencil λa B to lie in the interior of the half-plane Reλ < α (α < ), is that, for any uniformly positive operator U 1, there exists an operator W such that A W B + B W A 2αA W A = 2U, (2.1) and a sufficient condition is that α + 1 is a regular value of the pencil λa B and there exists an operator W such that A W B + B W A 2αA W A. (2.2) Proof. Necessary condition. Suppose that σ(a, B) {λ C : Reλ < α}. Then, (α + 1) ρ(a, B) and the operator T = [(α 1)A B][(α + 1)A B] 1 is well defined and bounded. Now, using the conformal mapping z = ϕ(λ) = λ α+1 λ α 1 which transforms the vertical line Reλ = α into the unit circle z = 1, we obtain 1 It means that U = U and that Ux, x > for all x with x = 1.

18 M. BENABDALLAH, M. HARIRI zi T = ( λ α+1 λ α 1 )[(α + 1)A B][ (α + 1)A B] 1 = 1 (λ α 1){ (λ α + 1)[(α + 1)A B] [(α 1)A B][(α + 1)A B ] 1 (λ α 1)[(α 1)A B] } [(α + 1)A B] 1 = 2 (λ α 1) (λa B)[(α + 1)A B] 1. So, the operator zi T is invertible if and only if the pencil λa B is also invertible. Therefore, ρ(t ) = ρ(i, T ) = ϕ(ρ(a, B)). Passing to the complement, we conclude that σ(t ) = σ(i, T ) = ϕ(σ(a, B)). Consequently σ(t ) is in the unit disk. Using [2, Theorem 2], we conclude that there exists an operator W such that T W T W = G for all G, (2.3) which is equivalent to {[(α 1)A B][(α + 1)A B] 1 } W {[(α 1)A B][(α + 1)A B] 1 } W = G [(α + 1)A B ] 1 [(α 1)A B ]W [(α 1)A B][(α + 1)A B] 1 W = G [(α 1)A B ]W [(α 1)A B] [(α + 1)A B ]W [(α + 1)A B] = [(α + 1)A B ]G[(α + 1)A B] (2A W B + 2B W A 4αA W A) = [(α + 1)A B ]G[(α + 1)A B] A W B + B W A 2αA W A = 1[(α + 2 1)A B ]G[(α + 1)A B] A W B + B W A 2αA W A = 2U,

where In fact, and for each x H, we have EXPONENTIAL STABILITY 19 U = 1 4 [(α + 1)A B ]G[(α + 1)A B]. U = 1 4 [(α + 1)A B ]G[(α + 1)A B] = U, Ux, x = 1 4 [(α + 1)A B ]G[(α + 1)A B]x, x = 1 G[(α + 1)A B]x, [(α + 1)A B]x 4 = 1 4 Gy, y k 4 y 2, y = [(α + 1)A B]x, where k is a positive constant. But, Therefore Thus x 2 = [(α + 1)A B] 1 y 2 [(α + 1)A B] 1 2 y 2. y 2 x 2 [(α+1)a B] 1 2. Ux, x 1 k 4 [(α+1)a B] 1 2 x 2 >. Consequently U, and (2.2) holds. Sufficient condition. If α +1 ρ(a, B) is a regular value for the pencil λa B, then the operator T = [(α 1)A B][(α + 1)A B] 1 is bounded and (2.2) becomes A W B + B W A 2αA W A = 1 2 [(α + 1)A B ]G[(α + 1)A B]. Therefore, G = W T W T (see (2.3)). Using again [2, Theorem 2], the spectrum σ(t ) will be inside the unit disk. We conclude that σ(a, B) = ϕ 1 (σ(t )) {λ : Reλ < α}, where λ = ϕ 1 (z) = α+ z+1 is a conformal mapping and Theorem z 1 2.4 is proved. Theorem 2.5. If (2.1) is satisfied for the pair of the positive uniform operators (W, U), then λ = α + 1 is not an eigenvalue for the pencil λa B. Proof. Suppose that λ = α + 1 is an eigenvalue. We denote by v the corresponding eigenvector. Then, [(α + 1)A B]v = or (α + 1)Av = Bv, and in the two cases the scalar product becomes Uv, v = 1 2 (A W B + B W A 2αA W A)v, v = 1 2 A W Bv, v 1 2 B W Av, v + αa W Av, v = 1 W Bv, Av 1 W Av, Bv + α W Av, Av 2 2 = W Av, Av <.

11 M. BENABDALLAH, M. HARIRI We obtain a contradiction, with the hypothesis U, since W. Consequently Theorem 2.5 is proved. Corollary 2.6. In the case of a finite-dimensional space H, the following statements are equivalent: (a) The equation (1.2) is exponentially stable; (b) σ(a, B) = σ p (A, B) {λ : Reλ < α}; (c) There exists a positive definite matrix W such that A W B + B W A 2αA W A. 3. Quasi-linear implicit equation In this section, we give some stability conditions of the quasi-linear implicit equation of the form (1.1), using the following variation of constants method (Lemma 3.1) and the Gronwall Bellman inequality (Lemma 3.2). Remember that D = {x()} denotes the initial manifold subspace of H for the stationary equation (1.2). Lemma 3.1. Suppose that (i) the restriction operator A = A D on D is invertible; 2 (ii) for any τ, the space θ(τ, x(τ)) is in the domain of A and the function S(t τ)a 1 θ(τ, x(τ)) is integrable (with respect to τ), where {S(t)} t is the semigroup of the operators for (1.2). Then the quasi-linear equation (1.1) is equivalent to the integral equation x(t) = S(t)x + Lemma 3.2 (Gronwall Bellman ). (see [1]). If g(t) c + S(t τ)a 1 θ(τ, x(τ))dτ. (3.1) g(τ)h(τ)dτ for all t, where h is a continuous positive real function and c > is an arbitrary constant, then [ ] g(t) c exp h(τ)dτ. For the quasi-linear equation (1.1), we have the next result. Theorem 3.3. Suppose that (i) the equation (1.2) is well-posed; (ii) the quasi-linear operator θ(t, x(t)), for all t, transforms D into AD such that θ(t, x(t)) dt <. Then the quasi-linear equation (1.1) is exponentially stable. 2 In particular, if (1.2) is well-posed, then A is invertible.

EXPONENTIAL STABILITY 111 Proof. Thanks to Lemma 3.1, equation (1.1) is equivalent to (3.1). According to the hypothesis (i), we have S(t)x Me αt x, and S(t τ)a 1 θ(τ, x(τ)) Me α(t τ) θ(τ, x(τ)). Considering (i) and (ii), we have A 1 θ(τ, x(τ)) D. Using (3.1), we obtain or x(t) Me αt x + M e αt x(t) M x + M e α(t τ) θ(τ, x(τ)) x(τ) dτ e ατ θ(τ, x(τ)) x(τ) dτ. Applying Lemma 3.2 with g(t) = e αt x(t), h(τ) = M θ(τ, x(τ)), and c = M x, we obtain [ ] Thus, where e αt x(t) M x exp [ M 1 = M exp M M [ M x exp M x(t) M 1 e αt x, θ(τ, x(τ)) dτ ] θ(τ, x(τ)) dτ. ] θ(τ, x(τ)) dτ <. Corollary 3.4. If the conditions (i) and (ii) of Theorem 3.3 are fulfilled and (1.2) is exponentially stable, then the quasi-linear equation (1.1) is also exponentially stable. Remark 3.5. Theorem 3.3 represents the generalization of the Dini Hukuhara theorem [1], where A I, B T, θ(t, x(t)) T (t){x(t)}, and α =. Finally we provide the following example to illustrate our main result. Example 3.6. Consider (1.1) in the finite-dimensional spaces: ( ) ( ) ( ) 1 1 1 1 A =, B =, θ(t, x(t)) e 1 1 1 t, t. 1 In our case D = {(a, b) R 2 : b = }, ( ) λ + 1 λa B = λ + 1 1 AD = {(a, b) R 2 : a = b}, ( 1 λ 1 λ + 1, (λa B) 1 = 1 λ + 1 It is clear that θ(t, x(t)) : D AD, t, and A is invertible. ).

112 M. BENABDALLAH, M. HARIRI Since σ(a, B) = σ p (A, B) = { 1}, then (1.2) is exponentially stable (see Corollary 2.6). From Corollary 3.4, we conclude that the corresponding quasilinear equation (1.1) is also exponentially stable as far as, θ(t, x(t)) dt = = <. θ(t, x(t)) dt e t dt Acknowledgement. The authors would like to express their gratitude to the anonymous referees for their comments and suggestions that improve the last version of the manuscript. References 1. R. Bellman and K.L. Cooke, Differential-Difference Equations, Vol 6, Academic Press., London, 1963. 2. M. Benabdallah, A.G. Rutkas and A.A. Soloviev, On the stability of degenerate difference systems in Banach spaces, J. Sov. Math. 57 (1991) 3435 3439. 3. P.L. Butzer and H. Berens, Semigroups of Operators and Approximations, Springer-Verlag, New York-Berlin, 1967. 4. J.L. Daleckiĭ and M.G. Krein, Stability of Solutions of Differential Equations in Banach Space, Amer. Math. Soci. Providence, RI, 1975. 5. A. Favini and A. Yagi, Degenerate Differential Equations in Banach Spaces, Marcel Dekker Inc. New York-Basel-Hong Kong, 1999. 6. A.G. Rutkas, Spectral methods for studying degenerate differential-operator equations, J. Math. Sci. 144 (27), no. 4, 4246 4263. 7. L.A. Vlasenko, Evolutionary models with implicit and degenerate differential equations, Dniepropetrovsk, no. 4885, (Russian), 26. 1 Department of Mathematics, Faculty of Math and Computer, USTOran, 31, Algeria E-mail address: mehdibufarid@yahoo.fr; mehdi.benabdallah@univ-usto.dz 2 Department of Mathematics, University of Djillali Liabes, Sidi Bel-Abbes, 22, Algeria. E-mail address: haririmohamed22@yahoo.fr; mohamed.hariri@univ-sba.dz