Next Article in Journal
Convolution Properties of Certain Classes of Analytic Functions Defined by Jackson q-Derivative
Next Article in Special Issue
Remote Geotechnical Monitoring of a Buried Oil Pipeline
Previous Article in Journal
Graded Derived Equivalences
Previous Article in Special Issue
Fluid-Flow Approximation in the Analysis of Vast Energy-Aware Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Queueing-Inventory with One Essential and m Optional Items with Environment Change Process Forming Correlated Renewal Process (MEP) †

by
Jaison Jacob
1,
Dhanya Shajin
2,
Achyutha Krishnamoorthy
3,*,
Vladimir Vishnevsky
4,* and
Dmitry Kozyrev
4,5
1
Department of Mathematics, St. Aloysius College, Elthuruth, Thrissur 680611, India
2
Department of Mathematics, Sree Narayana College, Chempazhanthy, Thiruvanathapuram 695587, India
3
Centre for Research in Mathematics, CMS College, Kottayam 686001, India
4
V.A. Trapeznikov Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya Street, 117997 Moscow, Russia
5
Applied Probability and Informatics Department, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya Street, 117198 Moscow, Russia
*
Authors to whom correspondence should be addressed.
This paper is an extended version of our paper published in Ann. Oper. Res. (February 2021). https://doi.org/10.1007/s10479-021-03975-8. Accepted on 30 January 2021, Published on 22 February 2021.
Mathematics 2022, 10(1), 104; https://doi.org/10.3390/math10010104
Submission received: 12 November 2021 / Revised: 15 December 2021 / Accepted: 23 December 2021 / Published: 29 December 2021
(This article belongs to the Special Issue Distributed Computer and Communication Networks)

Abstract

:
We consider a queueing inventory with one essential and m optional items for sale. The system evolves in environments that change randomly. There are n environments that appear in a random fashion governed by a Marked Markovian Environment change process. Customers demand the main item plus none, one, or more of the optional items, but were restricted to at most one unit of each optional item. Service time of the main item is phase type distributed and that of optional items have exponential distributions with parameters that depend on the type of the item, as well as the environment under consideration. If the essential item is not available, service will not be provided. The lead times of optional and main items have exponential distributions having parameters that depend on the type of the item. The condition for stability of the system is analyzed by considering a multi-dimensional continuous time Markov chain that represent the evolution of the system. Under this condition, various performance characteristics of the system are derived. In terms of these, a cost function is constructed and optimal control policies of the different types of commodities are investigated. Numerical results are provided to give a glimpse of the system performance.

1. Introduction

Queueing inventory models have been extensively analyzed since 1992. Very few of these discuss multi-commodity systems in randomly changing environments. Queueing systems which evolve under influences from external sources have, for a long time, inspired interest. In real life situations, inventory systems are often subject to randomly changing exogenous environment conditions that affect the demand for the product, the supply, and the cost structure. The area of queues in random environments is today a field of active research in applied probability. Queueing systems with correlated arrival flow of customers give adequate mathematical models for different real world systems including computer and telecommunication systems, and network protocols [1]. The following papers are relevant to the present paper only in that the authors consider multi-commodity inventory systems without any specified main and optional items. FaizAl-Khayyal et al. [2] consider a multi-commodity network model in maritime routing and scheduling. They tried to optimize the cost and the quantity of each commodity with constrained production rates, consumption rates, and storage capacities in each port. This paper addresses the common problem faced in the maritime transportation of petrochemical products. The problem is formulated as a mixed-integer non-linear programming problem.
Two-echelon, multi-commodity supply chain network design problem is considered by Hannan et al. [3]. The authors formulate the problem as a mixed-integer programming model in deterministic, single-period, and multi-commodity contexts. They develop a heuristic solution procedure based on Lagrangian relaxation due to the complexity of the problem. The problem deals with developing an optimum strategy in locating and sizing factories and warehouses and minimizing the total cost of the system which includes the costs associated with production, storage, transportation, and lead-times of commodities.
Jin QIN et al. [4] discuss an optimization problem in multi-commodity logistic network design. The problem is formulated as a non-linear mixed-integer programming model based on assumed normally distributed stochastic demands of the retailers. In this work, the strategic decisions regarding inventory controls and facility locations are incorporated simultaneously. The authors also developed a combined simulated annealing (CSA) algorithm to solve the problem. Optimization of the cost function involving costs associated with location, inventory, and transportation is also done.
Ronald et al. [5] examine a multi-commodity logistics network design problem with simultaneous emphasis on establishing ideal location facilities and distribution of commodities to ensure minimization of costs together with improved services. Factors such as the location of facilities and warehouses, storage capacity of the warehouses, and transportation routes involving all these are considered in the optimization problem. They later designed, tested, and compared a genetic algorithm and a specific problem heuristic on several realistic scenarios.
Claudio et al. [6] discuss multi-commodity inventory location models where inventory control policies are reviewed periodically and continuously under modular stochastic capacity constraints. The model is formulated as a mixed-integer non-linear programming model. The logistic problem of supplying certain commodities to the warehouses, which acts as customer service centers, from a single factory is under consideration. An objective function with factors associated with the selection of warehouses, type, and quantity of commodities to be assigned, type of customers to be served, etc., is optimized. They have developed a Lagrangian heuristic to obtain a feasible integer solution at each iteration of the subgradient method.
Ali et al. [7] examine a dynamic multi-commodity inventory and facility location problem in steel supply chain networks. Demand is assumed to be stochastic with normal distribution. In this model, the authors suggested a potential production capacity with emergency and shared safety stocks. The authors have presented a mixed-integer non-linear programming model and a mixed-integer linear programming model in this paper. The paper focuses on the strategic and tactical design of steel supply chain networks.
In Shajin et al. [8], the authors discuss a single server multi-commodity queueing inventory system with one essential and m optional items. They were the first to introduce the concept of optional items for sales/service. Customer arrival follows the Markovian arrival process, service completion with respect to the essential inventory follows Phase-type distribution and that with respect to optional inventories follows an exponential distribution. In this model, immediately after the service of an essential item, the customer either leaves the system with probability p or with probability 1 p the customer goes for optional item(s). The system is assumed to be idle either in the absence of an essential item or when there is no customer in the system. Each customer is allowed to purchase only one unit of the essential item, whereas more than one type of optional items can be purchased with an imposed restriction of, at most, one unit per item. The stability condition is obtained by using the well-known fact that the left drift rate should be less than that of the right drift. Optimization of the control variables with respect to the cost function is also done numerically.
The following papers deal with queueing inventory systems influenced by randomly changing environments. In Song et al. [9], the authors consider an inventory model where the rate of demand is dependent on the environment variables. These variables can be anything, such as different stages in the life cycle of the particular inventory or changes in various factors linked with the economy, etc. They not only derived basic characteristics of the optimal policies but also observed the influence of various patterns in problem data on optimal policies and developed algorithms for computing optimal policies
Özekici et al. [10] describe inventory models with unreliable suppliers in randomly changing environments. The environment change follows a Markov chain. The dependence of the stock-flow equations of the system on random environments is represented by a two-dimensional stochastic process. Under specified conditions, they have derived an optimality condition for the base-stock policy and (s, S) policy. Computational issues and some extensions are also determined.
A single item inventory model which is observed periodically in a randomly changing environment is considered in Erdem et al. [11]. All the model parameters are dependent on a time-homogenous Markov chain environment. The replenishment quantity is minimum{ Order quantity, Vendors capacity}. The problem is analyzed in single, multiple, and infinite periods. In all these cases, the authors prove that the optimal base-stock level depends on the state of the environment. Comparisons of the results with the case when the replenishment quantity equals the quantity ordered is also done.
Perry et al. [12] discuss production-inventory models with an unreliable facility operating in a two-state random environment. The system is characterized by a production machine. The production can even be stopped purposefully when there is a limited stocking capacity. When the machine is in ON period, the input into the buffer is assumed to be continuous and uniform until the threshold is reached, whereas the output from the buffer follows a compound Poisson process during OFF periods. Two different models are discussed and the factors controlling OFF periods are determined.
A continuous review (s, S) inventory system in a randomly changing environment is discussed in Feldman et al. [13] and its steady-state distribution obtained. The demand process is an environment dependent compound Poisson process when the environment is in a fixed state during an interval of time. The environmental process follows a continuous-time Markov process.
Kalpakam et al. [14] consider a lost sales (s, S) inventory system in a random environment. No backlog is allowed. The demand and supply rates are influenced by the environment process which is a finite irreducible Markov chain in continuous time. They have obtained the transform solution of the inventory level distribution and also an efficient algorithm to evaluate the long run system state is provided. Moreover, transient and limiting values of the mean reorder and shortage rates are also obtained. Goh et al. [15] discuss price-dependent inventory models with discount offers at random times. The offer is accepted when the inventory position is lower than a threshold level. Three different pricing policies are considered in which demand is induced by the retailer’s price variation. They have obtained expressions for optimal order quantities, prices, and profits under the assumptions of constant demand rates.
Highlights of this paper are:
  • It considers multi-commodity inventory with positive service time [16] in finite number of randomly changing environments;
  • The first paper to introduce optional items for service in random environments;
  • Except for one item (essential), all others are optional;
  • The customer demand process follows Markovian arrival process (MAP);
  • The environment change process follows marked Markovian environment arrival process (MMEAP)[n] of order n;
  • Service time of customers, being served with the essential inventory follows phase type distribution and that w.r.t optional item(s) follows exponential distribution (depending on the environment). The latter has a parameter, depending on the specific item(s) demanded by the customer.
The rest of the paper is organized as follows. The mathematical formulation of the model including the stability condition and the steady-state probability vector is described in detail in Section 2. Section 3 deals with some system performance measures and in Section 4, the construction of the cost function for optimizing the system control variables is discussed. Numerical illustrations and the numerical analysis of the cost function are discussed in Section 5. Section 6 gives the conclusion followed by references.
Notations and abbreviations used:
  • ( s , S ) ordering policy: An inventory policy which says that when the inventory level falls below a certain minimum number s, the order for replenishment is made to restore the inventory to a maximum number S;
  • e = Column vector of appropriate order with all its entries as 1 s;
  • 0 ¯ = Matrix of appropriate order with all its entries as 0;
  • I k = Identity matrix of order k;
  • [ G ] p q = ( p , q ) t h element of the matrix G;
  • C T M C : Continuous time Markov chain;
  • L I Q B D : Level Independent Quasi-Birth and Death process;
  • M A P = Markovian arrival process;
  • ( M M E A P ) [ n ] = Marked Markovian environment arrival process (with n distinct environments);
  • C D = The Kronecker product of two given matrices C m × n and D p × q , given by ( [ C ] p q D ) of order m p × n q ;
  • C D = The Kronecker sum of two square matrices B and C of orders m and n, respectively, given by C I n + I m D ;
  • C c o r = The Correlation coefficient;
  • M A P p = Markovian arrival process with positive C c o r ;
  • M A P n = Markovian arrival process with negative C c o r ;
  • P H ( γ , T ) : Phase type distribution with the initial probability vector γ and the transition generator matrix T.

