







































Copyright © CC-BY-NC 2019, BJMSR 

 

Bangladesh Journal of Multidisciplinary Scientific Research; Vol. 1, No. 2 
ISSN 2687-850X      E-ISSN 2687-8518  

July-September; 2019 
Published by Centre for Research on Islamic Banking & Finance and Business  

 

     46 

Application of Single-Server Queue System in Performance Analysis of Shuttle 
Bus Operation: A Case Study of Federal University of Technology Akure 

 
 
 

Kanyio Olufunto Adedotun 
Department of Transport Management Technology 

The Federal University of Technology,Akure  
Ondo State, Nigeria 

E-mail: oobembe97@gmail.com 
 
 
Abstract 
This study has examined the performance of University transport bus shuttle based on utilization using a Single-server queue 
system which occur if arrival and service rate is Poisson distributed (single queue) (M/M/1) queue. In the methodology, Single-
server queue system was modelled based on Poisson Process with the introduction of Laplace Transform. Also, PASTA was 
introduced in queuing systems with Poisson arrivals. It is concluded that the performance of University transport bus shuttle is 
93 percent which indicates a very good performance such that the supply of shuttle bus in FUTA is capable of meeting the 
demand. This study can be improved upon by examining the peak and off-peak period of traffic in the two major corridors 
(North gate and South gate) of FUTA, the economic cost of operating bus shuttle services can also be examined. 
 
Keywords: Single-Server Queue System; Transit Shuttle Buses; FUTA  
 
 
1.Introduction 
The issues arising from transportation has continually subjected to various debates in the urban societies. Globally, several 
attempts have been made to tackle the challenges, although the situation is not getting much better (Aderamo, 2012; Adanikin, 
Olutaiwo, and Obafemi, 2017; Sidiq, 2019). Managing transport infrastructures is crucial to facilitate accessible, affordable, 
reliable, safe, and efficient that movement of people and goods, which can be achieved by continous assessment of transport 
performance indicators.  Among the various modes of transport, the road transport is highly predominant which offers door-to-
door service as in the case of the Federal University of Technology Akure (FUTA). Road transport infrastructure in the 
University is commonly plied by shuttle buses for the movement of the students, staff, members and non-members of the 
Institution to and from the University on a daily basis. 

Among the noticeable transportation problems in the University are traffic congestion; longer commuting; public 
transport inadequacy; difficulties for tricycles to have access to routes being plied by shuttle buses, challenge of freight 
distribution from one end of the University to another end, and other challenges which all impacts the performance of the 
University transport shuttle. This study concentrates on the performance of University transport bus shuttle with the aim of 
examining the bus shuttle efficiency based on utilization. 

Adeniran and Kanyio (2019) have laid a foundation of model on single-server queue system which this study will 
absolutely rely on. A similar study was conducted by Adanikin, Olutaiwo, and Obafemi (2017) on the performance study of 
University of Ado Ekiti (UNAD) transit shuttle buses. They adopted traffic volume, speed, density and revenue as main 
parameters of performance of transport shuttles, and find that the morning peak period (8.00am to 9.00am) has 234 
vehicles/hr, evening peak period (2.00pm to 3.00pm) has 284 vehicles/hr, while the off-peak period (11.00am to 12.00pm) 
has 156 vehicles/hr. Also, The average stopping time was 6.55 minutes, average interval between arrivals of motorists was 16.40 
seconds, the average queue length was 14.23 people, and the average waiting time at the bus-stop 4.17 minutes. These values 
were obtained using the queuing theory and shows much commuters time is lost on transit queues. This study focuses on the 
performance of bus terminal in FUTA, and does not factor in other parameters such as peak period, traffic volume, traffic speed, 
density, and others. 

A very close study was examined by Sidiq (2019) on single-server queue system of shuttle bus performance in the 
Federal University of Technology Akure. In his study, Single server queue system was modelled based on Poisson Process with 
the introduction of Laplace Transform which was culled from Adeniran and Kanyio (2019). The study finds that the 
performance of University transport bus shuttle is 96.6 percent which indicates a very good performance such that the supply of 

mailto:oobembe97@gmail.com


 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

47 
               

