 Open Access
 Authors : Toky Basilide Ravaliminoarimalalason , Falimanana Randimbindrainibe
 Paper ID : IJERTV10IS100072
 Volume & Issue : Volume 10, Issue 10 (October 2021)
 Published (First Online): 23102021
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Occupied Resources in A System With Resource Sharing Between Servicebased Customers
Toky Basilide Ravaliminoarimalalason
Ecole Doctorale en Sciences et Techniques de lIngÃ©nierie et de lInnovation,
University of Antananarivo, Antananarivo, Madagascar
Falimanana Randimbindrainibe
Ecole Doctorale en Sciences et Techniques de lIngÃ©nierie et de lInnovation,
University of Antananarivo Antananarivo, Madagascar
Abstract This paper presents analytical and generalized expressions of the amounts of occupied resources in a system providing several services to its customers. This system can be assimilated to a multiservice queue. We considered that the required amount of resources for each service are deterministic, then random variables amount of resources is studied in second part. Analysis carried out on the resource occupations have made possible the dimensioning of the total amount of necessary resources that must be deployed on the queue server. They are to serve the customers up to a certain quality of service in terms of the probability of congestion. Reported to a singleservice case, we found from the analytical expressions established in this work the probability of congestion in the FryMolina traffic model. The obtained results are compared to digital simulations of case studies similar to a real system.
Keywords Queueing; Dimensioning; Multiservice; Resource occupation; Resource sharing

INTRODUCTION
a cFIFO multiservice queuing system while allowing variable rate based on the amounts of unused resources in the server.
In all these literatures, whether it is in the sharing of resources available to all the customers present in the system, or either with variable rate, or elastic traffic, the problem arises on the total amount of resources that the multiservice server queue has for its customers. This article proposes a generalization of a multiservice queue, with fixed or random amount of resources required for each service class, in order to determine the capacity required in the server to provide a certain quality of service. Our study is based on the analytical determinations of resource occupations in such a generalized system. It is also an extension of the work done in [14] and [15].

MODEL DESCRIPTION
Let be a system sharing its resources through customers requesting services. Each customer can request a service i
The dimensioning of the amount of resources necessary for
requiring a fixed resource Ri
from the system for a random
a system guaranteeing a certain quality of service began in the era of Erlang [1]. He established a famous Erlang B formula to determine the number required channels for a telephone network [1][2]. This formula only considers a single service (voice service) provided by a system (telephone network) to users who request it. Its formula is derived from the study of a finite capacity M/M/n/n queue. The results of his studies have been developed and extended by many researchers like Palm [3], Vaulot & Chaveau [4], Joys [5], Iversen [6] and many others.
average duration i . The arrival rate of customers requesting service i from the system is described by a Poisson process of intensity i .
1 N i1 i
1 N i1 i
By denoting by N the total number of services that this system can provide to his customers, assuming that the arrivals of customers of each service are independent, we can say that the arrival to this system forms a Poisson process whose intensity is expressed by … N . The
Then Iversen [7] studied the case of a multiservice system
rate pi of customers requesting service i is equal to
by evaluating the congestion rate of such a system with access
pi i /
for each i. These two expressions derive from the
control. Elastic traffic cases have been developed on multi service queueing systems by Hanczewski et al. [8][9]. In [8], they proposed a multiservice queue model with SDFIFO (State Dependent FIFO) discipline. This system allocates its resources to each class of service in a balanced fairness way. They further proposed a generalization of multiservice queues used in elastic traffic [10].
superposition property of independent Poisson processes [16]. Let Si be the distribution of the service i duration requested by a customer. It is a strictly positive real random variable with mean i . The duration of service requested by any customer
will be denoted by S. From the total probability formula, we can express that the cumulative distribution function FS (s) of
Other studies on a fairness resource sharing have been put forward by Bonald et al. [11] where all the resources are
S is equal to
FS (s)
N
i1
pi FSi