2. Mathematical Formulation

Consider a single server multi-commodity queueing inventory system with one essential and m optional inventories in n random environments. Only one environment will be in operation at any given time. The arrival of customers follows Markovian arrival process ( M A P ) with representation ( H 0 , H 1 ), where each H i for 0 i 1 is of order m 3 . The generator matrix of the underlying C T M C ( δ ( t ) , t 0 ) on the state space { 1 , 2 , 3 , , m 3 } is given by H = H 0 + H 1 . These matrices, H 0 and H 1 , are of the form
H 0 = h 11 ( 0 ) h 12 ( 0 ) h 1 m 3 ( 0 ) h 21 ( 0 ) h 22 ( 0 ) h 2 m 3 ( 0 ) h m 3 1 ( 0 ) h m 3 2 ( 0 ) h m 3 m 3 ( 0 ) , H 1 = h 11 ( 1 ) h 12 ( 1 ) h 1 m 3 ( 1 ) h 21 ( 1 ) h 22 ( 1 ) h 2 m 3 ( 1 ) h m 3 1 ( 1 ) h m 3 2 ( 1 ) h m 3 m 3 ( 1 )
where, h i i ( 0 ) = j = 1 , j i m 3 h i j ( 0 ) + j = 1 m 3 h i j ( 1 ) for 1 i m 3 . Thus, h i j ( 1 ) , 1 i , j m 3 gives the transition rate from i t h state to j t h state through an arrival, while h i j ( 0 ) , 1 i , j m 3 , gives the transition from i t h state to j t h state without an arrival. Note that the transition rate between the i t h states, given by h i i ( 1 ) , 1 i m 3 occurs only with an arrival. Let η be the steady-state probability vector of H. Then, η satisfy η H = 0 and η e = 1. The fundamental rate λ A of this M A P is given by λ A = η H 1 e which gives the expected number of arrivals per unit of time. The coefficient of variation C v a r of intervals between arrivals is calculated as C v a r = 2 λ A η ( H 0 ) 1 e 1 and coefficient of correlation C c o r of intervals between successive arrivals is given as C c o r = ( λ A η ( H 0 ) 1 H 1 ( H 0 ) 1 e 1 ) / C v a r .
There are n environments that occur randomly and the occurrence of the environments follows marked Markovian environment arrival process ( M M E A P [ n ] ) with representation ( D 0 , D 1 , D 2 ,⋯, D n ), where each D i for 0 i n is of order m 2 . As stated earlier, only one environment will be in operation at any given time. The change in environment is directed by the stochastic process { V ( t ) ; t 0 } which is an irreducible continuous time Markov chain with the state space { 1 , 2 , , m 2 } . The sojourn time of this chain in the state v , 1 v m 2 , is exponentially distributed with parameter λ ( v ) . When the sojourn time in the state v expires, the process { V ( t ) ; t 0 } jumps to the state v without any change in the environment with probability p ( v , v ) ( 0 ) where v , v { 1 , 2 , , m 2 } , v v . On the other hand, the process { V ( t ) ; t 0 } jumps to the state v with the arrival of l t h environment with probability p ( v , v ) ( l ) where v , v { 1 , 2 , , m 2 } , v v and l { 1 , 2 , , n } .
The behavior of the M M E A P is completely characterized by the matrices D l , l = 0 , 1 , 2 , n defined by
D 0 = λ ( 1 ) λ ( 1 ) p 1 , 2 ( 0 ) λ ( 1 ) p 1 , m 2 ( 0 ) λ ( 2 ) p 2 , 1 ( 0 ) λ ( 2 ) λ ( 2 ) p 2 , m 2 ( 0 ) λ ( m 2 ) p m 2 , 1 ( 0 ) λ ( m 2 ) p m 2 , 2 ( 0 ) λ ( m 2 ) ,
D l = λ ( 1 ) p 1 , 1 ( l ) λ ( 1 ) p 1 , 2 ( l ) λ ( 1 ) p 1 , m 2 ( l ) λ ( 2 ) p 2 , 1 ( l ) λ ( 2 ) p 2 , 2 ( l ) λ ( 2 ) p 2 , m 2 ( l ) λ ( m 2 ) p m 2 , 1 ( l ) λ ( m 2 ) p m 2 , 2 ( l ) λ ( m 2 ) p m 2 , m 2 ( l ) ,
l = { 1 , 2 , , n } The matrix D = l = 0 n D l represents the generator of the process { V ( t ) ; t 0 } .
Service time of those customers who are served with the essential item is phase type distributed with representation P H ( γ , T ) of order m 1 . This service time is the time until the undergoing Markov chain ( ζ ( t ) , t 0 ) with a finite state space { 1 , 2 , 3 , , m 1 + 1 } reaches the absorbing state m 1 + 1 . γ = ( γ 1 , γ 2 , , γ m 1 ) gives the initial probability of starting in any of the m 1 states. T is the generator matrix that gives transition rates within the states { 1 , 2 , 3 , , m 1 }. The absorption rates from the individual transient states { 1 , 2 , 3 , , m 1 } to the absorption state m 1 + 1 is given by T 0 = T e . μ = γ T 1 e gives the mean service of the customer.
Service time of those customers who are served with the optional items are environment dependent and they are exponentially distributed with parameter μ i k , where 1 k n and i ∈ { i 1 , ( i 1 i 2 ) , ( i 1 i 2 i 3 ) , , ( i 1 i 2 i 3 i m ) }. It is important to note that no order preference has been given to any element, i.e., ( i j i l ) = ( i l i j ) , and so on where each i l ∈ { 1 , 2 , 3 , , m } with 1 j l m .
In this model, a customer is allowed to demand exactly one unit of the essential inventory where as more than one type of optional inventory can be demanded by a customer with an imposed restriction of, at most, one item from each optional inventories. Service rates of the optional items are assumed to be environment dependent. The i t h optional item is served in the k t h environment with probability p i k , similarly the l t h and r t h optional inventories are served in the k t h environment with probability p l r k , and so on. If the demanded optional inventory is not available, the customer is expected to quit the system after acquiring the essential item together with those available optional inventories. The server is assumed to be in the idle state in the absence of customers, as well as essential inventories. Essential and optional inventories have exponentially distributed positive lead time with parameters β and β j for 1 j m , respectively. The essential inventories are under the ( s , S ) control policy whereas the environment dependent optional inventories are under the ( s i k , S i ) for i = 1 , 2 , 3 , , m and k = 1 , 2 , 3 , , n control policies in the k t h environment.
At any given time t, let N ( t ) , S ( t ) , E ( t ) , O l ( t ) , J 1 ( t ) , J 2 ( t ) , and J 3 ( t ) denote, respectively, number of customers in the system, status of environment, number of essential inventory items, number of l t h optional inventory items for 1 l m , service phase, environment phase and arrival phase of the customers. The status of the server at any given time t is defined as,
C ( t ) = 0 * Idle server 0 Essential inventory service l l t h optional inventory service l r combined service of l t h and r t h optional inventories l r u combined service of l t h , r t h and u t h optional inventories 123 , . . . , m combined service of all optional inventories
Let Δ be the collection of all the permitted combinations of different optional inventories and let C u denotes the server status, in general, for the combined service of u optional items, for 1 u m . Thus, the process Γ = { ( N ( t ) , S ( t ) , E ( t ) , C ( t ) , O 1 ( t ) , O 2 ( t ) , , O m ( t ) , J 1 ( t ) , J 2 ( t ) , J 3 ( t ) ) , t 0 } is a C T M C which is a level independent quasi-birth and death process ( L I Q B D ) with state space as follows
{ ( 0 , j , i , 0 * , i 1 , i 2 , , i m , j 2 , j 3 ) , 0 i S , 1 j n , 0 i r S r j for 1 r m , 1 j 2 m 2 , 1 j 3 m 3 }
{ ( n ¯ , j , 0 , 0 * , i 1 , i 2 , , i m , j 2 , j 3 ) , n ¯ 1 , 1 j n , 0 i r S r j for 1 r m , 1 j 2 m 2 , 1 j 3 m 3 }
{ ( n ¯ , j , i , 0 , i 1 , i 2 , , i m , j 1 , j 2 , j 3 ) , n ¯ 1 , 1 j n , 1 i S , 0 i r S r j for 1 r m , 1 j 1 m 1 , 1 j 2 m 2 , 1 j 3 m 3 }
{ ( n ¯ , j , i , C 1 , i 1 , i 2 , , i m , j 2 , j 3 ) , n ¯ 1 , 1 j n , 1 i S , C 1 Δ , 1 j 2 m 2 , 1 j 3 m 3 for if C 1 = l where 1 l m then 0 i k S k j for k { 1 , 2 , , m } { l } and 1 i l S l j }
{ ( n ¯ , j , i , C 2 , i 1 , i 2 , , i m , j 2 , j 3 ) , n ¯ 1 , 1 j n , 1 i S , C 2 Δ , 1 j 2 m 2 , 1 j 3 m 3 for if C 2 = l q where l j , 1 l , q m then 0 i k S k j for k { 1 , 2 , , m } { l , q } and 1 i k S k j for k { l , q } }
{ ( n ¯ , j , i , C 3 , i 1 , i 2 , , i m , j 2 , j 3 ) , n ¯ 1 , 1 j n , 1 i S , C 3 Δ , 1 j 2 m 2 , 1 j 3 m 3 for if C 3 = h q l where h q l , 1 h , q , l m then 0 i k S k j for k { 1 , 2 , , m } { h , q , l } and 1 i k S k j for k { h , q , l } }   { ( n ¯ , j , i , 12 m , i 1 , i 2 , , i m , j 2 ) , n ¯ 1 , 1 j n , 1 i S , 12 m Δ , 1 j 2 m 2 , 1 j 3 m 3 , 1 i k S k j for k { 1 , 2 , , m } } .
The infinitesimal generator Q of the system is of the form
Q = A 00 A 01 A 10 A 1 A 0 A 2 A 1 A 0
Matrices A 01 and A 0 are of order a × b and a × a , respectively, their entries are due to the arrival of customers following MAP with representation ( H 0 , H 1 ) . Matrices A 10 and A 2 are of order c × a and c × c , respectively, their entries are due to the service of essential inventories following phase type distribution with representation P H ( γ , T ) and also due to the environment dependent, exponentially distributed service of optional inventories. Matrices A 00 and A 1 are square matrices of order a and c, respectively, their entries includes the replenishment rates of the inventories in addition to the negative sign of sum of other entries of the same row found in A 01 , A 0 , A 10 , and A 2 , where a = n ( S + 1 ) Π k = 1 m ( S k + 1 ) m 2 m 3 , b = n ( S + 1 ) Π k = 1 m ( S k + 1 ) m 1 m 2 m 3 and c = n ( Π k = 1 m ( S k + 1 ) m 2 m 3 + S u Δ l u ) .
When u = 0 , then l 0 = π k = 1 m ( S k + 1 ) m 1 m 2 m 3 .
When u = i , then l i = π k = 1 , k i m ( S k + 1 ) S i m 2 m 3 for 1 i m .
When u = i j , then l i j = π k = 1 , k i , j m ( S k + 1 ) Π k { i , j } S k m 2 m 3 for 1 i , j m
When u = h i j , then l h i j = π k = 1 , k h , i , j m ( S k + 1 ) Π k { h , i , j } S k m 2 m 3 for 1 h , i , j m and so on.
When u = 12 m , then l 12 m = π k = 1 m S k m 2 m 3 .
In order to have a better understanding of the system, a detailed illustration of the model has been provided in Appendix A by fixing the number of optional items m = 2 and the number of environments n = 2 . All the transitions and resultant component matrices are shown clearly in the Appendix A.

