Next Article in Journal
Modeling and Experimental Verification of Plasma Jet Electromagnetic Signals
Previous Article in Journal
Modelling Key Performance Indicators in a Gamified Waste Management Tool
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Derivation of the Optimal Solution for the Economic Production Quantity Model with Planned Shortages without Derivatives

Department of Strategic Management & Operations, Woodbury School of Business, Utah Valley University, Orem, UT 84058, USA
Modelling 2022, 3(1), 54-69; https://doi.org/10.3390/modelling3010004
Submission received: 22 October 2021 / Revised: 14 December 2021 / Accepted: 18 December 2021 / Published: 5 January 2022

Abstract

:
In this paper, we study a reformulation of the Economic Production Quantity (EPQ) problem. We study a more general version of the problem first and derive the conditions for an optimal solution, as well as the optimal solution itself, all without using derivatives. Then, we apply the approach to the reformulated EPQ problem. This version of the EPQ problem has been tackled by a number of researchers, wherein they have derived the conditions for the optimal solution and proposed algebraic derivations. However, their derivations for the conditions, as well as the optimal solution, have been shown to be questionable. Other than being questionable, the existing approaches are so complicated that they defeat the purpose of simplifying the optimization by using a derivative-free approach. We propose a correct and more succinct, much less complicated approach to derive the conditions and the optimal solution without using derivatives.

1. Introduction

The basic EOQ model was developed by Harris [1]. In this pioneering work, he states that “the solution to this problem requires higher mathematics.” Grübbstrom and Erdem [2] show that the EOQ problem (and its backordering extension) can be solved algebraically. Cárdenas-Barrón [3] extends the same approach to the Economic Production Quantity (EPQ) model with backorders, and Wee et al. [4] extend it to the EOQ problem with a one-time discount offer. Huang [5] extends it to planned backorders and defective items. Ronald et al. [6] use analytic geometry, and Chang et al. [7] simplify the algebraic approach for both the EOQ and EPQ models in the former articles.
Sphicas [8] diverges from the former approaches and develops what he calls the “completing the perfect square” approach in solving the EOQ problem with linear and fixed backordering costs. Omar [9] generalizes this new approach. Cárdenas-Barrón [10] combines it with analytic geometry. Wee and Chung [11] extend it to a two-echelon inventory ordering problem. Chung [12] improves the approach of Wee and Chung [11].
Teng [13] introduces the arithmetic mean–geometric mean inequality approach. Lung [14] and Cárdenas-Barrón [15] show that this approach is not general. Cárdenas-Barrón [15] lists three conditions for this approach to work correctly, and Cárdenas-Barrón [16] extends it to the EOQ and EPQ problems with backorders. Ouyang et al. [17] use it for the EOQ problem with defective items and partially permissible delay in payments. Cárdenas-Barrón [18] uses it for a two-echelon problem similar to the one in Wee and Chung [11]. Cárdenas-Barrón [19] uses analytic geometry in conjunction with this approach for the EOQ and EPQ problems. Lin [20] modifies the approach in Cárdenas-Barrón [16]. Teng et al. [21] combine the completing the perfect square and the arithmetic mean–geometric mean inequality approaches and apply them to the vendor–buyer inventory problem. Sphicas [22] algebraically determines the solution to the EOQ problem with fixed and linear backordering costs and calls it the generalized EOQ formula. Luo [23] points out some issues in Sphicas [22].
Minner [24] proposes the cost comparisons approach to find the optimal order interval in the EOQ and EPQ problems (and in their backordering extensions). Wee et al. [25] optimize the order quantity as opposed to the order interval using a similar approach. Cárdenas-Barrón [26] compares the earlier algebraic approaches and the cost comparison approach and derives the optimal backorder quantity in the EOQ and EPQ models. Chung [27] extends the cost comparisons approach to the vendor–buyer production–inventory model. Widyadana et al. [28] extend it to the deteriorating items model.
Chang et al. [7] pose an open question in the context of their derivation for the EPQ problem with planned backorders: solving a reformulation of the problem without using derivatives. The reformulated problem has the backorder quantity as the only decision variable because the order quantity has been eliminated based on the relationship between the optimal order and backorder quantities. The literature studying this open question is summarized in Table 1. Lau et al. [29] determine the conditions for the existence and uniqueness of the optimal solution for the posed problem. Chiu et al. [30] analyze the approach of Lau et al. [29] and point out questionable results and make corrections and improvements. Luo and Chou [31] point out questionable results in both Lau et al. [29] and Chiu et al. [30], complete their analyses, and propose a purely algebraic approach, as opposed to the former, which uses differential calculus. Çalışkan [32] shows that the algebraic derivation in Luo and Chou [31] is also questionable because important steps of the process have been left out and there are some incorrect conclusions along the process, so the open question has not been answered yet. We propose a correct and much simpler, easy to understand, and succinct approach for both the derivation of the optimal solution and the derivation of the necessary conditions for the problem. Our derivation is based on the method proposed in Çalışkan [33].

2. Description of the Problem

