A MARKOVIAN TWO COMMODITY QUEUEING-INVENTORY SYSTEM WITH COMPLIMENT ITEM AND CLASSICAL RETRIAL FACILITY Текст научной статьи по специальности «Медицинские технологии»

URAL MATHEMATICAL JOURNAL, Vol. 8, No. 1, 2022, pp. 90-116

DOI: 10.15826/umj.2022.1.009

A MARKOVIAN TWO COMMODITY QUEUEING-INVENTORY SYSTEM WITH COMPLIMENT ITEM AND CLASSICAL RETRIAL FACILITY

M. Nithyat, C. Sugapriyat, S. Selvakumartt, K. Jeganathantt'1, T. Harikrishnanttt

tDepartment of Mathematics, Queen Mary's College, Tamil Nadu, Chennai - 600004, India

ttRamanujan Institute for Advanced Study in Mathematics, University of Madras, Chepauk, Tamil Nadu, Chennai - 600005, India

tttDepartment of Mathematics, Guru Nanak College (Autonomous), Tamil Nadu, Chennai - 600042, India

1kjeganathan84@gmail.com

Abstract: This paper explores the two-commodity (TC) inventory system in which commodities are classified as major and complementary items. The system allows a customer who has purchased a free product to conduct Bernoulli trials at will. Under the Bernoulli schedule, any entering customer will quickly enter an orbit of infinite capability during the stock-out time of the major item. The arrival of a retrial customer in the system follows a classical retrial policy. These two products' re-ordering process occurs under the (s, Q) and instantaneous ordering policies for the major and complimentary items, respectively. A comprehensive analysis of the retrial queue, including the system's stability and the steady-state distribution of the retrial queue with the stock levels of two commodities, is carried out. The various system operations are measured under the stability condition. Finally, numerical evidence has shown the benefits of the proposed model under different random situations.

Keywords: Markov process, Compliment item, Infinite orbit, Waiting time.

1. Introduction

In an inventory business, substitution schemes play a key role in reducing customer losses. During the stock-out period of a demanded item, the flexible item can be used. For example, any product from dairy inventories has stochastic demand, and some of them can be used as substitutions for others. As sales and demand are directly proportional to the complement items, it is convenient to calculate the relevant product's future sales and demand. Besides, the business plans to launch a new product as a supplement to another product that will help the new product understand consumers' needs. Also, the company initially offers a supplement that eventually encourages a customer to make a purchase and continue to purchase the same product without a compliment in the future due to the market's product rivalry. Everyone in a mobile store, supermarket, car showroom, electrical equipment shop, etc., will experience these kinds of situations. We are motivated by these products to evaluate a TC inventory system, including complimentary items. The different varieties of this complimentary technique shall be defined as follows:

• Discounting is the crucial factor used to market tomorrow's cash flow and how the seller uses pricing to accomplish a particular market goal. More precisely, it deals with the customer's psychological reaction to those types of prices. For example, a textile showroom offers a 20%

discount or 30% discount. Online shopping platforms like Amazon, Flipkart and Myntra, etc., offer a discount of 20% on all domestic products, etc.

• The company increases the quantity of product to a certain percentage as an offer rather than giving a discount in price. For example, the Colgate company increases its quantity by 20% for the same price to attract more customers. A detergent, soap, or powder company increases a considerable percentage to sell the product effectively.

• Similarly, few firms provide buy one get one offer (may be the same item or a different item) to increase sales and profit. For example, buy one shirt and get the same shirt free, a mobile phone with a memory card, a laptop with a pen drive, etc.

So far, several researchers have studied two product inventory structures with a finite orbit. Since multi-commodity is more challenging than the single commodity system, most researchers have studied the single commodity queuing inventory model. For more details about a single commodity, one may refer [6, 11-15, 17, 20, 21, 26].

Nowadays, the modeling of multi-product inventory systems receives more attention. Many businesses and firms have gradually begun to use multi-commodity inventory systems in modern computer technologies. Kalpakam and Arivarignan [16] analyzed a joint reordering policy of a multi-item inventory such that the replenishment duration of every new order is zero. In [10], Goyal and Satir explored inventory control models where a range of items are jointly replenished. Under deterministic and stochastic demand conditions, the approaches available for determining the economic operating strategy for jointly replenished items are reviewed. Sivazlian [27] looks at the stationary properties of a multi-commodity inventory analysis periodically. The stochastic model assumes a dyadic replenishment strategy with proportional costs and a single set-up cost. The optimal operating costs for individual and joint ordering policies were compared in numerical scenarios.

Yadavalli et al. [36] assumed a Markovian arrival process (MAP) for a TC continuous review inventory system with three categories of arriving customers. Sivakumar et al. [28] suggested a TC perishable stochastic inventory method under continuous study at a service facility with a finite waiting hall. TC's are supposed to be interchangeable. That is, if either of the stocks is empty instantly, the other commodity may be used to meet the demand. Sivakumar et al. [29] assumed an inventory system with a loss in sales where each commodity's lifetime and lead time of a joint reorder of two commodities are all independent exponential distributions. Serife Ozkar and Umay Uzunoglu Kocerciteser investigated a TC queueing-inventory model with an individual ordering policy in which there are two types of customers: priority (Type-1) and ordinary (Type-2). Customers of Type-1 request commodity-1 only, while customers of Type-2 demand commodity-2. Each customer's arrival is subjected to an independent Poisson process of varying rates. Senthil Kumar [25] studied a TC inventory system that was subjected to discrete-time review, with each commodity's demand determined by an independent Bernoulli process. Anbazhagan and Arivarig-nan [2] elaborately studied independent Poisson demand processes of the TC continuous review inventory system. They made coordinated reorders whenever each commodity's stock level is less than or equal to its reorder level. TC inventory schemes were studied under separate ordering policies by Anbazhagan and Arivarignan [3, 4].

Sivakumar [31] initiated a retrial policy on the TC inventory system where the demand for any commodity is substitutable with others if the demanded commodity is not available. However, Yadavalli et al. [33] made the ordering quantity of each up to its maximum stock level. A TC inventory system with a single server was considered by Binitha Benny et al. [8]. The buffer capacity in this paper will be finite. Customers arrive in a Poisson process, and the demand for each type of commodity, or both types of commodities, is defined using specific probabilities.

Krishnamoorthy and Merlymol Joseph [19] discussed a continuous review of TC inventory problems with bulk demand. Assume that the model has both the commodities' probability of a demand equal to zero. Instantaneous replenishment and no shortages are permitted. Sivakumar et al. [30] studied a TC continuous review inventory system with a renewal of demand and ordering policy, a combination of policies referred to as the ordering of individual commodities and the ordering of both commodities jointly. Krishnamoorthy et al. [18] dealt with a TC inventory system with zero lead time in which the Poisson arrival of any customer is satisfied either with one or both commodities under a prefixed probability distribution. Yadavalli et al. [34] studied two commodity inventory systems in which the customers' demand patterns for each commodity follow Poisson and the demand for each commodity is fulfilled with another commodity with distinct probabilities.

Artaljeo et al. [7] began researching a classical retrial strategy in an inventory system. The reality of the classical nature of retrial policy in an inventory system has also been determined by Ushakumari [32]. A classical retrial queue with a single server with phase-type service facilities was addressed by Krishnamoorthy and Dhanya Shajin [22]. A classical retrial policy on an (s,Q) inventory system with a finite customer source in which multi-servers provide homogeneous services was extensively studied by Yadavalli et al. [35]. Srinivas R. Chakravarthy et al. [9] studied a multiserver, infinite-orbit retrial method that considers a classic customer retrial pattern.