2.1. Stability Condition

Let π = ( π 0 , π 1 , π 2 , , π S ) be the steady-state probability vector of A = A 0 + A 1 + A 2 , where π k = ( π k ( 1 ) , π k ( 2 ) , , π k ( n ) ) for 0 k S .
Then
π A = 0 , π e = 1 .
Refer to Appendix A for the component matrix representations From (1),
π 0 N 0 + π 1 M 0 = 0
π i N 1 + π i + 1 M = 0 ,   1 i s
π i N 2 + π i + 1 M = 0 ,   s + 1 i S 1
π 0 Z 0 + i = 1 s π i Z + π S N 2 = 0
where,
N 0 = [ Z 0 i j ] ; M 0 = [ M 0 i j ] for 1 i , j n
N k = [ L i j + Z k i j + M ¯ i j ] ,   1 k 2 for 1 i , j n
M = [ M i j ] ; Z 0 = [ Z 0 i j ] ; Z = [ Z i j ] for 1 i , j n
Solving Equations (2)–(5), we get
π i = π S W 0 i = 0 π S W i 1 i s π S W ^ i s + 1 i S
where,
W 0 = ( 1 ) S ( M N 2 1 ) S s 1 ( M N 1 1 ) s ( M 0 N 0 1 )
W i = ( 1 ) S i ( M N 2 1 ) S s 1 ( M N 1 1 ) s + 1 i
W ^ i = ( 1 ) S i ( M N 2 1 ) S i
The only unknown probability vector π S is obtained from the normalizing condition
π S W 0 + i = 1 s W i + i = s + 1 S W ^ i e = 1 .
Theorem 2.1.
The necessary and sufficient condition for the stability of queuing inventory system under study is
π S H 0 V 0 < π S H 1 V 1 + H 2 V 2
Proof. 
The queueing system with the generator Q under study is stable if, and only if,
π A 0 e < π A 2 e
Refer to Appendix A for the component matrix representations.
Using Equations (6)–(9) together with the matrices A 0 and A 2 we get
π A 0 e = π S i = 1 s W i + i = s + 1 S W ^ i L . e
where L = [ L i j ] for 1 i , j n ,
π A 2 e = π S W 1 ( M 0 + M ¯ ) + i = 2 s W i + i = s + 1 S W ^ i ( M ¯ + M ) . e
where M ¯ = [ M ¯ i j ] for 1 i , j n . Let
H 0 = i = 1 s W i + i = s + 1 S W ^ i , V 0 = L . e , H 1 = W 1 , V 1 = ( M 0 + M ¯ ) . e ,
H 2 = i = 2 s W i + i = s + 1 S W ^ i & V 2 = ( M ¯ + M ) . e
Then, by (10) we get the stated result. □

2.2. Steady State Probability Vector

Let x denote the steady state probability vector of the generator Q . Then, we have
x Q = 0 , x e = 1 .
Partitioning x as x = ( x 0 , x 1 , x 2 , . . . ) , from (23) we get
x 0 A 00 + x 1 A 10 = 0
x 0 A 01 + x 1 A 1 + x 2 A 2 = 0
x n 1 A 0 + x n A 1 + x n + 1 A 2 = 0 ; n 2
By assuming the stability condition, we see that x is obtained as (see [17])
x n = x 1 R n 1 ; n 2 ,
where R is the minimal non-negative solution of the matrix quadratic equation
R 2 A 2 + R A 1 + A 0 = 0
The boundary conditions are given by
x 0 A 00 + x 1 A 10 = 0
x 0 A 01 + x 1 [ A 1 + R A 2 ] = 0
From Equation (23) we get,
x 1 = x 0 K
and by the normalizing condition in ( 23 ) , we get
[ x 0 + x 0 K ( I R ) 1 ] e = 1
where
K = ( A 01 ) ( A 1 + R A 2 ) 1

3. System Performance Measures

  • Expected re-ordering rate of the essential inventory
    E R E = μ r = 1 t = 1 n k = 1 m i k = 1 S k j 1 = 1 m 1 j 2 = 1 m 2 j 3 = 1 m 3 x r ( t , s + 1 , 0 , i 1 , , i m , j 1 , j 2 , j 3 )
  • Expected re-ordering rate of l t h optional item in the a t h environment, 1 l m , 1 a n
    E R O ( l ) a = μ l a r = 1 i = 1 S ( u Δ l u i k = 0 , k u , 1 k m S k i k = 1 , k u , 1 k m S k ) j 2 = 1 m 2 j 3 = 1 m 3 x r ( a , i , u , i 1 , , i h 1 , s l a + 1 , i h + 1 , , j 2 , j 3 )
  • Expected number of customers in the system E C = i = 1 i . x i . e .
  • Expected number of essential inventories in the system
    E E I = t = 1 n i = 1 S k = 1 m i k = 1 S k j 2 = 1 m 2 j 3 = 1 m 3 i . x 0 ( t , i , 0 * , i 1 , , i m , j 2 , j 3 ) + r = 1 t = 1 n i = 1 S k = 1 m i k = 1 S k j 1 = 1 m 1 j 2 = 1 m 2 j 3 = 1 m 3 i . x r ( t , i , 0 , i 1 , , i m , j 1 , j 2 , j 3 ) + r = 1 t = 1 n i = 1 S ( u Λ i k = 0 , k u , 1 k m S k i k = 1 , k u , 1 k m S k ) j 2 = 1 m 2 j 3 = 1 m 3 i . x r ( t , i , u , i 1 , , i m , j 2 , j 3 )
  • Expected number of l t h optional inventories in the system for 1 l m .
    E O I ( l ) = t = 1 n i = 1 S k = 1 m i k = 1 S k j 2 = 1 m 2 j 3 = 1 m 3 i l . x 0 ( t , i , 0 * , i 1 , , i m , j 2 , j 3 ) + r = 1 t = 1 n i = 1 S k = 1 m i k = 1 S k j 1 = 1 m 1 j 2 = 1 m 2 j 3 = 1 m 3 i l . x r ( t , i , 0 , i 1 , , i m , j 1 , j 2 , j 3 ) + r = 1 t = 1 n i = 1 S ( u Λ i k = 0 , k l , k u , 1 k m S k i k = 1 , k = l , k u , 1 k m S k ) j 2 = 1 m 2 j 3 = 1 m 3 i l . x r ( t , i , u , i 1 , , i m , j 2 , j 3 )
  • Expected loss rate of customers in the absence of essential item
    E L = λ A r = 1 t = 1 n k = 1 m i k = 1 S k j 2 = 1 m 2 j 3 = 1 m 3 x r ( t , 0 , 0 * , i 1 , , i m , j 2 , j 3 )

4. Cost Function

In-order to optimize the inventory levels s , S , s i k , and S i for 1 i m , 1 k n , we construct the following cost function,
K ( s , s 1 1 , , s 1 n , , s m 1 , , s m n , S , S 1 , , S m ) = C 0 E R E + i = 1 m j = 1 n C i j E R O ( i ) j + C E I E E I + i = 1 m C O I ( i ) E O I ( i ) + C 1 E C + C 2 E L ,
where
  • C 0 = Fixed set up cost per unit of the essential item;
  • C i j = Fixed set up cost per unit of the i t h optional item in the j t h environment;
  • C E I = Carrying cost per unit of the essential item;
  • C O I ( i ) = Carrying cost per unit of the i t h optional item;
  • C 1 = Customer holding cost per unit time;
  • C 2 = Cost due to loss of goodwill per unit time, in the absence of the essential item.

5. Numerical Illustration

In this section, we provide the numerical illustration of the system performance measures with varied values of the underlying parameters. The model with one essential and two optional inventories is considered here.
The phase-type service process of the customer is characterized by
T = 8 4 1 4 ,   T 0 = 4 3 ,   γ = 0.3 0.7
where the mean service time is obtained as μ = 0.3107 .
To exhibit the correlation effect, we introduce two M A P arrival of customers coded ass M A P p and M A P n for the customer arrivals with positive correlation coefficient and negative correlation coefficient, respectively.
  • M A P p is defined by
    H 0 = 3.45 0.7 0.85 5 , H 1 = 2.3 0.45 0.5 3.65
    The coefficient of correlation C c o r = + 0.13 .
  • M A P n is defined by
    H 0 = 3.35 0.6 0.85 4.105 , H 1 = 2.3 0.45 0.005 3.25
    The coefficient of correlation C c o r = 0.23 .