We use the same notation as in Cárdenas-Barrón [3], Chang et al. [7], and Luo and Chou [31] to represent the parameters and variables:
D the demand per unit time K the cos t of setup per production batch h the cos t of inventory holding per unit per unit time b the cos t of backordering per unit per unit time c the cos t of production per unit of the item Q the batch size to use in each production cycle B the number of units to backorder in each production cycle
The total cost function for the EPQ problem with backorders that is used by Cárdenas-Barrón [3] and slightly revised later by Ronald et al. [6] and Chang et al. [7] is as follows:
C ( Q , B ) = b + h 2 ρ Q B 2 h B + h ρ 2 Q + K D Q + c D
where ρ = 1 D P . Chang et al. [7] reformulate Equation (1) by considering the following relationship between Q and B for the optimal solution:
Q ( B ) = 2 h ρ K D + b + h 2 ρ B 2
and by substituting Equation (2) into Equation (1), we obtain the following total cost function in only B:
C ( B ) = c D + h ( 1 + α ) B 2 + β B
where α = b h , and β = 2 ρ K D h . The question is to minimize Equation (3) without using derivatives, under the conditions that α > 0 , β > 0 and 0 < B < , which is equivalent to minimizing
f ( B ) = ( 1 + α ) B 2 + β B .
Lau et al. [29] extend the problem of minimizing Equation (4) to a more general problem of minimizing
f ( x ) = a x 2 + b x + c x
Both Equations (4) and (5) contain square roots but we have to consider the positive square roots. Since the underlying problem is the total cost for an inventory model, negative total cost is not possible. Furthermore, f ( B ) and f ( x ) technically would not be considered functions if both negative and positive square roots are considered, resulting in a one-to-two mapping between x and f ( x ) .

3. The Optimal Solution and Its Conditions