Anbazhagan and Jeganathan [5] addressed a two-commodity system with a complimenting item, where primary demand for the first commodity enters a finite orbit size N. Lakshmanan et al. [23] examine a two-commodity situation with a complement and frequent working holidays. Both commodities are independent of their ordering policies, and each customer orders service at a positive time. When the requested item is out of stock or the server is busy, each consumer is allowed to a given finite retrial orbit.

From the above studies, it is noticed that there is no work on TC inventory systems with complement items and an infinite orbit under the classical retrial policy. We consider that this concept is a gap in the inventory system until now. To fill such a gap in this field, we propose the model as a two-commodity inventory system involving an infinite orbit in which the orbital customers approach the system under the classical retrial policy. Sections 2, 3, 4, 5, 6, and 7 of the paper presented a model description, system analysis, waiting time analysis, measures of different system performances, cost analysis, and numerical illustration and conclusion, respectively.

1.1. Notation

We will use the following notation. Let the symbol 0 denotes the matrix with zero entries, let e be a convenient-sized column vector with one in each of the co-ordinates, a In be an nth order identity matrix. Let

dij :

1 if j = i,

0, otherwise,

5ij := 1 — 5,

H]

H(x) : =

[Bn]a,b : =

i ]i

1 if x > 0, 0, otherwise;

1 if a = 2, ■ ■ ■

n,

b = a — 1,

[CS2 ]a,b :

H:

0, otherwise;

1 if a = S2, b = S2,

0, otherwise;

= {1 < u < L, 0 < v < Si, 1 < w < S2}.

Perishable items for commodity-2 Perishable items for commodity-1

A

A

Y2

Commodity-1 available

No Yes

(1 - r)A

rA

Commodity-1(S1 ) Commodity-2(S2)

(1 - p) A

-<r

Loss

Decide to orbit

No Yes

pA

->•( Infinite orbit

Figure 1. Graphical illustration of the model

2. Explanation of model

A system of two commodities is considered, in which one commodity is labelled as a main commodity and the other one as a complementary product. Any customer's timely arrival follows an exponential distribution at a rate of A; S1 and S2 refer to the maximum stock level of main and complimentary products, respectively. Any order quantity of the main commodity follows (s,Q) whereas a (0, S2) ordering policy is used for the complimentary product. The inter-arrival times between successive reorders of the main commodity follow an exponential distribution with rate fi. However, the system receives any reorder of the complimentary product instantly when it is stocked out. Also, the main commodity and complimentary product are both independently perishable with rates and jy2 (1 < i < S1,1 < j < S2), respectively, and their deteriorating times follow exponential distributions where i and j denote their corresponding stock levels at that time. If there is no main commodity in the stock, any arrival customer enters into an infinite orbit at the rate of A with probability p. According to a purchase study, any new customer demands either the main commodity only or both commodities, at the rate of A with probability 1 — r or r, respectively (0 < r < 1). According to a classical retrial policy, any retrial customer demands the main commodity only or both commodities at the rate of k0 with probability 1 — r or r where k (k > 0) is the size of the retrial queue at that time. The proposed model of the two commodity queueing-inventory system is given in Fig. 1.

3. Analysis of the system

Consider a triplet (A(t),B(t),C(t)), where A(t),B(t), and C(t) denote the size of orbital customers, first commodity level, and second commodity level, respectively. According to the Markov property and the assumptions of the given model, the continuous-time discrete-state random pro-

cess

X(t) = {(A(t),B(t),C(t)), t > 0} is said to be a Markov chain and its state space D is defined as

D = {(u,v,w)| u = 0,1, 2,...; v = 0,1, 2,...,Si; w = 1, 2,.

3.1. Construction of infinitesimal generator matrix

The rate matrix of stochastic queueing inventory system (SQIS) X(t) is given by

U =

where

/ Uoo Uoi 0 0 0 ... o 0 0 0 0 •• \

Uio Uii Uoi 0 0 ... o 0 0 0 0 ••

0 U2o U2i Uoi 0 ... 0 0 0 0 0 ••

0 0 0 0 0 ••• Uno U N i Uoi 0 0 ••

0 0 0 0 0 ... 0 U(N +i)o U(N+i)i Uoi 0 ••

V . /

Pitt f pA, v' = v, v = 0, w' = w, w = 1, 2, S2;

[Uo i] (vw)(vlwl) 0, otherwise

(3.1)

(3.2)

for u = 1, 2, 3,...,

[U«o](vw)(v/w/) =

for u = 0, 1, 2,

ur0, v' = v — 1, v = 1,2,

ur0 v' = v — 1, v = 1,2,

u(1 — r)0, v' = v — 1, v = 1,2,

0, otherwise

, Si, w' = S2, w = 1; , Si, w' = w — 1, w = 2,3, ...,S2; , Si, w' = w, w = 1,2, ...,S2;

[UM«](vw)(v/w/) =

& v' = v + Q, v = 0, 1, ... , s,

w' = w, w = 1, 2, .. S2;

72, v' = v, v = 0, 1, 2, ..., Si,

w' = S2, w = 1;

wY2, v' = v, v = 0, 1, 2, ..., Si,

w' = w — 1, w = 2, 3, .. ., S2 ;

rA, v' = v — 1, v = 1, 2, ... , Si ,

w' = S2, w = 1;

rA, v' = v — 1, v = 1, 2, ... , Si ,

w' = w — 1, w = 2, 3, .. ., S2 ;

vyi + (1 — r)A, v' = v — 1, v = 1, 2, ... , Si

w' = w, w = 1, 2, .. ., S2 ;

—(pA + vyi + wy2 + u6> + ,0), v' = v, v = 0,

w' = w, w = 1, 2, .. ., S2 ;

—(A + vyi + wY2 + u6> + ,0), v' = v, v = 1, 2, ... , s,

w' = w, w = 1, 2, .. ., S2 ;

—(A + vyi + wY2 + u0), v' = v, v = s + 1,s + 2,...,Si,

w' = w, w = 1, 2, .. ., S2 ;

(3.3)

0,

otherwise.

Explanation of the above matrix structure. By the assumption of the proposed model, the submatrices in (3.1) are square matrices of order (S1 + 1)S2. First, consider the matrix U01. It contains the submatrix pAIs2 whose entries are nothing but the transition rate of arrival A enter into the orbit under the Bernoulli schedule as follows:

(u, 0,w) —V (u + 1,0,w), u = 0,1,2,..., w = 1,2, ...,S2, we get equation (3.2).

Next let us discuss the matrix Uu0, where u = 1,2,.... This matrix is equal to [B(Sl+1)(s2)]A, where A is the submatrix of dimension S2, diagonal entries of A are the transition rate of retrial arrivals demand for the main commodity d with probability (1 — r) and the remaining positive entries of A defined as [BS2 + CS2]0 with probability r means transition rate of retrial arrivals demand for both commodity, each retrials based on the classical retrial policy enter into getting service, it gives equation (3.1). Then transition rates are as follows:

(u,v, 1) -V (u — 1, v — 1, S2), u = 1,2,..., v = 1,2,..., S1, (u, v,w) —V (u — 1,v — 1,w — 1), u = 1,2,..., v = 0,1,2,...S1, w = 2, ...,S2,

u(1—r)d

(u, v,w) -> (u — 1,v — 1,w), u = 1,2,..., v = 1,2, ...S1, w = 1,2,...,S2.

Then consider the diagonal matrix Uuu, where u = 0,1,2,..., whose elements are of the transition rates as follows:

(1) ft denotes the rate of reorder transition which follows the (s, Q) ordering policy,

(u, v,w) —V (u, v + Q,w), u = 0,1,2,..., v = 0,1,2,...,s, w = 1,2, ...,S2;

(2) A is the rate of arrival transition enter into the service for demanding both commodity with probability r,

(u, v, 1) —V (u, v — 1, S2), u = 0,1,2,..., v = 1,2,..., S1;

(u, v,w) —V (u, v — 1,w — 1), u = 0,1,2,..., v = 1,2,...,S1, w = 2, ...,S2;

(3) A is the rate of arrival transition enter into the service for demanding the main commodity with probability (1 — r),

(u, v,w) (1 r)A> (u, v — 1,w), u = 0,1,2,..., v = 1,2, ...S1, w = 1,2, ...,S2;

(4) y1 indicates the perishable transition rates for the first commodity which depends on the number of present inventory level for the first commodity,

(u, v,w) —V (u, v — 1,w), u = 0,1,2,..., v = 1,2,...S1, w = 1,2, ...,S2;

(5) y2 indicates the perishable transition rates for the second commodity which depends on the number of present inventory level for the second commodity,

(u, v,w) —V (u, v,w — 1), u = 0,1,2,..., v = 0,1,2, ...,S1, w = 2,...,S2; (u, v, 1) (u,v,S2), u = 0,1,2,..., v = 0,1,2,..., S1;

(6) then the diagonal element is filled by the sum of all the entries in the corresponding rows with the negative sign to satisfy the sum of all entries in each row yield zero; we obtain equation (3.3). Hence all the submatrices obtained through the respective transitions give the infinitesimal generator matrix U as in equation (3.1).

3.2. Matrix geometric approximation

In this section, we find the steady-state probability vector $ and the system's stability condition.

3.2.1. Steady state analysis

Consider K to be the truncation process's cutoff point for the matrix-geometric approximation. To find the steady-state of the considered system using Neuts-Rao truncation method, we assume that Uuo = UKo and Uui = UKi for all u > K. The truncated system X(t)'s modified generator matrix is

/ Uoo U01 0 0 0 ■ ■ 0 0 0 0 0 ■ " \

U10 U11 U01 0 0 ■ ■ 0 0 0 0 0 ■

0 U20 U21 U01 0 ■ ■ 0 0 0 0 0 ■

0 0 0 0 0 ■ ■ UK0 Uk 1 U01 0 0 ■

0 0 0 0 0 ■ ■ 0 Uk 0 Uk1 U01 0 ■

V . ..

Theorem 1. The steady-state probability vector $ corresponds to the generator matrix UK, where UK = UK0 + UK1 + Uoi is given by

$(v) = $(Q)ev, v = 0,1,..., S1,

(3.4)

where

f (-1)Q-vFqEQ-iFQ-I ...Fv+iE-1,

ev =

S-v

(-1)2Q-v+^ [(FqE--iFq-1 ...Fs+i-v'E--1v/)x

v'=0

GE--v/(FS-v/E--v,-1FS-v/-1 . . . Fv+1E,

v-1)

v = 0,1,...,Q - 1, v = Q,

v = Q + 1, Q + 2,...,S1,

and $(Q) is obtained by solving

$(Q)

s-1

(-1)Q X] [(FQEQ-1FQ-1 . . . Fs+1-v/Es-v/)GE--v/(FS-v/E-- v/-1FS-v/-1 . . . Fv+1Ev 1)]

v'=0

Fq+1 + EQ + (-1)qFQE--1FQ-1 ... F1Eo-1G

and

Proof. We have where

Si

]T$(v)e = 1.

v'= 1

$UK = 0 and $e = 1,

Ev, v' = v,

v = 0,1, 2,... ,51;

[UK ]vv/ =

Fv, v' = v - 1, v = 1, 2,..., 51; G, v' = v + Q, v = 0,1,2,..., s; 0, otherwise.

0

The first equation of the above framework yields the following set of equations:

+ $v Ev = 0, v = 0,- 1, + $v Ev + $v-QG = 0, v = + 1,...,S1 - 1, (3.5)

$vEv + $v-QG = 0, v = S1.

We get equation (3.4) by recursively solving the set of equations (3.5) and using the normalising condition. □

Next, the stability condition in which the framework is stable is then determined. Theorem 2. The system's stability condition at the truncation point K is given by

r1pAe < r2K0e,

where

S2 Si S2

r1 = ^ $(°>w), r2 = ¿ ¿ .

w=1 v=1w=1

Proof. From the well known-result of Neuts [24] on the positive recurrence of Uk, we have

)U°1e < $(K)UK°e,

and, by exploiting the structure of the matrices U01 and UK°, we get, for v = 0,1,2,..., S1 and w = 1,2,... ,S2,

(v,w)U01e < (v,w)UK°e.

First,

(0), (1),..., (S1)]U°1e < (0), (1),... , (S1)]UK°e, where (i) = (v,w) and

