Next Article in Journal
Forecasting the Volatility of Cryptocurrencies in the Presence of COVID-19 with the State Space Model and Kalman Filter
Previous Article in Journal
Adaptive Broad Echo State Network for Nonstationary Time Series Forecasting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maximum-Profit Inventory Model with Generalized Deterioration Rate

1
College of Teacher Education, Weifang University of Science and Technical, Weifang 262799, China
2
School of Intelligent Manufacturing, Weifang University of Science and Technical, Weifang 262799, China
3
School of General Studies, Weifang University of Science and Technical, Weifang 262799, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3189; https://doi.org/10.3390/math10173189
Submission received: 31 July 2022 / Revised: 23 August 2022 / Accepted: 31 August 2022 / Published: 3 September 2022

Abstract

:
We developed a maximum profit inventory model with a generalized deterioration rate where the selling rate is dependent on the inventory level that is an extension of two published papers. A complete solution structure is provided to decide the optimal solution with reasonable conditions supported by numerical examples, and then we prove that the optimal solution is independent of the demand pattern. Numerical examples are provided to illustrate our findings. In a previously published paper, three examples had symmetric conditions to decide the local maximum solution. Our approach provides a reasonable explanation for this symmetric phenomenon. Our findings will help researchers develop new inventory models in the future.

1. Introduction

Hill [1] is a pioneer, creating an inventory model with ramp-type demands where the demand is separated into two phases. In the first phase, the demand is an increasing function, and then in the second phase, the demand is a constant function. There has been a trend of considering these kinds of inventory models with minimum-cost and maximum-profit problems. For example, we list several articles by Mandal and Pal [2] on decay products; Wu et al. [3] on the backlogging rate related to the waiting time; Wu and Ouyang [4] on two replenishment policies starting with a shortage or without a shortage; Wu [5] on the deteriorated items following the Weibull distribution; and Giri et al. [6] on three-parameter Weibull distribution. Wu and Ouyang [4] examined the cyclic deterioration of items. Deng [7], to improve upon the work of Wu et al. [3], showed the existence of an optimal solution. Manna and Chaudhuri [8] developed inventory models where the deterioration is dependent on time. Deng et al. [9] improved upon the work of Mandal and Pal [2] to prove that an optimal solution exists. Wu et al. [10] examined the maximum-profit problem in inventory models with stock-dependent selling rates and ramp-type demand. They constructed two inventory systems which were dependent on the changing point of the ramp-type demand, and then examined the optimal solution for each case. In Wu et al. [10], the demand is a ramp type, and Chuang [11] generalized it to arbitrary demand where the selling rate is stock-dependent, and then Chuang [11] derived a complete solution procedure. Several related papers include: Lin et al. [12] with net present value, Mahajan and Diatha [13] with continuous compounding, Alım and Beullens [14] with a flexible delivery option, Khakzad and Gholamian [15] with inspection and advanced payment, Khan et al. [16] with an all-units discount environment and price-sensitive demand, Rezagholifam et al. [17] with non-instantaneous deteriorating items and capacity constraint, Sarkar et al. [18] with the time value of money, and Tiwari et al. [19] with trade credit. These papers are worth mentioning to point out the current research trend for this kind of inventory system. Several related papers that studied inventory models with deteriorating items are worthwhile to mention. Padiyar et al. [20] solved an inventory model of decay products under a fuzzy environment and with shortages where the demand is related to the selling price. Acevedo-Chedid et al. [21] developed a four-level inventory model, consisting of the supplier, production plants, distribution, and retailers, to examine collaborative planning for decay items. Mishra et al. [22] solved a sustainable supply chain inventory system with carbon emissions and controllable deteriorating products in a greenhouse environment. Mashud et al. [23] constructed a sustainable inventory system for imperfect items with controllable emissions and decay products to reduce carbon dioxide through green technology investments. Mishra et al. [24] examined an EOQ inventory system where the demand is dependent on the selling price and inventory level with preservation technology investment, shortages, and a controllable deterioration rate. Tiwari et al. [25] developed a two-echelon inventory system with partial trade credits for decay products with an expiration date to provide solution algorithms and sensitivity analysis. Khan et al. [26] constructed inventory models of deteriorating items with an expiration date where the demand is related to the selling price. They presented a solution algorithm to find the optimal solution. Padiyar et al. [27] considered an integrated EPQ inventory system of producer and buyer, with an imperfect production process for decay items that have two different demands: the producer with exponential demand, and the buyer with triangular demand, where the production rate is related to the demand. We also point out the following two other important papers related to inventory systems. Yadav et al. [28] examined a two-level integrated inventory system of manufacturer and retailer where the lead time is crashable with a service-level constraint, and the setup cost has a learning and forgetting effect. Navarro et al. [29] studied a maximum-profit EPQ inventory model with a supply chain of three levels, manufacturers, distribution centers, and retailers, to consider the effect of the marketing effort on the demand under multiple-item environments. We provide Table 1 below to compare our paper with the above-mentioned articles to indicate the current trend of research directions. Our findings will help researchers to realize maximum-profit inventory models.
Based on Table 1, the selling price and purchasing cost have been studied by many papers. It shows that maximum-profit models are more prevalent than minimum-cost models. Many models consider the deterioration cost, revealing that decision-makers pay attention to decay problems. Sensitivity analysis and shortages are the third and fourth most common issues in Table 1. Sensitivity analysis examines the influence of decision variables by the variation in constant parameters. Investigation of shortages can prevent stocking too many products to provide a balance between inventory and backorders. The solution method is another important research issue, showing that besides developing new inventory systems, researchers also focus on finding the optimal solution. Our paper reviews the above-mentioned six issues to indicate that it follows the main research direction.

2. Notation and Assumptions

To develop an inventory system where the demand is dependent on the inventory level, to be compatible with Wu et al. [10] and Chuang [11], we use their assumptions and notations, except the deterioration function is generated from a constant function to an arbitrary function. Moreover, we constructed an auxiliary function, m ( t 1 ) . The notations are given below:
(1)
The lead time is zero because replenishments are instantaneous with an infinite replenishment rate.
(2)
A is the ordering cost per order.
(3)
c is the purchasing cost per unit.
(4)
C d is the deterioration cost per unit.
(5)
C h is the inventory holding cost per unit, per unit of time.
(6)
C s is the shortage cost per unit, per unit of time.
(7)
s is the selling price per unit.
(8)
t 1 is the time when the inventory level drops to zero.
(9)
t 1 * is the optimal solution for the variable t 1 .
(10)
T is the finite time horizon under consideration.
(11)
α is the portion of extra selling related to inventory level. We define those sales related to inventory level as “extra sells”. For a technical reason, we assume that α ( s + C d ) < C h , which will be explained in Section 3.
(12)
μ is the time that the ramp-type demand changes from an increasing phase to a constant phase.
(13)
θ ( t ) is the deteriorated fraction of the on-hand inventory deterioration per unit of time.
(14)
I ( t ) is the on-hand inventory level at the time t over the ordering cycle [ 0 , T ] .
(15)
The shortage is allowed and fully backlogged.
(16)
Without referring to extra sales, the demand rate R(t) is a positive continuous function.
(17)
D ( t ) is the consumption sale rate (the actual rate), which is stock-dependent with D ( t ) = { R ( t ) + α I ( t )   , I ( t ) > 0 , R ( t )   , I ( t ) 0 ,   where α is a positive constant.
(18)
Z ( t 1 ) is the total profit that consists of sale revenue, order cost, purchase cost, holding cost, deterioration cost, and shortage cost.
(19)
m ( t 1 ) is an auxiliary function such that m ( t 1 ) = 0 and d d t 1 Z ( t 1 ) = 0 have the same solution.
We recall that (a) t 1 is the time when the inventory level drops to zero, and (b) μ is the time when the ramp-type demand changes from an increasing phase to a constant phase. Depending on the relation between t 1 and μ , Wu et al. [10] divided their model into two cases: Case 1, t 1 μ , and Case 2, t 1 μ . The deterioration rate is a constant function in Wu et al. [10]. Hence, there were two objective functions in Wu et al. [10] with different domains. Chuang [11] studied the same inventory model as Wu et al. [10], and he also considered the deterioration cost. Chuang [11] did not write down the explicit formula for the ramp-type demand, such that there is only one objective function in Chuang [11]. Our paper extended the work of Wu et al. [10] and Chuang [11] considering a deterioration rate from a constant function to a non-negative function. We followed the solution procedure of Chuang [11]. We focus on (i) the demand function and (ii) the deterioration function, both in an implicit setting, such that there is only one objective function in our paper.

3. Our proposed Inventory Model

We recall the inventory systems with ramp-type demand. These kinds of inventories were first proposed by Hill [1], and then further investigated by Mandal and Pal [2], Deng et al. [9], Wu and Ouyang [4], and Wu et al. [10]. Replenishment occurs at the time t = 0 when the inventory level begins at its maximum level. From t = 0 to t 1 , the inventory level decreases owing to deterioration and demand, R ( t ) . At t 1 , the inventory level drops to zero, after which shortages occurred for the period ( t 1 , T ) , and all of the unsatisfied demand during the shortage period ( t 1 , T ) is completely backlogged. Since the actual selling rate is dependent on the inventory level, the analytical model is described by the following two differential equations:
d d t I ( t ) + ( θ ( t ) + α )   I ( t ) = R ( t ) , 0 < t < t 1
and
d d t I ( t ) = R ( t ) , t 1 < t < T
We try to solve the differential system of Equations (1) and (2), then
I ( t ) = t t 1 R ( y ) e t y ( θ ( x ) + α ) d x d y
for 0 t t 1 , and
I ( t ) = t t 1 R ( x ) d x
for t 1 t T .
The purchased quantity at t = 0 plus the backorder quantity at t = T is the sum of I ( 0 ) and I ( T ) ; then, the purchased quantity will be
I ( 0 ) I ( T ) = 0 t 1 R ( y ) e 0 y ( θ ( x ) + α ) d x d y + t 1 T R ( x ) d x
The holding cost during the period [ 0 , t 1 ] is evaluated by changing the order of integration
C h 0 t 1 I ( t ) d t = C h 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y
The number of deteriorated items during the period, [ 0 , t 1 ] , is evaluated
I ( 0 ) 0 t 1 D ( x ) d x = 0 t 1 R ( y ) ( e 0 y ( θ ( x ) + α ) d x 1 ) d y α 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y
The shortage cost during the period [ t 1 , T ] is computed by changing the order of integration
C s t 1 T I ( t ) d t = C s t 1 T ( T x ) R ( x ) d x
We evaluate the sale revenue per replenishment to derive
s ( 0 t 1 D ( t ) d t I ( T ) ) = s α 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y + s 0 T R ( x ) d x
We obtain the total profit per replenishment cycle to denote it as the sale revenue minus the total cost, where the total cost is consistent with deterioration cost, shortage cost, holding cost, purchasing cost, and ordering cost.
Z ( t 1 ) = 1 T { s α 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y + s 0 T R ( x ) d x A
c [ 0 t 1 R ( y ) e 0 y ( θ ( x ) + α ) d x d y + t 1 T R ( x ) d x ] C h 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y C s t 1 T ( T x ) R ( x ) d x C d [ 0 t 1 R ( y ) ( e 0 y ( θ ( x ) + α ) d x 1 ) d y α 0 t 1 R ( y ) 0 y e t y ( θ ( x ) + α ) d x d t d y ] }
Based on Equation (10), we derive
d d t 1 Z ( t 1 ) = R ( t 1 ) T { ( s α C h + α C d ) 0 t 1 e t t 1 ( θ ( x ) + α ) d x d t ( c + C d ) ( e 0 t 1 ( θ ( x ) + α ) d x 1 ) + C s ( T t 1 ) }
Based on Equation (11), we assume an auxiliary function, say m ( t 1 ) , to simplify the later discussion, as follows
m ( t 1 ) = ( s α C h + α C d ) 0 t 1 e t t 1 ( θ ( x ) + α ) d x d t ( c + C d ) ( e 0 t 1 ( θ ( x ) + α ) d x 1 ) + C s ( T t 1 )
such that
d d t 1 Z ( t 1 ) = R ( t 1 ) T m ( t 1 )
for 0 t 1 T . Hence, d d t 1 Z ( t 1 ) and m ( t 1 ) have the same roots.
From Equation (12), it yields
m ( t 1 ) = ( s α C h + α C d ) ( 1 + ( θ ( t 1 ) + α ) 0 t 1 e t t 1 ( θ ( x ) + α ) d x d t ) ( c + C d ) ( θ ( t 1 ) + α ) e 0 t 1 ( θ ( x ) + α ) d x C s
We rewrite Equation (14) as
m ( t 1 ) = α ( s + C d ) C h C s + ( θ ( t 1 ) + α ) { ( s α + α C d C h ) 0 t 1 e t t 1 ( θ ( x ) + α ) d x d t ( c + C d ) e 0 t 1 ( θ ( x ) + α ) d x }
We also compute
m ( 0 ) = C s T > 0
to
m ( T ) = ( α ( s + C d ) C h ) 0 t 1 e t t 1 ( θ ( x ) + α ) d x d t ( c + C d ) ( e 0 t 1 ( θ ( x ) + α ) d x 1 ) < 0
From Equations (15) and (17), we notice that without giving an explicit expression of the deterioration function θ ( t ) , researchers cannot compute those integrations related to d d t m ( t ) . The alternative approach is to consider numerical support to face the monotonic problem of d d t m ( t ) and m ( T ) < 0 .
We can say that the most feasible way to prove that m ( t 1 ) = 0 has a unique solution for 0 < t 1 < T is to show that m ( t 1 ) is a strictly decreasing function. Hence, we try to find conditions to guarantee d d t 1 m ( t 1 ) < 0 and m ( T ) < 0 .
We check the coefficients in Equations (15) and (17) to find that
α ( s + C d ) C h < 0
and
α ( s + C d ) C h C s < 0
Owing to the inequality in Equation (18), the inequality in Equation (19) is implied, so we will try to verify that the inequality in Equation (18) is valid.
Consequently, we can derive our desired result of the existence and uniqueness of only one root for m ( t 1 ) , which is the optimal solution for Z ( t 1 ) , which is the unique solution for d d t 1 Z ( t 1 ) = 0 .
Since our proposed model is new, in the literature, there are no inventory models that have considered those parameters altogether. Hence, we have to break our condition into two parts:
α s < 0.5 C h
and
α C d < 0.5 C h
The reason we separate Equation (18) into two parts is explained below.
We recall that in the numerical examples of Wu et al. [10] and Chuang [11], α = 0.01 , s = 20 C h = 3 , such that
α s = 0.2 < 0.5 C h = 1.5
Therefore, our assertion of Equation (21) is supported by the numerical examples of Wu et al. [10] and Chuang [11].
On the other hand, Wu et al. [10] did not include the deterioration cost per unit, C d , in their inventory models. Chuang [11] has already developed an inventory model with the deterioration cost. However, Chuang [11] did not inform us of his deterioration cost. In comparison with Wu et al. [10], in the numerical examples of Chuang [11], he still assumed that C d = 0. Hence, we must refer to another paper, Lin [30], and his numerical example, C d = 3 and C h = 10 , to imply that
α C d = 0.03 < 0.5 C h = 5
Thus, our claim of Equation (22) is confirmed by the numerical example of Lin [30].
Hence, our extra condition of Equation (18) is supported by the numerical examples of Wu et al. [10], Lin [30], and Chuang [11].
Since exponential functions are positive, we imply that those two integrations of exponential functions are also positive. The two coefficients before integrations are negative. We combine our observations to derive
m ( t 1 ) < 0
We show that m ( t 1 ) is a decreasing function from m ( 0 ) > 0 to m ( T ) < 0 .
Hence, there is a unique point, say t 1 * , that is satisfying m ( t 1 * ) = 0 . We recall that m ( t 1 ) and Z ( t 1 ) have the same sign, such that
Z ( t 1 ) > 0 ,   for   0 < t 1 < t 1 *
and
Z ( t 1 ) < 0 ,   for   t 1 * < t 1 < T
We prove that t 1 * is the optimal solution for our maximum-profit inventory model. Equation (12) reveals that the demand rate, R ( t 1 ) , is not included in the expression of m ( t 1 ) , so the optimal solution is independent of the demand rate. This finding is consistent with the results of Lin [30] and Hung [31], which are minimum-cost inventory models, that is, the optimal solution is independent of the demand pattern.

4. Numerical Examples

We recall numerical examples in Chuang [11] with the following data: the setup (the ordering) cost A = 50 , the holding cost C h = 3 , the shortage cost C s = 5 , the purchasing cost c = 15 , the actual sale rate D 0 = 400 , the selling price s = 20 , the planning horizon T = 1 , the ratio of extra selling related to the inventory level α = 0.01 , and the deterioration rate θ = 0.05 . In Chuang [11], his inventory model did not include the deterioration cost. Hence, we refer to Lin [30] to adopt C d = 3 . To check the influence of the variation for parameters, we examine a detailed sensitivity analysis to alter the values from the decreasing percentages of 30%, 20%, and 10%, to increasing percentages of 10%, 20%, and 30% for the setup cost A, the holding cost C h , the shortage cost C s , the purchasing cost c, the selling price s, the planning horizon T, the ratio of extra selling related to the inventory level α , and the deterioration function θ . We list the results for variation in the setup cost in Table 2.
Our findings in Table 2 show that the optimal solution of t 1 * is independent of the setup cost. Our findings for variation in the holding cost are listed in Table 3.
Our results in Table 3 show that the optimal solution of t 1 * has a negative relationship with the holding cost. Our results for variation in the shortage cost are presented in Table 4.
Our findings in Table 4 show the optimal solution t 1 * has a positive relationship with the shortage cost. Our findings for variation in the purchasing cost are shown in Table 5.
Our results in Table 5 they assert that the optimal solution t 1 * has a negative relationship with the purchasing cost. Our results for variation in the selling price are listed in Table 6.
Our results in Table 6 denote that the optimal solution t 1 * has a positive relationship with the selling price. Our findings for variation in the planning horizon are expressed in Table 7.
Our findings in Table 7 imply that the optimal solution t 1 * has a positive relationship with the planning horizon. Our results for variation in the ratio of extra selling are denoted in Table 8.
Our results in Table 8 claim that the optimal solution t 1 * has a positive relationship with the ratio of extra selling. Our findings for the results for variation in the deterioration function are shown in Table 9.
Our findings in Table 9 point out that the optimal solution t 1 * has a negative relationship with the deterioration function. Our results for variation in the deterioration cost are expressed in Table 10.
Our results in Table 10 claim that the optimal solution t 1 * has a negative relation with the deterioration cost.

5. Application of Our Approach

When θ ( t ) is degenerated to a constant function, denoted as θ , we can rewrite Equation (12) as
m ( t 1 ) = [ α ( s c ) C h θ ( c + C d ) θ + α ] ( e ( θ + α ) t 1 1 ) + C s ( T t 1 )
which is the same result as the criterion of Chuang [11] for the optimal solution, revealing that Chuang [11] is a special case of our model.
To be comparable with Wu et al. [10], we consider the same numerical examples as their Example 1,   A = 50 , C h = 3 , C s = 5 , D 0 = 400 , T = 1 , c = 15 , s = 20 , α = 0.01 , and θ = 0.05 .
Wu et al. [10] did not consider the deterioration cost. In comparing our work with that of Wu et al. [10], we suppose C d = 0 . in the following discussion.
Wu et al. [10] studied a ramp-type demand as
D ( t ) = { D 0 t ,   0 t μ , D 0 μ ,   μ t T .
Depending on the relationship between μ and t 1 , Wu et al. [10] partitioned their model into two cases: Case A: t 1 μ , and Case B: μ t 1 . Consequently, they constructed their profit function Z 1 ( t 1 ) and Z 2 ( t 1 ) for Cases A, and B, respectively.
For Case A, Wu et al. [10] obtained
α ( s c ) ( C h + θ c ) θ + α [ e ( θ + α ) t 1 1 ] + C s ( T t 1 ) = 0
to decide their optimal solution for μ t 1 T .
For Case B, Wu et al. [10] derived the same result as Equation (29) to decide their optimal solution for 0 t 1 μ .
Wu et al. [10] did not explain why their results for Cases A and B are identical under two different domains. They examined three examples for their inventory model with (i) Example 1, μ = 0.4 , (ii) Example 2, μ = 0.5953 , and (iii) Example 3, μ = 0.6 .
For Example 1, Wu et al. [10] claimed that the local maximum solution for Z 1 ( t 1 ) is μ and the local maximum solution for Z 2 ( t 1 ) is t 1 * = 0.5953 , with
Z 2 ( t 1 * = 0.5953 ) > Z 2 ( μ = 0.4 ) = Z 1 ( μ = 0.4 )
For Example 2, Wu et al. [10] found that the local maximum solution for Z 1 ( t 1 ) is μ and the local maximum solution for Z 2 ( t 1 ) is μ , with
Z 2 ( μ = 0.5953 ) = Z 1 ( μ = 0.5953 )
For Example 3, Wu et al. [10] found that the local maximum solution for Z 1 ( t 1 ) is t 1 * = 0.5953 and the local maximum solution for Z 2 ( t 1 ) is μ , with
Z 2 ( μ = 0.6 ) = Z 1 ( μ = 0.6 ) < Z 1 ( t 1 * = 0.5953 )
Wu et al. [10] were not aware that their three examples had the same global maximum point, t 1 * = 0.5953 . Hence, Wu et al. [10] did not know that the maximum solution for these kinds of inventory models is independent of the demand type such that there are two divisions: the first one is to construct Z 1 ( t 1 ) and Z 2 ( t 1 ) for Cases A and B, respectively, concerning the relation between t 1 and μ , which is unnecessary. The second one is to develop three examples, based on the relation t 1 * and μ , which is redundant. Hence, applying our approach, researchers can solve this kind of inventory model without those lengthy derivations proposed by Wu et al. [10].
At last, but not least, Wu et al. [10] and Chuang [11] mentioned that
t 1 * = 0.5953
In the following, we list four related values to clearly indicate the values of m ( t 1 ) decreasing from positive to negative. Hence, our solution is the right result.
On the other hand, we derive that m ( 0.570511 ) = 9.922 × 10 6 , m ( 0.570512 ) = 1.093 × 10 6 , m ( 0.570513 ) = 7.736 × 10 6 , and m ( 0.5953 ) = 0.219 to indicate that the maximum solution should be revised as
t 1 * = 0.570512
to show that the work of Wu et al. [10] and Chuang [11] needs revisions.
We can also apply our analytical approach to consider the deterioration rate function as a linear function such that θ ( t ) = a + b t , and then rewrite Equation (12) as follows:
m ( t 1 ) = ( s α C h + α C d ) 0 t 1 e t t 1 ( a + α + b x ) d x d t ( c + C d ) ( e 0 t 1 ( a + α + b x ) d x 1 ) + C s ( T t 1 ) = ( s α C h + α C d ) e ( a + α ) t 1 + ( b t 1 2 / 2 ) 0 t 1 e ( a + α ) t ( b t 2 / 2 ) d t ( c + C d ) ( e ( a + α ) t 1 + ( b t 1 2 / 2 ) 1 ) + C s ( T t 1 )
Under the condition of θ ( t ) = a + b t , solving m ( t 1 ) = 0 will be a challenging problem because, in the integration of Equation (35), the following quadratic form denoted as E, with
E = 0 t 1 e ( a + α ) t ( b t 2 / 2 ) d t
appeared.
For a linear term in the exponential function, researchers can derive that
0 t 1 e ( a + α ) t d t = 1 a + α [ 1 e ( a + α ) t 1 ]
to indicate that when the deterioration function is a constant term as   θ ( t ) = a , then the integration can be directly solved as Equation (37).
When we extend our inventory model for the deterioration from a constant term to a linear function as θ ( t ) = a + b t , the integration problem of Equation (36) cannot be directly solved to an exact result.
Therefore, for a quadratic form in the exponential function of Equation (36), we change the expression of Equation (36) to
0 t 1 e ( a + α ) t ( b t 2 / 2 ) d t = 2 b e ( a + α ) 2 / 2 b ( a + α ) / 2 b ( b / 2 ) ( t 1 + ( ( a + α ) / b ) ) e y 2 dy
with
y = b 2 ( t + a + α b )
and
b 2 t 2 ( a + α ) t = ( b 2 ( t + a + α b ) ) 2 + ( a + α ) 2 2 b
For the improper integration of 0 e y 2 dy , researchers know that
0 e y 2 dy = π 2
However, for another definite integration, for example, that denoted as F, with
F = 3 5 e y 2 dy ,  
researchers cannot find the exact solution for Equation (42). To the best of our knowledge, many practitioners have applied Taylor’s series expansion of the exponential function to derive an infinite series for the definite integration of Equation (42)
3 5 e y 2 dy = 3 5 k = 0 ( y 2 ) k k ! dy = k = 0 ( 1 ) k k ! ( 5 2 k + 1 3 2 k + 1 2 k + 1 )
Depending on the accuracy, practitioners selected a number, denoted as “m”, to accept that
3 5 e y 2 dy k = 0 m ( 1 ) k k ! ( 5 2 k + 1 3 2 k + 1 2 k + 1 )
where the estimation error is less than 5 2 m + 3 3 2 m + 3 ( m + 1 ) !   ( 2 m + 3 ) , because the infinite series on the right-hand side of Equation (44) is an alternative series.
The second way to solve the value of Equation (42) is to consider the well-known fact that finite integration is the limit of the Riemann sum. It shows that
3 5 e y 2 dy = lim n k = 1 n e ( 3 + 2 k n ) 2 2 n .  
Based on the above discussion, in the beginning, we will adopt the second approach to solve m ( t 1 ) = 0 , when θ ( t ) = a + b t .
We construct an example for a linear deterioration function with a = 0.05 and b = 0.1 to extend the deterioration function to a linear form such that A = 50 , C d = 3 , C h = 3 , C s = 5 , D 0 = 400 , T = 1 , c = 15 , s = 20 , α = 0.01 , θ = 0.05 , and θ ( t ) = 0.05 + 0.1 t .
For a given value of t 1 , we uniformly partition the integration interval [ 0 , t 1 ] into twenty subintervals, and then we estimate m ( t 1 ) using the following approximation:
m ( t 1 ) = ( s α C h + α C d ) e ( a + α ) t 1 + ( b t 1 2 / 2 ) t 1 20 k = 1 20 e ( a + α ) k 20 t 1 b 2 ( k 20 t 1 ) 2 ( c + C d ) ( e ( a + α ) t 1 + ( b t 1 2 / 2 ) 1 ) + C s ( T t 1 )
We begin to search for the optimal solution, t 1 * , which is the solution of m ( t 1 ) = 0 , with the knowledge that m ( 0 ) > 0 , m ( T = 1 ) < 0 , and m ( t 1 ) is a decreasing function.
Many methods can locate the value of the optimal solution, t 1 * . In the following, we will use a direct numerical method to list values of m ( t 1 ) that will clearly reveal that m ( t 1 ) is a decreasing function from a positive value to a negative value, and then we can find the solution of m ( t 1 ) = 0 .
For t 1 = 0 , 0.1 , , 1 , we list the estimation result of m ( t 1 ) in the following Table 11.
Based on Table 11, we know that the optimal solution, t 1 * , satisfies
0.5 < t 1 * < 0.6
Hence, for t 1 = 0.50 ,   0.51 , ,   0.6 , we list the estimation result of m ( t 1 ) in Table 12.
Based on Table 12, we know that the optimal solution, t 1 * , satisfies
0.53 < t 1 * < 0.54
Hence, for t 1 = 0.530 ,   0.531 , ,   0.540 , we list the estimation result of m ( t 1 ) in Table 13.
Based on Table 13, we know that the optimal solution, t 1 * , satisfies
0.530 < t 1 * < 0.531
Hence, for t 1 = 0.5300 ,   0.5301 , ,   0.5310 , we list the estimation result of m ( t 1 ) in Table 14.
Based on Table 14, we know that the optimal solution, t 1 * , satisfies
0.5301 < t 1 * < 0.5302
We compare
m ( 0.5301 ) = 5.113 × 10 4
and
m ( 0.5302 ) = 4.949 × 10 4 ,  
to select
t 1 * = 0.5302

6. Conclusions

In this paper, we constructed a maximum-profit inventory model to show that the optimal solution is independent of the demand such that the deterioration rate is generalized to the general expression. We show that without writing the implicit expression of demand function and deterioration function, through our approach, complicated solution procedures with several different objective functions can be avoided. Our paper is not only an extension of the work of Wu et al. [10] and Chuang [11] but also provides a complete solution structure to find the optimal solution.

Author Contributions

Conceptualization, Y.-L.W.; methodology, Y.-L.W.; software, Y.-L.W. and M.-L.C.; validation, P.J.; formal analysis, Y.-L.W.; investigation, Y.-L.W.; resources, M.-L.C.; data curation, M.-L.C.; writing—original draft preparation, Y.-L.W.; writing—review and editing, Y.-L.W.; visualization, M.-L.C.; supervision, Y.-L.W.; project administration, Y.-L.W.; funding acquisition, Y.-L.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hill, R.M. Inventory models for increasing demand followed by level demand. J. Oper. Res. Soc. 1995, 46, 1250–1259. [Google Scholar] [CrossRef]
  2. Mandal, B.; Pal, A.K. Order level inventory system with ramp type demand rate for deteriorating items. J. Interdiscipl. Math. 1998, 1, 49–66. [Google Scholar] [CrossRef]
  3. Wu, J.W.; Lin, C.; Tan, B.; Lee, W.C. An EOQ model with ramp type demand rate for items with Weibull deterioration. Int. J. Info. Manag. Sci. 1999, 10, 41–51. [Google Scholar]
  4. Wu, K.S.; Ouyang, L.Y. A replenishment policy for deteriorating items with ramp type demand rate. Proc.-Natl. Sci. Counc. Repub. China Part A Phys. Sci. Eng. 2000, 24, 279–286. [Google Scholar]
  5. Wu, K.S. An EOQ inventory model for items with Weibull distribution deterioration, ramp type demand rate and partial backlogging. Prod. Plan. Control 2001, 12, 787–793. [Google Scholar] [CrossRef]
  6. Giri, B.C.; Jalan, A.K.; Chaudhuri, K.S. Economic Order Quantity model with Weibull deterioration distribution, shortage and ramp-type demand. Int. J. Syst. Sci. 2003, 34, 237–243. [Google Scholar] [CrossRef]
  7. Deng, P.S. Improved inventory models with ramp type demand and Weibull deterioration. Int. J. Inf. Manag. Sci. 2005, 16, 79–86. [Google Scholar]
  8. Manna, S.K.; Chaudhuri, K.S. An EOQ model with ramp type demand rate, time dependent deterioration rate, unit production cost and shortages. Eur. J. Oper. Res. 2006, 171, 557–566. [Google Scholar] [CrossRef]
  9. Deng, P.S.; Lin, R.; Chu, P. A note on the inventory models for deteriorating items with ramp type demand rate. Eur. J. Oper. Res. 2007, 178, 112–120. [Google Scholar] [CrossRef]
  10. Wu, K.S.; Ouyang, L.Y.; Yang, C.T. Retailer’s optimal ordering policy for deteriorating items with ramp-type demand under stock-dependent consumption rate. Int. J. Inf. Manag. Sci. 2008, 19, 245–262. [Google Scholar]
  11. Chuang, J.P.C. Inventory model with stock dependent selling rate. Afr. J. Bus. Manag. 2012, 6, 780–786. [Google Scholar]
  12. Lin, R.; Lin, J.S.J.; Chen, K.; Julian, P.C. Note on inventory model with net present value. J. Interdiscip. Math. 2007, 10, 587–592. [Google Scholar] [CrossRef]
  13. Mahajan, S.; Diatha, K. Minimizing the Discounted Average Cost Under Continuous Compounding in the EOQ Models with a Regular Product and a Perishable Product. Am. J. Oper. Manag. Inf. Syst. 2018, 3, 52–60. [Google Scholar] [CrossRef]
  14. Alım, M.; Beullens, P. Joint inventory and distribution strategy for online sales with a flexible delivery option. Int. J. Prod. Econ. 2020, 222, 107487. [Google Scholar] [CrossRef]
  15. Khakzad, A.; Gholamian, M.R. The effect of inspection on deterioration rate: An inventory model for deteriorating items with advanced payment. J. Clean. Prod. 2020, 254, 120117. [Google Scholar] [CrossRef]
  16. Khan, M.A.A.; Ahmed, S.; Babu, M.S.; Sultana, N. Optimal lot-size decision for deteriorating items with price-sensitive demand, linearly time-dependent holding cost under all-units discount environment. Int. J. Syst. Sci. Oper. Logist. 2020, 9, 61–74. [Google Scholar] [CrossRef]
  17. Rezagholifam, M.; Sadjadi, S.J.; Heydari, M.; Karimi, M. Optimal pricing and ordering strategy for non-instantaneous deteriorating items with price and stock sensitive demand and capacity constraint. Int. J. Syst. Sci. Oper. Logist. 2020, 9, 121–132. [Google Scholar] [CrossRef]
  18. Sarkar, S.; Giri, B.C.; Sarkar, A.K. A vendor-buyer inventory model with lot-size and production rate dependent lead time under time value of money. RAIRO-Oper. Res. 2020, 54, 961–979. [Google Scholar] [CrossRef]
  19. Tiwari, S.; Cárdenas-Barrón, L.E.; Shaikh, A.A.; Goh, M. Retailer’s optimal ordering policy for deteriorating items under order-size dependent trade credit and complete backlogging. Comput. Ind. Eng. 2020, 139, 105559. [Google Scholar] [CrossRef]
  20. Padiyar, S.V.S.; Singh, S.R.; Punetha, N. Inventory system with price dependent consumption for deteriorating items with shortages under fuzzy environment. Int. J. Sust. Agric. Manag. Inf. 2021, 7, 218–231. [Google Scholar] [CrossRef]
  21. Acevedo-Chedid, J.; Salas-Navarro, K.; Ospina-Mateus, H.; Villalobo, A.; Sana, S.S. Production system in a collaborative supply chain considering deterioration. Int. J. Appl. Comput. Math. 2021, 7, 69. [Google Scholar] [CrossRef]
  22. Mishra, U.; Mashud, A.H.M.; Tseng, M.-L.; Wu, J.-Z. Optimizing a sustainable supply chain inventory model for controllable deterioration and emission rates in a greenhouse farm. Mathematics 2021, 9, 495. [Google Scholar] [CrossRef]
  23. Mashud, A.H.M.; Roy, D.; Daryanto, Y.; Ali, M.H. A sustainable inventory model with imperfect products, deterioration, and controllable emissions. Mathematics 2020, 8, 2049. [Google Scholar] [CrossRef]
  24. Mishra, U.; Cárdenas-Barrón, L.E.; Tiwari, S.; Shaikh, A.A.; Graza, G.T. An inventory model under price and stock dependent demand for controllable deterioration rate with shortages and preservation technology investment. Ann. Oper. Res. 2017, 254, 165–190. [Google Scholar] [CrossRef]
  25. Tiwari, S.; Cárdenas-Barrón, L.E.; Goh, M.; Shaikh, A.A. Joint pricing and inventory model for deteriorating items with expiration dates and partial backlogging under two-level partial trade credits in supply chain. Int. J. Prod. Econ. 2018, 200, 16–36. [Google Scholar] [CrossRef]
  26. Khan, A.A.; Shaikh, A.A.; Panda, G.C.; Konstantaras, I.; Taleizadeh, A.A. Inventory system with expiration date: Pricing and replenishment decisions. Comput. Ind. Eng. 2019, 132, 232–247. [Google Scholar] [CrossRef]
  27. Padiyar, S.V.S.; Kuraie, V.C.; Bhagat, N.; Singh, S.R.; Chaudhary, R. An integrated inventory model for imperfect production process having preservation facilities under fuzzy and inflationary environment. Int. J. Math. Model. Numer. Optim. 2022, 12, 252–286. [Google Scholar] [CrossRef]
  28. Yadav, D.; Singh, S.R.; Kumar, S.; Cárdenas-Barrón, L.E. Manufacturer-retailer integrated inventory model with controllable lead time and service level constraint under the effect of learning-forgetting in setup cost. Sci. Iran. 2022, 29, 800–815. [Google Scholar]
  29. Navarro, K.S.; Chedid, J.A.; Florez, W.F.; Mateus, H.O.; Cárdenas-Barrón, L.E.; Sana, S.S. A collaborative EPQ inventory model for a three-echelon supply chain with multiple products considering the effect of marketing effort on demand. J. Ind. Manag. Optim. 2020, 16, 1613–1633. [Google Scholar] [CrossRef]
  30. Lin, S.W. Inventory models with managerial policy independent of demand. Eur. J. Oper. Res. 2011, 211, 520–524. [Google Scholar] [CrossRef]
  31. Hung, K.C. An inventory model with generalized type demand, deterioration and backorder rates. Eur. J. Oper. Res. 2011, 208, 239–242. [Google Scholar] [CrossRef]
Table 1. Comparison of our present work and some related papers.
Table 1. Comparison of our present work and some related papers.
PaperDIIPDDPCSPSASSMSDDSC
[10]X XX XXX
[11]X XX XXX
[20]X X XXX X
[21]XX XX X X
[22]X X XXXX X
[23]XXXXXX X
[24]X X XXXXX
[25]X XXXXXX X
[26]X XXXXXX
[27]XX X XX X
[28] XX X X
[29] X XXX X
oursX XXXXXX
D deterioration; II imperfect items; PDD price-dependent demand; PC purchasing cost; SP selling price; SA sensitivity analysis; S shortages; SM solution method; SDD stock-dependent demand; SC supply chain.
Table 2. Variation in the setup cost A.
Table 2. Variation in the setup cost A.
A−30%−20%−10%base+10%+20%+30%
t 1 * 0.56080.56080.56080.56080.56080.56080.5608
Table 3. Variation in the holding cost C h .
Table 3. Variation in the holding cost C h .
C h −30%−20%−10%base+10%+20%+30%
t 1 * 0.62450.60170.58050.56080.54240.52520.5090
Table 4. Variation in the shortage cost C s .
Table 4. Variation in the shortage cost C s .
C s −30%−20%−10%base+10%+20%+30%
t 1 * 0.47260.50580.53490.56080.58400.60480.6236
Table 5. Variation in the purchasing cost c.
Table 5. Variation in the purchasing cost c.
c −30%−20%−10%base+10%+20%+30%
t 1 * 0.57850.57250.56660.56080.55520.54970.5442
Table 6. Variation in the selling price s.
Table 6. Variation in the selling price s.
s −30%−20%−10%base+10%+20%+30%
t 1 * 0.55700.55830.55960.56080.56210.56340.5647
Table 7. Variation in the planning horizon T.
Table 7. Variation in the planning horizon T.
T−30%−20%−10%base+10%+20%+30%
t 1 * 0.39340.44930.50510.56080.61640.67200.7275
Table 8. Variation in the ratio of extra selling α .
Table 8. Variation in the ratio of extra selling α .
θ −30%−20%−10%base+10%+20%+30%
t 1 * 0.56010.56030.56060.56080.56110.56130.5616
Table 9. Variation in the deterioration function θ .
Table 9. Variation in the deterioration function θ .
θ −30%−20%−10%base+10%+20%+30%
t 1 * 0.57960.57320.56690.56080.55480.54890.5432
Table 10. Variation in the deterioration cost C d .
Table 10. Variation in the deterioration cost C d .
C d −30.−20%−10%base+10%+20%+30%
t 1 * 0.56370.56270.56180.56080.55990.55890.5580
Table 11. The estimation result for m ( t 1 ) with t 1 = 0 , 0.1 , , 1 .
Table 11. The estimation result for m ( t 1 ) with t 1 = 0 , 0.1 , , 1 .
t 1 00.10.20.30.40.5
m ( t 1 ) 5.0004.1053.1882.2501.2880.302
t 1 0.60.70.80.91
m ( t 1 ) −0.709−1.747−2.814−3.910−5.038
Table 12. The estimation result for m ( t 1 ) , with t 1 = 0.50 ,   0.51 , ,   0.6 .
Table 12. The estimation result for m ( t 1 ) , with t 1 = 0.50 ,   0.51 , ,   0.6 .
t 1 0.500.510.520.530.540.55
m ( t 1 ) 0.3020.2020.102 1.518 × 10 3 −0.099−0.200
t 1 0.560.570.580.590.60
m ( t 1 ) −0.302−0.403−0.505−0.607−0.709
Table 13. The estimation result for m ( t 1 ) , with t 1 = 0.530 ,   0.531 , ,   0.540 .
Table 13. The estimation result for m ( t 1 ) , with t 1 = 0.530 ,   0.531 , ,   0.540 .
t 1 0.5300.5310.5320.5330.5340.535
m ( t 1 ) 1.518 × 10 3 8.54 × 10 3 −0.019−0.029−0.039−0.049
t 1 0.5360.5370.5380.5390.540
m ( t 1 ) −0.059−0.069−0.079−0.089−0.099
Table 14. The estimation result for m ( t 1 ) , with t 1 = 0.5300 ,   0.5301 , ,   0.5310 .
Table 14. The estimation result for m ( t 1 ) , with t 1 = 0.5300 ,   0.5301 , ,   0.5310 .
t 1 0.53000.53010.53020.53030.53040.5305
m ( t 1 ) 1.518 × 10 3 5.11 × 10 4 4.95 × 10 4 1.50 × 10 3 2.50 × 10 3 3.51 × 10 3
t 1 0.53060.53070.53080.53090.5310
m ( t 1 ) 4.52 × 10 3 5.53 × 10 3 6.53 × 10 3 7.54 × 10 3 8.54 × 10 3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.-L.; Chen, M.-L.; Julian, P. Maximum-Profit Inventory Model with Generalized Deterioration Rate. Mathematics 2022, 10, 3189. https://doi.org/10.3390/math10173189

AMA Style

Wang Y-L, Chen M-L, Julian P. Maximum-Profit Inventory Model with Generalized Deterioration Rate. Mathematics. 2022; 10(17):3189. https://doi.org/10.3390/math10173189

Chicago/Turabian Style

Wang, Yu-Lan, Ming-Li Chen, and Peterson Julian. 2022. "Maximum-Profit Inventory Model with Generalized Deterioration Rate" Mathematics 10, no. 17: 3189. https://doi.org/10.3390/math10173189

APA Style

Wang, Y. -L., Chen, M. -L., & Julian, P. (2022). Maximum-Profit Inventory Model with Generalized Deterioration Rate. Mathematics, 10(17), 3189. https://doi.org/10.3390/math10173189

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

Article Metrics

Back to TopTop