, where
FSi
(s)
is the
S
S
available to all service classes and to all the users who present
cumulative distribution function of
Si . We will denote by
themselves in front of the server. An extension of Haddad and Mazumdar [12] assesses a certain congestion in such a system.
Fc (s)
its complementary cumulative distribution function.
Hanczewski et al. [13] have proposed another type of discipline called cFIFO (continuous FIFO) for the multi service queue. They approximated the convolutional model of
The system is then equivalent to a queue, with a single server with a total of R resources. Customers arrive at this queue according to a Poisson process of intensity . The
service duration of any customer is a random variable S with
customer arriving uniformly between the interval (0, t)
and
an average . These customers will use an amount Ri of
still being served at time t is equal to p 1 t F (s)ds .
c
c
resources from the queue server depending on the type i of service requested. If customers are being served by the system, using resources r, the sum of all the amounts of resources used by these customers, another arriving customer wanting to use
an amount of resources less than or equal to R r will be
t t 0 S
Conditionally on At m , the number N(t) of customers still served at time t is therefore equal to n with a probability:
m pn 1 p mn
served immediately. If the amount of resources requested by the arriving customer is greater than R r , then it will be queued waiting for a served customer to complete and release resources.
If the initial amount R of resources available in the server is infinite, there is a M/GI/ queue (in terms of customers). Customers who arrive will be immediately served.
t t
n
n
This is due to the fact that n customers are still being served (probability pt ) at instant t, and therefore m n customers have finished their services before instant t (probability 1 pt ).
Then, by applying the formula of the total probabilities for all the possible values of m, we obtain:
n (t) P(N (t) n) P( At m).P(N (t) n  At m)
mn
t m
.et . m . pn (1 p )mn
mn m!
n t t
t m
.et . m! . pn (1 p )mn
mn m!
pn
n!(m n)! t t
t m pn .t n t mn
t et
.(1 p )mn t et
.(1 p )mn
n! (m n)! t n! (m n)! t
mn mn
tp n
t k
tp n
t tp k
Fig. 1. System description.


ANALYTICAL EXPRESSIONS
t
n!
tp n
et
k 0
.(1 p )k
t
t
k !
tp n
t et t
n! k 0 k !

Sharing of infinite resources
We will express the total amount of resources used by the
t etet tpt
n!
t etpt n!
customers described in paragraph 2. Note by n
the
n
S
S
0 S
c
c
probability that the system serve n customers.
Theorem
0 .exp
n!
t
F (s)ds
n follows a Poisson distribution with parameter
Equilibrium is obtained when t tends towards infinity, i.e.:
Fc (s).ds , i.e. :
n lim n (t)
0 S
n
.
Fc (s)ds
Fc (s)ds
n
.exp .
Fc (s)ds
(1)
t
n
0 S
.exp
Fc (s)ds
Proof
n n! 0 S 0 S
n! 0 S
Let N(t) be the number of customers being served at instant
t, assuming that N (0) 0 , and n (t) P(N(t) n) the probability that the system serves n customers at instant t.
By replacing by the expressions of and FS within our system providing N services to its customers, we therefore have:
On the one hand, we denote by ( A )
the process of
N n
customer arrivals to the queue, and (T )
t t0
their arrival times,
i
N n
Si
Si
n n0
i.e. the instants of the process. Conditionally on the counting
i1
n n!
. 0
Fc (s)ds
measure At m of the process, the random variables T1,…,Tm
i1
(2)
are arranged in ascending order, and uniformly over the
N
.exp .
N
Fc (s)ds
interval (0, t) . On the other hand, the probability that a
i1
i 0
Si
i1
customer will have a service duration greater than s is Conditionally on the number k of customers that the server is
S S
S S
P(S s) 1 F (s) Fc (s) . Then, the probability that a
serving, we will note rk
the amount of resources used by these
customers. Let
P(R r)
be the cumulative distribution
function of the probability that a total amount r of resources of this system is used by the customers it serves.
Thus, we have:

Sharing of finite resources
For the case of finite resources, it is obvious to think of a M/GI/C queue where C denotes a certain finite capacity of the system. But the analytical resolution of this queue is still an
P(R r) k .P(rk r)
k 0
(3)
open problem, and the solutions proposed so far are only
where P(rk r) denotes the cumulative distribution function of the probability that the amount r of resources occupied by the k customers is rk .
Now we are assuming that the system serves k customers. The cumulative distribution function of the probability that the total amount of occupied resources is r is P(rk r) .
A customer requesting a service i ( i 1,…, N ) uses a
approximations [17][18][19]. In addition, in our case, we have a variable capacity depending on the customers served by the system and their requirements in terms of amount of resources, so this will further worsen the complexity of our problem.
In our research, we focus on dimensioning the amount of required resources to satisfy customers in their uses of the services provided by the system.
Being given that the system has a capacity C in terms of
resource
Ri . There are Ntuples of coefficients
(k1,…, kN )
total amount of available resources it can allocate to customers.
verifying
N k .R r
such that
k … k k . k
If the requirement of the next customer is greater than the
i1 i i
1 N i
amount of available resources in the system, we can say that
indicates the number of customers using the service i among k
the customers served by the system.
k
there is a blocking phenomenon. Two scenarios can arise: either the customer in question is rejected, or it is put on hold. In both these cases, the capacity in terms of amount of
P(r r) . pk1 … pkN
1
N
1
N
(4)
resources is insufficient.
k
k1 ,…,kN
k1 ,…, kN
On the one hand, if we are going to deploy a lot of resources.
where
k1R1 …kN RN r
pi denotes the probability that a randomly chosen
Some of them may not be used. It is therefore a waste of resources. In some cases, the availability of resources is
customer, from among k served customers, uses the service i,
i.e. pi i / .
chargeable and it will therefore be a loss in terms of investment. On the other hand, the deployment of few resources can lose some customers as we described in the
and k
k ! .
previous blocking phenomenon. It is therefore necessary to
1 N 1 2 N
1 N 1 2 N
k ,…, k k !k !…k !
The expression of the cumulative distribution function of the occupied resources in the system is therefore given by:
determine the right total amount of resources to put in place. For this, a certain tolerance threshold must be adopted such that the probability of rejection or of placing customers on hold is less than this tolerance.
Let's assume that we have deployed an infinite amount of
P(R r) .
k pk1 …pkN
(5)
resources. For this, the probability of using an amount of at
k
k ,…, k
. 1 N
most r resources is given by equation (5) recalled below:
k 0
k1 ,…,kN 1 N
i
i
k1R1…kN RN r
k
where:
P(R r) k . k
. pki
N n
k 0
ki
ki Ri r
i i
i
i1
N n
c
c
. F (s)ds
0
0
In other words, the probability of using an amount of resources greater than C is given by:
n n!
S
i1 i
k
N
.exp .
N
Fc (s)ds
B P(R C) .
. pki
i S
k k i
i1
0 i1 i
k 0
ki
ki Ri C
i i
In the following, we will simplify the writing of the equation (5) by:
We can say that this is the probability of blocking the system of capacity C because since from this amount of
k
resources C, it can no longer serve more customers. It will be
i
i
P(R r) k . k
. pki
called probability of congestion.
k 0
ki
ki Ri r
i i
For dimensioning with a tolerance , the objective is therefore to determine a capacity C such that:
By adopting the notations ki k1 ,…, kN , Ri R1 ,…, RN and
k
k R k R … k R .
B P(R C) . . pki
(6)
i i 1 1 N N
k k i
k 0
ki
ki Ri C
i i