shuttle bus in FUTA is capable of meeting the demand. This present study was conducted barely two weeks after the study of 
Sidiq was conducted, and will replicate the methodology of Sidiq as culled from Adeniran and Kanyio (2019). 

 
2. Methodology 
2.1 Concept of Queuing  
The concept of queue was first used for the analysis of telephone call traffic in 1913 (Copper, 1981; Gross and Harris, 1985; 
Bastani, 2009). In a system that deals with the rate of arrival and service rate, waiting time is inevitable and it is always 
influenced by queue length. It is therefore crucial to minimize the waiting time to the lowest level in the bus terminal (Jain, 
Mohanty and Bohm, 2007). This is referred to as queuing system (Adeniran and Kanyio, 2019).  
The basic application of queue is shown in Figure 1, also the basic quantities are: 

i. Number of customers in queue L (for length); 
ii. Time spent in queue W for (wait) 

 

Figure1: Basic application of queue 

Source: Adeniran and Kanyio (2019) 

Examples of queue system are: 

1. Single-server queue system: This is also referred to as single queue, single server. It is simple if arrivals and services are 
Poisson distributed (M/M/1) queue. It has limited number of spots and not difficult. Figure 2 depicts single-server 
queue system. 

 

Figure 2: Single-server queue system 

Source: Adeniran and Kanyio (2019) 

2. Multi-server queue system: This is comprises of single queue, many servers (M/M/c) queue. The c is referred to as 
Poisson servers. Figure 3 depicts multiple-server queue system. 

 

Figure 3: Multiple-server queue system 

Source: Adeniran and Kanyio (2019) 



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

48 
               

In single-server queue system, arrival and service processes are Poisson such that 

a. Customers arrive at an average rate of λ per unit time; 

b. Customers are serviced at an average rate of µ per unit time; 

c. Interarrival and inter-service time are exponential and independent; 

d. Hypothesis of Poisson arrivals is reasonable; and 

e. Hypothesis of exponential service times are not so reasonable (Adeniran and Kanyio, 2019) 

In order to explain how the queuing system works, there is need to first introduce the Poisson Process (PP). It has exceptional 
properties and is a very important process in queuing theory. To simplify the model, we often assume customer arrivals follow a 
PP. The Laplace Transform (LT) is also a very powerful tool that was adopted in the analysis (Trani, 2011). Apart from PP 
and LT, there is focus on the queue model itself (Adeniran and Kanyio, 2019). 

2.2 Modelling of Single Queue System 
2.2.1 Queuing system from Poisson Process and “PASTA” 
The Poisson Process (PP) is important in queue theory due to its outstanding properties. According to Adan and Resing 

(2015), queuing system is achieved as “let N(t) be the number of arrivals in [0, t] for a PP with rate λ, i.e. the time between 

successive arrivals is exponentially distributed with parameter λ and independent of the past. Then N(t) has a Poisson 

distribution with parameter λt. This is culled from Adeniran and Kanyio (2019). 

P(N(t) = k) = 
(λt)k

k!
 e−λt, for k = 0, 1, 2………………… Equation 1 

The mean, and coefficient of variation of N(t) are  

Mean: E(N(t)) = λt;  

Coefficient of Variation: c2N(t) = 
1

λt
 ………………… Equation 2 

By the memoryless property of Poisson distribution, it can be verified that 

P(arrival in (t, t + ∆t]) = λ∆t + 0(∆t) ………………… Equation 3 

Hence, when ∆t is small, 

P(arrival in (t, t + ∆t)) ≈ λ∆t ………………… Equation 4 

In each small time interval of length ∆t the occurrence of an arrival is equally likely. In other words, Poisson arrivals occur 
completely randomly in time. The Poisson Process is an extremely useful process for modelling purposes in many practical 
applications. An important property of the Poisson Process is called “PASTA” (Poisson Arrivals See Time Averages). 

PASTA is meant for queuing systems with Poisson arrivals, (M/./. systems), arriving vehicles find on average the same situation 
in the queuing system as an outside observer looking at the system at an arbitrary point in time. More precisely, the fraction of 
vehicles finding on arrival the system in some state A is exactly the same as the fraction of time the system is in state A. 

2.2.2 Laplace Transform 
The Laplace transform LX(s) of a nonnegative random variable X with distribution function f(x) is define as: 

LX(s) = E(e−sX) = ∫ e−sX∞

x=0
 f(x)dx  ………………… Equation 5 

It can be noted that  

LX(0) = E(e−X .0) = E(1) = 1  ………………… Equation 6 

and 

L1
X(0) = E((e−sX)1)|s=0 

 = E(−Xe−sX)|s=0 



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

49 
               

 = −E(X) ………………… Equation 7 

Correspondingly, 

L (k)
x

 (0) = (−1)kE(Xk) ………………… Equation 8 

There are many useful properties of Laplace Transform. These properties can make calculations easier when dealing with 
probability. For instance, let X, Y, Z be three random variables with  

