- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Journal of Applied Mathematics

Volume 2012 (2012), Article ID 793801, 19 pages

http://dx.doi.org/10.1155/2012/793801

## Random -Policy Geo/G/1 Queue with Startup and Closedown Times

Department of Accounting Information, National Taichung University of Science and Technology, Taichung 404, Taiwan

Received 23 July 2012; Revised 6 October 2012; Accepted 7 October 2012

Academic Editor: Han H. Choi

Copyright © 2012 Tsung-Yin Wang. 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.

#### Abstract

This study investigates a random -policy Geo/G/1 queue with startup and closedown times. is newly determined every time a new cycle begins. When random customers are accumulated, the server is immediately turned on but is temporarily unavailable to the waiting customers. It needs a startup time before starting providing service. After all customers in the system are served exhaustively, the server is shut down by a closedown time. Using the generating function and supplementary variable technique, analytic solutions of system size, lengths of state periods, and sojourn time are derived.

#### 1. Introduction

This paper deals with a random *-*policy Geo/G/1 queueing system in which the random variable , the startup time, and the closedown time obey the general distributions, respectively. The system of turning on and turning off the server depends on the number of customers in the queue. is newly determined every time a new cycle begins. When the queue length reaches a random threshold (), the server is instantly turned on but is temporarily unavailable to the waiting customers. The server needs the startup time before starting each of his service periods. Once the startup is over, the server immediately starts serving the waiting customers until the system is empty. As soon as the system becomes empty, the server also needs a closedown time to be shut down.

It is assumed that the random threshold , service time , startup time ,* and *closedown time are all general distributions with probability mass functions, means, variances and probability generating functions are as follows:

Interarrival time is a geometric distribution with parameter rate and the possible value of is on the integers 1,2,…. In startup and closedown periods, the system allows the customers to enter the system to be served. At the instant of the end of closedown, if there are customers in the system, the service is immediately started without startup time. Arriving customers form a single waiting line based on first-come, first-served discipline. All customers arriving to the system are assumed to be served exhaustively. Furthermore, various stochastic processes are independent of each other.

The controllable queueing systems possess its applications in wide fields such as manufacturing/production systems, communication networks, and computer systems. Because the -policy is analytically easier than other policies, many researchers concentrated on this type of service policy. The -policy M/G/1 queueing system has been well studied by many queueing researchers (see [1–4] and others). -policy G//M/1 queues with the finite and infinite capacity have been studied by Ke and Wang [5], and Zhang and Tian [6], respectively. Some related studies on M/G/1 queueing system with a startup time have been reported (see [7–11], etc). Bisdikian [12] applied the decomposition property to derived queue size for a random -policy M/G/1 queue in which is a random variable. He also investigated the analytic solutions of waiting time for both the FIFO and LIFO service disciplines. On the other hand, various authors analysed queueing models under several combinations of server vacations and -policy. The investigations of this type can be seen in [13–16]. Extending the combinations of server vacations and -policy to server with startup can be found in [17, 18]. Recently, Arumuganathan and Jeyakumar [19] considered a bulk queue with multiple vacations, setup times with -policy, and closedown times. After completing a service, if the queue length, , is less than “a”, then the server performs closedown work and then takes vacations. The server returning from a vacation, if is still less than “a”, then the server leaves for another vacation and so on, until is greater than “b”. Ke [20] derived the distribution of various system characteristics for two different kinds of policy M/G/1 queueing system with breakdown, startup and closedown time. Recently, Choudhury et al. [21] investigated bulk arrival queue with -policy by introducing a delay time for commencement of service after a breakdown. Each time a service is to start there should be at least customers in the system. Kuo et al. [22] dealt with the optimal operation of the -policy M/G/1 queue with server breakdowns, general startup and repair times. When the number of customers in the system reaches , turn the server on with probability and leave the server off with probability (). Such a system policy is called -policy. They developed system performance measures and provided an efficient procedure to determine the optimal threshold of () that minimized the total expected cost. Ke et al. [23] study the operating characteristics of an queueing system with*-*policy and at most vacations.

While many continuous time queueing systems with -policy have been studied, their discrete time counterparts have received little attention in the literature. Takagi [24] derived the queue size and waiting time under the -policy Geo/G/1 queue with batch arrival. Böhm and Mohanty [25] investigated -policy for the Geo/Geo/1 queue involving batch arrival and batch service, respectively. Moreno [26] extended a modified -policy issue, where the first customers of each consecutive service period are served together and the rest of customers are served singly. She gave detailed derivations of system characteristics for a discrete time Geo/G/1 queue and developed a cost function to search the optimal operating -policy at a minimum cost. Furthermore, Moreno [27] analyzed a discrete time single serve queue with a generalized -policy and setup-closedown times. In [27], the author derived the formulae for various system performance measures, such as queue and system lengths, the expected length of the vacation, setup, and busy and closedown periods, and performed a numerical investigation on the expected cost function. Recently, Wang and Ke [28] discussed a discrete-time Geo/G/1 queue, in which if the customers are accumulated to , the server operates a () policy.

However, only very few works in the literature concerned with *-*policy queueing systems with startup, and closedown time have been done. Especially, the researches of the discrete time *-*policy queueing models with startup, and closedown time are really very rare. Moreover, in the past works of the *-*policy, is a fixed number except for Bisdikian [12]. To my best knowledge, the discrete time case of *-*policy where is a random variable has never been studied. In the real situation, the server with startup and closedown times is a natural abstraction and the number may vary depending on different instances of the operation of the system. For example, Streaming technique is used for compression of the audio/video files so they can be retrieved and played by remote viewers in real time. When a user wants to play video file and activates the streaming player, the streaming player will download and store the uncertain size of data (random , is newly determined every time for a new cycle beginning, depending the network bandwidth) in the buffer in advance before playing the data. The user can watch the video file before the entire video file has been downloaded. However, when the streaming player is activated for playing the video, it needs a short time to start up. It also needs a shutdown time to be closed as all received data has been played.

The paper is structured as follows. In the next section, we formulate the system as an embedded three-dimensional Markov chain and provide the stationary joint distribution of system size and the server’s status. In Section 3, the idle,startup, busy, and closedown periods are derived. In Section 4 the explicit forms for the mean waiting time of a customer in the system conditioned on the various states are obtained. The mean waiting time of a customer in the system is also obtained and this result confirms the Little’s formula. Section 5 gives the numerical aspects to illustrate the effect of the varying parameters on the expected length in the system.

#### 2. Model Formulation and Stationary Distribution

In continuous time queues an arrival and a departure never happen simultaneously (i.e., the probability of an arrival and a departure occurring simultaneously is zero) hence the order of an arrival and a departure can be easily distinguishable. However, in discrete-time system, time is treated as a discrete variable (slot), and an arrival and a departure can only occur at boundary epochs of time slots (i.e., an arrival and a departure may occur concurrently in a slot). If we want to compute the number of customers in the system at time slot and let it has a precise meaning, the order of the arrivals and departures must be stated. Whether an arrival or a departure is recorded into the number of customers in the system at time slot , there are two different agreements: one is called the late arrival system (LAS) if a potential arrival occurs within and a potential departure occurs within ; the other is called the early arrival system (EAS) if a potential departure occurs within and a potential arrival occurs within . These concepts and other related ones can be found in Takagi [24] and Hunter [29]. The event occurring in denotes the event occurring immediately before slot boundary and the event occurring in denotes that it is occurring immediately after slot boundary. Figure 1 depicts time epochs for the LAS and the EAS. LAS has two variants: LAS with immediate access and LAS with delayed access. The difference between them is when a customer arrives late in the th slot during the server which is free, the service is started in the th slot (LAS with immediate access) or the service is started in the ()st slot (LAS with delayed access). Because the management policy of LAS with immediate access has no obvious applications to computer and communication systems, we adopt the LAS policy with delayed access in the presented model and for any real number , we denote . Figure 2 depicts time epochs under a natural extension of the LAS for the present model.

Let denote the server state at time :

And let be a supplementary random variable defined as follows:

Let the random variable indicate the number of customers in the system at . The sequence of is a Markov chain whose state spaceis as follows:

Let us define the following limiting probabilities:

The steady-state Kolmogorov equations are given by where

The left-hand sides in (2.5)–(2.8) represent the steady-state probabilities that states observed immediately after the current slot boundary changes to states observed immediately after the next slot boundary. For example, the left-hand side, , of (2.5) denotes the probability from the current state transiting to the next state where the server is idle, the threshold is and there are customers in the system. It should be noted that the next state depends only on the current state. From the current state transiting to the next state, there are three cases. (a) Given the current state that the server is idle, the threshold is , and there are customers in the system, the probability of no customer arriving from the current state to the next state is . The probability that the server is idle, the threshold is , and there are customers in the system at current state is . Hence the joint probability is . (See the first term of right-hand side for (2.5)). (b) Given the current state that the server is during closedown period, there are nocustomers in the system, and the remaining closedown time is one slot at current, the probability that no customers arrive at the next state and the threshold is is . The probability that the server is during closedown period, there are no customers in the system, and the remaining closedown time is one slot is . Therefore, the joint probability is . (See the second term of right-hand side for (2.5)). (c) Given the current state that the server is idle, there are customers in the system, and the threshold is , the probability that a customer arrives from the current state to the next state is . Hence the unconditional probability is . (See the third term of right-hand side for (2.5)). By using the similar approach, we can obtain (2.6)–(2.8) for denoting the next state that the server is startup, the next state that the server is busy, and the next state that the server is closedown, respectively.

Define the following generating functions: where and .

From (2.5), when and , we have

For and in (2.5) it yields that

From (2.11) and (2.12), we get

Hence, the generation function of an idle server is given by

Multiplying (2.6) by and summing over after multiplying (2.6) by and summing over , we obtain

Applying the same procedure to (2.7) and (2.8), respectively, we obtain

Inserting in (2.15), (2.16), and (2.17) yields

Setting in (2.20) and yields

Substituting (2.18), (2.20), and (2.24) into (2.22) and (2.23) yields

Let denote the system size. Following (2.14), (2.21), (2.25) and (2.26), the probability generating function (p.g.f.) is given by

##### 2.1. The Derivation of

Using the normalization condition, , yields where .

Substituting (2.28) into (2.27) gives where is the p.g.f. of the system size in the classical Geo/G/1 queue and

##### 2.2. Stationary Distribution of the Server State

Let us define the following: the probability that the server is idle; the probability that the server is startup; the probability that the server is turned on (working); the probability that the server is shut down; the probability that the system is empty.

Substituting (2.28) into (2.14), (2.21), (2.25), and (2.26) yields

Setting and in and using (2.26) yield

Hence, the probability that the system is empty is given by

##### 2.3. The Expected System Size

Differentiating in (2.29) and setting , we note that the numerator and denominator are both 0. Based on L’Hospital’s rule twice, the expected system size is given by where is the expected system size in the system for the classical Geo/G/1 queue.

#### 3. The Idle, Startup, Busy, and Closedown Periods

This section studies the idle, startup, busy, and closedown periods. An idle period starts at the departure instant of a customer which leaves the system empty and terminates when accumulated customers reach , where may be 1,2,… with probability . A startup period begins at the end of an idle period and terminates at the beginning of a service. A busy period starts at the beginning of a service and terminates when a service is completed and the system is empty. A closedown period begins at the end of a busy period and ends at the completion of closedown time.

Define the following: the length of an idle period; the length of a startup period; the length of a busy period; the length of a closedown period.

##### 3.1. The Idle Period

According to the definition, the p.g.f. of is given by and the mean length of idle period is

##### 3.2. The Startup Period

Next, we derive the startup period. A startup period is the amount of time to start up the server and ends when the startup expires. The p.g.f. and mean length of are given by and

##### 3.3. The Busy Period

Let be the p.g.f. of busy period of classical Geo/G/1 with late arrival delay access. From Takagi [24], we have . To derive the p.g.f. of the length of the busy period for the present model, we consider the following.

*Case 1. *The probability that there are arrivals at the end of idle period and the server is ready to serve is , . The p.g.f. for the length of the busy period extended by these arrivals is

*Case 2. *The p.g.f. for the number of arrivals that arrive during startup time is . Therefore, the p.g.f. for the length of the busy period extended by these arrivals is .

*Case 3. *The p.g.f. for the number of arrivals that arrive during a closedown period with arrivals is given by and the number of closedowns in a cycle obeys the geometric distribution with parameter . Therefore, we can obtain the p.g.f. for the number of customers that arrive during the closedown period with arrivals:

Consequently, the p.g.f. for the length of the busy period extended by these arrivals is given by

Because the arrivals in idle period, startup period, busy period and closedown period are independent. Hence, the p.g.f. for the length of the busy period is given by and the mean length is

##### 3.4. The Closedown Period

The probabilities of no customer and at least one customer in the system at the completion of a closedown are and , respectively. This process is Bernoulli process. Hence, the distribution of the number of shutdowns in a cycle is a geometric distribution with parameter . Besides, the p.g.f. for the length of a shutdown time is . Consequently, the p.g.f. of is given by and the mean length of is given by

From (3.2), (3.3), (3.8), and (3.10), the mean length of the service cycle, , is given by

#### 4. Sojourn Time in the System

When a customer enters a service system, one of his most concerning issues is his sojourn time in the system. The sojourn time in the system of an arrival depends on the server status at the arrival epoch. Let us define the following p.g.f.s: the p.g.f. of the sojourn time in the system of a customer (that arrives at an arbitrary slot) conditioning that the server state is idle; the p.g.f. of the sojourn time in the system of a customer (that arrives at an arbitrary slot) conditioning that the server state is startup; the p.g.f. of the sojourn time in the system of a customer (that arrives at an arbitrary slot) conditioning that the server state is busy; the p.g.f. of the sojourn time in the system of a customer (that arrives at an arbitrary slot) conditioning that the server state is closedown; the unconditional p.g.f. of the sojourn time in the system of a customer (that arrives at an arbitrary slot).

There are four situations of the server at the epoch of a tagged arrival.

*Case 1. *Suppose the server is idle will go into startup state if customers are accumulated in the queue. A tagged customer arrives and finds customers () in the system at that moment the server is idle. The sojourn time in the system of this tagged customer consists of (i) the time of waiting arrivals for the server going into startup state; (ii) the startup time; (iii) the total service time of the preceding customers; (iv) the service time of the tagged customer:

*Case 2. *A tagged customer arrives during startup period and finds () customers in the system. The tagged customer’s sojourn time in the system is the remaining startup time plus the service time of the (including himself) customers:

*Case 3. *A tagged customer arrives while the server is busy and finds customers in the system. In this case, his sojourn time in the system consists of (i) the remaining service time of the customer being served; (ii) the service time of the customers (including himself). Hence, we have

*Case 4. *A tagged customer arrives while the server is during the closedown period and finds () customers in the system. In this case, the sojourn time in the system of the customer consists of the remaining closedown time, and the service time of the customers:

By unconditioning (4.1)–(4.4), the sojourn time in the system of an arbitrary customer is given by

Finally, we obtain where is the p.g.f. of the customer’s waiting time in the queue.

From (4.6), the expected sojourn time in the system is given by which is in accordance with and this result confirms Little’s formula.

#### 5. Numerical Examples

In this section, we present some numerical examples to compare the mean length in the system of the presented model with respect to no startup, no closedown, and both no startup and no closedown. From (2.37), we can easily obtain the mean lengths in the system with no startup, no closedown, and both no startup and no closedown as follows:(i)no startup: (ii)no closedown: (iii)no startup and no closedown:

We assume that the startup time and the closedown time are geometric distributions with parameter rates 0.4 and 0.5, respectively, and consider the following four cases. For convenience, if a random has a geometric distribution with parameter , we denote . Similarly, if a random has a discrete uniform distribution on the integers from to , we denote .

*Case 1. *Let , , and the values of the parameter vary from 0.52 to 096 by the increment 0.04.

*Case 2. *Let , , and the values of the parameter vary from 0.52 to 096 by the increment 0.04.

*Case 3. *Let , , and vary the values of from 0.3 to 0.5 by the increment 0.02.

*Case 4. *Let , , and vary the values of from 0.3 to 0.5 by the increment 0.02.

Cases 1–4 are depicted in Figures 3–6, respectively. Figures 3 and 4 display that the mean length in the system decreases in service rate and Figures 5–6 display that the mean length in the system increases in arrival rate. In Figure 3 (also see Figure 5), the mean length in the system with startup and closedown is always greater than the mean length in the system with no startup and no closedown. But in Figure 4 (also see Figure 6), the mean length in the system with no startup and no closedown is always greater than the mean length in the system with no startup. This is because we choose the different threshold distributions. Because the server needs a startup time before starting to serve customers, the startup time will increase number of customers in the system. Moreover, the system allows the customers to enter the system in closedown period and the server immediately serves the customers in the system without startup time. Comparing it with the system with no startup time and no closedown, it is efficient in decreasing the customers in the system. Hence, we observe the phenomenon that the mean length in the system without closedown is greater than the others and the mean length in the system without startup is smaller than the others.

#### 6. Conclusions

In this work, we investigated a random -policy Geo/G/1 queue with startup and closedown times. The analytical results of mean lengths of the system size, idle period, startup period, busy period, and closedown period were derived. The most careful thing of a customer entering a service system, his sojourn time in the system, is also derived. We obtained the steady-state distributions of the waiting time of a customer in the system conditioned on the various states. The analytic solution of unconditional mean waiting time was also obtained and the result showed that Little’s formula still holds.

#### Acknowledgments

The author thanks the anonymous reviewers whose comments and suggestions were very helpful in improving this paper. This research was supported in part by the National Science Council, Taiwan, under Contract no. NSC 101–2221-E-025-005-.

#### References

- M. Yadin and P. Naor, “Queueing systems with a removable service station,”
*Operational Research Quarterly*, vol. 14, pp. 393–405, 1963. View at Google Scholar - D. P. Heyman, “Optimal operating policies for M/G/1 queueing system,”
*Operations Research*, vol. 16, pp. 362–382, 1968. View at Google Scholar - C. E. Bell, “Characterization and computation of optimal policies for operating an M/G/1 queueing system with removable server,”
*Operations Research*, vol. 19, pp. 208–218, 1971. View at Google Scholar - K. H. Wang and J. C. Ke, “A recursive method to the optimal control of an M/G/1 queueing system with finite capacity and infinite capacity,”
*Applied Mathematical Modelling*, vol. 24, no. 12, pp. 899–914, 2000. View at Publisher · View at Google Scholar · View at Scopus - J.-C. Ke and K.-H. Wang, “A recursive method for the $N$ policy $G/M/1$ queueing system with finite capacity,”
*European Journal of Operational Research*, vol. 142, no. 3, pp. 577–594, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - Z. G. Zhang and N. Tian, “The $N$ threshold policy for the GI/M/1 queue,”
*Operations Research Letters*, vol. 32, no. 1, pp. 77–84, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J. Medhi and J. G. C. Templeton, “A poisson input queue under
*N*-policy and with a general start up time,”*Computers and Operations Research*, vol. 19, no. 1, pp. 35–41, 1992. View at Google Scholar · View at Scopus - H. Takagi, “A M/G/1/K queues with
*N*-policy and setup times,”*Queueing Systems. Theory and Applications*, vol. 14, no. 1-2, pp. 79–98, 1993. View at Publisher · View at Google Scholar - H. W. Lee and J. O. Park, “Optimal strategy in
*N*-policy production system with early set-up,”*Journal of the Operational Research Society*, vol. 48, no. 3, pp. 306–313, 1997. View at Google Scholar · View at Scopus - S. Hur and S. J. Paik, “The effect of different arrival rates on the
*N*-policy of M/G/1 with server setup,”*Applied Mathematical Modelling*, vol. 23, no. 4, pp. 289–299, 1999. View at Google Scholar · View at Scopus - K. H. Wang, T. Y. Wang, and W. L. Pearn, “Optimal control of the N policy M/G/1 queueing system with server breakdowns and general startup times,”
*Applied Mathematical Modelling*, vol. 31, no. 10, pp. 2199–2212, 2007. View at Publisher · View at Google Scholar · View at Scopus - C. C. Bisdikian, “Random
*N*-policy,” in*Proceedings of IEEE 13th Networking for Global Communications (INFOCOM '94)*, pp. 1178–1182, June 1994. View at Scopus - O. Kella, “The threshold policy in the M/G/1 queue with server vacations,”
*Naval Research Logistics*, vol. 36, no. 1, pp. 111–123, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - H. W. Lee, S. S. Lee, J. O. Park, and K. C. Chae, “Analysis of the ${M}^{X}/G/1$ queue with $N$-policy and multiple vacations,”
*Journal of Applied Probability*, vol. 31, no. 2, pp. 476–496, 1994. View at Publisher · View at Google Scholar - K. C. Chae and H. W. Lee, “${M}^{X}/G/1$ vacation models with
*N*-policy: heuristic interpretation of the mean waiting time,”*Journal of the Operational Research Society*, vol. 46, no. 2, pp. 258–264, 1995. View at Google Scholar · View at Scopus - S. S. Lee, H. W. Lee, S. H. Yoon, and K. C. Chae, “Batch arrival queue with
*N*-policy and single vacation,”*Computers and Operations Research*, vol. 22, no. 2, pp. 173–189, 1995. View at Google Scholar · View at Scopus - G. V. Krishna Reddy, R. Nadarajan, and R. Arumuganathan, “Analysis of a bulk queue with $N$-policy multiple vacations and setup times,”
*Computers & Operations Research*, vol. 25, no. 11, pp. 957–967, 1998. View at Publisher · View at Google Scholar - J. C. Ke, “The optimal control of an M/G/1 queueing system with server startup and two vacation types,”
*Applied Mathematical Modelling*, vol. 27, no. 6, pp. 437–450, 2003. View at Publisher · View at Google Scholar · View at Scopus - R. Arumuganathan and S. Jeyakumar, “Steady state analysis of a bulk queue with multiple vacations, setup times with
*N*-policy and closedown times,”*Applied Mathematical Modelling*, vol. 29, no. 10, pp. 972–986, 2005. View at Publisher · View at Google Scholar · View at Scopus - J. C. Ke, “On M/G/1 system under
*NT*policies with breakdowns, startup and closedown,”*Applied Mathematical Modelling*, vol. 30, no. 1, pp. 49–66, 2006. View at Publisher · View at Google Scholar · View at Scopus - G. Choudhury, J.-C. Ke, and L. Tadj, “The $N$-policy for an unreliable server with delaying repair and two phases of service,”
*Journal of Computational and Applied Mathematics*, vol. 231, no. 1, pp. 349–364, 2009. View at Publisher · View at Google Scholar - C.-C. Kuo, K.-H. Wang, and S.-L. Lee, “Optimal control of the $\langle p,N\rangle $-policy M/G/1queue with server breakdowns and general startup Time,”
*International Journal of Information and Management Sciences*, vol. 20, no. 4, pp. 565–577, 2009. View at Google Scholar - J.-C. Ke, H.-I. Huang, and Y.-K. Chu, “Batch arrival queue with $N$-policy and at most $J$ vacations,”
*Applied Mathematical Modelling*, vol. 34, no. 2, pp. 451–466, 2010. View at Publisher · View at Google Scholar - H. Takagi,
*Queueing Analysis: A Foundation of Performance Evaluation: Discrete-time Systems*, vol. 3, North-Holland, Amsterdam, The Netherland, 1993. - W. Böhm and S. G. Mohanty, “On discrete time Markovian $N$-policy queues involving batches,”
*Sankhyā A*, vol. 56, no. 1, pp. 144–163, 1994. View at Google Scholar - P. Moreno, “A discrete-time single-server queue with a modified $N$-policy,”
*International Journal of Systems Science*, vol. 38, no. 6, pp. 483–492, 2007. View at Publisher · View at Google Scholar - P. Moreno, “Analysis of a Geo/G/1 queueing system with a generalized $N$-policy and setup-closedown times,”
*Quality Technology & Quantitative Management*, vol. 5, no. 2, pp. 111–128, 2008. View at Google Scholar - T. Y. Wang and J. C. Ke, “The randomized threshold for the discrete-time Geo/G/1 queue,”
*Applied Mathematical Modelling*, vol. 33, no. 7, pp. 3178–3185, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. J. Hunter,
*Mathematical Techniques of Applied Probability, Discrete Time Models: Techniques and Applications*, vol. 2, Academic Press, New York, NY, USA, 1983.