# An Inventory Model with Finite Replenishment Rate, Trade Credit Policy and Price-Discount Offer

Authors: Biswajit Sarkar; Shib Sankar Sana; Kripasindhu Chaudhuri
Journal: Journal of Industrial Engineering
 (2013)

---

## Abstract

When some suppliers offer trade credit periods and price discounts to retailers in order to increase the demand of their products, retailers have to face different types of discount offers and credits within which they have to take a decision which is the best offer for them to make more profit. The retailers try to buy perfect-quality items at a reasonable price, and also they try to invest returns obtained by selling those items in such a manner that their business is not hampered. In this point of view, we consider aneconomic order quantity (EOQ) model for various types of time-dependent demand when delay in payment and price discount are permitted by suppliers to retailers. The models of various demand patterns are discussed analytically. Some numerical examples and graphical representations are considered to illustrate the model.


---

## Body

## 1. Introduction

Many classical inventory models assume that demand is constant. In present marketing environment, few items follow constant demand. Many product’s demands follow variable time-varying demand. The recent trend of the marketing system is to provide more buy opportunities to the retailer by the supplier by offering different discounts. To take the discount opportunity, retailers prefer to buy more beyond their capacity of buying. As a result, the supplier has the opportunity to sell more for better earning. This is the benefit of the supplier. The classical inventory model does not consider the delay time concept or variable demand. The proposed model considers time-varying demand and delay in payments along with finite replenishment rate.The basic well-knownsquare root formula for the EOQ of the item was formulated by Harris [1] based on constant demand. Donaldson [2] extended the constant demand to linear time-dependent demand model analytically with finite time horizon. Following Donaldson [2], significant contribution in this direction came out from researchers like Goyal [3], Goswami and Chaudhuri [4], Goyal et al. [5], and others. Hariga and Benkherouf [6] discussed an optimal and heuristic replenishment model for deteriorating items with an exponentially time-varying demand. Wee [7] studied a deterministic lot size inventory model for deteriorating items with shortage and decline market. Khanra and Chaudhuri [8] extended an inventory model with quadratic increasing demand over a finite time horizon and shortages. Sana and Chaudhuri [9] studied an inventory model with linear trend demand incorporating shortages. Cárdenas-Barrón [10] discussed the derivation of inventory models by using analytic geometry and algebra. Sarkar et al. [11] explained an inventory model with quadratic time-varying demand by considering Euler-Lagrange method.It is common to all that every customer prefers to buy more at reduced price. Some researchers like Abad [12], and Kim and Hwang [13] developed the traditional quantity discount model. In the traditional EOQ model, it was assumed that the retailer pays the purchasing cost when he received the items from a supplier. In trade-credit policy, the supplier allows a certain fixed period to pay the purchasing cost. This fixed period which is settled by the supplier is called the credit period to the retailer. During this credit period, the supplier sells items to the retailer with different types of discounts to obtain more profit as early as possible during the credit period. Depending on this policy, Goyal [14] first developed an inventory model with permissible delay in payments. Aggarwal and Jaggi [15] developed an inventory model with an exponentially deteriorating rate by considering permissible delay in payments. Chu et al. [16] extended Goyal’s [14] model by considering the case of deterioration. Jamal et al. [17] extended an inventory model with shortages.Teng [18] developed an EOQ model for a retailer to order small lot size in order to take the benefit of permissible delay in payments. Arcelus et al. [19] developed an inventory model by considering the retailer’s maximizing profit and inventory policies for vendor’s trade promotion offer of price/credit on the purchase of perishable items. Huang [20] extended an inventory model of retailer’s inventory system as a cost minimization model to determine the retailer’s optimum inventory cycle time and optimal order quantity. Huang [21] developed an economic production quantity (EPQ) model of retailer’s inventory system to investigate the optimal retailer’s decisions under two levels of trade credit policy. Cárdenas-Barrón [22] extended optimal ordering policies in response to a discount offer. Teng et al. [23] explained optimal ordering decisions with returns and excess inventory. Sarkar [24] discussed an inventory model with delay in payments in the presence of imperfect production. Sarkar [25] developed an inventory model with delay in payments and time-varying deterioration rate. Forghani et al. [26] explained an inventory model in the single period inventory system with price adjustment.This paper considers an inventory model for credit periods and price-discount offers with different types of time varying demand and constant supply rateK up to time t=t1. During [0,t1], inventory piles up by adjusting the demand. The accumulated inventory level at time t1 depletes gradually to meet the demand, and the level reaches zero level at time T(t1≤T). The agreement between the supplier and the retailer is such that total purchasing cost of whole amount (Kt1) would be paid within the time R(R>t1) with purchasing cost at discount rate ρ. The different delay periods with different discount rates on the purchasing cost are permitted by the supplier to the retailer. During the credit period, the retailer can earn interest by selling items whereas interest of purchasing cost is charged against the delay of excess time of credit of payment period by the retailer to the supplier.The rest of the paper is designed as follows. The mathematical model is presented in Section2. In Section 3, numerical examples are given. Finally, concluding remarks are explained in Section 4. See Table 1 for the comparison of this model with previous works.Table 1
Comparisons of this model with previous works.
Author/authors
Linear/quadratic demand
Exponential  demand
Constant demand
Tread credit policy
Harris [1]
✓
Donaldson [2]
✓
Goyal [14]
✓
✓
Goyal [3]
✓
Goswami and Chaudhuri [4]
✓
Hariga and Benkherouf [6]
✓
✓
Teng [18]
✓
✓
Khanra and Chaudhuri [8]
✓
Arcelus et al. [19]
✓
✓
Sana and Chaudhuri [9]
✓
Huang [20]
✓
✓
Huang [21]
✓
✓
Cárdenas-Barrón [10]
✓
✓
Teng et al. [23]
✓
✓
Sarkar [24]
✓
✓
This paper
✓
✓
✓
✓


## 2. Mathematical Model

We consider the following notation to develop the model.T*:
the optimal length of inventory cycle (decision variable)t1*:
the optimal duration of replenishment (decision variable)Q1(t):
on-hand inventory at timet(0≤t≤t1)Q2(t):
on-hand inventory at timet(t1≤t≤T)D(t):
time-varying demand rateK:
constant replenishment/supply rateR:
variable delay periodRi:
ith permissible delay periodρi:
discount rate on purchasing cost atith permissible delay periodC1:
ordering cost per orderC2:
unit holding cost per unit time, excluding interest chargeC3:
purchasing cost per unitC4:
maximum retail price per unitP:
selling price per unitIc:
rate of interest gaining due to the credit balanceIf:
rate of interest due to financing inventoryT:
length of the inventory cyclet1:
duration of the replenishmentZ1i:
average profit of the system whenT≥RiZ2i:
average profit of the system whenT≤Ri.The following assumptions are considered to develop this model.(1)
The inventory system involves only single type of product.(2)
The demand rate is constant or time dependent (quadratic, linear, and exponential).(3)
Different discount rates on the purchasing cost for different delay periods are considered.(4)
Replenishment rate is instantaneously infinite, but its size is finite.(5)
Time horizon is infinite.(6)
Lead time is negligible.(7)
Neither shortage nor backlogging is considered.The cycle starts with zero inventory at supply rateK. The replenishment or supply continues up to time t1. During the time span [0,t1], inventory piles up by adjusting the demand in the market. This accumulated inventory level at time t1 depletes gradually to meet the demand and it reaches zero level at time t=T(T>t1). Generally, the supplier offers delay period R(R>t1) to the retailer to pay the total purchasing cost (C3Kt1) of items. For different delay periods Ri(i=1,2,3), different discounts ρi(i=1,2,3) of purchasing cost are offered to the retailer by the supplier. In this direction, we consider the purchasing cost of different delay periods as follows:
(1)C3={C4(1-ρ1)whenR=R1,C4(1-ρ2)whenR=R2,C4(1-ρ3)whenR=R3,∞whenR>R3,
where Ri′s are the ith permissible delay to settle the purchasing cost at which the discount rate to the retailer is ρi. Also C3 tends to be ∞ at R>R3. That is, at infinite purchasing cost, the retailer never purchases any item from the supplier. Indirectly, the supplier would not supply the product to the retailer while delay period R exceeds R3. In our model two cases may arise.Case 1 (T≥R).
That is, inventory cycle lengthT is larger or equal to the credit period R (see Figure 1). When T≥R, there are some profits based on credit balance during the delay period and there is some interest charged due to financing inventory during [R,T].Figure 1
Inventory versus time (Case1  T≥R).Case 2 (T≤R).
That is, inventory cycle lengthT is smaller or equal to the credit period R (see Figure 2). When T≤R, there are some profits based on credit balance during the delay period and there is no interest charged due to financing inventory.Figure 2
Inventory versus time (Case2  T≥R).The governing differential equations of this model are(2)dQ1dt=K-D(t)withQ1(0)=0,0≤t≤t1,dQ2dt=-D(t)withQ2(T)=0,t1≤t≤T.From (2), we obtain
(3)Q1(t)=Kt-∫0tD(t)dt,0≤t≤t1,(4)Q2(t)=∫tT‍D(t)dt,t1≤t≤T.By utilizing the continuity att1, Q1(t1)=Q2(t1), we have
(5)Kt1-∫0t1‍D(t)dt=∫t1T‍D(t)dt,that is,(6)t1=1K∫0TD(t)dt.Now we formulate Cases1 and 2 using (3) to (6).Case 1 (T≥R).
When the inventory cycle(T) is larger or equal to the credit period R, the holding cost excluding the interest charges is C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}.
The profit gains due to credit balance during the delay period[0,R] are IcP∫0R‍(R-t)D(t)dt.
The interest charged for financing inventory during[R,T] is IfC3∫RT‍Q2(t)dt.
Therefore, the total profit is (see Figure1)
(7)P1i=[∫t1Ri‍(P-C3)Kt1P1i=55.+IcP{∫0t1‍(Ri-t)D(t)dt+∫t1Ri‍(Ri-t)D(t)dt}P1i=55.-C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}P1i=55.-IfC3∫RiT‍Q2(t)dt-C1]fori∈{1,2,3}.
Hence, the average profit is(8)Z1i=P1iT=1T[∫RiT‍(P-C3)Kt1=1T55.+IcP{∫0t1‍(Ri-t)D(t)dt+∫t1Ri‍(Ri-t)D(t)dt}=1T55.-C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}=1T55.-IfC3∫RiT‍Q2(t)dt-C1]fori∈{1,2,3}.Case 2 (T≤R).
When the inventory cycle (T) is smaller or equal to the credit period R, the holding cost excluding the interest charges is C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}.
The profit gains due to credit balance during the time[0,R] are IcP[∫0T‍(T-t)D(t)dt+Kt1(R-T)].
Therefore, the total profit is (see Figure2)
(9)P2i=[∫0t1∫t1T‍‍(P-C3)Kt155..+IcP{∫t1T‍∫0t1‍(T-t)D(t)dt55..+icp551+∫t1T‍(T-t)D(t)dt+Kt1(Ri-T)}55..-C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}-C1]5.5-C2{∫0t1‍Q1(t)dt+∫t1T‍}fori∈{1,2,3}.
Hence, the average profit is(10)Z2i=P2iT=1T[{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}{∫t1T‍∫0t1‍(T-t)D(t)dt(P-C3)Kt1=1T55+IcP{∫0t1‍(T-t)D(t)dt=1T55+IcP55+∫t1T‍(T-t)D(t)dt+Kt1(Ri-T)}=1T55-C2{∫0t1‍Q1(t)dt+∫t1T‍Q2(t)dt}-C1{∫0t1‍(T-t)D(t)dt]C2=1T55-C1{∫0t1‍(T-t)D(t)dt55{∫0t1‍Q1(t)dt}fori∈{1,2,3}.
Now our objective is to maximizeZ=Sup{Z1i,Z2i} for i∈{1,2,3} and obtain the optimal replenishment period t1* and inventory cycle length T*. We discuss various time-dependent or constant demands by utilizing this general formula.


### 2.1. Quadratic Demand Pattern

The quadratic time-dependent demand is of the formD(t)=a+bt+ct2 where a>0, b>0, and c>0. This trend of demand is applied to the products like essential commodities and seasonal goods.From (2), we have
(11)Q1(t)=Kt-at-bt22-ct33,Q2(t)=a(T-t)+b2(T2-t2)+c3(T3-t3).
Using these Q1(t) and Q2(t) in (8) and (10), we obtain
(12)Z1i=1T[(T3Ri-Ri44)(P-C3)Kt1+IcP{aRi22+bRi36+cRi412}-C2{Kt122+aT22+bT33+cT44-aTt1-bT2t12-cT3t13}-IfC3{aT22+bT33+cT44-a(TRi-Ri22)-b2(T2Ri-Ri33)-c3(T3Ri-Ri44)}-C1],Z2i=1T[+cT44-aTt1-bT2t12-cT3t13}-C1(P-C3)Kt1+IcP{aT22+bT36+cT412+Kt1(Ri-T)}-C2{Kt122+aT22+bT33+cT44-aTt1-bT2t12-cT3t13}-C1].
Using (6) in the above objective functions, we obtain
(13)Z1i=[(P-C3)(a+bT2+cT23)+IcPT{aRi22+bRi36+cRi412}-C2{12K(a2T+b2T34+c2T59+abT2+bcT43+2acT33)+aT2+bT23+cT34-a2TK-abT2K-2acT33K-b2T34K-bcT43K-c2T59K}-IfC3{aT2+bT23+cT34-aRi+aRi22T-bTRi2+bRi36T-cT2Ri3+cRi412T}-C1T],Z2i=[(P-C3)(a+bT2+cT23)+IcP{aRi+bTRi2+cT2Ri3-aT2-bT23-cT34}-C2{12K(a2T+b2T34+c2T59+abT2+bcT43+2acT33)+aT2+bT23+cT34-a2TK-abT2K-2acT33K-b2T34K-bcT43K-c2T59K}-C1T].To obtain the optimal cycle length, we construct two lemmas as follows.Lemma 1.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if [4ξ1(c)(T*)6+ξ2(b,c)(T*)5+ξ3(a,b,c)(T*)4+ξ4(a,b,c)(T*)3-2ξ6(a,b,c)]<0. Otherwise, Z1i has a maximum value at unique T*={T∣(c/3)T3+(b/2)T2+aT-KRi=0,T>0} if [ξ1(c)Ri6+ξ2(b,c)Ri5+ξ3(a,b,c)Ri4+ξ4(a,b,c)Ri3+ξ5(a,b)Ri2+ξ6(a,b,c)]>0 holds.Proof.
See AppendixA.Lemma 2.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if [5(T*)6ξ1(c)+(T*)5ξ2(b,c)+2(T*)4ξ7(a,b,c)+(T*)3ξ8(a,b,c)-2C1]<0. Otherwise, Z2i has a maximum value at T=Ri if [(Ri)6ξ1(c)+(Ri)5ξ2(b,c)+(Ri)4ξ7(a,b,c)+(Ri)3ξ8(a,b,c)+(Ri)2ξ9(a,b,c)+C1]>0.Proof.
See AppendixA.


### 2.2. Linear Time-Dependent Demand Pattern

We consider linear time-dependent demand which can be obtained if we substitutec=0 in the form of quadratic demand. We get D(t)=a+bt, a>0, b>0 (see for instances Donaldson [2], Goyal [3], and Goswami and Chaudhuri [4]). Generally, in some computer games, computer android applications, we obtain this type of linear time-dependent demand.Now from (2), we have
(14)Q1(t)=Kt-at-bt22,
and Q2(t)=a(T-t)+(b/2)(T2-t2), respectively.UsingQ1(t) and Q2(t) in (8) and (10), we obtain
(15)Z1i=1T[(T2Ri-Ri33)(P-C3)Kt1+IcP{aRi22+bRi36}-C2{Kt122+aT22+bT33-aTt1-bT2t12}-IfC3{aT22+bT33-a(TRi-Ri22)-b2(T2Ri-Ri33)}-C1],Z2i=1T[{Kt122+aT22+bT33-aTt1-bT2t12}(P-C3)Kt1+IcP{aT22+bT36+Kt1(Ri-T)}-C2{Kt122+aT22+bT33-aTt1-bT2t12}-C1].Using (6) in the above objective functions, we have
(16)Z1i=[(P-C3)(a+bT2)+IcPT{aRi22+bRi36}-C2{12K(a2T+b2T34+abT2)+aT2+bT23-a2TK-abT2K-b2T34K}-IfC3{aT2+bT23-aRi+aRi22T-bTRi2+bRi36T}-C1T],Z2i=[b2T34K(P-C3)(a+bT2)+IcP{aRi+bTRi2-aT2-bT23}-C2{12K(a2T+b2T34+abT2)+aT2+bT23-a2TK-abT2K-b2T34K}-C1T].To obtain the optimal cycle length, we construct the following lemmas.Lemma 3.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if [2(T*)4ξ3(a,b,0)+(T*)3ξ4(a,b,0)-2ξ6(a,b,0)]<0. Otherwise Z1i has a global maximum value at T=((-a)+(a2+2KbRi))/b if Ri≤2(K-a)/b holds.Proof.
See AppendixB.Lemma 4.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if [2(T*)4ξ3(a,b,0)+(T*)3ξ8(a,b,0)-2C1]<0. Otherwise Z2i has a maximum value at T=Ri if [(Ri)4ξ3(a,b,0)+(Ri)3ξ8(a,b,0)+(Ri)2ξ5(a,b)+C1]>0.Proof.
See AppendixB.


### 2.3. Constant Demand Pattern

We consider the demand as constant which can be found by substitutingb=0 in the linear demand. We obtain D=a, a>0. See for instances Harris [1] and Cárdenas-Barrón [10, 22]. This type of demand is usually found in product’s life cycle.Now from (2), we have Q1(t)=(K-a)t and Q2(t)=a(T-t), respectively.UsingQ1(t) and Q2(t) in (8) and (10), we get
(17)Z1i=1T[(P-C3)Kt1+IcP(aRi22)-C2{Kt122+aT22-aTt1}-IfC3{aT22-a(TRi-Ri22)}-C1],(18)Z2i=1T[(P-C3)Kt1+IcP{aT22+Kt1(Ri-T)}-C2{Kt122+aT22-aTt1}-C1],
respectively. Using (6) in the above functions, we obtain
(19)Z1i=[(P-C3)a+IcPaRi22T-C2{aT2-a2T2K}-IfC3{aT2-aRi+aRi22T}-C1T],Z2i=[{aT2-a2T2K}-C1T(P-C3)(a)+IcP{aRi-aT2}-C2{aT2-a2T2K}-C1T].
To obtain optimal cycle length, we formulate the following lemmas.Lemma 5.
If the conditionsRi2<K{2C1+aRi2(IfC3-IcP)}/(C2a(K-a)+IfC3aK)<K2Ri2/a2 and K{2C1+aRi2(IfC3-IcP)}>0 hold, then Z1i has a global maximum at T=T*=[K{2C1+aRi2(IfC3-IcP)}/(C2a(K-a)+IfC3aK)]1/2.Proof.
See AppendixC.Lemma 6.
If0<2C1K/(C2a(K-a)-IcPaK)<Ri2 hold, then Z2i has a global maximum at T=T*=[2C1K/(C2a(K-a)-IcPaK)]1/2.Proof.
See AppendixC.


### 2.4. Exponential Demand Pattern

We consider another important type of demand as exponential type, that is,D(t) varies exponentially with time t. In this case, the demand rate increases very fast as in seasonal goods, new computer parts like RAM and data storage device. For this type of demand rate, we consider D(t)=aexp(bt) where a>0, b>0.From (2), we have
(20)Q1(t)=Kt+ab{1-exp(bt)},Q2(t)=ab{exp(bT)-exp(bt)}.
Using these Q1(t) and Q2(t) in (8) and (10), we obtain
(21)Z1i=1T[-exp(bT)b-t1exp(bT))(P-C3)Kt1-IcP{aRib+ab2(1-exp(bRi))}-C2{-exp(bT)b-t1exp(bT))Kt122+ab(-exp(bT)bt1+1b+Texp(bT)-exp(bT)b-t1exp(bT))-exp(bT)b-t1exp(bT)Kt122)}-IfC3ab{-exp(bT)-exp(bRi)b}-C1(T-Ri)exp(bT)-exp(bT)-exp(bRi)b}-C1-exp(bT)b-t1exp(bT))],Z2i=1T[{Kt122+ab(t1+1b+Texp(bT)-exp(bT)b(P-C3)Kt1-IcP{aTb+ab2(1-exp(bT))-Kt1(Ri-T)(aTb)}-C2{Kt122+ab(t1+1b+Texp(bT)-exp(bT)b-t1exp(bT)+1b)Kt122}-C1].
From (6), we have
(22)t1=1K∫0T‍aexp(bt)dt=aKb[exp(bT)-1].Substituting this inZ1i and Z2i, we have
(23)Z1i=abT[exp(bT)-exp(bRi)b(P-C3)(exp(bT)-1)-IcP{Ri+1b(1-exp(bRi))}-C2{-aexp(2bT)2Kb-a2Kb1b+Texp(bT)-exp(bT)b+aexp(bT)Kb-aexp(2bT)2Kb-a2Kb}-IfC3{-exp(bT)-exp(bRi)b(T-Ri)exp(bT)-exp(bT)-exp(bRi)b}-bC1a+1bexp(bT)-exp(bRi)b],Z2i=abT[1-exp(bT)b+(P-C3)(exp(bT)-1)-IcP{1-exp(bT)b+Ri+(T-Ri)exp(bT)}-C2{-aexp(2bT)2Kb-a2Kb)1b+Texp(bT)-exp(bT)b+aexp(bT)Kb-aexp(2bT)2Kb-a2Kb)}-bC1a].Lemma 7.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if (d2Z1i/dT2)|T=T*<0. Otherwise, Z1i has a unique maximum value T=T*=(1/b)[ln(1+KbRi/a)], if (dZ1i/dT)|T=Ri>0.Proof.
See AppendixD.Lemma 8.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if (d2Z2i/dT2)|T=T*<0. Otherwise, Z2i has a maximum value T=Ri, if (dZ2i/dT)|T=Ri>0Proof.
See AppendixD.


## 2.1. Quadratic Demand Pattern

The quadratic time-dependent demand is of the formD(t)=a+bt+ct2 where a>0, b>0, and c>0. This trend of demand is applied to the products like essential commodities and seasonal goods.From (2), we have
(11)Q1(t)=Kt-at-bt22-ct33,Q2(t)=a(T-t)+b2(T2-t2)+c3(T3-t3).
Using these Q1(t) and Q2(t) in (8) and (10), we obtain
(12)Z1i=1T[(T3Ri-Ri44)(P-C3)Kt1+IcP{aRi22+bRi36+cRi412}-C2{Kt122+aT22+bT33+cT44-aTt1-bT2t12-cT3t13}-IfC3{aT22+bT33+cT44-a(TRi-Ri22)-b2(T2Ri-Ri33)-c3(T3Ri-Ri44)}-C1],Z2i=1T[+cT44-aTt1-bT2t12-cT3t13}-C1(P-C3)Kt1+IcP{aT22+bT36+cT412+Kt1(Ri-T)}-C2{Kt122+aT22+bT33+cT44-aTt1-bT2t12-cT3t13}-C1].
Using (6) in the above objective functions, we obtain
(13)Z1i=[(P-C3)(a+bT2+cT23)+IcPT{aRi22+bRi36+cRi412}-C2{12K(a2T+b2T34+c2T59+abT2+bcT43+2acT33)+aT2+bT23+cT34-a2TK-abT2K-2acT33K-b2T34K-bcT43K-c2T59K}-IfC3{aT2+bT23+cT34-aRi+aRi22T-bTRi2+bRi36T-cT2Ri3+cRi412T}-C1T],Z2i=[(P-C3)(a+bT2+cT23)+IcP{aRi+bTRi2+cT2Ri3-aT2-bT23-cT34}-C2{12K(a2T+b2T34+c2T59+abT2+bcT43+2acT33)+aT2+bT23+cT34-a2TK-abT2K-2acT33K-b2T34K-bcT43K-c2T59K}-C1T].To obtain the optimal cycle length, we construct two lemmas as follows.Lemma 1.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if [4ξ1(c)(T*)6+ξ2(b,c)(T*)5+ξ3(a,b,c)(T*)4+ξ4(a,b,c)(T*)3-2ξ6(a,b,c)]<0. Otherwise, Z1i has a maximum value at unique T*={T∣(c/3)T3+(b/2)T2+aT-KRi=0,T>0} if [ξ1(c)Ri6+ξ2(b,c)Ri5+ξ3(a,b,c)Ri4+ξ4(a,b,c)Ri3+ξ5(a,b)Ri2+ξ6(a,b,c)]>0 holds.Proof.
See AppendixA.Lemma 2.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if [5(T*)6ξ1(c)+(T*)5ξ2(b,c)+2(T*)4ξ7(a,b,c)+(T*)3ξ8(a,b,c)-2C1]<0. Otherwise, Z2i has a maximum value at T=Ri if [(Ri)6ξ1(c)+(Ri)5ξ2(b,c)+(Ri)4ξ7(a,b,c)+(Ri)3ξ8(a,b,c)+(Ri)2ξ9(a,b,c)+C1]>0.Proof.
See AppendixA.


## 2.2. Linear Time-Dependent Demand Pattern

We consider linear time-dependent demand which can be obtained if we substitutec=0 in the form of quadratic demand. We get D(t)=a+bt, a>0, b>0 (see for instances Donaldson [2], Goyal [3], and Goswami and Chaudhuri [4]). Generally, in some computer games, computer android applications, we obtain this type of linear time-dependent demand.Now from (2), we have
(14)Q1(t)=Kt-at-bt22,
and Q2(t)=a(T-t)+(b/2)(T2-t2), respectively.UsingQ1(t) and Q2(t) in (8) and (10), we obtain
(15)Z1i=1T[(T2Ri-Ri33)(P-C3)Kt1+IcP{aRi22+bRi36}-C2{Kt122+aT22+bT33-aTt1-bT2t12}-IfC3{aT22+bT33-a(TRi-Ri22)-b2(T2Ri-Ri33)}-C1],Z2i=1T[{Kt122+aT22+bT33-aTt1-bT2t12}(P-C3)Kt1+IcP{aT22+bT36+Kt1(Ri-T)}-C2{Kt122+aT22+bT33-aTt1-bT2t12}-C1].Using (6) in the above objective functions, we have
(16)Z1i=[(P-C3)(a+bT2)+IcPT{aRi22+bRi36}-C2{12K(a2T+b2T34+abT2)+aT2+bT23-a2TK-abT2K-b2T34K}-IfC3{aT2+bT23-aRi+aRi22T-bTRi2+bRi36T}-C1T],Z2i=[b2T34K(P-C3)(a+bT2)+IcP{aRi+bTRi2-aT2-bT23}-C2{12K(a2T+b2T34+abT2)+aT2+bT23-a2TK-abT2K-b2T34K}-C1T].To obtain the optimal cycle length, we construct the following lemmas.Lemma 3.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if [2(T*)4ξ3(a,b,0)+(T*)3ξ4(a,b,0)-2ξ6(a,b,0)]<0. Otherwise Z1i has a global maximum value at T=((-a)+(a2+2KbRi))/b if Ri≤2(K-a)/b holds.Proof.
See AppendixB.Lemma 4.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if [2(T*)4ξ3(a,b,0)+(T*)3ξ8(a,b,0)-2C1]<0. Otherwise Z2i has a maximum value at T=Ri if [(Ri)4ξ3(a,b,0)+(Ri)3ξ8(a,b,0)+(Ri)2ξ5(a,b)+C1]>0.Proof.
See AppendixB.


## 2.3. Constant Demand Pattern

We consider the demand as constant which can be found by substitutingb=0 in the linear demand. We obtain D=a, a>0. See for instances Harris [1] and Cárdenas-Barrón [10, 22]. This type of demand is usually found in product’s life cycle.Now from (2), we have Q1(t)=(K-a)t and Q2(t)=a(T-t), respectively.UsingQ1(t) and Q2(t) in (8) and (10), we get
(17)Z1i=1T[(P-C3)Kt1+IcP(aRi22)-C2{Kt122+aT22-aTt1}-IfC3{aT22-a(TRi-Ri22)}-C1],(18)Z2i=1T[(P-C3)Kt1+IcP{aT22+Kt1(Ri-T)}-C2{Kt122+aT22-aTt1}-C1],
respectively. Using (6) in the above functions, we obtain
(19)Z1i=[(P-C3)a+IcPaRi22T-C2{aT2-a2T2K}-IfC3{aT2-aRi+aRi22T}-C1T],Z2i=[{aT2-a2T2K}-C1T(P-C3)(a)+IcP{aRi-aT2}-C2{aT2-a2T2K}-C1T].
To obtain optimal cycle length, we formulate the following lemmas.Lemma 5.
If the conditionsRi2<K{2C1+aRi2(IfC3-IcP)}/(C2a(K-a)+IfC3aK)<K2Ri2/a2 and K{2C1+aRi2(IfC3-IcP)}>0 hold, then Z1i has a global maximum at T=T*=[K{2C1+aRi2(IfC3-IcP)}/(C2a(K-a)+IfC3aK)]1/2.Proof.
See AppendixC.Lemma 6.
If0<2C1K/(C2a(K-a)-IcPaK)<Ri2 hold, then Z2i has a global maximum at T=T*=[2C1K/(C2a(K-a)-IcPaK)]1/2.Proof.
See AppendixC.


## 2.4. Exponential Demand Pattern

We consider another important type of demand as exponential type, that is,D(t) varies exponentially with time t. In this case, the demand rate increases very fast as in seasonal goods, new computer parts like RAM and data storage device. For this type of demand rate, we consider D(t)=aexp(bt) where a>0, b>0.From (2), we have
(20)Q1(t)=Kt+ab{1-exp(bt)},Q2(t)=ab{exp(bT)-exp(bt)}.
Using these Q1(t) and Q2(t) in (8) and (10), we obtain
(21)Z1i=1T[-exp(bT)b-t1exp(bT))(P-C3)Kt1-IcP{aRib+ab2(1-exp(bRi))}-C2{-exp(bT)b-t1exp(bT))Kt122+ab(-exp(bT)bt1+1b+Texp(bT)-exp(bT)b-t1exp(bT))-exp(bT)b-t1exp(bT)Kt122)}-IfC3ab{-exp(bT)-exp(bRi)b}-C1(T-Ri)exp(bT)-exp(bT)-exp(bRi)b}-C1-exp(bT)b-t1exp(bT))],Z2i=1T[{Kt122+ab(t1+1b+Texp(bT)-exp(bT)b(P-C3)Kt1-IcP{aTb+ab2(1-exp(bT))-Kt1(Ri-T)(aTb)}-C2{Kt122+ab(t1+1b+Texp(bT)-exp(bT)b-t1exp(bT)+1b)Kt122}-C1].
From (6), we have
(22)t1=1K∫0T‍aexp(bt)dt=aKb[exp(bT)-1].Substituting this inZ1i and Z2i, we have
(23)Z1i=abT[exp(bT)-exp(bRi)b(P-C3)(exp(bT)-1)-IcP{Ri+1b(1-exp(bRi))}-C2{-aexp(2bT)2Kb-a2Kb1b+Texp(bT)-exp(bT)b+aexp(bT)Kb-aexp(2bT)2Kb-a2Kb}-IfC3{-exp(bT)-exp(bRi)b(T-Ri)exp(bT)-exp(bT)-exp(bRi)b}-bC1a+1bexp(bT)-exp(bRi)b],Z2i=abT[1-exp(bT)b+(P-C3)(exp(bT)-1)-IcP{1-exp(bT)b+Ri+(T-Ri)exp(bT)}-C2{-aexp(2bT)2Kb-a2Kb)1b+Texp(bT)-exp(bT)b+aexp(bT)Kb-aexp(2bT)2Kb-a2Kb)}-bC1a].Lemma 7.
When(dZ1i/dT)|T=T*=0 exists for T*∈[Ri,∞), then Z1i has a maximum value at T=T* if (d2Z1i/dT2)|T=T*<0. Otherwise, Z1i has a unique maximum value T=T*=(1/b)[ln(1+KbRi/a)], if (dZ1i/dT)|T=Ri>0.Proof.
See AppendixD.Lemma 8.
When(dZ2i/dT)|T=T*=0 exists for T*∈[0,Ri], then Z2i has a maximum value at T=T* if (d2Z2i/dT2)|T=T*<0. Otherwise, Z2i has a maximum value T=Ri, if (dZ2i/dT)|T=Ri>0Proof.
See AppendixD.


## 3. Numerical Examples

Example  A. We consider the following parametric values in appropriate units: C1=$250 per order, C2=$20/unit/year, R1=120/365 year, R2=150/365 year, R3=180/365 year, ρ1=10%, ρ2=5%, ρ3=2%, a=100 units, b=50 units, c=10 units, D(t)=a+bt+ct2, C4=$120/units, P=$180/units, If=0.16/$/year, Ic=0.13/$/year, and K=500 units.Then the optimal solutions are{Z11=$7522.33, T*=0.85 year, t1*=0.21year}, {Z12=$6980.84, T*=0.72 year, t1*=0.17year},  {Z13=$6760.68, T*=0.66 year, t1*=0.15year}, {Z21=$7174.54, T*=0.33 year, t1*=0.07year}, {Z22=$6852.25, T*=0.41 year, t1*=0.91year}, {Z23=$6725.72, T*=0.49 year, t1*=0.11year}.Among the above optimal solutions, the optimal solution is{Z11=$7522.33, T*=0.85 year, t1*=0.21 year} global maximum. The average profit function is highly nonlinear. Figure 3 shows the concavity of the functions and comparisons of the cost functions.Figure 3
The figures of the average profit functions are in descending order asZ11, Z21, Z12, Z22, Z13, and Z23.Example  B. We consider the following parametric values in appropriate units: C1=$250 per order, C2=$20/unit/year, R1=120/365 year, R2=150/365 year, R3=180/365 year, ρ1=10%, ρ2=5%, ρ3=2%, a=100 units, b=50 units, D(t)=a+bt, C4=$120/units, P=$180/units, If=0.16/$/year, Ic=0.13/$/year, and K=500 units.Then the optimal solutions are{Z11=$7404.27, T*=0.70 year, t1*=0.17year}, {Z12=$6895.63, T*=0.62 year, t1*=0.14year}, {Z13=$6689.02, T*=0.58 year, t1*=0.13year}, {Z21=$7149.18, T*=0.33 year, t1*=0.07year}, {Z22=$6816.19, T*=0.41 year, t1*=0.09year}, and {Z23=$6677, T*=0.49 year, t1*=0.11year}.Among the above optimal solutions, the optimal solution is{Z11=$7404.27, T*=0.70 year, t1*=0.17year} global maximum. The average profit function is highly nonlinear. Figure 4 shows the concavity of the functions and comparisons of the cost functions.Figure 4
The figures of the average profit functions are in descending order asZ11, Z21, Z12, Z22, Z13, and Z23.Example  C. We consider the following parametric values in appropriate units: C1=$250 per order, C2=$20/unit/year, R1=120/365year, R2=150/365year, R3=180/365year, ρ1=10%, ρ2=5%, ρ3=2%, a=100 units, D(t)=a, C4=$120/unit,  P=$180/unit,  If=0.16/$/year,  Ic=0.13/$/year, and K=500.Then the optimal solutions are{Z11=$6566.5, T*=0.36year,  t1*=0.07year}, {Z12=$6160.63, T*=0.35year,  t1*=0.07 year},  {Z13=$6004.87, T*=0.33 year,  t1*=0.07year}, {Z21=$6561.20, T*=0.33 year,  t1*=0.07year}, {Z22=$6143.73, T*=0.41year,  t1*=0.08year}, and {Z23=$5915.57, T*=0.49year,  t1*=0.10year}.From the above optimal solutions, the optimal solution is{Z11=$6566.5, T*=0.36year,  t1*=0.07year} global maximum. We obtain a closed type formula for T*, and Figure 5 shows the concavity of the functions and comparisons of the cost functions.Figure 5
The figures of the average profit functions are in descending order asZ11, Z21, Z12, Z22, Z13, and Z23.Example  D. We consider the following parametric values in appropriate units: C1=$250 per order, C2=$10/unit/month, R1=1 month, R2=2  months, R3=3  months, ρ1=20%, ρ2=15%, ρ3=10%, a=50 units, b=0.16 units, D(t)=aexp(bt), C4=$200/unit, P=$265/unit, If=0.16$/month, Ic=0.13$/month, and K=500 units.Then the optimal solutions are{Z11=$5840.84, T*=5.86  months, t1*=0.97 month}, {Z12=$5227.55, T*=4.51  months, t1*=0.66 month}, {Z13=$4712.47, T*=3.32  months, t1*=0.44 month}, {Z21=$5270.15, T*=1 month, t1*=0.11 month}, {Z22=$5081.96, T*=2  months, t1*=0.24 month}, and {Z23=$4710.04, T*=3  months, t1*=0.39 month}.From the above optimal solutions, the optimal solution is{Z11=$5840.84, T*=5.86  months, t1*=0.97 month} global maximum. The average profit function is highly non-linear. Figure 6 shows the concavity of the functions and comparisons of the cost functions.Figure 6
The figures of the average profit functions are in descending order asZ11, Z21, Z12, Z22, Z13, and Z23.Sensitivity Analysis. The sensitivity analysis of the key parameters is given in Table 2 for quadratic demand, Table 3 for linear demand, Table 4 for constant demand, and Table 5 for exponential demand. (i)
If ordering cost increases, then material handlling cost, shipping cost, and placing order’s cost increase; as a result the total relevant profit decreases.(ii)
If the unit holding cost per unit item increases, the total profit of the system decreases.(iii)
Increasing value of maximum retailer price increases the total purchasing cost of the whole system which decreases the total profit.(iv)
Increasing value of selling price per item increases the total selling price per lot which reduces the total cost of the whole system. Therefore, the total profit of the system increases.(v)
If we increase the replenishment rate then the holding cost increases which results the decreasing value of the total profit.Table 2
Sensitivity analysis for the quadratic demand pattern.
Parameters
%
Z
1
Z
2
Z
3
Quadratic demand Case1
C
0
-
50
%
2.14
2.80
3.24
-
25
%
1.02
1.31
1.49
+
25
%
-
0.95
-
1.20
-
1.35
+
50
%
-
1.86
-
2.32
-
2.59
C
H
-
50
%
8.17
6.89
6.27
-
25
%
3.34
2.89
2.66
+
25
%
-
2.50
-
2.25
-
2.11
+
50
%
-
4.50
-
4.11
-
3.87
C
R
-
50
%
18.02
25.83
29.23
-
25
%
55.72
58.33
59.46
+
25
%
-
42.46
-
47.02
-
49.56
+
50
%
-
82.40
-
92.16
-
97.71
P
-
50
%
-
137.21
-
148.03
-
154.16
-
25
%
-
70.15
-
75.12
-
77.87
+
25
%
78.03
166.01
221.24
+
50
%
54.07
223.77
221.87
K
-
50
%
4.76
3.76
3.34
-
25
%
1.24
1.05
0.95
+
25
%
-
0.66
-
0.57
0.52
+
50
%
-
1.06
-
0.92
0.85
Quadratic demand Case2
C
0
-
50
%
5.30
4.44
3.77
-
25
%
2.65
2.22
1.88
+
25
%
-
2.65
-
2.22
-
1.88
+
50
%
-
5.30
-
4.44
-
3.77
C
H
-
50
%
2.01
2.70
3.37
-
25
%
1.01
1.35
1.69
+
25
%
-
1.01
-
1.35
-
1.69
+
50
%
-
2.01
-
2.70
-
3.37
C
R
-
50
%
81.72
92.20
98.91
-
25
%
40.86
46.00
49.45
+
25
%
-
40.86
-
46.00
-
49.45
+
50
%
-
81.72
-
92.20
-
98.91
P
-
50
%
-
139.04
-
149.34
-
156.05
-
25
%
-
69.52
-
74.67
-
78.03
+
25
%
69.52
74.67
78.03
+
50
%
139.04
149.34
156.05
K
-
50
%
1.08
1.47
1.88
-
25
%
0.36
0.49
0.62
+
25
%
-
0.22
-
0.29
-
0.37
+
50
%
-
0.36
-
0.49
-
0.62Table 3
Sensitivity analysis for the linear demand pattern.
Parameters
%
Z
1
Z
2
Z
3
Linear demand Case1
C
0
-
50
%
2.67
3.32
3.75
-
25
%
1.26
1.54
1.72
+
25
%
-
1.16
-
1.40
-
1.54
+
50
%
-
2.24
-
2.69
-
2.96
C
H
-
50
%
5.84
5.29
4.97
-
25
%
2.54
2.34
2.21
+
25
%
-
2.07
-
1.95
-
1.86
+
50
%
-
3.82
-
3.63
-
3.47
C
R
-
50
%
112.63
121.07
124.66
-
25
%
47.64
51.69
53.57
+
25
%
-
41.94
-
46.70
-
49.32
+
50
%
-
82.21
-
92.11
-
97.72
P
-
50
%
-
137.74
-
148.48
-
154.55
-
25
%
-
69.87
-
74.97
-
77.78
+
25
%
77.20
80.96
82.68
+
50
%
169.90
176.81
179.62
K
-
50
%
3.14
2.80
2.60
-
25
%
0.92
0.84
0.78
+
25
%
-
0.51
-
0.47
-
0.44
+
50
%
-
0.83
-
0.77
-
0.73
Linear demand Case2
C
0
-
50
%
5.32
4.46
3.80
-
25
%
2.66
2.23
1.90
+
25
%
-
2.66
-
2.23
-
1.90
+
50
%
-
5.32
-
4.46
-
3.80
C
H
-
50
%
2.38
3.24
4.13
-
25
%
1.01
1.35
1.68
+
25
%
-
1.01
-
1.35
-
1.68
+
50
%
-
2.01
-
2.69
-
3.37
C
R
-
50
%
81.74
92.22
98.92
-
25
%
40.87
46.11
49.46
+
25
%
-
40.87
-
46.11
-
49.46
+
50
%
-
81.74
-
92.22
-
98.92
P
-
50
%
-
139.07
-
149.37
-
78.08
-
25
%
-
69.54
-
74.69
-
78.04
+
25
%
69.54
74.69
78.04
+
50
%
139.07
149.37
156.08
K
-
50
%
1.08
1.47
1.86
-
25
%
0.36
0.49
0.62
+
25
%
-
0.22
-
0.29
-
0.37
+
50
%
1.08
1.47
1.86Table 4
Sensitivity analysis for the constant demand pattern.
Parameters
%
Z
1
Z
2
Z
3
Constant demand Case1
C
0
-
50
%
6.39
7.18
7.80
-
25
%
2.86
3.18
3.42
+
25
%
−2.74
−2.73
−2.90
+
50
%
−4.68
−5.15
−5.46
C
H
-
50
%
2.35
2.40
2.37
-
25
%
1.13
1.16
3.79
+
25
%
−1.07
−1.10
−3.60
+
50
%
−2.08
−2.14
−4.60
C
R
-
50
%
82.26
92.65
93.61
-
25
%
41.13
46.31
45.47
+
25
%
−41.13
−46.30
−50.49
+
50
%
−82.25
−92.59
−98.38
P
-
50
%
−139.56
−150.26
−156.00
-
25
%
−69.82
−75.23
−78.76
+
25
%
69.92
75.51
74.24
+
50
%
139.96
151.45
159.23
K
-
50
%
1.13
1.16
1.15
-
25
%
0.37
0.38
0.37
+
25
%
−0.22
−0.22
−0.22
+
50
%
−0.36
−0.37
−0.37
Constant demand Case2
C
0
-
50
%
5.80
4.95
4.28
-
25
%
2.90
2.48
2.14
+
25
%
−2.90
−2.48
−2.14
+
50
%
−5.80
−4.95
−4.28
C
H
-
50
%
2.00
2.68
3.33
-
25
%
1.00
1.34
1.67
+
25
%
−1.00
−1.34
−1.67
+
50
%
−2.00
−2.68
−3.33
C
R
-
50
%
82.30
92.78
99.40
-
25
%
41.15
46.39
49.70
+
25
%
−41.15
−46.39
−49.70
+
50
%
82.30
−92.78
−99.40
P
-
50
%
−140.10
−150.40
−157.02
-
25
%
−70.05
−75.20
−78.51
+
25
%
70.05
75.20
78.51
+
50
%
140.10
150.40
157.02
K
-
50
%
1.00
1.34
1.67
-
25
%
0.33
0.44
0.56
+
25
%
−0.20
−0.27
−0.33
+
50
%
−0.33
−0.44
−0.56Table 5
Sensitivity analysis for the exponential demand pattern.
Parameters
%
Z
1
Z
2
Z
3
Exponential demand Case1
C
0
-
50
%
0.37
0.54
0.83
-
25
%
0.18
−0.27
0.41
+
25
%
−0.18
−0.26
−0.39
+
50
%
−0.36
−0.52
−0.77
C
H
-
50
%
—
—
—
-
25
%
18.16
13.01
9.39
+
25
%
−6.93
−5.23
−3.78
+
50
%
−10.30
−7.89
−5.80
C
R
-
50
%
—
—
—
-
25
%
—
—
—
+
25
%
−46.01
−46.01
−50.15
+
50
%
−84.03
−84.03
−95.06
P
-
50
%
—
—
—
-
25
%
−71.44
−76.39
−83.42
+
25
%
—
—
—
+
50
%
—
—
—
K
-
50
%
—
—
5.03
-
25
%
2.82
1.80
1.16
+
25
%
−1.25
−0.86
−0.57
+
50
%
−1.97
−1.37
−0.92
Exponential demand Case2
C
0
-
50
%
2.37
1.23
0.88
-
25
%
1.86
0.61
0.44
+
25
%
-
1.18
-
0.61
-
0.44
+
50
%
-
2.37
-
1.23
-
0.88
C
H
-
50
%
2.36
5.42
9.72
-
25
%
1.18
2.71
4.86
+
25
%
-
1.18
-
2.71
-
4.86
+
50
%
-
2.36
-
2.36
-
9.72
C
R
-
50
%
82.31
98.56
122.63
-
25
%
41.15
49.28
61.31
+
25
%
-
41.15
-
49.28
-
61.31
+
50
%
-
82.31
-
98.56
—
P
-
50
%
—
—
—
-
25
%
-
68.52
-
77.61
-
91.62
+
25
%
68.52
77.61
91.62
+
50
%
137.04
155.21
183.23
K
-
50
%
0.56
1.37
2.62
-
25
%
0.18
0.46
0.87
+
25
%
-
0.11
-
0.27
-
0.52
+
50
%
-
0.18
-
0.46
-
0.87
—  indicates the value of the parameter is either imaginary or the value does not exist at that point.