Z = X +Y and X, Y are independent.  

Then the Laplace Transform of Z can be found as: 

LZ(s) = LX(s) · LY (s) ………………… Equation 9 

Moreover, when Z with probability P equals X, with probability 1 − P equals Y, then  

LZ(s) = PLX(s) + (1 − P)LY (s) ………………… Equation 10 

Laplace Transforms of some useful distributions can now be introduced. 

a. Suppose X is a random variable which follows an exponential distribution with rate λ. The Laplace Transform of X is  

LX(s) = 
(λ)

λ + s
  ………………… Equation 11 

b. Suppose X is a random variable which follows an Erlang − r distribution with rate λ. Then X can be written as: 

X = X1 + X2 + · · · + Xr  ………………… Equation 12 

where Xi are i.i.d. exponential with rate λ. Therefore, we have 

LX(s) = LX1(s) · LX2(s)… LXr(s) 

 = (
λ

λ + s
)

n

  ………………… Equation 13 

c. Suppose X is a constant real number c, then 

LX(s) = E(e−sX) 

 = E(e−sc) 

 = e−sc   ………………… Equation 14 (culled from Adeniran and Kanyio, 2019) 
 

2.2.3 Basic queuing systems 
Kendall’s notation shall be used to describe a queuing system as denoted by:  
A/B/m/K/n/D   …………………. Equation 15 (Adan and Resing, 2016) 
Where 
 A: distribution of the interarrival times 
 B: distribution of the service times 
 m: number of servers 
 K: capacity of the system, the maximum number of passengers in the system including        the one being 
serviced 
  n: population size of sources of passengers 
 D: service discipline 
G shall be used to denote general distribution, M used for exponential distribution (M stands for Memoryless), D be used for 
deterministic times (Sztrik, 2016).  
A/B/m is also used to describe a queuing system, where:  
 A stands for distribution of interarrival times,  
 B stands for distribution of service times and  
 m stands for number of servers.  
Hence M/M/1 denotes a system with Poisson arrivals, exponentially distributed service times and a single server.  
M/G/m denotes an m- server system with Poisson arrivals and generally distributed service times, and so on. 
In this section, the basic queuing models (M/M/1 system), which is a system with Poisson arrivals, exponentially distributed 
service times and a single server. The following part is retrieved from Queuing Systems (Adan and Resing, 2016). 

Firstly, it is assumed that inter-arrivals follows an exponential distribution with rate λ, and service time follows the exponential 
distribution with rate µ. Further, in the single service model, to avoid queue length instability, it is assume that: 
According to Adanikin, Olutaiwo and Obafemi (2017),  



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

50 
               

 

Utilization (R) = 
Average Arrival Rate (λ)

Average service rate (μ)
 < 1 …………………….. Equation 16   

 
Here R is the fraction of time the server is working (called the utility factor). Time-dependent behaviour of this system will be 
considered firstly, then the limiting behaviour. Let Rn(t) denote the probability that at time t there are n passengers in the 
system.  

Then by equation 3, when ∆t → 0, 

R0(t + ∆t) = (1 − λ∆t)R0(t) + µ∆tR1(t) + 0(∆t) …………………….. Equation 17 

Rn(t + ∆t) = λ∆tRn−1(t) + (1 − (λ + µ)∆t)Rn(t) + µ∆tRn+1(t) + 0(∆t) ……………….. Equation 18 
where n = 1,2, ... 

Hence, by tending ∆t → 0, the following infinite set of differential equations for Rn(t) will be obtained. 

R1
0 (t) = −λR0(t) + µR1(t) …………………….. Equation 19  

R1
 n(t) = λpn−1(t) − (λ + µ)Rn(t) + µRn+1(t), n = 1,2, …………………….. Equation 20 

It is very difficult to solve these differential equations. However, when we focus on the limiting or equilibrium behaviour of this 
system, it is much easier. 

It was revealed by (Sztrik, 2016) that when t →∞, R1
n(t) → 0 and Rn(t) → Rn. It follows that the limiting probabilities Rn 

satisfy equations 

0 = −λR0 + µR1………………………….. Equation 21 

0 = λRn−1 − (λ + µ)Rn + µRn+1,  , n = 1,2, ……………………… Equation 22 
Moreover, Rn also satisfy 

∑ Rn = 1∞
n=0  ………………………….. Equation 23 

which is called the normalization equation. We can also use a flow diagram to derive the normalization equations directly. For 
an M/M/1 system, the flow diagram is shown in figure 4: 

 
Figure 4: Process diagram for M/M/1 Queue, k=1,2,3,... (Ademoh and Anosike, 2014; Adeniran and Kanyio, 2019) 
The rate matrix of the system is: 
 