Table 1 and Table 2 show the effect of μ 1 1 and μ 1 2 , respectively, the environment dependent and exponentially distributed service completion rates of the first optional inventory on different performance measures. From the table, it is clear that the values of the performance measures show similar behavior in both the cases, with M A P p and M A P n , respectively. Values of E R O I 1 1 , E R O I 1 2 , E R O I 2 1 , E R O I 2 2 , E E I , E O I ( 1 ) , E O I ( 2 ) , and E L seems to increase, respectively, with increased values of μ 1 1 and μ 1 2 whereas, the values of E R E and E C seems to decrease, respectively, with the increased values of μ 1 1 and μ 1 2 .
Table 3 and Table 4 show the effect of μ 2 1 and μ 2 2 , respectively, the environment dependent and exponentially distributed service completion rates of the second optional inventory on different performance measures. From the table, it is clear that the values of the performance measures show similar behavior in both the cases, with M A P p and M A P n , respectively. The values of E R E , E R O I 1 1 , E R O I 1 2 , E R O I 2 1 , E R O I 2 2 , E E I , E O I ( 1 ) , E O I ( 2 ) , and E L seems to increase, respectively, with increased values of μ 2 1 and μ 2 2 whereas the value of E C seems a decrease, respectively, with the increased values of μ 2 1 and μ 2 2 .
Table 5 and Table 6 show the effect of μ 12 1 and μ 12 2 , respectively, the environment dependent and exponentially distributed service completion rates of the combined optional inventories on different performance measures. From the table it is clear that the values of the performance measures show similar behavior in both the cases with M A P p and M A P n , respectively. The values of E R E , E R O I 1 1 , E R O I 1 2 , E R O I 2 1 , E R O I 2 2 , E E I , E O I ( 1 ) , E O I ( 2 ) , and E L seems to increase, respectively, with increased values of μ 12 1 and μ 12 2 whereas the value of E C seems to decrease, respectively, with the increased values of μ 12 1 and μ 12 2 .

Numerical Optimization of the Cost Function

Numerical optimization of the cost function constructed in Section 4 is shown in this section. Table 7 and Table 8 show the effect of values of the ( s , S ) control policy with respect to the essential inventory in the cost incurred on the system with M A P p and M A P n , respectively. Both the tables give a concrete idea of the optimality of the variables S and s of the cost function. As expected, it can be observed that the cost incurred increases with the increased value of S , the maximum level of the main inventory. The factors, such as holding cost, contribute to this behavior. The cost is seen minimum when the values S and s are close enough, i.e., when the re-order level is close to the maximum inventory level. Additionally, the cost incurred on the system is seen to increase with the increase in the difference between the values S and s, i.e., when the re-order level is not closer to the maximum inventory level. This shows that it is not profitable to hold more essential inventories in the system at any given time and that the re-ordering level of the essential inventory is to be kept very close to the value of S.
Table 9 and Table 10 show the effect of the pairs ( s 1 1 , S 1 ) and ( s 1 2 , S 1 ) concerning the first optional item in the cost incurred on the system under consideration with M A P p and M A P n , respectively. The table shows that the value of the cost incurred tends to vary in proportional to the value of S 1 . The optimum minimum cost value is seen to be obtained with the decrease in the values of S 1 , s 1 1 , and s 1 2 , respectively.
Table 11 and Table 12 show the effect of the pairs ( s 2 1 , S 2 ) and ( s 2 2 , S 2 ) concerning the second optional inventory in the cost incurred on the system with M A P p and M A P n , respectively. Both the tables gives an idea on the condition for optimal values of the variables S 2 , s 2 1 , and s 2 2 . It can be seen that the cost function show a decreasing tendency with decreased values of S 2 , s 2 1 , and s 2 2 .

6. Conclusions

We studied a single server multi-commodity queueing inventory system with one essential and m optional items in n random environments. The condition for stability of the system is obtained. Under this condition, different performance measures of the system are derived. A cost function involving these measures and inventory control variables is constructed. Optimization of the cost function along with the control variables is also done numerically. The obtained numerical results showed huge resemblances with what we see and experience around us. A very familiar example is a car/truck (Heavy automobiles) showroom, where one can see only very few items (main item) displayed. In this example, one can think of the additional accessories as the optional inventories. Here, booking of the essential inventory with other optional inventories has to be done as per the requirement of the customer.

Author Contributions

Conceptualization, J.J. and A.K.; methodology, J.J. and D.S.; software, D.S.; validation, J.J., D.S., and A.K.; formal analysis, A.K., V.V., and D.K.; writing—original draft preparation, J.J.; writing—review and editing, J.J., D.S., V.V., and D.K.; supervision, A.K.; project administration, A.K.; funding acquisition, V.V. and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported by the RUDN University Strategic Academic Leadership Program (recipient D.K., formal analysis) and funded by Russian Foundation for Basic Research according to the research project number 19-29-06043 (recipients V.V. and D.K., writing—review and editing).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their gratitude to the referees for valuable suggestions that improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Illustration : ( 1 ) : Fix the number of optional items m = 2 and the number of environments n = 2 , the state space of the system is precisely,
{ ( 0 , j , i , 0 * , i 1 , i 2 , j 2 , j 3 ) , 0 i S , 1 j 2 , 0 i 1 S 1 , 0 i 2 S 2 , 1 j 2 m 2 , 1 j 3 m 3 }
{ ( n ¯ , j , 0 , 0 * , i 1 , i 2 , j 2 , j 3 ) , n ¯ 1 , 1 j 2 , 0 i r S r j for 1 r 2 , 1 j 2 m 2 , 1 j 3 m 3 }   { ( n , j , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) , n ¯ 1 , 1 j 2 , 1 i S , 0 i r S r j for 1 r 2 , 1 j 1 m 1 , 1 j 2 m 2 , 1 j 3 m 3 }   { ( n , j , i , 1 , i 1 , i 2 , j 2 , j 3 ) , n ¯ 1 , 1 j 2 , 1 i S , 1 i 1 S 1 j ; 0 i 2 S 2 j , 1 j 2 m 2 , 1 j 3 m 3 }   { ( n , j , i , 2 , i 1 , i 2 , j 2 , j 3 ) , n ¯ 1 , 1 j 2 , 1 i S , 0 i 1 S 1 j ; 1 i 2 S 2 j , 1 j 2 m 2 , 1 j 3 m 3 }   { ( n , j , i , 12 , i 1 , i 2 , j 2 , j 3 ) , n ¯ 1 , 1 j 2 , 1 i S , 1 j 2 m 2 , 1 j 3 m 3 ; 1 i k S k j for k { 1 , 2 } } .
The transitions rates are:
  • Transition rates due to customer arrival.
    (a)
    ( 0 , l 1 , i , 0 * , i 1 , i 2 , j 2 , j 3 ) ( 1 , l 2 , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 )
    at the rate γ j 1 [ H 1 ] j 3 , j 3 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 , 1 j 2 m 2 , and 1 j 3 , j 3 m 3
    where 1 i S ; 0 i r S r for 1 r 2 .
    (b)
    ( n ¯ , l 1 , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) ( n ¯ + 1 , l 2 , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) at the rate [ H 1 ] j 3 , j 3 when l 1 = l 2
    for 1 l 1 , l 2 2 and 1 j 1 m 1 , 1 j 2 m 2 and 1 j 3 , j 3 m 3 where n ¯ 1 ; 1 i S ; 0 i r S r for 1 r 2 .
    (c)
    ( n ¯ , l 1 , i , C u , i 1 , i 2 , j 2 , j 3 ) ( n ¯ + 1 , l 2 , i , C u , i 1 , i 2 , j 2 , j 3 )
    at the rate [ H 1 ] j 3 , j 3 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 2 m 2 , 1 j 3 , j 3 m 3 where n ¯ 1 ; 1 i S and,
    when
    C u = 1 1 i 1 S 1 ; 0 i 2 S 2 ; C u = 2 0 i 1 S 1 ; 1 i 2 S 2 ; C u = 12 1 i 1 S 1 ; 1 i 2 S 2
  • Transition rates due to the service completion of essential and optional items.
    (a)
    Transition rates from level 1 to level 0 due to the service completion of optional inventory
    i.
    ( 1 , l 1 , i , 1 , i 1 , i 2 , j 2 , j 3 ) ( 0 , l 2 , i , 0 * , i 1 1 , i 2 , j 2 , j 3 )
    at the rate μ 1 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 1 i 1 S 1 ; 0 i 2 S 2 .
    ii.
    ( 1 , l 1 , i , 2 , i 1 , i 2 , j 2 , j 3 ) ( 0 , l 2 , i , 0 * , i 1 , i 2 2 , j 2 , j 3 )
    at the rate μ 2 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 0 i 1 S 1 ; 1 i 2 S 2 .
    iii.
    ( 1 , l 1 , i , 12 , i 1 , i 2 , j 2 , j 3 ) ( 0 , l 2 , i , 0 * , i 1 1 , i 2 1 , j 2 , j 3 )
    at the rate μ 12 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 1 i 1 S 1 ; 1 i 2 S 2 .
    (b)
    Transition rates from level 1 to level 0 due to the service completion of essential inventory
    i.
    ( 1 , l 1 , i , 0 , 0 , 0 , j 1 , j 2 , j 3 ) ( 0 , l 2 , i 1 , 0 * , 0 , 0 , j 2 , j 3 )
    at the rate T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where 1 i S
    ii.
    ( 1 , l 1 , i , 0 , 0 , i 2 , j 1 , j 2 , j 3 ) ( 0 , l 2 , i 1 , 0 * , 0 , i 2 , j 2 , j 3 )
    at the rate η 1 T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3 where 1 i S ; 1 i 2 S 2 and η 1 = p + ( 1 p ) ( p 1 + p 12 ) .
    iii.
    ( 1 , l 1 , i , 0 , i 1 , 0 , j 1 , j 2 , j 3 ) ( 0 , l 2 , i 1 , 0 * , i 1 , 0 , j 2 , j 3 )
    at the rate η 2 T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3 where 1 i S ; 1 i 1 S 1 and η 2 = p + ( 1 p ) ( p 2 + p 12 ) .
    iv.
    ( 1 , l 1 , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) ( 0 , l 2 , i 1 , 0 * , i 1 * , i 2 * , j 2 , j 3 )
    at the rate p T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where 1 i S ; 1 i k = i k * S k for k { 1 , 2 } .
    (c)
    Transition rates from level n ¯ to level n ¯ 1 for n ¯ 2 due to the service completion of optional inventory
    i.
    ( n ¯ , l 1 , i , 1 , i 1 , i 2 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i , 1 , i 1 1 , i 2 , j 2 , j 3 )
    at the rate μ 1 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , n ¯ 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 1 i 1 S 1 ; 0 i 2 S 2 .
    ii.
    ( n ¯ , l 1 , i , 2 , i 1 , i 2 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i , 2 , i 1 , i 2 1 , j 2 , j 3 )
    at the rate μ 2 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , n ¯ 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 0 i 1 S 1 ; 1 i 2 S 2 .
    iii.
    ( n ¯ , l 1 , i , 12 , i 1 , i 2 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i , 12 , i 1 1 , i 2 1 , j 2 , j 3 )
    at the rate μ 12 l 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , n ¯ 2 , 1 j 2 m 2 , 1 j 3 m 3
    where 1 i S ; 1 i 1 S 1 ; 1 i 2 S 2 .
    (d)
    Transition rates from level n ¯ to level n ¯ 1 for n ¯ 2 due to the service completion of essential inventory
    i.
    ( n ¯ , l 1 , 1 , 0 , 0 , 0 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , 0 , 0 * , 0 , 0 , j 2 , j 3 )
    at the rate T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 i S .
    ii.
    ( n ¯ , l 1 , 1 , 0 , 0 , i 2 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , 0 , 0 * , 0 , i 2 , j 2 , j 3 )
    at the rate η 1 T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 i 2 S 2
    and η 1 = p + ( 1 p ) ( p 1 + p 12 ) .
    iii.
    ( n ¯ , l 1 , 1 , 0 , i 1 , 0 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , 0 , 0 * , i 1 , 0 , j 2 , j 3 )
    at the rate η 2 T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 , 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 i 1 S 1
    and η 2 = p + ( 1 p ) ( p 2 + p 12 ) .
    iv.
    ( n ¯ , l 1 , 1 , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , 0 , 0 * , i 1 , i 2 , j 2 , j 3 )
    at the rate p T j 1 when l 1 = l 2
    for 1 l 1 , l 2 2 ; 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 i k S k for k { 1 , 2 } .
    v.
    ( n ¯ , l 1 , i , 0 , 0 , 0 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i 1 , 0 , 0 , 0 , j 1 , j 2 , j 3 )
    at the rate T j 1 when l 1 = l 2
    for 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 l 1 , l 2 2 ; 1 i S .
    vi.
    ( n ¯ , l 1 , i , 0 , 0 , i 2 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i 1 , 0 , 0 , i 2 , j 1 , j 2 , j 3 )
    at the rate η 1 T j 1 when l 1 = l 2
    for 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 l 1 , l 2 2 , 2 i S ; 1 i 2 S 2 and η 1 = p + ( 1 p ) ( p 1 + p 12 ) .
    vii.
    ( n ¯ , l 1 , i , 0 , i 1 , 0 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i 1 , 0 , i 1 , 0 , j 1 , j 2 , j 3 )
    at the rate η 2 T j 1 when l 1 = l 2
    for 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3
    where n ¯ 2 ; 1 l 1 , l 2 2 , 2 i S ; 1 i 1 S 1 and η 2 = p + ( 1 p ) ( p 2 + p 12 ) .
    viii.
    ( n ¯ , l 1 , i , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) ( n ¯ 1 , l 2 , i 1 , 0 , i 1 , i 2 , j 1 , j 2 , j 3 ) at the rate p T j 1 when l 1 = l 2
    for 1 j 1 m 1 ; 1 j 2 m 2 ; 1 j 3 m 3 where 1 l 1 , l 2 2 , n ¯ 2 ; 2 i S ; 1 i k S k for k { 1 , 2 } .