(0)A/s2 , (1)0,..., (S)0]e < (0)A, (1)A,..., (S)A]e.

The left-hand side becomes

(0)A/s2 = (0, w)pA. On the other hand, due to the structure of A, the right-hand side becomes

(v)A = (v, w), (v, w),..., (v, w)]K0.

Therefore, the last inequality becomes

S2 Si S2

^ (0, w)pAe < ^ ^ (v, w)K0e.

w=1 v=1w=1

Hence,

r1pAe < r2K0e,

where

S2 Si S2 r1 = ^ $(°>w), r2 = ¿¿ $(v'w),

w=1 v=1w=1

as desired. □

3.3. Stationary probability vector

The regularity of the Markov process X(t) with the state space D can be seen from the structure of the rate matrix U and Theorem 2. Henceforth, the limiting probability distribution defined as

X(ra) = lim Pr [A(t) = u, B(t) = v, C(t) = w | A(0),B(0),C(0)1

t—y^o

exists and is independent of the initial state. Let x = (x(o),x(i),.. .) satisfy

XU = 0, xe = 1.

We can partition the vector x(u) as

x(u) = (x(u'o), x(u>i),..., X(u'Sl)) , u > 0,

and

X(u'v) = (xKv,1),XKv'2),...,XKv'S2)) , u > 0, 0 < v < S1.

3.3.1. Computation of the matrix R

Theorem 3. Utilizing the vector x and the specific structure of U, R can be determined by

R2Uko + RUki + Uoi = 0, where R is the minimal nonnegative solution of the matrix quadratic equation (MNSMQE).

P r o o f. Since the Markov process is a regular, the stationary probability distribution exists and is given by

xU = 0, xe = 1.

In order to express the solution in a recursive form, we assume that

x(u) = x(K)Ru, u > K, where the spectrum of R is less than 1, which is ensured by the stability condition. Then, we get

x(m)(R2Uxo + RUki + Uoi) = 0, m = K,K + 1,K + 2,... . Since the above equation is true for all m = K, K + 1, N + 2,..., we get

(R2Uko + RUki + Uoi) =0, m = K, K + 1, K + 2,... . Then R is the MNSMQE, and let us assume that the matrix R is of the form

R =

( R00 R01 00 00

R0Si 0

0 0

This matrix R has only S2 nonzero rows of dimension (Si + 1)(S2). The structure of the block matrix R0v/, where v' € {0,1, ■ ■ ■ , S1}, is of the form

Rqv —

''Qv' /21 ''Qv' /31

1 Qv'

/12

1 Qv' l22 1 Qv' 32 1 Qv'

/13

1 Qv' 23 1 Qv' 33 1 Qv'

114 1 Qv' 124 1 Qv' 134 1 Qv'

115 1 Qv' 25 1 Qv' 35 1 Qv'

1 1 S21 2 1 S23 1^24 1 S2 5 ;

\ 1 Qv' 1 Qv' 1 Qv' 1 Qv' 1 Qv' • • • 1 Qv' /

11S2 \ 1 Qv' \ 12S2 1 Qv' 13S2 1 Qv'

This is also a square matrix of dimension S2. Now, exploiting the coefficient matrices UKo,1, and U01 with R2 and R equal to 0, we obtain a system of S2-dimensional equations as follows (for v = 0):

• for v' = 0,1,2,... ,Si - 1, w = 1,2,..., S2, and w' = 1,2,... ,S2 - 1,

( E CC+1) K (1 - r)0 + E CCÎv^+i^Kr^ + CvT 'Cif + Ci/" /+1)(w' +1)72

v=1

v=1

v(v.

+1 W(WW''+1) ((v' + 1)71 + (1 - r)A) + 1 W((vT+1)1)rA + ¿v'Q <WA + H (s - v') 1W(WW'-s) ß) — 0; for v' — 0,1,2, • • • ,S1 - 1, w — 1,2, • • • , S2, and w' — S2,

( E Ov7v''+1) K(1 - r)0 + E /WvX/vX(v'+1) Kr* + CV 'Civ0 + Cv1'72 + C(Wm) ((v' + 1)Y1

v=1 v=1

+(1 - r)A) + C(1v'+1) rA + ¿v'Q^ii'PA + H (s - v') qv-s) ß) — 0; for v' — S1, w — 1,2, • • • , S2, and w' — 1,2, • • • , S2 - 1,

( 1 Vi' oiP + 1 wir+1)(w' + 1)72 + Cswft) — 0; for v' — S1, w — 1,2, • • • , S2, and w' — S2,