Q =  

−λ λ 0 0 0 . . .
µ −(µ +  λ) λ 0 0 . . .

0 µ −(µ +  λ) λ 0 . . .

0 0 µ −(µ +  λ) λ . . .
. . . . . . . .  
. . . . . . . .  
. . . . . . . .

  ………………….. Equation 24 

 
 
 
Notice that the sum of each row equals 0. 
In order to determine the equations from the flow diagram, a global balance principle was adopted. Global balance principle 
states that for each set of states under the equilibrium condition, the flow out of set is equal to the flow into that set. Based on 
figure 1,   
  



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

51 
               

   

      𝑆𝑡𝑎𝑡𝑒 𝑅𝑎𝑡𝑒 𝐼𝑛 = 𝑅𝑎𝑡𝑒 𝑂𝑢𝑡
0 µ𝑅1 = 𝜆𝑅0

1 𝜆𝑅0 + µ𝑅2 = (𝜆 + µ)𝑅1

2 𝜆𝑅1 + µ𝑅3 = (𝜆 + µ)𝑅2 

 

 
 
This is exactly the normalization equation. To solve the equation, firstly, it was assume that 

(R) = 
𝐴𝑣𝑒𝑟𝑎𝑔𝑒 𝐴𝑟𝑟𝑖𝑣𝑎𝑙 𝑅𝑎𝑡𝑒 (𝜆)

𝐴𝑣𝑒𝑟𝑎𝑔𝑒 𝑠𝑒𝑟𝑣𝑖𝑐𝑒 𝑟𝑎𝑡𝑒 (𝜇)
  which is known as the utilization factor. From the equilibrium equation of state 0, we have: 

R1 = Rp0 …………………………. Equation 25 
When equation 25 was substituted into the equilibrium equation of state 1, then: 

λp0 + µp2 = (λ + µ)Rp0 

      = 
𝜆2

𝜇
 p0 …………………………….. Equation 26 

That is 

µp2 =  
𝜆2

𝜇
p0 ……………………………. Equation 27 

Therefore, 
 p2 = R2p0 …………………………. Equation 28 
Generally, 
 pk = Rkp0 …………………………. Equation 29 
Since 

∑∞
𝑛=0 pn = 1 …………………………………. Equation 30 

Using (1.29), we can replace pk by R0. Then 

∑∞
𝑛=0 Rnp0 = 1…………………………….. Equation 31 

That is 
1

1−𝑅
 p0 = 1 

p0 = 1 – R………………………………………. Equation 32 
Moreover, for any k, 

pk = Rk(1 −R) …………………………………. Equation 33 
Finally is the limiting probability pk in the M/M/1 system. The expected queue length L is given by 

E(L)  = ∑ i. pi

∞

i=0

 

= ∑ i. Ri 

∞

i=0

(1 − R) 

= R(1 − R) ∑ i. Ri 

∞

i=0

 

= R(1 − R) (∑〖i. Ri)i

∞

i=1

〗 

 = R(1-R) (
1

1−R
)i 

 = 
R

1−R
..............................Equation 34 (Culled from Adeniran and Kanyio, 2019) 

 
3. Results and Discussion 
3.1. Traffic Survey 
3.1.1 Stopping time of shuttle bus 
Stopping time refers to the total time duration the shuttle bus spends at the bus stop. The stopping time is made up of: 

a) The boarding stop time “A”: This is also made up of the time taken to close the door = 12 seconds; time taken by the 
driver to check the traffic before take-off = 5 seconds; and time taken to park the bus and open the bus for commuters 
= 5 seconds. 

 A = (12+5+5) = 22 seconds 
b) The average boarding time per passenger = “B1” = 17 seconds 
c) Number of passengers boarding = n1 = 18 passengers.  



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

52 
               

Mathematically, stopping time T = A + B1* n1  
Stopping Time = (22 + (17*18)) = (22+306) = 328 seconds 
Stopping Time = 5.47 minutes 
 
3.2 Waiting Time  
This is the length of time spent by the passengers at the bus stop before boarding a bus. It is also referred to as Delay. The 
queuing theory is employed in this study. 

a) Average arrival rate (λ) = 226 passengers/hour 

λ =
226

3600
  = 0.0628 

b) Average service rate (μ) = 243 passengers/hour 

μ =
243

3600
 = 0.0675 

c) Average interval between arrival =
1

λ
 