We will first derive the conditions for Equation (5) to have a unique interior point minimum x * , i.e., where f ( x ) has zero slope, as opposed to a boundary solution where f ( x ) has nonzero slope. First, we observe that Equation (5) can take on different shapes, it could be convex or concave, and its range may be entirely in R or it may be partially or entirely in C . f ( x ) does not have to be convex to have a unique interior minimum. However, if f ( x ) is convex and has a finite minimum that corresponds to a finite x, then it will be unique. We will derive the conditions for a unique interior minimum by studying f ( x ) and its shape with regard to when it is convex vs. when it is concave and when it has no real values.
Lemma 1.
Equation (5) has a unique interior minimum x * < only if g ( x ) = a x 2 + b x + c has at most one real root.
Proof. 
To be convex, f ( x ) has to satisfy Jensen’s inequality (Jensen [34]):
λ f ( x 1 ) + ( 1 λ ) f ( x 2 ) f ( λ x 1 + ( 1 λ ) x 2 )
for x 1 , x 2 R , x 1 x 2 , and 0 < λ < 1 . Because f ( x ) is complex-valued when g ( x ) < 0 , we assume that both g ( x 1 ) 0 and g ( x 2 ) 0 . Therefore, the following should hold:
λ a x 1 2 + b x 1 + c λ x 1 + ( 1 λ ) a x 2 2 + b x 2 + c ( 1 λ ) x 2 a ( λ x 1 + ( 1 λ ) x 2 ) 2 + b ( λ x 1 + ( 1 λ ) x 2 ) + c λ x 1 ( 1 λ ) x 2
After cancelling the identical terms on both sides, we obtain:
λ a x 1 2 + b x 1 + c + ( 1 λ ) a x 2 2 + b x 2 + c a ( λ x 1 + ( 1 λ ) x 2 ) 2 + b ( λ x 1 + ( 1 λ ) x 2 ) + c
We can square both sides and the equality will still hold because both sides are non-negative. We then obtain:
λ 2 ( a x 1 2 + b x 1 + c ) + ( 1 λ ) 2 ( a x 2 2 + b x 2 + c ) + 2 λ ( 1 λ ) ( a x 1 2 + b x 1 + c ) ( a x 2 2 + b x 2 + c ) a ( λ 2 x 1 2 + 2 λ ( 1 λ ) x 1 x 2 + ( 1 λ ) 2 x 2 2 ) + λ b x 1 + ( 1 λ ) b x 2 + c
This can further be simplified as:
λ ( λ 1 ) b x 1 + λ ( λ 1 ) b x 2 + 2 λ ( λ 1 ) c 2 λ ( λ 1 ) ( a x 1 2 + b x 1 + c ) ( a x 2 2 + b x 2 + c ) 2 λ ( λ 1 ) a x 1 x 2
Dividing both sides by λ ( λ 1 ) < 0 will change the direction of the inequality. Rearranging the terms after cancelling λ ( λ 1 ) , we obtain:
2 a x 1 x 2 + b ( x 1 + x 2 ) + 2 c 2 ( a x 1 2 + b x 1 + c ) ( a x 2 2 + b x 2 + c )
If the left-hand side of Equation (7), call it h ( x 1 , x 2 ) , is non-negative, we can again square both sides and the inequality will still hold. We assume for now that h ( x 1 , x 2 ) 0 and square both sides of Equation (7) to obtain:
4 a 2 x 1 2 x 2 2 + 4 a x 1 x 2 [ b ( x 1 + x 2 ) + 2 c ] + [ b ( x 1 + x 2 ) + 2 c ] 2 4 ( a x 1 2 + b x 1 + c ) ( a x 2 2 + b x 2 + c )
This can further be simplified as follows:
2 ( 4 a c b 2 ) x 1 x 2 + ( b 2 4 a c ) x 1 2 + ( b 2 4 a c ) x 2 2 0 ( b 2 4 a c ) ( x 1 x 2 ) 2 0 b 2 4 a c 0
We will now show that h ( x 1 , x 2 ) 0 for every x 1 and x 2 for which f ( x 1 ) , f ( x 2 ) R and f ( x 1 ) , f ( x 2 ) 0 . We know that b x 1 + c a x 1 2 and b x 2 + c a x 2 2 . Then, the following inequality should hold in general:
h ( x 1 , x 2 ) = 2 a x 1 x 2 + ( b x 1 + c ) + ( b x 2 + c ) 2 a x 1 x 2 a x 1 2 a x 2 2 = a ( x 1 x 2 ) 2
1.
When a 0 , h ( x 1 , x 2 ) 0 and therefore Equation (8) is applicable. Then, the following is true for f ( x ) :
(a)
If a < 0 and b 2 4 a c < 0 , f ( x ) C for all x.
(b)
If a < 0 and b 2 4 a c = 0 , f ( x ) consists of a single point at x = b 2 a ; f ( x ) C elsewhere.
(c)
If a < 0 and b 2 4 a c > 0 , then f ( x ) is concave between the roots of g ( x ) ; f ( x ) C elsewhere.
(d)
If a = 0 , f ( x ) is concave where g ( x ) 0 and f ( x ) C elsewhere.
2.
When a > 0 and b 2 4 a c > 0 , g ( x ) has two real roots and f ( x ) C between the roots. Therefore, we simply need to show that h ( x 1 , x 2 ) 0 for x 2 x 1 and x 2 x 1 . h ( x 1 , x 2 ) 0 for x 2 = x 1 . Without loss of generality, we fix x 1 at the smaller root. We will now show that h ( x 1 , x 2 ) 0 for every x 2 < x 1 . Substituting the smaller root for x 1 in h ( x 1 , x 2 ) results in:
h ( x 1 , x 2 ) = 2 a x 2 b b 2 4 a c 2 a + b b b 2 4 a c 2 a + b x 2 + 2 c = b 2 4 a c x 2 b 2 a b 2 2 a + 2 c
Because x 2 < x 1 = b 2 a b 2 4 a c 2 a , the following holds:
h ( x 1 , x 2 ) > b 2 4 a c b 2 a + b 2 4 a c 2 a b 2 a b 2 2 a + 2 c = 0
Therefore, h ( x 1 , x 2 ) 0 for x 2 < x 1 and Equation (8) is applicable. Similarly, we can fix x 1 at the greater root and, using a similar process, we can prove that h ( x 1 , x 2 ) 0 for x 2 > x 1 and Equation (8) is again applicable. Consequently, when a > 0 and b 2 4 a c > 0 , f ( x ) C between the roots and concave outside of the roots.
3.
When a > 0 and b 2 4 a c 0 , g ( x ) has, at most, one real root and f ( x ) R for all x. When b 2 4 a c = 0 , x = b 2 a is the only root. Again, without loss of generality, we fix x 1 at this root. Substituting the root for x 1 in h ( x 1 , x 2 ) results in:
h ( x 1 , x 2 ) = 2 a x 2 b 2 a b 2 2 a + b x 2 + 2 c = b 2 2 a + 2 c
Because b 2 4 a c , the following holds:
h ( x 1 , x 2 ) > 4 a c 2 a + 2 c = 0
Therefore, h ( x 1 , x 2 ) 0 for x 1 x 2 and the simplification is valid, and thus, f ( x ) is convex per Equation (8).
We observe that g ( x ) having at most one real root is necessary for a unique interior solution that is finite valued, but not sufficient. The following lemma further establishes the sufficiency conditions.
Lemma 2.
Equation (5) has a unique interior minimum > x * > 0 if and only if:
1.
b 2 4 a c 0
2.
a > 1
Proof. 
(1) is already established by Lemma 1. In order to prove (2), we will derive the closed-form euation for x * . Let
x u = x * + Δ x
x l = x * Δ x
for some Δ x > 0 . Then, the following will hold for x * and f ( x ) :
f ( x u ) f ( x * ) 0
f ( x l ) f ( x * ) 0
Substituting Equation (5) into Equation (16), we obtain:
a x u 2 + b x u + c x u a ( x * ) 2 + b x * + c + x * 0 a x u 2 + b x u + c a ( x * ) 2 + b x * + c ( x u x * ) 0 a x u 2 + b x u + c a ( x * ) 2 + b x * + c ( x u x * ) 1 0
Expanding the fraction by a x u 2 + b x u + c + a ( x * ) 2 + b x * + c and simplifying, we obtain:
a ( x u 2 ( x * ) 2 ) + b ( x u x * ) ( x u x * ) ( a x u 2 + b x u + c + a ( x * ) 2 + b x * + c ) 1 0 a ( x u + x * ) + b a x u 2 + b x u + c + a ( x * ) 2 + b x * + c 1
A similar analysis of Equation (17) will result in the following inequality:
a ( x l + x * ) + b a x l 2 + b x l + c + a ( x * ) 2 + b x * + c 1
Combining Equations (18) and (19), we obtain:
a ( x u + x * ) + b a x u 2 + b x u + c + a ( x * ) 2 + b x * + c 1 a ( x l + x * ) + b a x l 2 + b x l + c + a ( x * ) 2 + b x * + c
Let Δ x become gradually smaller and approach zero. Then, both sides of Equation (20) will approach one another and therefore they will both approach 1 as well. It then follows that:
2 a x * + b 2 a ( x * ) 2 + b x * + c = 1
Note that Equation (21) is equivalent to
f ( x ) = 2 a x * + b 2 a ( x * ) 2 + b x * + c 1 = 0
However, we have not used knowledge of the derivative of any functional form to obtain this result. Equation (21) yields:
2 a x * + b = 2 a ( x * ) 2 + b x * + c
We will now take the square of both sides to obtain a quadratic equation. Both negative and positive square roots of a ( x * ) 2 + b x * + c satisfy the resulting Equation (24). Hence, we obtain two solutions, one that corresponds to a slope of 0 for f ( x ) , and another that corresponds to a slope of 2 for f ( x ) . The one that corresponds to zero slope is the interior minimum point and it is required to be positive. We know that f ( x ) is required to be convex for a unique minimum, which means that its slope is negative to the left of the minimum point. Therefore, the larger of the two solutions of Equation (24) is the one that corresponds to the interior minimum point x * . Taking the square of both sides and rearranging the terms, we obtain:
4 ( a a 2 ) ( x * ) 2 + 4 ( b a b ) x * + ( 4 c b 2 ) = 0
The discriminant of Equation (24) could be calculated as follows:
Δ = [ 4 ( b a b ) ] 2 16 ( a a 2 ) ( 4 c b 2 ) = 16 ( 1 a ) ( b 2 4 a c )
Then, the solution to Equation (24) will be the larger of the two roots:
x * = 4 b ( 1 a ) + 4 ( 1 a ) ( b 2 4 a c ) 2 [ 4 a ( 1 a ) ] = 1 2 a b + b 2 4 a c 1 a
From (1), we know that 4 a c b 2 . Equation (26) requires a > 1 in order to have a positive value inside the square root so that x * R , and it requires a 1 to ensure that x * < . Therefore, we need a > 1 . It then follows that c 0 because 4 a c b 2 (this is not a separate condition, but a consequence of a > 1 and 4 a c b 2 ). □
We also need to establish conditions for x * > 0 and f ( x ) > 0 because the underlying problem is the EPQ problem with backorders, where the problem is defined on x > 0 and f ( x ) > 0 for all x > 0 .
Lemma 3.
f ( x ) > 0 for all x R if and only if 4 ( a 1 ) c > b 2 and a > 1 .
Proof. 
The condition f ( x ) > 0 can be simplified as follows:
a x 2 + b x + c x > 0 a x 2 + b x + c > x 2 ( a 1 ) x 2 + b x + c > 0
Equation (27) will be satisfied for all x R if and only if its graph entirely lies above the x axis—in other words, if it has no real roots and a 1 > 0 or a > 1 . Therefore, Δ = b 2 4 ( a 1 ) c < 0 , or 4 ( a 1 ) c > b 2 . □
Lemma 4.
If b 0 , then x * > 0 if and only if 4 c > b 2 .
Proof. 
If b 0 , x * > 0 requires:
b 2 4 a c 1 a > b b 2 4 a c 1 a > b 2 4 c > b 2
and if b 0 , there is no additional condition for guaranteeing that x * > 0 . □
We can now combine and summarize all of the necessary and sufficient conditions in the following theorem, as in Luo and Chou [31].
Theorem 1.
f ( x ) is positive, convex, and has a unique positive finite valued non-boundary minimum if and only if a > 1 and one of the following two conditions is satisfied:
1.
b < 0 and 4 ( a 1 ) c > b 2
2.
b 0 and 4 c > b 2
Proof. 
Immediate from Lemmas 1–4: if a > 1 , then 4 a c > 4 ( a 1 ) c > b 2 as well as 4 a c > 4 c > b 2 , which means, for case 1, we can drop b 2 4 a c 0 and keep 4 ( a 1 ) c > b 2 ; and for case 2, we can drop b 2 4 a c 0 and keep 4 c > b 2