## 4. Concluding Remarks

In most EOQ models, the retailer pays the purchasing cost of items to the supplier as soon as the items are received. In the competitive marketing environment, the supplier offers the retailer a delay period called thetrade credit period to stimulate the retailer to buy more items. Before the end of the trade credit period, the retailer can sell his products, accumulate revenue, and earn interest. Besides it, a higher interest is charged if the payment is not settled by the end of the trade credit period. Several papers discussing the topic have been mentioned in the literature that investigates inventory problems under various conditions. Many articles in this direction consider the replenishment rate to be instantaneously infinite, but its size is finite and the lot is placed in the beginning of the cycle. Some models consider finite replenishment/supply rate for constant demand rate. Our models consider the demand of the products as an increasing function of time. From the view point of the supplier, price discount on the purchasing cost of the items by retailer is given at a different delay period to motivate the retailer to buy more. Generally, suppliers consider maximum delay period, after which they would not take any risk of getting back money from retailers. As a result, the purchasing cost C3 is infinite when R is greater than R3; that is, retailer never purchases items at infinite cost. In addition, we establish eight effective and easy-to-use lemmas and figures to help the retailer to take optimal strategy for his marketing policy. We solved the model analytically. We found global optimal solutions for different types of demand patterns. We illustrate the numerical results graphically. As far as the knowledge of authors goes, such type of model has not yet been discussed in the EOQ literature. The application of our model in a different field of inventory models where the assumptions are valid. Further extension of the model may be generalized by considering shortages, deterioration, and stochastic demand. One immediate extension of this model is to consider partial backlogging.



---

*Source: 672504-2013-08-24.xml*