[Abstract] [PDF] [HTML] [Linked
References]
Comparative Study of Inventory Model for Duopolistic Market under Trade Credit
Deepa H Kandpal*, Khimya S Tinani#
Department of Statistics, Faculty of Science, The M. S. University of Baroda, Vadodara-390002, Gujarat, INDIA.
Corresponding Addresses:
*[email protected], #[email protected]
Abstract: This paper did some researches on stochastic inventory model in duopolistic market when permissible delay in payment was applied and the model was compared when the privilege is not offered. From the perspective of supplier’s and businessman mathematical model has been established. Spectral theory is used to derive explicit expression for the transition probabilities of a four state continuous time markov chain representing the status of the systems. These probabilities are used to compute the exact form of the average cost expression. We use concepts from renewal reward processes to develop average cost objective function. Optimal solution is obtained using Newton Rapson method in R programming. Finally sensitivity analysis for varying parameter values has been carried out.
Keywords: Future supply uncertainty, Permissible delay in payment, Duopolistic market, Two suppliers.
1. Introduction
This paper represent practical life situation by assuming that the supplier’s market is not monopolistic as competitive spirit in the business is increased especially after induction of multinational companies. In other words, it is assumed that the inventory manager may place his order with any one of two suppliers who are randomly available. Here we assume that the decision maker deals with two suppliers who may be ON or OFF. Here there are three states that correspond to the availability of at least one supplier that is states 0, 1 and 2 where as state 3 denotes the non availability of either of them. State 0 indicates that supplier 1 and supplier 2 both are available. Here it is assumed that one may place order to either one of the two suppliers or partly to both. State 1 represents that supplier 1 is available but supplier 2 is not available. State 2 represents that supplier 1 is not available but supplier 2 is available.
Silver [11] appears to be the first author to discuss the need for models that deal with supplier uncertainty. Articles by Parlar and Berkin [9] consider the supply uncertainty problem for a class of EOQ model with a single supplier where the availability and unavailability periods constitute an alternating Poisson process. Parlar and Berkin [9] assume that at any time the decision maker is aware of the availability status of the product although he does not know when the ON (available) and OFF (unavailable) periods will start and end. When the inventory level reaches the reorder point of zero and the status is ON, the order is received; otherwise the decision maker must wait until the product becomes available. Parlar and Perry [8] developed inventory model for non deteriorating items with future supply uncertainty considering demand rate d=1 for two suppliers. Kandpal and Gujarathi [4,5] has extended the model of Parlar & Perry [8] by considering demand rate greater than one and for deteriorating items for single supplier. Kandpal and Tinani [6] developed inventory model for deteriorating items with future supply uncertainty under inflation and permissible delay in payment for single supplier.
The traditional economic order quantity (EOQ) model is widely used as a decision tool for the control of inventory. The EOQ assumes that the retailer’s capitals are unrestricting and must be paid for the items as soon as the items are received. However, in practice, the supplier will offer the retailer a delay period which is the trade credit period, in paying for the amount of purchasing cost. Before the end of the trade credit period, the retailer can sell the goods and accumulate revenue and earn interest. A higher interest is charged if the payment is not settled by the end of the trade credit period. Permissible delay in payments is generally considered to be the equivalent form of discount contract. Goyal [3] has studied an EOQ system with deterministic demand and delay in payments is permissible which was reinvestigated by Chand and Ward [2]. Shah [10] developed model for deteriorating items when delay in payments is permissible by assuming deterministic demand. Aggarwal and Jaggi [1] developed a model to determine the optimum order quantity for deteriorating items under a permissible delay in payment.Kandpal and Tinani [7] developed perishable-inventory model under inflation and delay in payment allowing partial payment for single supplier.
Case 1: When permissible delay in payment is allowed. 2. Notations, Assumptions and Model
The inventory model here is developed on the basis of following assumptions.
(a) Demand rate d is deterministic and it is d>1.
(b) We define Xi and Yi be the random variables corresponding to the length of ON and OFF period respectively for ith supplier where i=1, 2. We specifically assume that Xi ~ exp (λi) and Yi ~ exp (µi). Further Xi and Yi are independently distributed.
(c) Ordering cost is Rs. k/order.
(d) Holding cost is Rs. h/unit/unit time.
(e) Shortage cost is Rs. p/unit.
(f) Time dependent part of the backorder cost is Rs. /unit/time.
(g) qi= order upto level i=0, 1, 2
(h) r=reorder upto level; qi and r are decision variables.
(i) Purchase cost is Rs. c/unit.
(j) T1i is a credit period allowed by ith supplier where i=1, 2 which is a known constant.
(k) T00 is cycle period.
(l) iei=Interest rate earned when purchase made from ith supplier where i=1, 2
ici=Interest rate charged by ith supplier where i=1, 2
(m) αi= Indicator variable for ith supplier i=1, 2
α1 = 0 if account is settled completely at T11
= 1 otherwise
α2 = 0 if account is settled completely at T12
= 1 otherwise
(n) =Interest earned over period (0 to ) =
(o) =Interest earned over period ( to ) upon interest earned previously.
(p) Interest charged by the ith supplier clearly (ici>iei) i=1, 2
In this paper, we assume that Supplier allows a fixed period ‘T1i’ to settle the account. During this fixed period no interest is charged by the ith supplier but beyond this period, interest is charged by the ith supplier under the terms and conditions agreed upon. Interest charged is usually higher than interest earned. The account is settled completely either at the end of the credit period or at the end of the cycle.During the fixed credit period T1i, revenue from sales is deposited in an interest bearing account.
The policy we have chosen is denoted by (q0, q1, q2, r). An order is placed for q i units i=0, 1, 2, whenever inventory drops to the reorder point r and the state found is i=0, 1, 2.When both suppliers are available, q0 is the total ordered from either one or both suppliers. If the process is found in state 3 that is both the suppliers are not available nothing can be ordered in which case the buffer stock of r units is reduced. If the process stays in state 3 for longer time then the shortages start accumulating at rate of d units/time. When the process leaves state 3 and supplier becomes available, enough units are ordered to increase the inventory to qi +r units where i=0, 1, 2.
=cost of ordering+ cost of holding inventory during a single interval that starts with an inventory of qi+r units and ends with r units.
= + i=0, 1, 2
=P (Being in state j at time t/starting in state i at time 0 ) i, j=0, 1, 2, 3
=long run probabilities i=0, 1, 2, 3
3. Optimal Policy Decision for the Model
For calculation of average cost objective function, we need to identify the cycles. Below given figure gives us the idea about cycles and their identification.
Fig 3.1: Inventory level and status process with two suppliers
Referring to Figure 3.1, we see that the cycles of this process start when the inventory goes up to a level of q0+r units. Once the cycle is identified, we construct the average cost objective function as a ratio of the expected cost per cycle to the expected cycle length.
i.e. Ac (q0, q1, q2, r) = where C00=E (cost per cycle) and T00=E (length of a cycle)
Analysis of the average cost function requires the exact determination of the transition probabilities , i, j=0, 1, 2, 3 for the four state CTMC. The solution is provided in the following lemma.
Lemma 3.1:Let t≥0, i, j=0, 1, 2, 3 be matrix of transition functions for the continuous time markov chain(CTMC). The exact transient solution is given as .where,
Proof: For proof refer Parlar and Perry[1996]
Corollary 3.2: The long run probabilities are
=
Proof: For proof refer Parlar and Perry[1996].
Define Ci0=E (cost incurred to the beginning of the next cycle from the time when inventory drops to r at state i=0, 1, 2, 3 and qi units are ordered if i=0, 1 or 2)
Lemma 3.3: is given by
i=0, 1, 2 (3.3.1)
(3.3.2)
Where with and
Proof:First consider i=0. Conditioning on the state of the supplier availability process when inventory drops to r, we obtain
The equation follows because being the initial inventory, when q0 units are used up we either observe state 0, 1, 2 or 3 with probabilities
, , and respectively. If we are in state 0 when r is reached, we must have incurred a cost of . On the other hand, if state j=1, 2, 3 is observed when inventory drops to r, then the expected cost will be with probability . The equation relating and are very similar but is obtained as
(3.3.1)
Here, is defined as the expected cost from the time inventory drops to r until either of the suppliers becomes available and it is computed as follows:
Now referring to fig 3.1., note that the cost incurred from the time when inventory drops to r and the state is OFF to the beginning of next cycle is equal to
Hence,
with as the rate of departure from state 3. This follows because if supplier availability process is in state 3 (OFF for both suppliers) when inventory drops to r, then the expected holding and backorder costs are equal to . If the process makes a transition to state 1, the total expected cost would then be . The probability of a transition from state 3 to state 1 is
P(Y1<Y2)= =
Multiplying this probability with the expected cost term above gives the first term of (3.3.1). The second term is obtained in a similar manner. Combining the results proves the lemma.
The following lemma provides a simpler means of expressing in an exact manner. To simplify the notation, we let , i=0, 1, 2 and i, j=0, 1, 2, 3.
Lemma 3.4: The exact expression for C00 is
(3.4.1)
where the pair [C10, C20] solves the system
= (3.4.2)
Proof: Rearranging the linear system of four equations in lemma(3.3) in matrix form gives
= (3.4.3)
We have from the last row of the system. Substituting this result in rows two and three and rearranging gives the system in (3.4.2). From the first row of (3.4.3) we obtain .
Hence above lemma is proved.
Define, Ti0=E [Time to the beginning of the next cycle from the time when inventory drops to r at state i=0, 1, 2, 3 and qi units are ordered if i=0, 1, 2]
Lemma 3.5: Expected cycle length is given by
where is the expected time from the time inventory drops to r until either supplier 1 or 2 becomes available.
Lemma 3.6:The exact expression for T 00 is
where the pair [T10, T20] solves the system.
=
The proof of the above two lemmas i.e. (3.5) and (3.6) are very similar to lemma (3.3) and (3.4).
Theorem 3.7: The Average cost objective function for two suppliers when delay in payment is considered is given by
Ac= =
Proof: Proof follows using Renewal reward theorem(RRT). The optimal solution for q0, q1, q2 and r is obtained by using Newton Rapson method in R programming.
4. Numerical Example
In this section we verify the results by a numerical example. We assume that
(i) k =Rs. 5/order, c=Rs.1/unit, d=20/units , h=Rs. 5/unit/time, π=Rs. 350/unit, =Rs.25/unit/time, α1=1, α2=1, ic1=0.11, ie1=0.02, ic2=0.13, ie2=0.04, T11=0.6, T12=0.8 , λ1=0.58, λ2=0.45, µ1=3.4, µ2=2.5
The last four parameters indicate that the expected lengths of the ON and OFF periods for first and second supplier are 1/λ1=1.72413794, 1/λ2=2.2222, 1/µ1=.2941176 and 1/µ2=.4 respectively. The long run probabilities are obtained as p0=0.7239588, p1=0.1303126, p2 =0.1234989 and p3=0.02222979. The optimal solution is obtained as
q0=4.032518, q1=27.63127, q2=27.43041, r=1.351202 and Ac= =6.499962
(ii) Keeping other parameters as it is, we consider α1=0 and α2=0. The optimal solution is obtained as
q0=5.72531, q1=30.10025, q2=30.0786,r=1.418232 and Ac= = 6.156267
(iii) Keeping other parameters as it is, we consider α1=1 and α2=0. The optimal solution is obtained as
q0=4.837963, q1=28.71599, q2=28.61197,r=1.407225 and Ac= = 6.315563
(iv) Keeping other parameters as it is, we consider α1=0 and α2=1. The optimal solution is obtained as
q0=4.68275, q1=28.45212, q2=28.32906,r=1.403514 and Ac= =6.362211
From the above numerical example, we conclude that the cost is minimum when account is settled at credit time given by the ith supplier.
Case 2: When permissible delay in payment is not allowed
The Average cost objective function for two suppliers when delay in payment is not considered reduces to
=
Numerical Example:
We assume thatk =Rs. 5/order, c=Rs.1/unit, d=20/units, h=Rs. 5/unit/time, π=Rs. 350/unit, =Rs.25/unit/time, λ1=0.58, λ2=0.45, µ1=3.4, µ2=2.5. The optimal solution is obtained as
q0=2.9845, q1=25.7849, q2=25.3637, r=1.2538 and Ac= =6.6247
On comparing the average cost of all the situations of case 1 we find that cost is maximum for case 2. This implies that the privilege offered proves to be bliss for entrepreneurs with minimum average cost than without privilege, so they are in business resulting in keen competition in the business.
5. Sensitivity Analysis for Case 1
(i) To observe the effect of varying parameter values on the optimal solution, we have conducted sensitivity analysis by varying the valueλ1 and keeping other parameter values fixed where α1=1 and α2=1. We resolve the problem to find optimal values of q0, q1, q2, r and Ac. The optimal values of q0, q1, q2 and Ac are plotted in Fig.5.1.
Table 5.1: Sensitivity Analysis Table by varying the parameter values of λ1 α1=1and α2=1
λ1 |
q0 |
q1 |
q2 |
r |
Ac |
0.5 |
3.717054 |
27.44038 |
27.2929 |
0.879031 |
6.551717 |
0.54 |
3.843987 |
27.48572 |
27.30287 |
1.123268 |
6.5296 |
0.58 |
4.032518 |
27.63127 |
27.43041 |
1.351202 |
6.499962 |
0.62 |
4.332144 |
27.97548 |
27.78156 |
1.563168 |
6.462281 |
0.66 |
4.89136 |
28.82827 |
28.68549 |
1.752149 |
6.414409 |
Fig 5.1: Sensitivity analysis graph for λ1(α1=1and α2=1)
We see that as λ1 increases i.e. expected length of ON period for 1st supplier decreases the value of q0, q1, q2 and r increases which result in decrease in average cost.
(ii) To observe the effect of varying parameter values on the optimal solution, we have conducted sensitivity analysis by varying the value λ1 and keeping other parameter values fixed where α1=0 and α2=0. We resolve the problem to find optimal values of q0, q1, q2, r and Ac. The optimal values of q0, q1, q2 and Ac are plotted in Fig.5.2.
Table 5.2: Sensitivity Analysis Table by varying the parameter values of λ1 α1=0and α2=0
λ1 |
q0 |
q1 |
q2 |
r |
Ac |
0.5 |
4.701216 |
28.5303 |
28.47871 |
0.99042 |
6.267463 |
0.54 |
5.102404 |
29.07288 |
29.02172 |
1.223672 |
6.219188 |
0.58 |
5.72531 |
30.10025 |
30.0786 |
1.418232 |
6.156267 |
0.62 |
6.714279 |
32.06434 |
32.09975 |
1.544664 |
6.073685 |
0.66 |
8.130382 |
35.44868 |
35.53845 |
1.56233 |
5.962785 |
Fig 5.2: Sensitivity analysis graph for λ1(α1=0and α2=0)
We see that as λ1 increases i.e. expected length of ON period for 1st supplier decreases the value of q0, q1, q2 and r increases which result in decrease in average cost.
(iii) To observe the effect of varying parameter values on the optimal solution, we have conducted sensitivity analysis by varying the value λ1 and keeping other parameter values fixed where α1=1 and α2=0. We resolve the problem to find optimal values of q0, q1, q2, r and Ac. The optimal values of q0, q1, q2 and Ac are plotted in Fig.5.3.
Table 5.3: Sensitivity Analysis Table by varying the parameter values of λ1 α1=1and α2=0
λ1 |
q0 |
q1 |
q2 |
r |
Ac |
0.5 |
4.17051 |
27.91604 |
27.81633 |
0.940845 |
6.406776 |
0.54 |
4.433138 |
28.17991 |
28.06476 |
1.185896 |
6.366733 |
0.58 |
4.837963 |
28.71599 |
28.61197 |
1.407225 |
6.315563 |
0.62 |
5.524972 |
29.85845 |
29.80357 |
1.59063 |
6.250326 |
.66 |
6.787964 |
32.46545 |
32.50053 |
1.683555 |
6.16309 |
Fig 5.3: Sensitivity analysis graph for λ1(α1=1and α2=0)
We see that as λ1 increases i.e. expected length of ON period for 1st supplier decreases the value of q0, q1, q2 and r increases which result in decrease in average cost.
(iv) To observe the effect of varying parameter values on the optimal solution, we have conducted sensitivity analysis by varying the value λ1 and keeping other parameter values fixed where α1=0 and α2=1. We resolve the problem to find optimal values of q0, q1, q2, r and Ac. The optimal values of q0, q1, q2 and Ac are plotted in Fig.5.4.
Table 5.4: Sensitivity Analysis Table by varying the parameter values of λ1 α1=0and α2=1
λ1 |
q0 |
q1 |
q2 |
r |
Ac |
0.5 |
4.13484 |
27.83933 |
27.73441 |
0.939671 |
6.425488 |
0.54 |
4.349138 |
28.03646 |
27.91042 |
1.181905 |
6.398667 |
0.58 |
4.68275 |
28.45212 |
28.32906 |
1.403514 |
6.362211 |
0.62 |
5.261809 |
29.36611 |
29.28288 |
1.595496 |
6.314107 |
0.66 |
6.438238 |
31.68256 |
31.69321 |
1.710866 |
6.247621 |
Fig 5.4: Sensitivity analysis graph for λ1(α1=0 and α2=1)
We see that as λ1 increases i.e. expected length of ON period for 1st supplier decreases the value of q0, q1, q2 and r increases which result in decrease in average cost.
Conclusion
From the above sensitivity analysis, we conclude that the cost is minimum when account is settled at credit time given by the ith supplier and is maximum when account is not settled at credit time given by the ith supplier. So it is advisable for the business man to settle the account at the credit time period given by the supplier, which means that we encourage the small business man by offering the privilege but simultaneously we discourage them for not clearing the account after the given credit time period.
References
- Aggarwal, S.P. and Jaggi, C.K. (1995), “Ordering policies of deteriorating items under permissible delay in payment”. Journal of the operational Research Society. 46,658-662.
- Chand, S. and Ward, J. (1987), “A note on Economic Order Quantity under conditions of Permissible Delay in Payments”. Journal of the Operational Research Society. 36, 83-84.
- Goyal, S.K. (1985), “Economic Order Quantity under Conditions of Permissible Delay in Payments”. Journal of the Operational Research Society. 34, 335-8.
- Kandpal, D.H. and Gujarathi, C.C. (2003), “Inventory Models of Future Supply Uncertainty”. Journal of Decision and Mathematical Science. 8, 79-91.
- Kandpal, D.H. and Gujarathi, C.C. (2006), ”Deteriorating items Inventory Model with Future Supply Uncertainty for Single Supplier”. Industrial Engineering Journal. 35(9), 23-28.
- Kandpal, D.H. and Tinani, K.S. (2009), ”Future supply uncertainty model for deteriorating items under inflation and permissible delay in payment for single supplier” Journal of Probability and Statistical Science. 7(2), 245-259.
- Kandpal, D.H. and Tinani, K.S. (2011), ”Perishable-Inventory Model under Inflation and Delay in Payment Allowing Partial Payment”. Journal of Indian Association for Productivity Quality and Reliability (IAPQR). Vol. 36(2), 111-131.
- Parlar, M. and Perry, D. (1996),“Inventory Models of Future Supply Uncertainty with Single and Multiple Suppliers”. Naval Research Logistic. 43, 191-210.
- Parlar, M. and Berkin, D. (1991),“Future Supply Uncertainty in EOQ Models”. Naval Research Logistics. 38,295-303.
- Shah, N.H. (1993), ”A Lot Size Model for Exponentially Decaying Inventory when Delay in payments is permissible”. Cahiers DUCERO. 35, 1-2,115-123.
- Silver, E.A. (1981), “Operations Research Inventory Management: A review and critique”. Operations Research. 29, 628-645.
Corresponding Author:
Deepa H Kandpal
Department of Statistics
Faculty of Science
The M.S. University of Baroda
Vadodara-390002, Gujarat, INDIA.
E-Mail ID: [email protected]
|