4. The Original Inventory Problem

In the original EPQ with the backordering problem of Equation (4), a = 1 + α ; b = 0 , and c = β . Clearly, a = 1 + b h > 1 ; therefore, the only remaining requirement is 4 β > 0 and this is also satisfied because β = 2 ρ K D h > 0 . Therefore, we can substitute these values of a, b, and c in Equation (26) and obtain the optimal solution:
B * = β α ( 1 + α ) , and f ( B * ) = α B
Contrary to the approach of Luo and Chou [31], our approach can directly be applied to Equation (4) as well. Let B * be the minimum point of f ( B ) that is not a boundary solution and let
B u = B * + Δ B
B l = B * Δ B
for some Δ B > 0 . Then, the following will hold due to B * being a non-boundary minimum point:
f ( B u ) f ( B * ) 0
f ( B l ) f ( B * ) 0
Equation (32) can further be simplified as follows:
( 1 + α ) B u 2 + β B u ( 1 + α ) ( B * ) 2 + β + B * 0 ( 1 + α ) B u 2 + β ( 1 + α ) ( B * ) 2 + β B u B * 1 0
Expanding the fraction by ( 1 + α ) B u 2 + β + ( 1 + α ) ( B * ) 2 + β results in
( 1 + α ) ( B u 2 ( B * ) 2 ) ( B u B * ) ( ( 1 + α ) B u 2 + β + ( 1 + α ) ( B * ) 2 + β ) 1 0 ( 1 + α ) ( B u + B * ) ( 1 + α ) B u 2 + β + ( 1 + α ) ( B * ) 2 + β 1
Applying the same approach to Equation (33) will result in:
( 1 + α ) ( B l + B * ) ( 1 + α ) B l 2 + β + ( 1 + α ) ( B * ) 2 + β 1
Again, combining Equations (34) and (35), we obtain:
( 1 + α ) ( B u + B * ) ( 1 + α ) B u 2 + β + ( 1 + α ) ( B * ) 2 + β 1 ( 1 + α ) ( B l + B * ) ( 1 + α ) B l 2 + β + ( 1 + α ) ( B * ) 2 + β
Let Δ B become gradually smaller and approach zero. Then, both sides of Equation (36) will approach one another and therefore they will both approach 1 as well. It then follows that:
2 ( 1 + α ) B * 2 ( 1 + α ) ( B * ) 2 + β = 1
Therefore, the solution to Equation (37) is:
B * = β α ( 1 + α )
f ( B ) is strictly convex if and only if the following holds for B 1 , B 2 R , B 1 B 2 and 0 < λ < 1 :
λ f ( B 1 ) + ( 1 λ ) f ( B 2 ) > f ( λ B 1 + ( 1 λ ) B 2 )
Equation (39) can be simplified further as follows:
λ ( 1 + α ) B 1 2 + β λ B 1 + ( 1 λ ) ( 1 + α ) B 2 2 + β ( 1 λ ) B 2 > ( 1 + α ) ( λ B 1 + ( 1 λ ) B 2 ) 2 + β λ B 1 ( 1 λ ) B 2 2 λ ( 1 λ ) β + 2 λ ( 1 λ ) [ ( 1 + α ) B 1 2 + β ] [ ( 1 + α ) B 2 2 + β ] > 2 λ ( 1 λ ) ( 1 + α ) B 1 B 2
Dividing both sides by the positive term 2 λ ( 1 λ ) , we obtain:
[ ( 1 + α ) B 1 2 + β ] [ ( 1 + α ) B 2 2 + β ] > ( 1 + α ) B 1 B 2 + β ( 1 + α ) β ( B 1 2 + B 2 2 2 B 1 B 2 ) > 0 ( B 1 B 2 ) 2 > 0
Therefore, B * > 0 is a unique non-boundary finite valued optimal solution that also satisfies f ( B * ) > 0 .