=
1

0.0628
 = 15.924 Seconds 

d) Average interval between service rate =
1

μ
   

=
1

0.0675
 = 14.815 Seconds 

e) Average queue length =
𝜆2

𝜇(𝜇−𝜆)
 

=
0.06282

0.0675(0.0675−0.0628)
 = =

0.0039

0.00032
= 12.19 Passengers 

f) Average waiting time in the queue =
λ 

𝜇(𝜇−𝜆)
 

=
0.0628

0.0675(0.0675−0.0628)
 = =

0.0628

0.00032
= 196.25 Seconds = 3.27 Minutes 

g) Average time spent in the system (bus stop) =
1

μ−λ 
 

=
1

0.0675 − 0.0628
= 212.76 Seconds = 3.55 Minutes 

h) Efficiency of shuttle bus operation based on bus stop utilization (R) = 
𝜆

𝜇
 

R =  
0.0628

0.0675
  = 0.93 

It is important to note that the Utilization factor is less than 1(R< 1), hence the performance of University transport bus 
shuttle is 93 percent which indicates a very good performance such that the supply of shuttle bus in FUTA is capable of meeting 
the demand.  
 
4. Conclusion and Recommendation 
This study has examined the performance of University transport bus shuttle based on utilization using a Single-server queue 
system which occur if arrival and service rate is Poisson distributed (single queue) (M/M/1) queue. It is concluded that the 
performance of University transport bus shuttle is 93 percent which indicates a very good performance such that the supply of 
shuttle bus in FUTA is capable of meeting the demand. This result is very close to that of Sidiq (2019) which finds that the 
performance of University transport bus shuttle is 96.6 percent. This study can be improved upon by examining the peak and 
off-peak period of traffic in the two major corridors (North gate and South gate) of FUTA, the economic cost of operating bus 
shuttle services can also be examined.  
 
References 
Adan, I., and Resing, J. (2015). Queuing Systems. Department of Mathematics and Computer 

Science, Eindhoven University of Technology, Netherlands, 2015. 
Adanikin, A., Olutaiwo, A., and Obafemi, T. (2017). Performance Study of University of Ado Ekiti (UNAD) Transit Shuttle 

Buses. American Journal of Traffic and Transportation Engineering, 2(5): 67-73 doi: 10.11648/j.ajtte.20170205.12. 
Ademoh, N. A., and Anosike, E. N. (2014). Queuing Modelling of Air Transport Passengers of Nnamdi Azikiwe International 

Airport Abuja, Nigeria Using Multi Server Approach. Middle East Journal of Scientific Research 21 (12): 2326-2338, 
2014 doi: 10.5829/idosi.mejsr.2014.21.12.21807. 

Adeniran, A. O., and Kanyio, O. A. (2019). Quantitative Model of Single-Server Queue System. Indian Journal of Engineering, 
16, 177-183. 

Aderamo, A. J. (2012). Urban transportation problems and challenges in Nigeria: A planner’s view. Prime Journals. 2(3), 198-
203. 



 Copyright © CC-BY-NC 2019, BJMSR 

 

 www.cribfb.com/journal/index.php/BJMSR       Bangladesh Journal of Multidisciplinary Scientific Research         Vol. 1, No. 2; 2019 

53 
               

Bastani P. (2009). A Queuing Model of Hospital Congestion. MSc. Thesis submitted to Department of Mathematics Simon 
Fraser University 

Copper, R. B. (1981). Introduction to Queuing Theory, 2nd Edition North Holland. 
Gross, D., and Harris, C. M. (1985). Fundamentals of Queuing Theory, 2nd  ed. John Wiley & Sons: New York. 
Jain, J. L., Mohanty, S. G., and Bohm, W. (2007). A Course on Queuing Models, Statistics: A series of Textbooks and 

Monographs, Chapman & Hall/CRC, Taylor & Francis Group. 
Sidiq, O. Ben. (2019). Single-Server Queue System of Shuttle Bus Performance: Federal University of Technology Akure as 

Case Study. American International Journal of Multidisciplinary Scientific Research, 5(3), 9-14 
Sztrik, J. (2012). Basic Queuing Theory. University of Debrecen, Faculty of Informatics. 
Trani, A., A. (2011). Introduction to Transportation Engineering, Introduction to Queuing Theory. Virginia Polytechnic 

Institute and State University, US, Pp: 2-46. 

 
 
 
Copyrights  
Copyright for this article is retained by the author(s), with first publication rights granted to the journal. This is an open-access 
article distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.org/licenses/by/4.0/). 
 