The infinitesimal generator Q of the system is obtained to be
Q = A 00 A 01 A 10 A 1 A 0 A 2 A 1 A 0
A 00 is a square matrix of order a. A 01 is a matrix of order a × b . A 10 is a matrix of order c × a . A 0 , A 1 , and A 2 are square matrices of order c, where a = 2 ( S + 1 ) ( S 1 + 1 ) ( S 2 + 1 ) m 2 m 3 , b = 2 ( S + 1 ) ( S 1 + 1 ) ( S 2 + 1 ) m 1 m 2 m 3 and c = 2 ( ( S 1 + 1 ) ( S 2 + 1 ) m 2 m 3 + S ( b 1 + b 2 + b 3 + b 4 ) ) where b 1 = ( S 1 + 1 ) ( S 2 + 1 ) m 1 m 2 m 3 , b 2 = S 1 ( S 2 + 1 ) m 2 m 3 , b 3 = ( S 1 + 1 ) S 2 m 2 m 3 and b 4 = S 1 S 2 m 2 m 3 .
The structure of A 00 , A 01 , A 10 , A 1 , A 0 , A 2 and their corresponding i j t h sub matrices A 00 i j , A 01 i j , A 10 i j , A 1 i j , A 0 i j , A 2 i j for 1 i , j 2 are obtained as given in Equations (A1)–(A12).
The structure of A 00 , A 01 , A 10 , A 0 , A 1 and A 2 are obtained as
A 00 =   1 2   1 2 A 00 11 A 00 12 A 00 21 A 00 22    
where each A 00 i j for 1 i , j 2 has the structure
A 00 i j = 0 1 s s + 1 S 0 1 s s + 1 S Z 0 i j Z ^ i j Z 1 i j ^ Z ^ i j Z 1 ^ i j Z ^ i j Z 2 ^ i j Z 2 ^ i j
A 00 =   1 2   1 2 A 01 11 A 01 12 A 01 21 A 01 22    
where each A 01 i j for 1 i , j 2 has the structure
A 01 i j = 0 1 S 0 1 S 0 ¯ L * i j L * i j
A 00 = 1 2 1 2 A 10 11 A 10 12 A 10 21 A 10 22
where each A 10 i j for 1 i , j 2 has the structure
A 10 i j = 0 1 2 S 0 1 2 S 0 ¯ M 0 i j M ^ i j M 0 i j M ^ i j M 0 i j M ^ i j
A 00 = 1 2 1 2 A 1 11 A 1 12 A 1 21 A 1 22
where each A 1 i j for 1 i , j 2 has the structure
A 1 i j = 0 1 s s + 1 S 0 1 s s + 1 S Z 0 i j Z 0 i j Z 1 i j Z i j Z 1 i j Z i j Z 2 i j Z 2 i j
A 00 = 1 2 1 2 A 0 11 A 0 12 A 0 21 A 0 22
where each A 0 i j for 1 i , j 2 has the structure
A 0 i j = 0 1 S 0 1 S 0 ¯ L i j L i j
A 00 = 1 2 1 2 A 2 11 A 2 12 A 2 21 A 2 22
where each A 2 i i for 1 i , j 2 has the structure
A 2 i j = 0 1 2 S 0 1 2 S 0 ¯ M 0 i j M ¯ i j M i j M ¯ i j M i j M ¯ i j
  • Submatrices of A 00 i j
    For i = j ,
    Z 0 i j = I s 1 i + 1 C 1 e s 1 i + 1 [ I ( S 2 + 1 ) β 1 i I m 2 m 3 ] 0 ¯ I S 1 s 1 i C 2 ,
    where
    C 1 = I s 2 i + 1 B 1 e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 2 ,
    C 2 = I s 2 i + 1 B 3 e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 4 ,
    B 1 = ( D 0 + D i ( β + β 1 i + β 2 i ) I m 2 ) H ,
    B 2 = ( D 0 + D i ( β + β 1 i ) I m 2 ) H ,
    B 3 = ( D 0 + D i ( β + β 2 i ) I m 2 ) H ,
    B 4 = ( D 0 + D i β I m 2 ) H .
    For i j , Z 0 i j = I ( S 1 + 1 ) ( S 2 + 1 ) [ D j I m 3 ]
    For i = j ,
    Z 1 i j ^ = I s 1 i + 1 C 1 ^ e s 1 i + 1 β 1 i I ( S 2 + 1 ) m 2 m 3 0 ¯ I S 1 s 1 i C 2 ^ , where
    C 1 ^ = I s 2 i + 1 B 1 ^ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 2 ^ ,
    C 2 ^ = I s 2 i + 1 B 3 ^ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 4 ^ ,
    B 1 ^ = ( D 0 + D i ( β + β 1 i + β 2 i ) I m 2 ) H 0 ,
    B 2 ^ = ( D 0 + D i ( β + β 1 i ) I m 2 ) H 0 ,
    B 3 ^ = ( D 0 + D i ( β + β 2 i ) I m 2 ) H 0 , B 4 ^ = ( D 0 + D i β I m 2 ) H 0 .
    For i j ,
    Z 1 i j ^ = I ( S 1 + 1 ) ( S 2 + 1 ) [ D j I m 3 ]
    For i = j ,
    Z 2 i j ^ = I s 1 i + 1 C 1 ¯ e s 1 i + 1 β 1 i I ( S 2 + 1 ) m 2 m 3 0 ¯ I S 1 s 1 i C 2 ¯ , where
    C 1 ¯ = I s 2 i + 1 B 1 ¯ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 2 ¯ ,
    C 2 ¯ = I s 2 i + 1 B 3 ¯ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 4 ¯ ,
    B 1 ¯ = ( D 0 + D i ( β 1 i + β 2 i ) I m 2 ) H 0 ,
    B 2 ¯ = ( D 0 + D i ( β 1 i ) I m 2 ) H 0 ,
    B 3 ¯ = ( D 0 + D i ( β 2 i ) I m 2 ) H 0 , B 4 ¯ = D 0 + D i H 0 .
    For i j ,
    Z 2 i j ^ = I ( S 1 + 1 ) ( S 2 + 1 ) [ D j I m 3 ]
    For i = j ,
    Z i j ^ = β I ( S 1 + 1 ) ( S 2 + 1 ) m 2 m 3 .
    For i j , Z i j ^ = 0 ¯
  • Submatrices of A 01 i j
    For i = j ,   L * i j = I ( S 1 + 1 ) ( S 2 + 1 ) [ γ ( I m 2 H 1 ) ]
    For i j , L * i j = 0 ¯
  • Submatrices of A 10 i j
    For i = j , M 0 i j = 0 1 2 S 1 0 1 2 S 1 m 0 m 1 m 1 m 1 , where
    m 0 = T 0 I m 2 m 3 0 ¯ 0 ¯ I S 2 η 1 T 0 I m 2 m 3 , where η 1 = p + ( 1 p ) ( p 1 + p 12 )
    m 1 = η 2 T 0 I m 1 m 2 0 ¯ I S 2 p T 0 I m 1 m 2 , where η 2 = p + ( 1 p ) ( p 2 + p 12 )
    For i j , M 0 i j = 0 ¯
    For i = j ,   M ^ i j = 0 1 2 12 0 1 2 12 0 ¯ μ 1 i I S 1 ( S 2 + 1 ) m 2 m 3 0 ¯ μ 2 i I ( S 1 + 1 ) S 2 m 2 m 3 0 ¯ μ 12 i I S 1 S 2 m 2 m 3 0 ¯ .
    For i j , M ^ i j = 0 ¯
  • Submatrices of A 1 i j
    For i = j ,   Z 0 i j = I ( S 1 + 1 ) ( S 2 + 1 ) [ γ β I m 2 m 3 ] .
    For i j , Z 0 i j = 0 ¯
    For i = j , Z i j = β I [ ( S 1 + 1 ) ( S 2 + 1 ) m 1 m 2 m 3 + S 1 ( S 2 + 1 ) m 2 m 3 + ( S 1 + 1 ) S 2 m 2 m 3 + S 1 S 2 m 2 m 3 ]
    For i j , Z i j = 0 ¯
    Z 1 i j = 0 1 2 12 0 1 2 12 G 1 G 12 G 13 G 14 G 2 G 3 G 4 where
    For i = j ,   G 1 = I s 1 i + 1 C 3 e s 1 i + 1 β 1 i I ( S 2 + 1 ) m 1 m 2 m 3 0 ¯ I S 1 s 1 i C 4 .
    C 3 = I s 2 i + 1 B 5 e s 2 i + 1 β 2 i I m 1 m 2 m 3 0 ¯ I S 2 s 2 i B 6 ,
    C 4 = I s 2 i + 1 B 7 e s 2 i + 1 β 2 i I m 1 m 2 m 3 0 ¯ I S 2 s 2 i B 8 ,
    B 5 = T [ ( D 0 + D i ) ( β + β 1 i + β 2 i ) I m 2 H 0 ] ,
    B 6 = T [ ( D 0 + D i ) ( β + β 1 i ) I m 2 H 0 ] ,
    B 7 = T [ ( D 0 + D i ) ( β + β 2 i ) I m 2 H 0 ] ,
    B 8 = T [ ( D 0 + D i ) β I m 2 H 0 ] .
    For i j ,   G 1 = I ( S 1 + 1 ) ( S 2 + 1 ) m 1 [ D j I m 3 ]
    For i = j , G 12 = 0 ¯ H 1 , H 1 = I S 1 ( S 2 + 1 ) [ ( 1 p ) p 1 T 0 I m 2 m 3 ]
    For i j , G 12 = 0 ¯
    For i = j , G 13 = I S 1 + 1 Z * , where
    Z * = 0 ¯ H 2 , H 2 = I S 2 [ ( 1 p ) p 2 T 0 I m 2 m 3 ] .
    For i j , G 13 = 0 ¯
    For i = j , G 14 = 0 ¯ I S 1 Z ¯ , where
    Z ¯ = 0 ¯ H 3 , H 3 = I S 2 [ ( 1 p ) p 12 T 0 I m 2 m 3 ] ,
    For i j , G 14 = 0 ¯
    For i = j ,   G 2 = I s 1 i C 5 e s 1 i β 1 i I ( S 2 + 1 ) m 2 m 3 0 ¯ I S 1 s 1 i C 6 , where
    C 5 = I s 2 i + 1 B 9 e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 10 ,
    C 6 = I s 2 i + 1 B 11 e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 12
    B 9 = [ ( D 0 + D i ) ( β + β 1 i + β 2 i + μ 1 i ) I m 2 ] H 0 ,
    B 10 = [ ( D 0 + D i ) ( β + β 1 i + μ 1 i ) I m 2 ] H 0 ,
    B 11 = [ ( D 0 + D i ) ( β + β 2 i + μ 1 i ) I m 2 ] H 0 ,
    B 12 = [ ( D 0 + D i ) ( β + μ 1 i ) I m 2 ] H 0 .
    For i j ,   G 2 = I S 1 ( S 2 + 1 ) [ D j I m 3 ]
    For i = j ,   G 3 = I s 1 i + 1 C 7 e s 1 i + 1 β 1 i I S 2 m 2 m 3 0 ¯ I S 1 s 1 i C 8 , where
    C 7 = I s 2 i B 13 e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 14 , C 8 = I s 2 i B 15 e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 16
    B 13 = [ ( D 0 + D i ) ( β + β 1 i + β 2 i + μ 2 i ) I m 2 ] H 0 ,
    B 14 = [ ( D 0 + D i ) ( β + β 1 i + μ 2 i ) I m 2 ] H 0 ,
    B 15 = [ ( D 0 + D i ) ( β + β 2 i + μ 2 i ) I m 2 ] H 0 ,
    B 16 = [ ( D 0 + D i ) ( β + μ 2 i ) I m 2 ] H 0 .
    For i j ,   G 3 = I ( S 1 + 1 ) S 2 [ D j I m 3 ]
    For i = j ,   G 4 = I s 1 i C 9 e s 1 i β 1 i I S 2 m 2 m 3 0 ¯ I S 1 s 1 i C 10 , where
    C 9 = I s 2 i B 17 e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 18 , C 10 = I s 2 i B 19 e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 20
    B 17 = [ ( D 0 + D i ) ( β + β 1 i + β 2 i + μ 12 i ) I m 2 ] H 0 ,
    B 18 = [ ( D 0 + D i ) ( β + β 1 i + μ 12 i ) I m 2 ] H 0 ,
    B 19 = [ ( D 0 + D i ) ( β + β 2 i + μ 12 i ) I m 2 ] H 0 ,
    B 20 = [ ( D 0 + D i ) ( β + μ 12 i ) I m 2 ] H 0 .
    For i j ,   G 4 = I S 1 S 2 [ D j I m 3 ]
    Z 2 i j = 0 1 2 12 0 1 2 12 G 1 ^ G 12 G 13 G 14 G 2 ^ G 3 ^ G 4 ^ ,
    where,
    For i = j ,   G ^ 1 = I s 1 i + 1 C ^ 3 e s 1 i + 1 β 1 i I ( S 2 + 1 ) m 1 m 2 m 3 0 ¯ I S 1 s 1 i C ^ 4 , where
    C 3 ^ = I s 2 i + 1 B 5 ^ e s 2 i + 1 β 2 i I m 1 m 2 m 3 0 ¯ I S 2 s 2 i B 6 ^ ,
    C 4 ^ = I s 2 i + 1 B 7 ^ e s 2 i + 1 β 2 i I m 1 m 2 m 3 0 ¯ I S 2 s 2 i B 8 ^ ,
    B 5 ^ = T [ ( D 0 + D i ) ( β 1 i + β 2 i ) I m 2 H 0 ] ,
    B 6 ^ = T [ ( D 0 + D i ) ( β 1 i ) I m 2 H 0 ] ,
    B 7 ^ = T [ ( D 0 + D i ) ( β 2 i ) I m 2 H 0 ] ,
    B 8 ^ = T [ ( D 0 + D i ) H 0 ] .
    For i j , G 1 ^ = G 1
    For i = j ,   G 2 ^ = I s 1 i C 5 ^ e s 1 i β 1 i I ( S 2 + 1 ) m 2 m 3 0 ¯ I S 1 s 1 i C 6 ^ , where
    C 5 ^ = I s 2 i + 1 B 9 ^ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 10 ^ ,
    C 6 ^ = I s 2 i + 1 B 11 ^ e s 2 i + 1 β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 12 ^
    B 9 ^ = [ ( D 0 + D i ) ( β 1 i + β 2 i + μ 1 i ) I m 2 ] H 0 ,
    B 10 ^ = [ ( D 0 + D i ) ( β 1 i + μ 1 i ) I m 2 ] H 0 ,
    B 11 ^ = [ ( D 0 + D i ) ( β 2 i + μ 1 i ) I m 2 ] H 0 ,
    B 12 ^ = [ ( D 0 + D i ) ( μ 1 i ) I m 2 ] H 0 .
    For i j , G 2 ^ = G 2 .
    For i = j ,   G 3 ^ = I s 1 i + 1 C 7 ^ e s 1 i + 1 β 1 i I S 2 m 2 m 3 0 ¯ I S 1 s 1 i C 8 ^ , where
    C 7 ^ = I s 2 i B 13 ^ e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 14 ^ , C 8 ^ = I s 2 i B 15 ^ e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 16 ^
    B 13 ^ = [ ( D 0 + D i ) ( β 1 i + β 2 i + μ 2 i ) I m 2 ] H 0 ,
    B 14 ^ = [ ( D 0 + D i ) ( β 1 i + μ 2 i ) I m 2 ] H 0 ,
    B 15 ^ = [ ( D 0 + D i ) ( β 2 i + μ 2 i ) I m 2 ] H 0 ,
    B 16 ^ = [ ( D 0 + D i ) ( μ 2 i ) I m 2 ] H 0 .
    For i j , G 3 ^ = G 3
    For i = j ,   G 4 ^ = I s 1 i C 9 ^ e s 1 i β 1 i I S 2 m 2 m 3 0 ¯ I S 1 s 1 i C 10 ^ , where
    C 9 ^ = I s 2 i B 17 ^ e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 18 ^ , C 10 ^ = I s 2 i B 19 ^ e s 2 i β 2 i I m 2 m 3 0 ¯ I S 2 s 2 i B 20 ^
    B 17 ^ = [ ( D 0 + D i ) ( β 1 i + β 2 i + μ 12 i ) I m 2 ] H 0 ,
    B 18 ^ = [ ( D 0 + D i ) ( β 1 i + μ 12 i ) I m 2 ] H 0 ,
    B 19 ^ = [ ( D 0 + D i ) ( β 2 i + μ 12 i ) I m 2 ] H 0 ,
    B 20 ^ = [ ( D 0 + D i ) ( μ 12 i ) I m 2 ] H 0 .
    For i j , G 4 ^ = G 4
  • Submatrices of A 0 i j
    For i = j ,
    L i j = 0 1 2 12 0 1 2 12 I ( S 1 + 1 ) ( S 2 + 1 ) m 1 m 2 H 1 I S 1 ( S 2 + 1 ) m 2 H 1 I ( S 1 + 1 ) S 2 m 2 H 1 I S 1 S 2 m 2 H 1 .
    For i j , L i j = 0 ¯
  • Submatrices of A 2 i j
    For i = j ,   M i j = 0 1 2 S 1 0 1 2 S 1 m 0 ^ m 1 ^ m 1 ^ m 1 ^ , where
    m 0 ^ = d i a g T 0 I m 2 m 3 0 ¯ 0 ¯ I S 2 [ η 1 d i a g T 0 I m 2 m 3 ] ,
    where η 1 = p + ( 1 p ) ( p 1 + p 12 )
    m 1 ^ = η 2 d i a g T 0 I m 1 m 2 0 ¯ 0 ¯ I S 2 p d i a g T 0 I m 1 m 2 ,
    where η 2 = p + ( 1 p ) ( p 2 + p 12 )
    For i j , M i j = 0 ¯
    M ¯ i j = 0 1 2 12 0 1 2 12 0 ¯ μ 1 i I S 1 ( S 2 + 1 ) m 2 m 3 μ 2 i I ( S 1 + 1 ) S 2 m 2 m 3 μ 12 i I S 1 S 2 m 2 m 3 .
    For i j , M ¯ i j = 0 ¯