(1vv' Cw' ) + 1vv'72 + 1 vs — 0-

After solving all such equations, one can obtain the elements of the matrix R. In this case,

C(v/) are the diagonal elements of the v'th diagonal submatrix of UK1. Theorem 4. The vector x can be determined by

□

(i+K-1) — x(K-1)r., i > 0,

X

due to the special structure of U, the fact that R is the MNSMQE

R2UkQ + RUK1 + UQ1 — 0,

and the vector x(i), i > 0,

K

X(i) —

ct!(q) n UjQ(-Uj-1), 0 < i < K - 1, j=i+1

ff!(Q) R(i-K), i > K,

(3.6)

where

k-i K

a = [1 + X(o) ^ [] Ujo(—Uj-i)e]-i, (3.7)

i=o j=i+i

and X(0) can be computed by using the normalising condition

X(o) (I — R)-ie = 1.

Proof. The subvector x(o),X(i),... ,X(K-i) and the block partitioned matrix of U give the set of equations (1 < i < K — 1)

X(o)U

oo + X(i)Uio = 0,

X(i-i)Uoi + X(i)Uii + X(i+i)U(i+i)o = 0, (3.8)

X(K-2)Uoi + X(K-i)(U(K-i)i + RUko) = 0.

Using (3.8) repeatedly, we find

X(o) = X(i)Uio(—Uo)-i

and

X(i) = X(2)U2o(—Ui)-i,

where

Ui = (Uii + Uio(—Uo)-iUoi)Uo = Uoo.

Next,

X(2) = X(3)Uso(—U2)-i,

where

U2 = (U2i + U2o(—Ui)-iUoi). On continuing this procedure up to K — 1 times, we get

X(i) = X(i+i)U(i+i)o(—Ui)-i, 0 < i < K — 1, (3.9)

where

U i Uio, i = 0,

Ui (Uii — Uio(—Ui-i)-iUoi), 1 < i < K.

For the next, we use the block Gaussian elimination method to find the vectors (X(K), X(K+i), X(K+2) ...). The nonboundary states subvector (x(k), X(K+i), X(K+2),...) satisfies the relation

(X(K),X(K+i),X(K+2)...)

( Uk Uoi 0 0 0 ■ ■ ■ \

Uko Uki Uoi 0 0 •••

0 Uko Uk i Uoi 0 •••

v . . . . .... J

(3.10)

Let us assume that

a ^ X(i)e, X(i) = a-iX(K+i), i > 0.

i=K

0

From (3.10), we get

X(K >U + x(K+1> U K0 = 0, x(K+i> = x(K+ 1>R, i > 1,

which implies that

X(0>UK + X(1>UK0 = 0 X(i> = X(i-1>R, i > 1,

that is,

X (0> [Uk + RUko] = 0. (3.11)

Since

we have

i>.

Ex (i>e = 1, i=0

X(o) (I — R)-ie = 1. (3.12)

As a result,

X(o)

is the only solution to equations (3.11) and (3.12). Hence,

X(i) = aX(o)R(i-K), i > K. (3.13)

Again, by (3.9) and (3.13), we get (3.6). Using

ro

EX(i) e = 1

i=o

and (3.6), we get

K-i K ro

aX(o)E II Ujo(—Uj-i)e + aX(o) ^ R(i-K)e = 1, i=o j=i+i K

which gives a as in (3.7). □

4. Waiting time analysis

Waiting time (WT) is the time interval between an epoch when a demand approaches the orbit and the moment when its time of operation completion occurs. Using the Laplace-Stieltjes transform (LST), we look at the WT of demand in orbit. To find the orbital demand's waiting period, we naturally limit the orbit size to a finite size. The continuous random variable Wo represents the waiting time distribution of an orbit customer.

4.1. WT of orbital customers

Theorem 5. The probability that an orbital demand will not wait in the orbit is determined as follows:

P {Wo = 0} = 1 — no, (4.1)

where

no = E E X(

L-1 S2

(u,0,w>

u=1w=1

Proof. Since the zero and positive waiting time probability sum is 1, we have

P{Wo = 0} + P{Wo > 0} = 1. (4.2)

Clearly, the probability of positive WT of orbital demand can be determined as

L-1 S2

P{Wo > 0} = E Ex(U'°'W)- (4.3)

u=1w=1

Equation (4.3) can be found easily using Theorem 4. Substituting it into equation (4.2), we get the stated result as desired in (4.1). □

To enable the distribution of Wo, we define some complimentary variables. Suppose that the queueing inventory system is at state (u, v,w), u > 0 at an arbitrary time t, and

(1) Wo(u, v, w) is the time until chosen demand becomes satisfied;

(2) the LST of Wo(u,v,w) is *Wo(u, v, w)(y) and we denote Wo by *Wo(y);

(3) *Wo(y) = E[eyWo] is the LST of unconditional waiting time (UWT);

(4) *Wo(u,v,w)(y) = E[eyWo(u'v'w)] is the LST of conditional waiting time (CWT).

Theorem 6. The LST

{*Wo(u,v,w)(y), (u,v,w) € H*, where H* = H U{*}}

satisfies the system

Zo(y)*Wo(y) = -0e(u, v, w), (u, v, w) € H, (4.4)

Zo(y) = (P - yi),

the matrix P is derived from U by deleting the state (0, v, w), 0 < v < S1, 1 < w < S2, {*} is the absorbing state, and the absorption appears if the orbital demand finds the positive commodities.

Proof. To analyse the CWT, we apply the first step analysis as follows:

for and

Next, for and

*Wo(u,0,w)(y) = ^ *Wo(u + l,0,w)(y)+6wl^ *Wo(u,0,S2)(y) *Wo(u,0,w - 1 )(y) + ^ *W0(u,Q,w)(y)

1 < u < L, v = 0, 1 < w < S2 a = (y + pA + ¿w1w72 + ¿w1w72 + ft). 1 < u < L, 1 < v < S1, 1 < w < S2

(4.5)

b=( y+^w1rA+(1-r)A + H (s-v)ft+v71+5w1w72+5w1w72+^w1(u-1)r0+5w1(u-1)r^,

we get

*Wo(u,v,w)(y) = *Wo(u,v - 1 ,S2)(y) + (1 br)X *Wo(u,v - 1 ,w)(y)

+ H{'\ 11)13 *W0(u,v + Q,w)(y) + ^*W0(u,v-l,w)(y)+ôwl^ *W0(u,v,S2)(y)

b b b (4.6)

*Wo(u,v,w-l)(y) + ôwl{U~b1)r0 *W0(u-l,v-l,S2)(y)

+ôwi(U~b1)rd *W0(u -l,v- l,w - 1 )(y) +

From equations (4.5) and (4.6), we attain a coefficient matrix of the unknowns as a block tridiagonal, which yields the stated result as in (4.4). □

Theorem 7. The nth moments of conditional waiting time is given by

dn+i dra+i

ZoM-jjrTT *wo(y) - (" + D^TFT *Wo(V) = 0

and

dn+i

P r o o f. Using linear equations obtained in Theorem 6, we get a recursive algorithm for finding a conditional and unconditional waiting times. Now, differentiating equations (4.5) and (4.6) (n + 1) times and setting y = 0, we obtain

for and

Next, for and

E[WZ+l(u,0,w)] = ^E[W:+l(u + l,0,w)]+öwl^E[W:+l(u,0,S2)]

+öwi^E[W:+1(u,0,w - 1)] + ^E[W?+1(u,Q,w)]

1 < u < L, v = 0, 1 < w < S2 a = (y + pA + ¿w1W72 + ¿w1W72 + ß). 1 < u < L, 1 < v < S1, 1 < w < S2

(4.7)

b=(y+&w1rA+(1-r)A+H (s-v)ß+v71+5w1W72+5w1W72+^w1(u-1)r0+5w1(u-1)r^,

we get

S[W0n+1(M,v,iy)] = - 1, S'a)] + {1~br)XE[WZ+1(u,v - l,w)}

+ H{S~V)ßE[W:+\u,v + Q,w)\ + ^E[W:+1(u,v - 1 ,w)]+öw1^E[W:+1(u,v,S2)}

b b b (4 8)

- 1)] + ~bi)r6E[W:+1(u -l,v- 1, S2)]

+4i ^ ~b)rd E[W>l+l{u -l,v- l,w - 1)] + 9-.

With reference to equations (4.7) and (4.8), one can determine the unknowns E[W™+1(u, v, w, x)] in terms of moments of one order less. Setting n = 0, we obtain the desired moments of particular order in an algorithmic way. □

Theorem 8. The LST of UWT of orbital demand is given by

* W0(y) = 1 - no + no * W0(u + 1,v, w)(y). (4.9)

Proof. Using Poisson arrival see time averages (PASTA) property, one can obtain the LST of Wo as follows:

*Wo(y) = x(i) *Wo(u,v,w)(y), 0 < u < L, 0 < v < S1, 0 < w < S2. (4.10)

Using the expressions (4.10), we get the stated result. Considering the Euler and Post-Widder algorithms in [1] for the numerical inversion of (4.9), we obtain the desired result. □

Theorem 9. The nth moment of UWT, by the above theorem, is given by

L—1 Si S2

E[W0n ] = ¿0„ + (1 - ¿0n )££ E X(u'v'w)E[Won(u + 1, v, w)]. (4.11)

rn

u=0 v=0 w=1

Proof. To determine the moments of Wo, we differentiate the equation in Theorem 8 n times and calculate at y = 0 to obtain the desired result, which gives the nth moment of UWT in terms of the CWT of the same order. □

Theorem 10. The expected waiting time of an orbital demand is defined by

L—1 Si S2

E[Wo] = ¿x(M'v,w)E[W0(u + 1,v,w)]. (4.12)

u=° v=0 w=1

Proof. Using equation (4.11) in Theorem 9 and substituting n = 1, we get the desired result as in (4.12). □

5. Measures of various performances of the system

In this section, the following measures of corresponding performance are used to obtain the expected total cost under the steady state probability vector.

1. EIC1 denotes the expected level of commodity 1. Using the steady state probability vector x, we define EIC1 as

ro S1 S2

eic 1 = ££ £vx(

u=0 v=1 w=1

2. EIC2 denotes the expected level of commodity 2 (the compliment item). Using the steady state probability vector x, we define EIC2 by

œ S1 S2

EIC £wx(

u=0 v=0 w=1

3. ERC1 denotes the expected reorder rate of commodity 1. Reorder for Q items is placed whenever the system reaches to s from s + 1 under the (s,Q) reordering policy. Therefore, ERC1 is given by

ro S2

ERC1 = £ £ M + (s + 1)Y1 + A] x(u's+1'w).

u=° w=1

4. ERC2 denotes the expected reorder rate of commodity 2. As instantaneous reordering policy is considered for the commodity 2, S2 items are replenished immediately whenever system drops from 1. Then ERC2 is defined by

ro S1 ro

ERC2 = £ £ [urtf + Y2 + rA] x(u'v>1) + £ Y2X(u'0,1).

u=°v=1 M=°

5. ECRO denotes the expected number of customers in the orbit. Therefore, it is defined by

ro S1 S2

ECRO = £ £ £ ux(u'v'w).

u=1v=° w=1

6. EORC denotes the overall rate of retrial customers. The customer from the orbit can try to buy the product in the system irrespective of the product's availability. Then we have

ro S1 S2

EORC =

u=1v=° w=1

7. ESRC denotes the successful rate of retrial customers. Whenever the orbit customer finds that there exists a positive commodity 1, then their retrial process will be successful. It is given by

ro S1 S2

ESRC =

u=1v=1 w=1

8. EPC1 denotes the expected number of perishable commodity 1. Due to the life time, commodity 1 can be perishable at anytime. The mean number of perishable commodity 1 is defined as

ro S1 S2

EPC 1 = ££ £vyiX(

«=0 v=1 w=1

9. The expected number of perishable commodity 2 denoted by EPC2 is given by

ro S1 S2

EPC2 = ££ £wY2X(u'v'w).

«=0 v=0 w=1

10. When the customer finds that commodity 1 is empty, they leave the system with a probability of (1 — p). Therefore, an expected customer lost in the system is defined as an ECL at any time by

œ S2

(u,0,w)

ECL = ££(1 - p)Ax(u'°'w).

«=0 w=1

6. Cost analysis and numerical illustration

Here, we discuss the feasibility of a proposed model through the system characteristics and sufficient economic illustrations. The expected total cost (ETC) is given by

ETC (S1, S 2) = Ch1EIC 1 + Ch2EIC 2 + C^ERC 1 + C^ERC 2 + CP1EPC1 +CP2EPC 2 + Cw E[W0] + CiECL.

To compute the ETC per unit time, the following costs are considered.

Ch1: Carrying cost of commodity 1/unit item.

Ch2: Carrying cost of commodity 2/unit item.

Cs1: Ordering cost of commodity 1/order.

Cs2: Ordering cost of commodity 2/order.

Cp1: Perishable cost of commodity 1/unit item.

Cp2: Perishable cost of commodity 2/unit item.

Cw: Waiting cost of an orbiting customer/unit customer.

Ci: Cost of a customer lost/unit customer.

6.1. Numerical illustration

Numerical analysis is an applied mathematical technique that allows a staggeringly large amount of data to be processed and analyzed for trends, thereby aiding in forming conclusions. To do such numerical illustrations, we fix the cost values as

Ch1 = 0.45, Ch2 = 0.15, CS1 = 20.5, CS2 = 5.5, Cp1 = 2, Cp2 = 0.8, Cw = 2.6, Ci = 6.6,

and the parameters are A = 1.6, p = 0.94, r = 0.92, ft = 1.23, and d = 0.65.

Case (i). Under the following assumptions on various costs and parameters, we illustrate the optimality and convexity of the cost function with independent ordering quantities of both commodities according to the following given range of S1 and S2 (Table 1). Let S1 = 15, 16, 17, 18, 19 and S2 = 2,3,4, 5,6. Using Table 1, under the range of S1, we obtain the minimum expected total cost for each S2 which is noted in bold script. Similarly, under the given range of S2, the minimum expected total cost for each S1 is noted with an underline. From these discussions, the least optimum expected total cost exists at S1 = 17 and S2 = 4, which is given in bold script with an underline (Table 1). Since the proposed model holds the convex properties on the expected

Table 1. ETC rate as a function of Si and S2

S\S'2 2 3 4 5 6

15 7.52903 7.12110 7.08464 7.11847 7.18228

16 7.34905 7.02337 6.99638 7.03800 7.10035

17 7.27243 7.00295 6.98645 7.02999 7.09535

18 7.28788 7.02966 7.02417 7.07014 7.13718

19 7.30836 7.09246 7.09167 7.14019 7.20857

total cost under the variation of pair of parameters (S1, S2), it can be apply to the real-life product sales business. This model gives the optimal total cost of the entire system for the fixed set of parameters. In a business, one can run it successfully if the entire business process is balanced. That

X

Figure 2. E[Wo] vs A and ß

X

Figure 3. E[Wo] vs A and yi

is, a balanced business means that it will maintain a good relationship between the customer and the system operator. To run a successful business, the customer-owner relationship is important, but at the same time, our business does not fall down. This aspect is determined by the optimal results of the system in a business. In such a way, the proposed model will give the assurance of providing an efficient business.

Case (ii). In this study, various parameters influencing the customers lost and the expected total cost rate (from Figs. 2-7) are discussed.

1. Under a Bernoulli's schedule, the rate of customers entering into the orbit increases as A increases, and so both the expected waiting time and expected total cost also increase (Fig. 2-Fig. 7).

2. According to any given increment of mean replenishment time of a commodity 1 (1/ft), expected waiting time and expected total cost are found to be more sensitive at the higher value of A. Nevertheless, we also notice that both measures are poorly sensitive at the optimum value of A (Fig. 5-Fig. 8).

3. According to any given increment of mean spoilage time of a commodity 1, the rate of the expected waiting time is not significant at any value of A (Fig. 4-Fig. 7).

4. Also d makes a significant effect on expected total cost, and further it is highly significant at the higher value of A, expected waiting time increases as A increases (Fig. 4 and Fig. 7).

According to this analysis, we can relate this numerical to real-life phenomena. In a business, the system manager will pay attention to controlling the average waiting time of a customer and the total cost of the system through the balanced mechanism. Here, the probability p plays such a role in the system. Then the reorder process of a business also plays an important role in reducing the actual waiting time of a customer and the expected total cost of the system. Whenever the system manager controls the average replenishment time, which does not exceed its limit as much as possible, this analysis can be made easier by the case (ii). Similarly, the arrival rate influences those system metrics.

Case (iii). The choice of a customer's demand is for either a single commodity (commodity 1) or both (commodity 1 and commodity 2). Hence, all the measures that are relevant to commodity 1 do not make any significant differences. But the rate of choice of demand for both commodities makes a significant change in the measures relevant to commodity 2. Fig. 8 to Fig. 13 show that some significant effects on expected total cost with the expected reorder rate of commodity 2 are noticed according to the independent decision to make the amount of commodity 1 and commodity 2.

1. Under Bernoulli's theory, if the probability of demand for both commodities (r) increases, then the expected reorder rate of commodity 2 and the expected total cost increase (Fig. 8 and Fig. 13).

2. Also, we determine that each measure of this study depends linearly on both commodities demand.

3. When comparing S1, s, and S2, the value of S1 is more sensitive on expected total cost with expected reorder rate of commodity 2 (Fig. 8-Fig. 10).

4. Also we notice that both the measures are highly sensitive at the higher value of s and the lower value of S1 and S2.

5. From this study, it is clear that both commodities' ordered quantities have a significant effect when the demand rate for both commodities increases.

The purchase of both commodities by a customer is decided by the probability r and its complement. This concept will be very helpful to analyze how many customers can buy both products and single products in the inventory sales business. Apart from that, an efficient businessman must have enough knowledge about the available stocking positions in the system, then only they can make a plan to place a reorder of the required items. If the requirement for both commodities increases, the system manager will pay attention to observing the reorder point of the commodity and the available stocks. This analysis will be the required analysis to have such knowledge.

X

Figure 4. E[Wo ] vs A and d

A

Figure 5. ETC vs A and ft

Case (iv). Suppose there is no stock in the system, the preference of a customer enters into the orbit under a Bernoulli's schedule. Using Tables 2-4, we determine the following merits of the proposed model with various measures according to the decision for making the number of orders/production.

1. For every cycle, if the production or order quantity of commodity 1 increases, the expected number of customers in orbit, the expected reorder rate of commodity 1, and the expected number of customers lost decrease significantly.

2. Also, we determine that the expected inventory of commodity 1 and its expected perishable quantity will significantly increase. But we can obtain the minimum expected total cost with a suitable order quantity.

a

Figure 6. ETC vs A and y1

a

Figure 7. ETC vs A and 9

3. Further, if p increases, then the rate of customers entering the orbit increases and the expected number of customers lost decreases. If p = 1, then the system cannot find any customers lost.

Also, we notice that the expected number of customers in orbit, the expected reorder rate, the expected inventory level of commodity 1, and its expected perishable quantity are all increasing at a slower rate due to the classical retrial policy.

Table 2. Response of Si vs p on various measures of the system

Si P EIC1 ERC 1 EPC1 ECRO ECL ETC

0.88 0.74143 0.16624 0.15570 1.20715 0.16461 9.97351

11 0.94 0.75841 0.16587 0.15926 1.28318 0.08190 9.64499

1.00 0.77503 0.16554 0.16275 1.35851 0.00000 9.32038

0.88 1.50279 0.00958 0.31558 1.20451 0.16425 7.27306

17 0.94 1.53739 0.01081 0.32285 1.27996 0.08169 6.98229

1.00 1.57109 0.01207 0.32992 1.35463 0.00000 6.69543

0.88 2.33520 0.00081 0.49039 1.20436 0.16423 7.79526

23 0.94 2.39262 0.00100 0.50245 1.27974 0.08168 7.49968

1.00 2.44854 0.00123 0.51419 1.35435 0.00000 7.20727

Table 3. Response of s vs p on various measures of the system

S P EIC 1 ERC 1 EPC 1 ECRO ECL ETC

0.88 1.77499 0.00136 0.37274 1.20443 0.16424 7.32728

3 0.94 1.81612 0.00166 0.38138 1.27984 0.08169 7.02110

1.00 1.85604 0.00199 0.38977 1.35449 0.00000 6.71855

0.88 1.50279 0.00958 0.31558 1.20451 0.16425 7.27306

5 0.94 1.53739 0.01081 0.32285 1.27996 0.08170 6.98229

1.00 1.57109 0.01207 0.32992 1.35463 0.00000 6.69543

0.88 1.26500 0.07046 0.26565 1.20469 0.16427 8.34847

7 0.94 1.29626 0.07324 0.27221 1.28018 0.08171 8.09013

1.00 1.32691 0.07585 0.27865 1.35492 0.00000 7.83216

Table 4. Response of S2 vs p on various measures of the system

P EIC2 ERC2 EPC2 ECRO ECL ETC

0.88 1.35030 0.23295 0.43209 1.20451 0.16425 7.54911

3 0.94 1.35122 0.23706 0.43239 1.27996 0.08170 7.27871

1.00 1.35215 0.24126 0.43268 1.35463 0.00000 7.01223

0.88 1.98000 0.13627 0.63360 1.20451 0.16425 7.27306

5 0.94 1.98326 0.13651 0.63464 1.27996 0.08170 6.98229

1.00 1.98652 0.13683 0.63568 1.35463 0.00000 6.69543

0.88 2.55762 0.11346 0.81844 1.20451 0.16425 7.38211

7 0.94 2.56362 0.11326 0.82036 1.27996 0.08170 7.09008

1.00 2.56954 0.11311 0.82225 1.35463 0.00000 6.80177

Case (v). Suppose that any corresponding cost rates of both commodities may be changed, then we notice some significant effects of the proposed model according to the given parameters' values (from Table 5).

1. For any holding cost and perishable cost, the expected total cost rate decreases as the setup cost decreases.

Table 5. Expected Total cost with different combinations of various cost

Cpi Cp2 Chi Chi Cs i — 5.5 Cs i = 20.5

Cs 2 — 5.5 Cs2 = 20.5 Cs 2 — 5.5 Cs 2 = 20.5

0.8 0.8 0.15 0.15 5.97149 8.01918 6.13365 8.18134

0.45 6.56647 8.61416 6.72863 8.77632

0.45 0.15 6.43271 8.48040 6.59486 8.64255

0.45 7.02769 9.07538 7.18984 9.23754

2.0 0.15 0.60 6.73306 8.78075 6.89522 8.94291

0.80 7.32804 9.37573 7.49020 9.53789

0.45 1.00 7.19428 9.24197 7.35644 8.64255

2.00 7.78926 9.83695 7.95142 9.99911

2.0 0.8 0.15 0.15 6.35891 8.40660 6.52107 8.56876

0.45 6.95389 9.00158 7.11605 9.16374

0.45 0.15 6.82013 8.86782 6.98229 9.02998

0.45 7.41511 9.46280 7.57727 9.62496

2.0 0.15 0.60 7.12049 9.16818 7.28264 9.33034

0.80 7.71547 9.76316 7.87762 9.92532

0.45 1.00 7.58170 9.62940 7.74386 9.79155

2.00 8.17668 10.22438 8.33884 10.38653

2. On the other hand, the expected total cost rate is maximum if the setup costs are equal and maximal.

3. If holding costs and perishability costs of commodity 2 increase, then the expected cost rate increases significantly.

4. The expected cost rate is highly sensitive when all costs increase.

7. Conclusion

We investigate the Markovian TC inventory system with a classical retrial facility. The system allows a customer who has purchased a free product to conduct Bernoulli trials at will. Also, it is assumed that the retrial process follows a classical retrial policy and an (s,Q) ordering policy for replenishment. The system's stability is derived through the matrix-geometric approximation; indeed, the optimum total cost is computed along with the stationary probability vector. Further, through the numerical illustration, the merits of the proposed model are explained. This model helps the producers to cope with the market as well as the public to achieve better success in the sale of their products and earn profits. Evaluating the proposed model in an economy is one of the most crucial decision-making variables that a business must analyze to survive and grow in a competitive market.

Acknowledgements

The authors would like to sincerely thank two anonymous reviewers for their valuable comments and constructive suggestions, which are very helpful in improving the presentation of this paper.

0.165

C\l

0.16

o

E

E 0.155

o

O

o

■ffi 0.15

ra

ce

r

0 ■g 0.145

o

0

ce

T3 0.14

0

tj

0

a. x 0.135

LU

0.13

—v- — S1 =11 .....

■■ S1 =17 .....■ V

- A S1 =23 ..... -V..... ......

--V

.-V

-V"

-A-

■ A -

0 -A- ■

■ -A-

■ A-

-A ■

.44

.52

.68 .76

r

.84

.92

Figure 8. ERC2 vs Si and r

0.146

CM

0.144

o E

E 0.142

o

O

o

■ffi ra ce

0.14

a 0.138

o 0 ce

0.136

o 0

£ 0.134T " LU

0.132

A- s=3 - — s=5 -■V— s=7

-V"

,-V

,-V

-V

O.....

..........«

..........0.....

.......o

.................

¿Y-

____A-

.--A----AT"

____A

____A--- ■

.44

.52

.68 .76

r

.84

.92 1

Figure 9. ERC2 vs s and r

REFERENCES

1. Abate J., Whitt W. Numerical inversion of Laplace transforms of probability distributions. ORSA J. Comput., 1995. Vol. 7, No. 1. P. 36-43. DOI: 10.1287/ijoc.7.1.36

2. Anbazhagan N., Arivarignan G. Two commodity continuous review inventory system with coordinated reorder policy. Int. J. Inf. Manag. Set., 2000. Vol. 11, No. 3. P. 19-30.

3. Anbazhagan N., Arivarignan G. Two commodity inventory system with individual and joint ordering policies. Int. J. Manag. Syst., 2003. Vol. 19, No. 2. P. 129-144.

4. Anbazhagan N., Arivarignan G. Analysis of two commodity Markovian inventory system with lead time. Korean J. Comput. & Appl. Math., 2001. Vol. 8, No. 2. Art. no. 427. P. 427-438. DOI: 10.1007/BF02941976

5. Anbazhagan N., Jeganathan K. Two-commodity Markovian inventory system with compliment and retrial demand. J. Adv. Math. Comput. Set., 2013. Vol. 3, No. 2. P. 115-134. DOI: 10.9734/BJMCS/2013/1862

C\l T3

o E E o O

o

■ffi ra a:

<D T3 O <D

a:

o 0 a. x LU

0.25

_______

_______V"

_______

__________

_____-v-

_________

_______V

0.2 -

0.15 -

! I-................'

-v— S2=2 -0--S2=4 -A - S2=6

i.................o................$................................>:<................o.........

-----A----A----A-----A----A----A----

0.1

.44

.52

.68 .76

r

.84

.92

Figure 10. ERC2 vs S2 and r

10

9.5

u> o O

JS 8.5 o

o 0

.

X LU

^t--------v---

6.5

------V--------V-

--V"-------v---

—v- — S1=11

0 ■■ S1=17

- A — S1=23

7.5/i-----A----A----A-----A----A----A----^

7™................$.................¡0:................0................................$................................"

.44 .52

.68 . 76 .84 .92

r

Figure 11. ETC vs Si and r

8

6. Arivarignan G., Keerthana M., Sivakumar B. A production inventory system with renewal and retrial demands. In: Applied Probability and Stochastic Processes. Joshua V., Varadhan S, Vishnevsky V. (eds). Infosys Sci. Found. Ser. Singapor: Springer, 2020. P. 129-142. DOI: 10.1007/978-981-15-5951-8.9

7. Artalejo J. R., Krishnamoorthy A., Lopez-Herrero M. J. Numerical analysis of (s, S) inventory systems with repeated attempts. Ann. Oper. Res., 2006. Vol. 141. P. 67-83. DOI: 10.1007/s10479-006-5294-8

8. Benny B., Chakravarthy S.R., Krishnamoorthy A. Queueing-inventory system with two commodities. J. Indian Soc. Probab. Stat., 2018. Vol. 19, No. 1. P. 437-454. DOI: 10.1007/s41096-018-0052-1

9. Chakravarthy S.R., Krishnamoorthy A., Joshua V. C. Analysis of a multi-server retrial queue with search of customers from the orbit. Perform. Eval., 2006. Vol. 63, No. 8. P. 776-798. DOI: 10.1016/j.peva.2005.09.002

10. Goyal S.K., Satir A. T. Joint replenishment inventory control: deterministic and stochastic models. European J. Oper. Res., 1989. Vol. 38, No. 1. P. 2-13. DOI: 10.1016/0377-2217(89)90463-3

8.2

m o O

~B 7.6 o

ü <u d X HI

7.4

7.2

6.8

-A— s=3

■O.....s=5

-Jv*-- s=7

7 -

.44

. -A----A----A-----A----A -

.................0................0................O................»

-A ■

■ -i\ ■I!

.52

.66 .74

r

.84

.92

Figure 12. ETC vs s and r

7.3

7.25

- 7.2

tn o

o

~B 7.15^-o

Î3 7.1 0

a. Ax

UU 7.05

6.95

-V"

,-v'

--V-- S2=2 S2=4

-V

____A----A

____A-----A----A

____A----A

.44

.52

.68

.76

.84

.92

Figure 13. ETC vs S2 and r

8

7

r

11. Jeganathan K., Reiyas M.A., Padmasekaran S., Lakshmanan K. An M/EK/1/N queueing-inventory system with two service rates based on queue lengths. Int. J. Appl. Comput. Math.., 2017. Vol. 3, No. 1. P. 357-386. DOI: 10.1007/s40819-017-0360-2

12. Jeganathan K., Melikov A. Z., Padmasekaran S. et al. A stochastic inventory model with two queues and a flexible server. Int. J. Appl. Comput. Math., 2019. Vol. 5, No. 1. Art. no. 21. P. 1-27. DOI: 10.1007/s40819-019-0605-3

13. Jeganathan K., Reiyas M.A. Two parallel heterogeneous servers Markovian inventory system with modified and delayed working vacations. Math. Comput. Simulation, 2020. Vol. 172. P. 273-304. DOI: 10.1016/j.matcom.2019.12.002

14. Jeganathan K., Harikrishnan T., Selvakumar S. et al. Analysis of interconnected arrivals on queueing-inventory system with two multi-server service channels and one retrial facility. Electronics, 2021. Vol. 10, No. 5. P. 1-35. DOI: 10.3390/electronics10050576

15. Jose K.P., Reshmi P.S. A production inventory model with deteriorating items and retrial demands. OPSEARCH, 2021. Vol. 58. P. 71-82. DOI: 10.1007/s12597-020-00471-8

16. Kalpakam S., Arivarignan G. A coordinated multicommodity (s, S) inventory system. Math. Comput. Model., 1993. Vol. 18, No. 11. P. 69-73. DOI: 10.1016/0895-7177(93)90206-E

17. Keerthana M., Saranya N., Sivakumar B. A stochastic queueing-inventory system with renewal demands and positive lead time. Eur. J. Ind. Eng., 2020. Vol. 14, No. 4. P. 443-484. DOI: 10.1504/ejie.2020.10029435

18. Krishnamoorthy A., Iqbal Basha R., Lakshmy B. Analysis of two commodity inventory problem. Int. J. Inf. Manag. Sen., 1994. Vol. 5, No. 1. P. 127-136.

19. Krishnamoorthy A., Merlymol J., Ravindranathan. Analysis of a bulk demand two commodity inventory problem. Calcutta Stat. Assoc. Bull, 1999. Vol. 49, No. 1-2. P. 193-194. DOI: 10.1177/0008068319990108

20. Krishnamoorthy A., Manikandan R., Lakshmy B. A revisit to queueing-inventory system with positive service time. Ann. Oper. Res., 2015. Vol. 233. P. 221-236. DOI: 10.1007/s10479-013-1437-x

21. Krishnamoorthy A., Shajin D., Lakshmy B. GI/M/1 type queueing-inventory systems with postponed work, reservation, cancellation and common life time. Indiana J. Pure Appl. Math., 2016. Vol. 47, No. 2. P. 357-388. DOI: 10.1007/s13226-016-0192-5

22. Krishnamoorthy A., Shajin D. MAP/PH/1 retrial queueing-inventory system with orbital search and reneging of customers. In: Analytical and Computational Methods in Probability Theory (ACMPT 2017). Rykov V., Singpurwalla N., Zubkov A. (eds). Lecture Notes in Comput. Sci., vol. 10684. Cham: Springer, 2017. P. 158-171. DOI: 10.1007/978-3-319-71504-9.15

23. Lakshmanan K., Padmasekaran S., Jeganathan K. Mathematical analysis of queueing-inventory model with compliment and multiple working vacations. Int. J. Eng. Adv. Teehnol., 2019. Vol. 8, No. 6. P. 42394240. DOI: 10.35940/ijeat.F9003.088619

24. Neuts M.F. Matrix-Geometrie Solutions in Stochastic Models: An Algorithmic Approach. New York: Dover Publication Inc., 1994. 332 p.

25. Senthil Kumar P. A discrete time two commodity inventory system. Int. J. Engrg. Res. Technol., 2016. Vol. 5, No. 1. P. 708-719. DOI: 10.17577/IJERTV5IS010602

26. Shajin D., Krishnamoorthy A., Dudin A. N. et al. On a queueing-inventory system with advanced reservation and cancellation for the next K time frames ahead: the case of overbooking. Queueing Syst., 2020. Vol. 94. P. 3-37. DOI: 10.1007/s11134-019-09631-0

27. Sivazlian B.D. Stationary analysis of a multicommodity inventory system with interacting set-up costs. SIAM J. App. Math., 1971. Vol. 20. Vol. 2. P. 264-278. https://www.jstor.org/stable/2099925?seq=1

28. Sivakumar B., Anbazhagan N., Arivarignan G. A two-commodity perishable inventory system. ORiON, 2005. Vol. 21, No. 2. P. 157-172. DOI: 10.5784/21-2-27

29. Sivakumar B., Anbazhagan N., Arivarignan G. Two commodity continuous review perishable inventory system. Int. J. Inf. Manag. Sa., 2006. Vol. 17, No. 3. P. 47-64.

30. Sivakumar B., Anbazhagan N., Arivarignan G. Two-commodity inventory system with individual and joint ordering policies and renewal demands. Stoch. Anal. Appl., 2007. Vol. 25, No. 6. P. 1217-1241. DOI: 10.1080/07362990701567314

31. Sivakumar B. Two-commodity inventory system with retrial demand. Eur. J. Oper. Res., 2008. Vol. 187, No. 1. P. 70-83. DOI: 10.1016/j.ejor.2007.02.036

32. Ushakumari P. V. On (s, S) inventory system with random lead time and repeated demands. J. Appl. Math. Stochastic Anal, 2006. Vol. 2006. Art. no. 81508. P. 1-22. DOI: 10.1155/JAMSA/2006/81508

33. Yadavalli V. S. S., Anbazhagan N., Arivarignan G. A two-commodity stochastic inventory system with lost sales. Stoch. Anal. Appl., 2004. Vol. 22, No. 2. P. 479-497. DOI: 10.1081/SAP-120028606

34. Yadavalli V. S.S., Arivarignan G., Anbazhagan N. Two commodity coordinated inventory system with Markovian demand. Asia-Pac. J. Oper. Res., 2006. Vol. 23, No. 4. P. 497-508. DOI: 10.1142/S0217595906001005

35. Yadavalli V. S. S., Sivakumar B., Arivarignan G., Olufemi Adetunji. A finite source multi-server inventory system with service facility. Comput. Ind.. Eng., 2012. Vol. 63. P. 739-753. DOI: 10.1016/j.cie.2012.04.014

36. Yadavalli V. S.S., Adetunji O., Sivakumar B., Arivarignan G. Two-commodity perishable inventory system with bulk demand for one commodity. S. Afr. J. Ind. Eng., 2010. Vol. 21, No. 1. P. 137-155. DOI: 10.7166/21-1-73