5. Limitations and Further Research

The proposed method does not explicitly use the derivative of any function and, instead, derives the optimal solution using only algebraic manipulations. It is simple, succinct, and intuitive, which makes it a useful pedagogical tool to teach inventory models to students with limited skills in differential calculus. On the other hand, it implicitly uses the first derivative of the total cost function, by obtaining it through algebraic operations, as shown in Çalışkan [33]. It implicitly assumes that the total cost function is continuous and continuously differentiable. When this is not the case, it will require further steps and will become more complicated. However, the same is true for any other approach when the total cost function is not continuous. Further research may apply the proposed method to total cost functions that are not continuous, and to other problems in inventory management and operations research. It may also be adapted to situations in which there are constraints and/or more decision variables in the problem.

6. Conclusions

In this paper, we answer the open question posed in Chang et al. [7]: deriving the optimal solution for a reformulation of the EPQ problem with backorders using an approach that does not require derivatives. The earlier approach proposed by Luo and Chou [31] is not entirely correct and complete, as shown in Çalışkan [32], so the open question is still outstanding. We also point out that the approach in Luo and Chou [31] is far too complicated that it defeats the purpose, which is simplifying the optimization by not using differential calculus so that undergraduate or even high school students, or practitioners without a working knowledge of differential calculus, can understand it. Our approach, on the other hand, is simple, succinct, and easy to follow. We consider the generalized version of the open question proposed by Lau et al. [29] and derive the necessary and sufficient conditions for the problem to have an acceptable optimal solution for the underlying EPQ problem, as well as the optimal solution itself. We also show that our approach could be applied directly to the original inventory problem with ease. Finally, we study the general problem to provide more insights and show that it can assume different characteristics depending on its parameters (Appendix A).

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. General Properties of f(x)

In Section 3, we established the conditions for f ( x ) to be positive and to have a unique interior point minimum that is positive so that f ( x ) makes sense as the total cost function of an inventory model. In its general form, f ( x ) takes on different shapes depending on the parameters a, b, and c. This section is not a part of the answer to the open question challenge of solving the EPQ problem with backordering in Chang et al. [7] without using derivatives, as we have already addressed this in Section 3 and Section 4. This section provides further insights into the function f ( x ) . Therefore, we freely use differential calculus in this section. We first analyze the limiting behavior of f ( x ) . In order to calculate the limits of f ( x ) as x or x , we can transform f ( x ) into the following form:
f ( x ) = a x 2 + b x + c x = x 2 a + b x + c x 2 x = | x | a + b x + c x 2 x f ( x ) = a + b x + c x 2 1 1 x , x < 0 a + b x + c x 2 1 1 x , x 0
Lemma A1.
The following limits are valid for f ( x ) :
1.
if a = 1 then lim x f ( x ) = and lim x f ( x ) = b 2
2.
if 1 > a > 0 then lim x f ( x ) = and lim x f ( x ) =
3.
if a > 1 then lim x f ( x ) = lim x f ( x ) =
4.
if a = 0 and b > 0 then lim x f ( x ) C and lim x f ( x ) =
5.
if a = 0 and b < 0 then lim x f ( x ) = and lim x f ( x ) C
6.
if a < 0 then lim x f ( x ) C and lim x f ( x ) C
Proof 
1.
For x , we simply evaluate Equation (A1) for x < 0 , which results in as both the numerator and denominator are negative and they approach −2 and 0, respectively. For x , we have to apply L’óspital’s rule because evaluating the integral results in 0 0 indeterminacy:
lim x a + b x c x 2 1 1 x = lim x b x 2 + 2 c x 3 2 a + b x + c x 2 1 x 2 = lim x b 2 c x 2 a + b x + c x 2 = b 2 a = b 2
2.
For x , evaluating Equation (A1) for x < 0 results in , as both the numerator and denominator are negative. As for x , evaluating Equation (A1) for x > 0 results in because the numerator is negative but the denominator is positive.
3.
For x , evaluating Equation (A1) for x < 0 results in , as both the numerator and denominator are negative. As for x , evaluating Equation (A1) for x > 0 results in because both the numerator and denominator are positive.
4.
If a = 0 and b > 0 , then lim x f ( x ) C because g ( x ) < 0 for x < c b . Additionally, lim x f ( x ) = 1 0 = because the denominator is positive.
5.
If a = 0 and b < 0 , then lim x f ( x ) = 1 0 = because the denominator is negative. Additionally, lim x f ( x ) C because g ( x ) < 0 for x > c b .
6.
If a < 0 , then g ( x ) < 0 as x as well as x and therefore lim x f ( x ) C and lim x f ( x ) C .
We can now summarize the general results regarding the shape and the minimum point x * of f ( x ) . The following theorem lists all of the possibilities for the shape, range, and the minimum point of f ( x ) .
Theorem A1.
Let the range of f ( x ) be Y. If g ( x ) has two real roots, let x o l and x o u be the roots. Then, the following statements hold for f(x):
1.
If b 2 4 a c < 0 and a > 0 , f(x) is strictly convex and Y R . Furthermore:
(a)
If a < 1 , x * =
(b)
If a = 1 , x * =
(c)
If a > 1 , x * = 1 2 a b ± b 2 4 a c 1 a
2.
If b 2 4 a c < 0 and a < 0 , f(x) is neither convex nor concave and Y C
3.
If b 2 4 a c = 0 and a > 0 , f(x) is convex but not strictly convex and Y R (piecewise linear). Furthermore:
(a)
If a < 1 , x * =
(b)
If a = 1 , x * = ( b 2 , )
(c)
If a > 1 , x * = b 2 a
4.
If b 2 4 a c = 0 and a < 0 , f(x) is both convex and concave because f ( x ) is a single point; x * = b 2 a , and Y = { b 2 a } for x = b 2 a , and Y C for x b 2 a
5.
If b 2 4 a c > 0 , f ( x ) is strictly concave and x * = a r g m i n x { f ( x o l ) , f ( x o u ) } and furthermore:
(a)
if a > 0 , Y C for x o u x x o l and Y R for x < x l and x > x u
(b)
if a < 0 , Y R for x o u x x o l , and Y C for x < x l and x > x u
6.
If a = 0 , f ( x ) is strictly concave and furthermore:
(a)
if b > 0 , x * = , and Y R for x x o and Y C for x < x o
(b)
if b < 0 , x * = c b , and Y R for x x o and Y C for x > x o
Proof 
1.
By Lemma 1, f ( x ) is strictly convex when b 2 4 a c < 0 and a > 0 . Furthermore, g ( x ) > 0 for all x R and therefore Y R .
(a)
If a < 1 , x * = because f ( x ) is strictly monotone decreasing and lim x f ( x ) = by Lemma A1 (2).
(b)
If a = 1 , x * = because f ( x ) is strictly monotone decreasing and lim x f ( x ) = b 2 by Lemma A1 (1).
(c)
If a > 1 , x * = 1 2 a b ± b 2 4 a c 1 a by Lemma 2.
2.
If b 2 4 a c < 0 and a < 0 , g ( x ) < 0 for all x R and has no real roots. Therefore, f ( x ) is neither convex nor concave and Y C
3.
If b 2 4 a c = 0 and a > 0 , f ( x ) is convex but not strictly convex by Lemma 1. Because Δ = b 2 4 a c = 0 for g ( x ) , c = b 2 4 a and g ( x ) can be expressed as a perfect square: g ( x ) = ( a x + b a 2 a ) 2 . Therefore, f ( x ) can be expressed as follows:
f ( x ) = a x + b a 2 a 2 x = a x + b a 2 a x = ( a + 1 ) x b a 2 a , x < b 2 a ( a 1 ) x + b a 2 a , x b 2 a
which is obviously piecewise linear with a corner point of x = b 2 a and Y R . Furthermore:
(a)
If a < 1 , x * = because f ( x ) is strictly monotone decreasing and lim x f ( x ) = .
(b)
If a = 1 , x * = b 2 , , which means infinitely many optimal solutions along the horizontal line of f ( x ) = b 2 .
(c)
If a > 1 , the unique minimum occurs at the corner point: x * = b 2 a .
4.
If b 2 4 a c = 0 , f ( x ) is convex but not strictly convex by Lemma 1; and if a < 0 , g ( x ) has a single real root x = b 2 a and g ( x ) < 0 for all x R , x b 2 a . Therefore, x * = b 2 a and Y = { b 2 a } for x = b 2 a , and Y C for x b 2 a
5.
If b 2 4 a c > 0 , f ( x ) is strictly concave by Lemma 1 and g ( x ) has two real roots. Therefore, x * = a r g m i n x { f ( x o l ) , f ( x o u ) } .
(a)
if a > 0 , g ( x ) < 0 between its roots. Therefore, Y C for x o u x x o l , and Y R for x < x l and x > x u .
(b)
if a < 0 , g ( x ) < 0 outside the interval between the two roots. Therefore, Y R for x o u x x o l , and Y C for x < x l and x > x u .
6.
If a = 0 , f ( x ) is strictly concave by Lemma 1 and g ( x ) has a single real root x = c b .
(a)
if b > 0 , g ( x ) < 0 for all x R x < c b . Therefore, Y R for x c b and Y C for x < c b . In addition, x * = because lim x f ( x ) = .
(b)
if b < 0 , g ( x ) < 0 , for all x R x > c b . Therefore, Y R for x c b and Y C for x > c b . In addition, x * = c b because f ( x ) is strictly monotone decreasing.
Figure A1Figure A3 demonstrate the various shapes of f ( x ) that correspond to the 12 cases that are listed in Theorem A1. As we can see from the graphs, only Case 1c satisfies the requirements of > f ( x * ) > 0 , f ( x ) > 0 for x R , and f ( x * ) = 0 , which essentially is the only case that satisfies all of the conditions in Theorem 1.
Figure A1. The graph of f ( x ) for Case 1(a) through Case 2 from Theorem A1.
Figure A1. The graph of f ( x ) for Case 1(a) through Case 2 from Theorem A1.
Modelling 03 00004 g0a1
Figure A2. The graph of f ( x ) for Case 3(a) through Case 4 from Theorem A1.
Figure A2. The graph of f ( x ) for Case 3(a) through Case 4 from Theorem A1.
Modelling 03 00004 g0a2
Figure A3. The graph of f ( x ) for Case 5(a) through Case 6(b) from Theorem A1.
Figure A3. The graph of f ( x ) for Case 5(a) through Case 6(b) from Theorem A1.
Modelling 03 00004 g0a3

References

  1. Harris, F. How Many Parts to Make At Once. Fact. Mag. Manag. 1913, 10, 135ߝ136, 152. [Google Scholar] [CrossRef]
  2. Grubbström, R.W.; Erdem, A. The EOQ with backlogging derived without derivatives. Int. J. Prod. Econ. 1999, 59, 529–530. [Google Scholar] [CrossRef]
  3. Cárdenas-Barrón, L.E. The economic production quantity (EPQ) with shortage derived algebraically. Int. J. Prod. Econ. 2001, 70, 289–292. [Google Scholar] [CrossRef]
  4. Wee, H.; Chung, S.; Yang, P. Technical Note A Modified EOQ Model with Temporary Sale Price Derived without Derivatives. Eng. Econ. 2003, 48, 190–195. [Google Scholar] [CrossRef]
  5. Huang, Y.F. The deterministic inventory models with shortage and defective items derived without derivatives. J. Stat. Manag. Syst. 2003, 6, 171–180. [Google Scholar] [CrossRef]
  6. Ronald, R.; Yang, G.K.; Chu, P. Technical note: The EOQ and EPQ models with shortages derived without derivatives. Int. J. Prod. Econ. 2004, 92, 197–200. [Google Scholar] [CrossRef]
  7. Chang, S.J.; Chuang, J.P.; Chen, H.J. Short comments on technical note—The EOQ and EPQ models with shortages derived without derivatives. Int. J. Prod. Econ. 2005, 97, 241–243. [Google Scholar] [CrossRef]
  8. Sphicas, G.P. EOQ and EPQ with linear and fixed backorder costs: Two cases identified and models analyzed without calculus. Int. J. Prod. Econ. 2006, 100, 59–64. [Google Scholar] [CrossRef]
  9. Omar, M.; Zubir, M.B.; Moin, N. An alternative approach to analyze economic ordering quantity and economic production quantity inventory problems using the completing the square method. Comput. Ind. Eng. 2010, 59, 362–364. [Google Scholar] [CrossRef]
  10. Cárdenas-Barrón, L.E. Optimal manufacturing batch size with rework in a single-stage production system–-A simple derivation. Comput. Ind. Eng. 2008, 55, 758–765. [Google Scholar] [CrossRef]
  11. Wee, H.M.; Chung, C.J. A note on the economic lot size of the integrated vendor–-Buyer inventory system derived without derivatives. Eur. J. Oper. Res. 2007, 177, 1289–1293. [Google Scholar] [CrossRef]
  12. Chung, K.J. “A note on the economic lot size of the integrated vendor–-buyer inventory system derived without derivatives”: A comment. Eur. J. Oper. Res. 2009, 198, 979–982. [Google Scholar] [CrossRef]
  13. Teng, J.T. A simple method to compute economic order quantities. Eur. J. Oper. Res. 2009, 198, 351–353. [Google Scholar] [CrossRef]
  14. Leung, K.N.F. Some comments on “A simple method to compute economic order quantities”. Eur. J. Oper. Res. 2010, 201, 960–961. [Google Scholar] [CrossRef]
  15. Cárdenas-Barrón, L.E. A simple method to compute economic order quantities: Some observations. Appl. Math. Model. 2010, 34, 1684–1688. [Google Scholar] [CrossRef] [Green Version]
  16. Cárdenas-Barrón, L.E. An easy method to derive EOQ and EPQ inventory models with backorders. Comput. Math. Appl. 2010, 59, 948–952. [Google Scholar] [CrossRef] [Green Version]
  17. Ouyang, L.Y.; Chang, C.T.; Shum, P. The EOQ with defective items and partially permissible delay in payments linked to order quantity derived algebraically. Cent. Eur. J. Oper. Res. 2012, 20, 141–160. [Google Scholar] [CrossRef]
  18. Cárdenas-Barrón, L.E.; Wee, H.M.; Blos, M.F. Solving the vendor–-Buyer integrated inventory system with arithmetic–-Geometric inequality. Math. Comput. Model. 2011, 53, 991–997. [Google Scholar] [CrossRef] [Green Version]
  19. Cárdenas-Barrón, L.E. The derivation of EOQ/EPQ inventory models with two backorders costs using analytic geometry and algebra. Appl. Math. Model. 2011, 35, 2394–2407. [Google Scholar] [CrossRef] [Green Version]
  20. Lin, S.S.C. Note on “The derivation of EOQ/EPQ inventory models with two backorders costs using analytic geometry and algebra”. Appl. Math. Model. 2019, 73, 378–386. [Google Scholar] [CrossRef]
  21. Teng, J.T.; Cárdenas-Barrón, L.E.; Lou, K.R. The economic lot size of the integrated vendor–-buyer inventory system derived without derivatives: A simple derivation. Appl. Math. Comput. 2011, 217, 5972–5977. [Google Scholar] [CrossRef]
  22. Sphicas, G.P. Generalized EOQ formula using a new parameter: Coefficient of backorder attractiveness. Int. J. Prod. Econ. 2014, 155, 143–147. [Google Scholar] [CrossRef]
  23. Luo, X.R. A Detailed Examination of Sphicas (2014), Generalized EOQ Formula Using a New Parameter: Coefficient of Backorder Attractiveness. Int. J. Prod. Econ. 2019, 11, 931. [Google Scholar] [CrossRef] [Green Version]
  24. Minner, S. A note on how to compute economic order quantities without derivatives by cost comparisons. Int. J. Prod. Econ. 2007, 105, 293–296. [Google Scholar] [CrossRef]
  25. Wee, H.M.; Wang, W.T.; Chung, C.J. A modified method to compute economic order quantities without derivatives by cost-difference comparisons. Eur. J. Oper. Res. 2009, 194, 336–338. [Google Scholar] [CrossRef]
  26. Cárdenas-Barrón, L. A note on how to compute economic order quantities without derivatives by cost comparisons: Some comments. Int. J. Appl. Manag. Sci. 2010, 2, 198–204. [Google Scholar] [CrossRef]
  27. Chung, C.J. An easy method to derive the integrated vendor—Buyer production—Inventory model with backordering using cost–difference rate comparison approach. Math. Comput. Model. 2013, 57, 632–640. [Google Scholar] [CrossRef]
  28. Widyadana, G.A.; Cárdenas-Barrón, L.E.; Wee, H.M. Economic order quantity model for deteriorating items with planned backorder level. Math. Comput. Model. 2011, 54, 1569–1575. [Google Scholar] [CrossRef]
  29. Lau, C.; Chou, E.; Dementia, J. Criterion to ensure uniqueness for minimum solution by algebraic method for inventory model. Int. J. Eng. Appl. Sci. 2016, 3, 71–73. [Google Scholar]
  30. Chiu, C.; Li, Y.; Julian, P. Improvement for Criterion for Minimum Solution of Inventory Model with Algebraic Approach. J. Bus. Manag. 2017, 19, 73–78. [Google Scholar] [CrossRef]
  31. Luo, X.R.; Chou, C.S. Technical note: Solving inventory models by algebraic method. Int. J. Prod. Econ. 2018, 200, 130–133. [Google Scholar] [CrossRef]
  32. Çalışkan, C. On “Technical note: Solving inventory models by algebraic method”. J. Stat. Manag. Syst. 2021, 24, 1533–1541. [Google Scholar] [CrossRef]
  33. Çalışkan, C. A general approach for the derivation of optimal solutions without derivatives. Int. J. Syst. Sci. Oper. Logist. 2021, 1–12. [Google Scholar] [CrossRef]
  34. Jensen, J.L.W.V. Sur les fonctions convexes et les inégalités entre les valeurs moyennes. Acta Math. 1906, 30, 175–193. [Google Scholar] [CrossRef]
Table 1. Summary of the previous literature most related to the present paper.
Table 1. Summary of the previous literature most related to the present paper.
ArticleContribution
Chang et al. [7]Open question: deriving the optimal solution for a reformulation of the EPQ problem with backorders
Lau et al. [29]The conditions for the existence and uniqueness of the optimal solution for the posed problem
Chiu et al. [30]Analysis of Lau et al. [29] and corrections and improvements of the former
Luo and Chou [31]Analysis of both Lau et al. [29] and Chiu et al. [30], a purely algebraic approach as opposed to the former
Çalışkan [32]Proof that Luo and Chou [31] is also partially incorrect and incomplete; correction and completion of the former
Çalışkan [33]A general approach to optimize objective functions without derivatives
This paperDefinitively answering the open question in Chang et al. [7]; correction of the former approaches; analysis of the general problem in Lau et al. [29]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Çalışkan, C. Derivation of the Optimal Solution for the Economic Production Quantity Model with Planned Shortages without Derivatives. Modelling 2022, 3, 54-69. https://doi.org/10.3390/modelling3010004

AMA Style

Çalışkan C. Derivation of the Optimal Solution for the Economic Production Quantity Model with Planned Shortages without Derivatives. Modelling. 2022; 3(1):54-69. https://doi.org/10.3390/modelling3010004

Chicago/Turabian Style

Çalışkan, Cenk. 2022. "Derivation of the Optimal Solution for the Economic Production Quantity Model with Planned Shortages without Derivatives" Modelling 3, no. 1: 54-69. https://doi.org/10.3390/modelling3010004

APA Style

Çalışkan, C. (2022). Derivation of the Optimal Solution for the Economic Production Quantity Model with Planned Shortages without Derivatives. Modelling, 3(1), 54-69. https://doi.org/10.3390/modelling3010004

Article Metrics

Back to TopTop