#### Abstract

We investigate the dynamics of a delayed stochastic mathematical model to understand the evolution of the alcohol consumption in Spain. Sufficient condition for stability in probability of the equilibrium point of the dynamic model with aftereffect and stochastic perturbations is obtained via Kolmanovskii and Shaikhet general method of Lyapunov functionals construction. We conclude that alcohol consumption in Spain will be constant (with stability) in time with around 36.47% of nonconsumers, 62.94% of nonrisk consumers, and 0.59% of risk consumers. This approach allows us to emphasize the possibilities of the dynamical models in order to study human behaviour.

#### 1. Introduction

In this paper, we propose a mathematical framework to model social epidemics. To be precise, we propose delayed and stochastic consideration on mathematical models to analyze human behaviors related to addictions.

Hereditary systems or systems with delays are very popular in researches (see, e.g., [1–4] and the references therein). In this paper, a nonlinear dynamic alcohol consumption model [5] is generalized via adding distributed delays. Sufficient conditions for existence of the positive equilibrium points of this system are obtained. It is also supposed that this nonlinear system is exposed to additive stochastic perturbations of white noise that are directly proportional to the deviation of the system current state from the equilibrium point. Such type of stochastic perturbations first was proposed in [6, 7] and successfully used later in [8–11]. One of the important point of this assumption is that the equilibrium point is the solution of the stochastic system too. In this case, the influence of the stochastic perturbations on the considered system is small enough in the neighbourhood of the equilibrium point and big enough if the system state is far enough from the equilibrium point.

The considered nonlinear system is linearized in the neighborhood of the positive point of equilibrium, and sufficient condition for asymptotic mean square stability of the zero solution of the constructed linear system is obtained via Kolmanovskii and Shaikhet general method of Lyapunov functionals construction (GMLFC) that is used for stability investigation of stochastic functional-differential and difference equations [12–22]. Since the order of nonlinearity more than 1, this condition is also sufficient one [23, 24] for stability in probability of the initial nonlinear system by stochastic perturbations.

This way of stability investigation was successfully used for investigation of different mathematical models of systems with delays: SIR epidemic model [6], predator-prey model [7, 11], Nicholson blowflies model [9], inverted pendulum [25–27].

The present paper is organized as follows. Section 2 presents the delayed model including the study of its equilibrium points. Stochastics characteristics of the delayed model are shown in Section 3. Section 4 is related to the study of stability of the equilibrium point of the delayed stochastics model proposed. The main conclusions derived from the study are presented in the last section.

#### 2. Delayed Mathematical Model

Taking into account the proposal presented by Rosenquist et al. in [28], we consider alcohol consumption habit as susceptible to be transmitted by peer pressure or social contact. This fact lead us to propose an epidemiological type mathematical model to study this social epidemic.

Let be nonconsumers, individuals that have never consumed alcohol or they infrequently have alcohol consumption, and nonrisk consumers, individuals with regular low consumption; to be precise, men who consume less than 50 cc of alcohol every day and women who consume less than 30 cc of alcohol every day. Let be risk consumers, individuals with regular high consumption, that is, men who consume more than 50 cc of alcohol every day and women who consume more than 30 cc of alcohol every day.

Considering homogeneous mixing [29], that is, each individual can contact with any other individual (peer pressure), dynamic alcohol consumption model is given by the following nonlinear system of ordinary differential equations with distributed delay:
where(i)*μ*: birth rate in Spain;(ii)*γ*: rate at which a risk consumer becomes a nonconsumer; (iii): death rate in Spain;(iv)*β*: transmission rate due to social pressure to increase alcohol consumption (family, friends, marketing, TV, etc.);(v)*d*: augmented death rate due to alcohol consumption. Accidents at work, traffic accidents and diseases derived by alcohol consumption are considered. The information available to calculate the augmented death rate due to alcohol consumption is aggregate. This fact does not allow us to consider the difference between the augmented death rate in nonrisk and risk consumers;(vi)*α*: rate at which a nonrisk consumer moves to the risk consumption subpopulation.

It is supposed that the parameters and the function are nonnegative and the following condition holds: In the particular case , where , is Dirac's function, system (2.1) is a system with discrete delay . The case of a system without delay () is considered in [5].

It is assumed that when a nonconsumer individual is *infected* by alcohol consumers, there is a time during which the alcohol consumption habit develops in nonconsumer and it is only after that time that nonconsumer individual becomes nonrisk consumer. It is also assumed that at any time alcohol consumers population is simply proportional to alcohol consumers population at time . Let be the fraction of nonconsumer population that takes time to become alcohol consumers (nonrisk consumer).

##### 2.1. Normalization of the Delayed Model

Put From (2.1) and (2.6), it follows that Adding the first three equations (2.1) by virtue of (2.7), we obtain

It is easy to see that and similarly Thus, putting from (2.1), (2.5), (2.7), (2.9), (2.10), we have

Via (2.7), the last equation can be rejected, and, as a result, we obtain the system of two differential equations

##### 2.2. Existence of the Equilibrium Point

Via (2.5), (2.7), (2.13), the point of equilibrium is defined by the system of the algebraic equations:

Lemma 2.1. *If , then system (2.14) has the unique positive solution if and only if
**
If , then system (2.14) does not have positive solutions.*

*Proof. *From two first equations of system (2.14), we have
or

Thus, via (2.14), the equilibrium point is defined by the system of the algebraic equations (2.19) and
It is easy to check that, by condition (or ), the existence of the solution of (2.19) in the interval is equivalent to the condition that is equivalent to (2.17). If , then and . So, (2.19) cannot have positive roots. The proof is completed.

*Example 2.2. *Consider the values of the parameters from [5]:

Then, condition (2.17) for these values of the parameters holds, and the solution of systems (2.19) and (2.20) is
or in the percents , , .

#### 3. Stochastic Perturbations, Centralization, and Linearization

Let us suppose that system (2.13) is exposed to stochastic perturbations type of white noise , which are directly proportional to the deviation of system (2.13) state from the point , that is, Here, are the mutually independent standard Wiener processes, the stochastic differential equations of system (3.1) are understanding in Ito sense [30].

To centralize system (3.1) in the equilibrium point, put now Then, from (3.1), it follows that or Via (2.14) from (2.16), it follows that where Note that, for , the numbers , , are positive.

Rejecting the nonlinear terms in (3.5), we obtain the linear part of system (3.5)

#### 4. Stability of the Equilibrium Point

Note that nonlinear system (3.5) has the order of nonlinearity more than 1. Thus, as it follows from [23, 24], sufficient conditions for asymptotic mean square stability of the zero solution of linear part (3.7) of nonlinear system (3.5) at the same time are sufficient conditions for stability in probability of the zero solution of nonlinear system (3.5) and therefore are sufficient conditions for stability in probability of the solution of system (3.1).

To get sufficient conditions for asymptotic mean square stability of the zero solution of system (3.7), rewrite this system in the form where

*Definition 4.1. *The trivial solution of (2.20) is called as follows.(i)*Mean square stable* if for any there exists a such that for any provided that the initial function , , satisfies the condition .(ii)*Asymptotically mean square stable* if it is mean square stable and satisfies the condition for each initial function .(iii)*Stable in probability* if for any and there exists such that the solution of (2.20) satisfies the condition for any initial function such that .

Following the GMLFC [12–14] for stability investigation of (4.1), consider the auxiliary equation without memory and the matrix equation where , the matrix is defined in (3.6), (4.2).

If matrix equation (4.4) has a positive definite solution , then the function is a Lyapunov function for (4.3).

Note that matrix equation (4.4) can be represented as the system of the equations with the solution where

Lemma 4.2. *If and , that is,
**
then the matrix with entries (4.7) is a positive definite one.*

*Proof. *Note that via (4.9) , . Thus, from (4.7), it follows that for arbitrary
Therefore, via (4.9), , . Let us show that . Really, from
it follows that by . Substituting into this inequality from (4.7), we have or
Since , , the obtained inequality holds for arbitrary . Thus, for arbitrary , the matrix with entries (4.7) is a positive definite one. The proof is completed.

Theorem 4.3. *If conditions (4.9) hold and for some the matrix entries (4.7) satisfy the condition
**
then the solution of system (3.1) is stable in probability.*

*Proof. *Note that the order of nonlinearity of system (3.1) is more than one. Therefore, via [23, 24] to get for this system conditions of stability in probability, it is enough to get conditions for asymptotic mean square stability of the zero solution of linear part (3.7) of this system. Following the GMLFC [13–15], we will construct a Lyapunov functional for system (3.7) in the form , where , , is a positive definite solution of system (4.6) with entries (4.7) and will be chosen below.

Let be the infinitesimal operator [30] of system (3.7). Then, via (3.7), (4.6),
Via (2.11), (2.5), we have and for some . Using these inequalities, we obtain
where

Putting
via (2.5), (2.11), we have . Therefore, via (4.15), (4.16) for the functional , we have
Thus, if
then the zero solution of system (3.7) is [3] asymptotically mean square stable. From (4.19), it follows that
Thus, if for some condition (4.13) holds, then there exists such that conditions (4.20) (or (4.19)) hold too, and therefore the zero solution of system (3.7) is asymptotically mean square stable. From here and [23, 24], it follows that the zero solution of system (3.5) and therefore the equilibrium point of system (3.1) is stable in probability. The proof is completed.

*Example 4.4. *Consider system (3.1) with the value of the parameters , , , , , and the equilibrium point given in (2.21) and (2.22) and the levels of noises , . We consider this value for and as an example.

From (3.6), it follows that the values of system (3.7) parameters are , , , and conditions (4.9), hold: , . Put . Then, via (4.7), , , and condition (4.13) holds: +. Thus, system (3.1) solution is stable in probability.

*Example 4.5. *Consider system (3.1) with the previous values of all parameters except for the levels of noises that are , . These values of and are selected taking into account sample errors of the monitoring of the epidemic (alcohol consumption) in Spain in the period 199720100332007. More details about these sample errors are shown in [31].

Then, conditions (4.9) hold: , . Put . Then, via (4.7), , , , and condition (4.13) holds: +. Thus, system (3.1) solution is stable in probability.

Let us get now two corollaries from Theorem 4.3 which simplify verification of stability condition (4.13). Put

*Remark 4.6. *From (4.7) and (4.21), it follows that and *. *Thus, condition (4.13) is equivalent to the condition *. *

Corollary 4.7. *If conditions (4.9) hold and , then the solution of system (3.1) is stable in probability.*

*Proof. *Note that from follows that . Putting , we obtain , that is, condition (4.13) holds. The proof is completed.

Corollary 4.8. *If conditions (4.9) hold and , then the solution of system (3.1) is stable in probability.*

*Proof. *Let us suppose that . Then, the minimum of the function is reached by
Substituting into the function , we obtain that the condition is equivalent to the condition . The proof is completed.

*Example 4.9. *Consider system (3.1) with the values of the parameters from Example 4.4. Calculating , , , we obtain, , . From Corollary 4.8, it follows that the solution of system (3.1) is stable in probability.

*Example 4.10. *Consider system (3.1) with the values of the parameters from Example 4.5. Calculating , , , we obtain , . From both Corollary 4.7 and Corollary 4.8 it follows that the solution of system (3.1), is stable in probability.

#### 5. Conclusions

In this work, a modelling approach based on delayed and stochastic differential equations is proposed to understand social behaviours and their evolutionary trends. Taking into account this approach, we can know how social habits can evolve in the future. Considering the study proposed related to the alcohol consumption habit, we note that around 63.5% of Spanish population is (and will be) alcohol consumers (risk or nonrisk consumers). The existence of a equilibrium point (stable in probability) in (*a*: nonconsumers), (*m*: nonrisk consumers), (*r*: riskconsumers), see (2.22), allows us to confirm it.

This mathematical approach is an example as to how delayed and stochastic models can also be an useful tool to model human behaviour. We consider that this approach can be an interesting framework for public health authorities and policy makers.

Note also that the system (2.1) with delay is considered here first. The conditions for the existence of the equilibrium point of the considered system by stochastic perturbations are new result too. The main result of the paper is the conditions for stability in probability of the equilibrium point by stochastic perturbations. The here proposed research method can be used for stability investigation of other important models.