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.
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
when the inventory level begins at its maximum level. From
to
, the inventory level decreases owing to deterioration and demand,
. At
, the inventory level drops to zero, after which shortages occurred for the period
, and all of the unsatisfied demand during the shortage period
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:
and
We try to solve the differential system of Equations (1) and (2), then
for
, and
for
.
The purchased quantity at
plus the backorder quantity at
is the sum of
and
; then, the purchased quantity will be
The holding cost during the period
is evaluated by changing the order of integration
The number of deteriorated items during the period,
, is evaluated
The shortage cost during the period
is computed by changing the order of integration
We evaluate the sale revenue per replenishment to derive
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.
Based on Equation (10), we derive
Based on Equation (11), we assume an auxiliary function, say
, to simplify the later discussion, as follows
such that
for
Hence,
and
have the same roots.
From Equation (12), it yields
We rewrite Equation (14) as
From Equations (15) and (17), we notice that without giving an explicit expression of the deterioration function , researchers cannot compute those integrations related to . The alternative approach is to consider numerical support to face the monotonic problem of and .
We can say that the most feasible way to prove that has a unique solution for is to show that is a strictly decreasing function. Hence, we try to find conditions to guarantee and .
We check the coefficients in Equations (15) and (17) to find that
and
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 , which is the optimal solution for , which is the unique solution for .
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:
and
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],
,
, such that
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,
, 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
= 0. Hence, we must refer to another paper, Lin [
30], and his numerical example,
and
, to imply that
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
We show that is a decreasing function from to .
Hence, there is a unique point, say
, that is satisfying
. We recall that
and
have the same sign, such that
and
We prove that
is the optimal solution for our maximum-profit inventory model. Equation (12) reveals that the demand rate,
, is not included in the expression of
, 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
, the holding cost
, the shortage cost
, the purchasing cost
, the actual sale rate
, the selling price
, the planning horizon
, the ratio of extra selling related to the inventory level
, and the deterioration rate
. In Chuang [
11], his inventory model did not include the deterioration cost. Hence, we refer to Lin [
30] to adopt
. 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
, the shortage cost
, 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
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
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
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
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
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
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
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
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
has a negative relation with the deterioration cost.
5. Application of Our Approach
When
is degenerated to a constant function, denoted as
, we can rewrite Equation (12) as
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,
,
,
,
,
,
,
,
, and
.
Wu et al. [
10] did not consider the deterioration cost. In comparing our work with that of Wu et al. [
10], we suppose
. in the following discussion.
Wu et al. [
10] studied a ramp-type demand as
Depending on the relationship between
and
, Wu et al. [
10] partitioned their model into two cases: Case A:
, and Case B:
. Consequently, they constructed their profit function
and
for Cases A, and B, respectively.
For Case A, Wu et al. [
10] obtained
to decide their optimal solution for
.
For Case B, Wu et al. [
10] derived the same result as Equation (29) to decide their optimal solution for
.
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,
, (ii) Example 2,
, and (iii) Example 3,
.
For Example 1, Wu et al. [
10] claimed that the local maximum solution for
is
and the local maximum solution for
is
, with
For Example 2, Wu et al. [
10] found that the local maximum solution for
is
and the local maximum solution for
is
, with
For Example 3, Wu et al. [
10] found that the local maximum solution for
is
and the local maximum solution for
is
, with
Wu et al. [
10] were not aware that their three examples had the same global maximum point,
. 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
and
for Cases A and B, respectively, concerning the relation between
and
, which is unnecessary. The second one is to develop three examples, based on the relation
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
In the following, we list four related values to clearly indicate the values of decreasing from positive to negative. Hence, our solution is the right result.
On the other hand, we derive that
,
,
, and
to indicate that the maximum solution should be revised as
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
, and then rewrite Equation (12) as follows:
Under the condition of
, solving
will be a challenging problem because, in the integration of Equation (35), the following quadratic form denoted as E, with
appeared.
For a linear term in the exponential function, researchers can derive that
to indicate that when the deterioration function is a constant term as
, 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 , 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
with
and
For the improper integration of
, researchers know that
However, for another definite integration, for example, that denoted as F, with
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)
Depending on the accuracy, practitioners selected a number, denoted as “m”, to accept that
where the estimation error is less than
, 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
Based on the above discussion, in the beginning, we will adopt the second approach to solve , when .
We construct an example for a linear deterioration function with and to extend the deterioration function to a linear form such that , , , , , , , , , , and .
For a given value of
, we uniformly partition the integration interval
into twenty subintervals, and then we estimate
using the following approximation:
We begin to search for the optimal solution, , which is the solution of , with the knowledge that , , and is a decreasing function.
Many methods can locate the value of the optimal solution, . In the following, we will use a direct numerical method to list values of that will clearly reveal that is a decreasing function from a positive value to a negative value, and then we can find the solution of .
For
, we list the estimation result of
in the following
Table 11.
Based on
Table 11, we know that the optimal solution,
, satisfies
Hence, for
, we list the estimation result of
in
Table 12.
Based on
Table 12, we know that the optimal solution,
, satisfies
Hence, for
, we list the estimation result of
in
Table 13.
Based on
Table 13, we know that the optimal solution,
, satisfies
Hence, for
, we list the estimation result of
in
Table 14.
Based on
Table 14, we know that the optimal solution,
, satisfies