References

  1. Dudin, A.N.; Klimenok, V.I.; Vishnevsky, V.M. The Theory of Queuing Systems with Correlated Flows; Springer: Cham, Switzerland, 2020; pp. 1–410. [Google Scholar] [CrossRef]
  2. Al-Khayyal, F.; Hwang, S.-J. Inventory constrained maritime routing and scheduling for multi-commodity liquid bulk, Part I: Applications and model. Eur. J. Oper. Res. 2007, 176, 106–130. [Google Scholar] [CrossRef]
  3. Sadjady, H.; Davoudpour, H. Two-echelon, multi-commodity supply chain network design with mode selection, lead-times and inventory costs. Comput. Oper. Res. 2012, 39, 1345–1354. [Google Scholar] [CrossRef]
  4. Qin, J.; Shi, F.; Miao, L.; Tan, G. Optimal Model and Algorithm for Multi-Commodity Logistics Network Design Considering Stochastic Demand and Inventory Control. Syst. Eng. Theory Pract. 2009, 29, 176–183. [Google Scholar] [CrossRef]
  5. Askin, R.G.; Baffo, I.; Xia, M. Multi-commodity warehouse location and distribution planning with inventory consideration. Int. J. Prod. Res. 2014, 52, 1897–1910. [Google Scholar] [CrossRef]
  6. Araya-Sassi, C.; Paredes-Belmar, G.; Gutiérrez-Jarpa, G. Multi-commodity inventory-location problem with two different review inventory control policies and modular stochastic capacity constraints. Comput. Ind. Eng. 2020, 143, 106410. [Google Scholar] [CrossRef]
  7. Sabzevari Zadeh, A.; Sahraeian, R.; Homayouni, S.M. A dynamic multi-commodity inventory and facility location problem in steel supply chain network design. Int. J. Adv. Manuf. Technol. 2014, 70, 1267–1282. [Google Scholar] [CrossRef]
  8. Shajin, D.; Jacob, J.; Krishnamoorthy, A. On a queueing inventory problem with necessary and optional inventories. Ann. Oper. Res. 2021, 1–12. [Google Scholar] [CrossRef]
  9. Song, J.S.; Zipkin, P. Inventory control in a fluctuating demand environment. Oper. Res. 1993, 41, 351–370. [Google Scholar] [CrossRef]
  10. Özekici, S.; Parlar, M. Inventory models with unreliable suppliers in a random environment. Ann. Oper. Res. 1999, 91, 123–136. [Google Scholar] [CrossRef]
  11. Erdem, A.S.; Özekici, S. Inventory models with random yield in a random environment. Int. J. Prod. Econ. 2002, 78, 239–253. [Google Scholar] [CrossRef]
  12. Perry, D.; Posner, M.J.M. Production-inventory models with an unreliable facility operating in a two-state random environment. Probab. Eng. Inf. Sci. 2002, 16, 325–338. [Google Scholar] [CrossRef]
  13. Feldman, R.M. A continuous review (s, S) inventory system in a random environment. J. Appl. Probab. 1978, 15, 654–659. [Google Scholar] [CrossRef]
  14. Kalpakam, S.; Arivariqnan, G. A lost sales inventory system in a random environment. Stoch. Anal. Appl. 1989, 7, 367–385. [Google Scholar] [CrossRef]
  15. Goh, M.; Sharafali, M. Price-dependent inventory models with discount offers at random times. Prod. Oper. Manag. 2002, 11, 139–156. [Google Scholar] [CrossRef]
  16. Krishnamoorthy, A.; Shajin, D.; Viswanath, C.N. Inventory with Positive Service Time:A Survey. In Advanced Trends in Queueing Theory: Series of Books “Mathematics and Statistics”, Sciences; Anisimov, V., Limnios, N., Eds.; ISTE & J. Wiley: London, UK, 2019. [Google Scholar]
  17. Neuts, M.F. Matrix-Geometric Solutions in Stochastic Models: An Algorithmic Approach; The Johns Hopkins University Press: Baltimore, MD, USA, 1981. [Google Scholar]