What happens in case of single service system and unit resource for each customer?
Thus, equation (6) becomes:
k
In this part, we will see the case of a system offering a single service in which each served customer uses only a single
P(R C) k . k . p
ki
i
i
amount of rsource. For the case of exponential service duration with mean 1/ , we will see that from equations (5)
k 0
ki
ki Ri C
i i
k
and (6), we end up with the BCH (BusyCalls Held) traffic model of FryMolina [20][21][22].
The number of services provided by the system is then
N 1 . The customer arrival rate follows a Poisson process of
intensity 1 . Since the service duration is exponential, the
exp
k !
k !
k C
k
By replacing exp with , we find the congestion
k 0 k !
k k
cumulative distribution function of the random variable S of service duration is given by FS (s) 1 exp s . The rate of
rate /
k ! k !
k ! k !
k C k 0
posed by FryMolina in their traffic model..
customers requesting service 1 is therefore
p1 1 .
We note that this model, which is an alternative to the Erlang
B formula [1] [2], is used in certain dimensioning of telephony in
The cumulative distribution function of the total amount of occupied resources, from equation (5), is given by:
North America. It is observed that this model reflects much more the reality of observed traffics compared to the Erlang model [23] [24].
P(R r) .
k pk1 … pkN

Generalization for random amounts of resources
k
k ,…, k
. 1 N
requested by customers
k 0 k1 ,…,kN 1 N
k1R1 …kN RN r
This time, we consider that the amounts of resources requested by service i customers are no longer deterministic as
k
r
k1
in the previous paragraphs. Lets assume that the service i
k . k . p1
k .1k r (k) k
customers request a random amount of resource following a
k 0
k1 1
k1r
k 0
k 0
probability distribution Li with mean i . We will denote by
where:
Gi (r) the cumulative distribution function of the amount of
N k
resources requested by a customer using service i.
We still assume an infinite amount of resources available
i
i
i1 .
N Fc (s)ds
from the system server. According to the methodology that we
k k !
0 Si
adopt in paragraph 3.1, to determine the cumulative
k
k
i1
distribution function P(R r) of the resources occupation of
N
N c
the system, it is necessary to determine a cdf
P(R
r)
.exp i . FS (s)ds k
i
i
k N
i1 0
k
i1
conditional on k customers served by the system.
We decompose this number k of served customers into the
k ! . 0 F (s)ds
c
Si
Si
number ki of customers per service such as k k1 … kN .
i1
k k !
.exp .
N
c
c
F (s)ds
We can have
k1,…, kN
k1 !…kN !
possible combinations of
0
Si
i1
these customers. For each customer, the probability that he uses an amount of resource less than or equal to x is equal to
k k
. exp(s)ds .exp . exp(s)ds
Gi (x) if he requested service i (or with a probability of
pi in
k ! 0 0
other words).
k 1 k
.
k !
.exp
The probability that all customers then use a total amount of resources less than or equal to x for a given vector
k
k
exp
k1 ,…, kN is expressed by:
k ! k
p G (x
)… p G (x )
k ,…, k
1 1 11 1 1 1k
By denoting / .
1 N x11 … x
NkN x
1
k1 users of service 1
(7)
Then,
… pN GN (xN1 )… pNGN (xNkN )
r r k
r k
k users of service N
P(R r) k
k 0
k !
k !
k 0
exp exp . N
k !
k !
k 0
where
xij
denotes the amount of resources used by jth
B P(R C)
customer of service i.
k N *N
. pki .G*ki (x)
k
k
Equation (7) can be rewritten as a discrete convolution product:
k 0
k1…kN k xC
k1,…, k
N i1
i i
i1
k
It is recalled that this probability must be less than a
k ,…, k
p1G1 (x11 )… p1G1 (x1k1 )
tolerance for the determination of the sufficient capacity of
1 N x11 … xNkN x
k1 users of service 1
… pN GN (xN1 )… pN GN (xNkN )
kN users of service N
the system.
Remark:
When the amount of resources requested by customers, regardless of the service used, is an independent and
k p G (x) *…* p G (x)*
identically distributed random variable with a cdf
G(r) ,
k ,…, k
1 1 1 1
equation (9) can be written as:
1 N
k1 consecutive convolution products
…* pN GN ( x) * …* pN GN (x)
P(R r)
k *k N k
kN consecutive convolution products
k .
G
(x). pi i
k
k 0
k1…kN k xr k1,…, kN
i1
pk1 … pkN .G*k1 (x) *…* G*kN (x)
1
N 1
N
1
N 1
N
k1 ,…, kN


RESULTS AND DISCUSSIONS
Lets take a case study of resources sharing between
where
f (t) * g(t) f (t1 )g(t2 ) f (s)g(t s)
for
customers arriving at a system of a capacity C which will be
t1t2 t s0
determined. The system can provide one of the N = 4 services
two functions
f (t) and g(t) with positive variable t, and
to each customer. The arrivals of these customers to this
f *n (t) f (t) *…* f (t) the n successive convolution products
system form Poisson processes of respective intensity
of the function f (t) by itself.
Thus, with all the possible values of the vector k1 ,…, kN
1 2, 2 3, 3 2.5, and 4 5 arrivals per unit of time for the service i = 1, 2, 3, and 4. By the principle of superposition of the Poisson processes, these arrivals form a
and all the values of x r , we obtain the cdf
P(Rk r)
single process with intensity 4
i1
i1
i
12.5 .
conditionally to k customers served by the system:
The amounts of resources required for these services are
P(Rk
r)
k pk1 … pkN .G*k1 (x) *…* G*kN (x)
respectively R1 4, R2 3, R3 2, et R4 1 resources.
The service durations are assumed to be random following
k ,…, k 1 N 1
N (8)
an exponential distribution with respective means
k1 …kN k xr 1 N
1 / 1 / 3, 1/ 1 / 4, 1/ 1 / 5, and 1/ 1 / 6 units of
k
N *N
ki
*ki
1 2 3 4
time for the services i = 1 to 4.
k ,…, k
pi .Gi
(x)
From equation (5), the cumulative distribution function of
k1 …kN k xr 1
N i1
i1
the occupied resources in the infinite capacity system is equal
*N
where
f (x) f (x) *…* f
(x) for a family of functions
<>to:
i 1 N
i1
fi (x) .
k k k
i1,…,N
P(R r) k . . p11 …p44
And the cdf of the amount of occupied resources is given by:
k 0 k1 ,…,k4 k1 ,…, k4
k1R1 …k4R4 r
P(R r)
Such that pi i / denotes the probability that the
k N
*N
ki *ki
(9)
customer arriving at the system will use the service i.
And:
k .
k ,…, k
pi .Gi
(x)
k 0
k1…kN k xr 1
N i1
i1
k
4
k
4
For the case of a finite capacity of resources C, we can
k
k ! . 0
e isds
.exp .0
e isds
deduct from equation (9) the probability of the system being
i1
i1
blocked:
k 4
1 k
4 1
.
k !
.exp .
i1
i
i1 i
Curves in Figures 2a and 2b represent this cumulative distribution function:
Fig. 2a. Cdf of the amount of occupied resources for an infinite capacity
Figure 3. Resource sharing model for a system having four service
classes for its users
The probability of congestion is calculated from the percentage of time during which the amount of occupied resources equals to the capacity of the system. We then use the parameter NbUsedResources in the below algorithm to determine this percentage.
count = 0;
capacity = 11;
t = 0:0.01:2000;
nbUsedResourcesResampled = resample(NbUsedResources,t,'zoh'); usedResRsmpl = nbUsedResourcesResampled.Data; [timeResRsmpl,y]=size(usedResRsmpl);
for i = 1:timeResRsmpl
if usedResRsmpl(i,1)>= capacity count = count + 1;
end
Fig. 2b. Cdf of the amount of occupied resources for an infinite capacity
(zoomed in)
It is obvious from Figure 2a that the deployment of more than 20 resources does not influence the system. We deduce that the probability that the amount of occupied resources is greater than or equal to 20 is very small (close to zero). It is therefore not necessary to provide more than 20 resources to the system since not all of them will be used.
From these curves and equation (6), we can conclude a probability of congestion of 10% when the system has 11 resources, and this probability becomes 1% with 17 resources.
We used MATLAB – Simulink to simulate these two values of total amount of resources to verify these congestion probabilities.
Model presented in Figure 3 has been used. The customers (users) are generated from the User of service i blocks following a respective rate i described above.
end
pCongestion = count/timeResRsmpl; disp('Capacity = 11');
fprintf('The percentage of time the capacity is reached : %0.3g \n', pCongestion);
We found the following percentages for a capacity of 11 resources, and for a capacity of 17 resources at 5000 simulation time units.
>>
Capacity = 11
The percentage of time the capacity is reached : 0.1073
>>
Capacity = 17
The percentage of time the capacity is reached : 0.0166
These are the probabilities that any user, arriving at a random moment, will find all the resources occupied.
We then found the theoretical results of our model explained by the curves of Figures 2a and 2b which are respectively 10% and 1% for capacities 11 and 17. These results validate the model thus established.
As a concrete example, the model can be used for the dimensioning of codes (resources) of a third generation mobile communication network (3G network) which share these codes through its customers [25] [26]. This type of network provides multiple services, such as:

AMR voice service: 2 codes needed

VP service: 8 codes needed

PS R99 DL service: 8 codes needed

HSDPA service: 16 codes needed
These results remain valid for the case of nondiscrete resources. In this case, the amount of resources initially available in the system is a strictly positive real number. The amounts of resources necessary for each service i can also be strictly positive real numbers (not necessarily a nonzero natural number). The random variable designating the total amount of resources occupied by the customers is then a positive real random variable. And the expressions of equations (5) and (6) are already carried with the cumulative distribution functions which characterize the real random variables.


CONCLUSION
A multiservice system shares its resources through its customers according to each customer's requirements for each service requested. They arrive on this system according to Poisson processes characteristic of each serviceclass. The system is assimilated with a multiservice M/GI/C queue in terms of the number of customers. Two cases of capacities in terms of amount of resources were considered: the infinite case which allowed us to evaluate the distribution of occupancy of the queue resources, and the finite case from which we derived a probability of congestion from the infinite case. The analytical expressions established have made it possible for us to dimension the system by determining the amount of resources necessary for the server with an admissible probability of congestion. From the theoretical result and the case study presented, we found that it is quite possible to determine this necessary amount since from a certain value of the capacity, the characteristic of the system in terms of the number of occupied resources does not show much variation. We have shown that the result of the FryMolina model is found from our model and from our result in the singleservice case.
REFERENCES

Brockmeyer, E.; Halstrom, H.; Jensen, A.; Erlang, A.K. The life and works of A.K. Erlang. Copenhagen Telephone Co, 1948

Erlang, A.K. LÃ¸sning af nogle Problemer fra Sandsynlighedsregningen af Betydning for de automatiske Telefoncentraler (Soltion of some Problems in the Theory of Probabilities of Significance in Automatic Telephone Exchanges). Elektroteknikeren. 1917, vol 13, pp 513.

Palm, C. Table of the Erlang Loss Formula. Telefonaktiebolaget LM Ericsson, Stockholm, 1947.

Joys, L.A. Comments on the Engset and Erlang Formulae for Telephone Traffic Losses. Thesis. Report TF No. 25,/71,
Research Establishment, The Norwegian Telecommunications Administration. 1971.

Iversen, V.B. The Aformula. Teleteknik, English ed. 1980, vol 23 nÂ°2, pp 6479.

Vaulot, E.; Chaveau, J. Extension de la formule dErlang au Cas oÃ¹ le Trafic est Fonction du Nombre dAbonnÃ©s OccupÃ©s (Extension of the Erlang formula to the case where the traffic depends on the number of busy subscribers). Annales de TÃ©lÃ©communications. 1949, vol 4 nÂ°8, pp 319324.

Iversen, V.B. The exact evaluation of multiservice loss system with access control. Teleteknik, English ed. 1987, vol 31 nÂ°2, pp 5661.

Hanczewski, S.; Stasiak, M.; Weissenberg, J. A queueing model of a multiservice system with statedependent distribution of resources for each class of calls. IEICE Trans. Commun. 2004, vol E97B nÂ°8, pp 15921605.

Hanczewski, S.; Stasiak, M.; Weissenberg, J. The queueing model of a multiservice system with dynamic resource sharing for each class of calls. Computer Networks, CN 2013, Communications in Computer and Information Science. 2013, vol 370, pp 436445

Hanczewski, S.; Kmiecik, D.; Stasiak, M.; Weissenberg, J. Multi service queuing system with elastic traffic. Proc. IEICE Gen.
Conf. 2016, BS318

Bonald, T.; ProutiÃ¨re, A.; Virtamo, J. A queueing analysis of max min fairness, proportional fairness and balanced fairness.
Queueing Systems. 2006, vol 53 nÂ°1, pp 6584

Haddad, J.P.; Mazumdar, R.R. Congestion in large balanced multirate networks. Queueing Systems. 2013, vol 74 nÂ°2, pp 333 368

Hanczewski, S.; Kaliszan, A.; Stasiak, M. Convolution Model of a Queueing System with the cFIFO Service Discipline. Hindawi Publishing Corporation, Mobile Information Systems. 2016, vol 2016, Article ID 2185714, pp 115

Ravaliminoarimalalason, T. B; Randimbindrainibe, F. Distribution of Occupied Resources on A Discrete Resources Sharing in A Queueing System. International Journal of Engineering Research and Technology. 2021, vol 10 nÂ°1, pp 638 643

Ravaliminoarimalalason, T. B; Randimbindrainibe, F.; Rakotomalala M., Distribution of occupied resources on a fractional resource sharing in a queueing system. Int. J. of Operational Research. 2021. DOI 10.1504/IJOR.2021.10039107

Crane, H.; McCullagh, P. Poisson superposition processes.
Journal of Applied Probability. 2015, vol 52 nÂ°4, pp 10131027

Hokstad, Per. Approximations for the M/G/m Queue. Operations Research. 1978, vol 26 nÂ°3, pp 510523

Tijms, H. C.; Van Hoorn, M. H.; Federgruen, A. Approximations for the SteadyState Probabilities in the M/G/c Queue. Advances in Applied Probability. 1981, vol 13 nÂ°1, pp 186206

Ma, B. N. W.; Mark, J. W. Approximation of the Mean Queue Length of an M/G/c Queueing System. Operations Research. 1995, vol 43 nÂ°1, pp 158165.

Molina, E.C. The Theory of Probability Applied to Telephone Trunking Problems. The Bell System Technical Journal. 1922, vol 1 nÂ°2, pp 6981

Fry, T.C. Probability and its Engineering Uses. New York, 1928

Molina, E.C. Application of the Theory of Probability to Telephone Trunking Problems. The Bell System Technical Journal. 1927, vol 6 nÂ°3, pp 461494

Iversen, V.B. Teletraffic engineering handbook, ITUD SG 2/16 & ITC, 2001

Iversen, V.B. Teletraffic engineering and network planning,
DTU Library, 2015

Kasapovi, S.; Sarajli, N. OVSF code assignment in UMTS networks. 20th International Crimean Conference Microwave & Telecommunication Technology. 2010, pp 429432

Bernhard, H.W.; Seidenberg, P.; Althoff, M.P. UMTS : The Fundamentals. John Wiley & Sons, 2003