Table 1. Effects of μ 1 1 (fix μ 1 2 = 4 ) and μ 1 2 (fix μ 1 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 1. Effects of μ 1 1 (fix μ 1 2 = 4 ) and μ 1 2 (fix μ 1 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
M A P p
345 μ 1 2 : μ 1 1 543
4.5 × 10 3 2.9 × 10 3 2.5 × 10 3 E R E 2.9 × 10 3 3.8 × 10 3 6.5 × 10 3
4.4 × 10 3 5.2 × 10 3 5.4 × 10 3 E R O I 1 1 5.2 × 10 3 4.2 × 10 3 2.4 × 10 3
3.8 × 10 3 6.0 × 10 3 7.9 × 10 3 E R O I 1 2 6.0 × 10 3 5.7 × 10 3 4.1 × 10 3
1.3 × 10 3 1.58 × 10 3 1.65 × 10 3 E R O I 2 1 1.6 × 10 3 1.4 × 10 3 9.4 × 10 4
1.9 × 10 3 2.2 × 10 3 2.248 × 10 3 E R O I 2 2 2.2 × 10 3 2.1 × 10 3 1.5 × 10 3
0.8010 0.9202 0.9568 E E I 0.9202 0.8524 0.6213
0.6550 0.5306 0.5220 E C 0.5306 0.5648 1.0821
0.2643 0.2790 0.2837 E O I ( 1 ) 0.2790 0.2691 0.2358
0.3040 0.3261 0.3331 E O I ( 2 ) 0.3261 0.3122 0.2649
3.11 × 10 5 3.76 × 10 5 4.01 × 10 5 E L 3.76 × 10 5 3.45 × 10 5 2.35 × 10 5
Table 2. Effects of μ 1 1 (fix μ 1 2 = 4 ) and μ 1 2 (fix μ 1 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 2. Effects of μ 1 1 (fix μ 1 2 = 4 ) and μ 1 2 (fix μ 1 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
M A P n
345 μ 1 2 : μ 1 1 543
3.8 × 10 3 3.0 × 10 3 2.7 × 10 3 E R E 3.0 × 10 3 3.5 × 10 3 4.7 × 10 3
5.2 × 10 3 5.6 × 10 3 5.7 × 10 3 E R O I 1 1 5.6 × 10 3 4.6 × 10 3 3.5 × 10 3
4.9 × 10 3 7.1 × 10 3 9.2 × 10 3 E R O I 1 2 7.1 × 10 3 6.9 × 10 3 6.4 × 10 3
1.6 × 10 3 1.7 × 10 3 1.8 × 10 3 E R O I 2 1 1.8 × 10 3 1.7 × 10 3 1.4 × 10 3
2.3 × 10 3 2.4 × 10 3 2.5 × 10 3 E R O I 2 2 2.5 × 10 3 2.4 × 10 3 2.29 × 10 2
1.0923 1.1719 1.2028 E E I 1.1719 1.1107 0.9824
0.3668 0.3316 0.3264 E C 0.3316 0.3455 0.4189
0.3078 0.3173 0.3213 E O I ( 1 ) 0.3173 0.3071 0.2886
0.3468 0.3596 0.3649 E O I ( 2 ) 0.3596 0.3470 0.3232
1.58 × 10 5 1.74 × 10 5 1.82 × 10 5 E L 1.74 × 10 5 1.65 × 10 5 1.43 × 10 5
Table 3. Effects of μ 2 1 (fix μ 2 2 = 4 ) and μ 2 2 (fix μ 2 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 3. Effects of μ 2 1 (fix μ 2 2 = 4 ) and μ 2 2 (fix μ 2 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
M A P p
345 μ 2 2 : μ 2 1 543
2.0 × 10 3 2.7 × 10 3 2.9 × 10 3 E R E 2.7 × 10 3 2.3 × 10 3 1.3 × 10 3
3.5 × 10 3 4.6 × 10 3 5.1 × 10 3 E R O I 1 1 4.6 × 10 3 3.8 × 10 3 2.0 × 10 3
4.2 × 10 3 5.7 × 10 3 6.3 × 10 3 E R O I 1 2 5.7 × 10 3 5.1 × 10 3 3.0 × 10 3
8.98 × 10 4 1.2 × 10 3 1.3 × 10 3 E R O I 2 1 1.2 × 10 3 8.28 × 10 4 3.50 × 10 4
1.2 × 10 3 2.0 × 10 3 2.7 × 10 3 E R O I 2 2 2.0 × 10 3 1.7 × 10 3 9.71 × 10 4
0.7109 0.8643 0.9238 E E I 0.8643 0.7737 0.5353
0.8829 0.5570 0.5190 E C 0.5570 0.6612 1.6951
0.2731 0.2943 0.3001 E O I ( 1 ) 0.2898 0.2818 0.2784
0.2993 0.3156 0.3229 E O I ( 2 ) 0.3156 0.3018 0.2681
2.50 × 10 5 3.47 × 10 5 3.91 × 10 5 E L 3.47 × 10 5 2.98 × 10 5 1.66 × 10 5
Table 4. Effects of μ 2 1 (fix μ 2 2 = 4 ) and μ 2 2 (fix μ 2 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 4. Effects of μ 2 1 (fix μ 2 2 = 4 ) and μ 2 2 (fix μ 2 1 = 5 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
M A P n
345 μ 2 2 : μ 2 1 543
2.4 × 10 3 2.8 × 10 3 2.9 × 10 3 E R E 2.8 × 10 3 2.5 × 10 3 1.9 × 10 3
4.4 × 10 3 5.1 × 10 3 5.4 × 10 3 E R O I 1 1 5.1 × 10 3 4.4 × 10 3 3.3 × 10 3
5.9 × 10 3 7.0 × 10 3 7.4 × 10 3 E R O I 1 2 7.0 × 10 3 6.5 × 10 3 5.4 × 10 3
1.2 × 10 3 1.4 × 10 3 1.5 × 10 3 E R O I 2 1 1.4 × 10 3 1.0 × 10 3 6.13 × 10 4
1.6 × 10 3 2.4 × 10 3 3.1 × 10 3 E R O I 2 2 2.4 × 10 3 2.1 × 10 3 1.7 × 10 3
0.9929 1.1144 1.1705 E E I 1.1144 1.0193 0.8552
0.4277 0.3423 0.3291 E C 0.3423 0.3857 0.5572
0.3112 0.3104 0.3130 E O I ( 1 ) 0.3104 0.3043 0.3003
0.3319 0.3477 0.3571 E O I ( 2 ) 0.3477 0.3305 0.3064
1.37 × 10 5 1.64 × 10 5 1.78 × 10 5 E L 1.64 × 10 5 1.47 × 10 5 1.17 × 10 5
Table 5. Effects of μ 12 1 (fix μ 12 2 = 5 ) and μ 12 2 (fix μ 12 1 = 4 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 5. Effects of μ 12 1 (fix μ 12 2 = 5 ) and μ 12 2 (fix μ 12 1 = 4 ): Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
M A P p
345 μ 12 2 : μ 12 1 543
6.4 × 10 4 2.9 × 10 3 3.5 × 10 3 E R E 3.9 × 10 3 3.5 × 10 3 2.9 × 10 3
1.1 × 10 3 5.1 × 10 3 6.1 × 10 3 E R O I 1 1 6.6 × 10 3 6.1 × 10 3 5.2 × 10 3
1.4 × 10 3 6.0 × 10 3 7.3 × 10 3 E R O I 1 2 8.2 × 10 3 7.3 × 10 3 6.0 × 10 3
3.4 × 10 4 1.5 × 10 3 1.8 × 10 3 E R O I 2 1 1.8 × 10 3 1.78 × 10 3 1.6 × 10 3
4.93 × 10 4 2.2 × 10 3 2.6 × 10 3 E R O I 2 2 2.9 × 10 3 2.6 × 10 3 2.2 × 10 3
0.4330 0.9077 1.027 E E I 1.097 1.027 0.9202
3.2931 0.5618 0.4106 E C 0.3745 0.4106 0.5306
0.1393 0.278 0.3115 E O I ( 1 ) 0.3337 0.3115 0.279
0.2273 0.3280 0.3461 E O I ( 2 ) 0.3615 0.3461 0.3261
6.77 × 10 6 3.28 × 10 5 4.18 × 10 5 E L 4.36 × 10 5 4.18 × 10 5 3.76 × 10 5
Table 6. Effects of μ 12 1 (fix μ 12 2 = 5 ) and μ 12 2 (fix μ 12 1 = 4 ) with M A P p and M A P n : Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
Table 6. Effects of μ 12 1 (fix μ 12 2 = 5 ) and μ 12 2 (fix μ 12 1 = 4 ) with M A P p and M A P n : Fix S = 4 , S 1 = 3 , S 2 = 3 , s = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 2 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 .
345 μ 12 2 : μ 12 1 543
1.9 × 10 3 2.9 × 10 3 3.3 × 10 3 E R E 3.5 × 10 3 3.3 × 10 3 3.0 × 10 3
3.7 × 10 3 5.5 × 10 3 6.1 × 10 3 E R O I 1 1 6.4 × 10 3 6.1 × 10 3 5.6 × 10 3
4.8 × 10 3 7.2 × 10 3 8.1 × 10 3 E R O I 1 2 8.8 × 10 3 8.1 × 10 3 7.1 × 10 3
1.1 × 10 3 1.7 × 10 3 1.8 × 10 3 E R O I 2 1 1.8 × 10 3 1.78 × 10 3 1.75 × 10 3
1.7 × 10 3 2.5 × 10 3 2.8 × 10 3 E R O I 2 2 3.1 × 10 3 2.8 × 10 3 2.5 × 10 3
0.888 1.1598 1.2702 E E I 1.3336 1.2702 1.1719
0.7408 0.3387 0.2854 E C 0.2708 0.2854 0.3316
0.2511 0.3174 0.3444 E O I ( 1 ) 0.3624 0.3444 0.3173
0.3187 0.3608 0.3800 E O I ( 2 ) 0.3949 0.3800 0.3596
9.72 × 10 6 1.56 × 10 5 1.82 × 10 5 E L 1.84 × 10 5 1.82 × 10 5 1.74 × 10 5
Table 7. Effect of ( s , S ) with M A P p on the Cost function: Fix S 1 = 2 , S 2 = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 7. Effect of ( s , S ) with M A P p on the Cost function: Fix S 1 = 2 , S 2 = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P p
S s 4321
8545.2348631.7666905.00761003.001
7403.4487711.095985.87751181.701
6292.7348435.4059677.8474914.2366
5295.0661297.7624612.2427933.1728
Table 8. Effect of ( s , S ) with M A P n on the Cost function: Fix S 1 = 2 , S 2 = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 8. Effect of ( s , S ) with M A P n on the Cost function: Fix S 1 = 2 , S 2 = 2 , s 1 1 = 1 , s 1 2 = 1 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P n
S s 4321
8633.9178667.4360841.9043879.9268
7547.5167709.7668899.4244995.1082
6406.4411492.9354647.1244786.6579
5398.9954397.9574586.8929781.3571
Table 9. Effect of ( s 1 1 , S 1 ) (fix s 1 2 = 1 ) and ( s 1 2 , S 1 ) (fix s 1 1 = 1 ) with M A P p in the cost incurred on the system: Fix S = 3 , s = 1 , S 2 = 2 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 9. Effect of ( s 1 1 , S 1 ) (fix s 1 2 = 1 ) and ( s 1 2 , S 1 ) (fix s 1 1 = 1 ) with M A P p in the cost incurred on the system: Fix S = 3 , s = 1 , S 2 = 2 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P p
123 s 1 2 : S 1 : s 1 1 321
935.5729894.1549863.96716922.4228935.3796935.5729
865.7215872.3495874.67485883.7009877.5753865.7215
809.9514824.6832809.56364796.8894827.6646809.9514
Table 10. Effect of ( s 1 1 , S 1 ) (fix s 1 2 = 1 ) and ( s 1 2 , S 1 ) (fix s 1 1 = 1 ) with M A P n in the cost incurred on the system: Fix S = 4 , S 2 = 2 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 10. Effect of ( s 1 1 , S 1 ) (fix s 1 2 = 1 ) and ( s 1 2 , S 1 ) (fix s 1 1 = 1 ) with M A P n in the cost incurred on the system: Fix S = 4 , S 2 = 2 , s 2 1 = 1 , s 2 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P n
123 s 1 2 : S 1 : s 1 1 321
509.8432494.74524.19616545.1854510.3812509.8432
524.6889520.1184501.43855562.9512543.2073524.6889
467.4061459.2036467.52964480.5882469.6073467.4061
Table 11. Effect of ( s 2 1 , S 2 ) (fix s 2 2 = 1 ) and ( s 2 2 , S 2 ) (fix s 2 1 = 1 ) with M A P p in the cost incurred on the system: Fix S = 3 , s = 1 , S 1 = 2 , s 1 1 = 1 , s 1 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 11. Effect of ( s 2 1 , S 2 ) (fix s 2 2 = 1 ) and ( s 2 2 , S 2 ) (fix s 2 1 = 1 ) with M A P p in the cost incurred on the system: Fix S = 3 , s = 1 , S 1 = 2 , s 1 1 = 1 , s 1 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P p
123 s 2 2 : S 2 : s 2 1 321
979.5430986.6257975.65256970.0710986.1526979.5430
885.1778854.4908851.08995840.7033844.198885.1778
692.9451696.7333693.51594699.5264701.7346692.9451
Table 12. Effect of ( s 2 1 , S 2 ) (fix s 2 2 = 1 ) and ( s 2 2 , S 2 ) (fix s 2 1 = 1 ) with M A P n in the cost incurred on the system: Fix S = 3 , S 1 = 2 , s 1 1 = 1 , s 1 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
Table 12. Effect of ( s 2 1 , S 2 ) (fix s 2 2 = 1 ) and ( s 2 2 , S 2 ) (fix s 2 1 = 1 ) with M A P n in the cost incurred on the system: Fix S = 3 , S 1 = 2 , s 1 1 = 1 , s 1 2 = 1 , m 1 = m 2 = m 3 = 2 , μ 1 1 = 5 , μ 1 2 = 4 , μ 2 1 = 6 , μ 2 2 = 4 , μ 12 1 = 3 , μ 12 2 = 5 , β = 6 , β 1 1 = 4 , β 1 2 = 5 , β 2 1 = 4 , β 2 2 = 3 , C 0 = 120 $ , C 1 1 = 35 $ , C 1 2 = 40 $ , C 2 1 = 45 $ , C 2 2 = 50 $ , C E I = 110 $ , C O I ( 1 ) = 500 $ , C O I ( 2 ) = 450 $ , C 1 = 70 $ , C 2 = 50 $ .
M A P n
123 s 2 2 : S 2 : s 2 1 321
850.5189853.9868843.16516842.5212854.2780850.5189
776.4571755.0671753.36825747.1692747.7253776.4571
631.4084634.4203634.02254640.1808639.8542631.4084
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jacob, J.; Shajin, D.; Krishnamoorthy, A.; Vishnevsky, V.; Kozyrev, D. Queueing-Inventory with One Essential and m Optional Items with Environment Change Process Forming Correlated Renewal Process (MEP). Mathematics 2022, 10, 104. https://doi.org/10.3390/math10010104

AMA Style

Jacob J, Shajin D, Krishnamoorthy A, Vishnevsky V, Kozyrev D. Queueing-Inventory with One Essential and m Optional Items with Environment Change Process Forming Correlated Renewal Process (MEP). Mathematics. 2022; 10(1):104. https://doi.org/10.3390/math10010104

Chicago/Turabian Style

Jacob, Jaison, Dhanya Shajin, Achyutha Krishnamoorthy, Vladimir Vishnevsky, and Dmitry Kozyrev. 2022. "Queueing-Inventory with One Essential and m Optional Items with Environment Change Process Forming Correlated Renewal Process (MEP)" Mathematics 10, no. 1: 104. https://doi.org/10.3390/math10010104

APA Style

Jacob, J., Shajin, D., Krishnamoorthy, A., Vishnevsky, V., & Kozyrev, D. (2022). Queueing-Inventory with One Essential and m Optional Items with Environment Change Process Forming Correlated Renewal Process (MEP). Mathematics, 10(1), 104. https://doi.org/10.3390/math10010104

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop