arxiv: v1 [physics.comp-ph] 17 Jan 2014

Transcription

1 An efficient method for soling a correlated multi-item inentory system Chang-Yong Lee and Dongu Lee The Department of Industrial & Systems Engineering, Kongu National Uniersity, Kongu South Korea (Dated: August 4, 208) arxi: [physics.comp-ph] 7 Jan 204 We propose an efficient method of finding an optimal solution for a multi-item continuous reiew inentory model in which a biariate Gaussian probability distribution represents a correlation between the demands of different items. By utilizing appropriate normalizations of the demands, we show that the normalized demands are uncorrelated. Furthermore, the set of equations coupled with different items can be decoupled in such a way that the order quantity and reorder point for each item can be ealuated independently from those of the other. As a result, in contrast to conentional methods, the solution procedure for the proposed method can be much simpler and more accurate without any approximation. To demonstrate the adantage of the proposed method, we present a solution scheme for a multi-item continuous reiew inentory model in which the demand of optional components depend on that of a anilla box, representing the customer s stochastic demand, under stochastic payment and budget constraints. We also perform a sensitiity analysis to inestigate the dependence of order quantities and reorder points on the correlation coefficient. I. INTRODUCTION In a competitie global market, a manufacturer usually needs to proide a wide ariety of products with a short amount of lead time to improe customer satisfaction and increase market share. When an order from a customer arries, differentiation to the tailored order is often postponed down the assembly line to achiee both product ariety and short lead time. In general, the more dierse a product, the longer lead time it requires. Modularization and postponement (or delayed differentiation) can be an effectie means to reduce the lead time while maintaining a wide ariety of products. Many modularization and postponement studies hae shown that these concepts offer an adantage in terms of reducing uncertainty and forecasting errors with regard to demand [8, 3, 4] in addition to creating product ariety and customization at low cost [5]. Thus, modularization and postponement hae become important concepts in the market to proide better serice to customers and make the business process more efficient. In a product line such as computer retailing or automobile assembly, concepts of modularization and postponement are realized by an assembly process that consists of a semi-finished product anilla box and optional components that are directly used in the final assembly. The anilla box consists of components, known as the commonality of parts, needed to assemble the final product with appropriate optional components, and the anilla box approach has been shown to be effectie under high ariance [3]. Applications of modularization and postponement to an inentory system require a multi-item model in which the demand of each of the seeral optional components depends on the demand or the presence of the anilla box. The anilla box represents the customer s stochastic demand, and optional components, in turn, depend on the demand of the anilla box for final assembly. Thus, the demands of the anilla box and optional components are stochastically correlated. As the inentory model becomes complicated because of the correlation, it is desirable to find an efficient and accurate method to sole the model system. A multi-item inentory model was proposed to comply with the concept of modularization and postponement [5]. The model consists of a anilla box and optional components in which the correlation between the two types of items is implemented as a biariate Gaussian probability distribution whereas the optional components are independent of each other. This model handles a continuous reiew inentory system in which an order quantity Q is placed wheneer an inentory leel reaches a certain reorder point r under the presence of serice leel and budget constraints. Subsequently, a stochastic payment is also included in the model in such a manner that the total inentory cost does not exceed a predetermined budget [6]. In this paper, we propose an efficient method to sole a correlated multi-item continuous reiew inentory model in which the correlation between the anilla box and an optional component is represented by a biariate Gaussian probability distribution. By using appropriate normalizations of the demands of items, we show that the set of Electronic address: clee@kongu.ac.kr

2 equations coupled with the anilla box and optional components can be reduced to sets of decoupled equations for each item. Furthermore, each set of decoupled equations is simplified in a closed form and soled without any approximation. Thus, the equations for each item can be soled independently of each other. The conentional method for soling such a model system is based on a heuristic of combining a Newton Raphson method and a Hadley Whitin iteratie procedure [6]. At each iteration, a candidate solution is found by using the Newton-Raphson method in which numerical integrations are carried out where required. The iteration proceeds until both Q and r sufficiently conerge. Briefly, the conentional method takes the set of simultaneous equations for Q s and r s as a whole and utilizes heuristic approximating procedures. Gien that the conentional method uses a rather complicated approximation and iteration, it requires heay computation time. In contrast, the proposed method does not rely on any approximation or heuristics. As a result, the solution procedure for the proposed method is much simpler, more accurate, and offers shorter computing time than the conentional method. We apply the proposed method to a correlated multi-item continuous reiew inentory model to demonstrate its usefulness. We also perform a sensitiity analysis in terms of the correlation to further characterize the behaior of the order quantity and reorder point of optional components. In addition, the proposed scheme can be used as a dependable method for a more generalized multi-item continuous reiew inentory model with much more complicated correlations among items. The rest of this paper is organized as follows. Section II discusses the normalization of the demands and introduces an illustratie model to show how the set of simultaneous equations can be decoupled and simplified. Section III describes the proposed method for soling the model system. Section IV presents experimental results and discusses the sensitiity analysis. Section V summarizes the study and gies our conclusions. 2 II. NORMALIZATION AND MULTI-ITEM INVENTORY MODEL A. Correlation and normalization Consider a multi-item inentory model that includes the correlation between a anilla box and an optional component. Furthermore, we allow multiple optional components and each optional component is dependent on the anilla box through a biariate Gaussian probability distribution, whereas optional components are independent of each other. A biariate Gaussian probability distribution function (PDF) of the random ariables X and X of the demand for the anilla box and the th optional component, respectiely, is gien by f(x,x ) = Σ exp { 2 ( x µ)t Σ ( x µ) }. () Here, the ariable ector x, the mean ector µ, the coariance matrix Σ, and the correlation coefficient ρ between X and X are expressed, respectiely, as [ ] [ ] [ ] x µ σ 2 x =, µ =, Σ = x µ σ 2, and ρ (2) σ with Σ being the determinant of the 2 2 matrix of Σ. Gien that X depends on X, we express the biariate PDF as a product of the marginal PDF of X and the conditional PDF of X gien X = x. In this way, the biariate PDF of Eq. () can be written as [ ( )] 2 f(x,x ) = exp { (x µ ) 2 } σ x µ +ρ σ 2σ 2 exp σ ρ 2 (x µ ) 2σ 2 ( ρ2 ) f X (x ) f X X (x x ). (3) In general, the demand for the anilla box is equal to the customer s demand, whereas the demand for each optional component depends on the safety stock of the anilla box. Gien that the safety stock of the anilla box is r µ, the demand for each optional component depends on the reorder point of the anilla box r. This implies that the conditional PDF of X is ealuated at X = r. Motiated by this characteristic, we define the normalized random ariables as Z X µ and Z X µ o σ o, (4)

3 3 where σ o σ ρ 2 and µ o µ +ρ σ (r µ ). (5) Note that the random ariable Z contains not only the demand X of the th optional component but also the reorder point r of the anilla box. With the normalization, the biariate PDF at X = r can be rewritten as f X (x ) f X X (x r ) = e z2 /2 e z2 /2. (6) We see from Eq. (6) that the normalization decomposes the biariate PDF into a product of two PDFs of Z and Z, both of which can be regarded as independent standard uniariate PDFs. This implies that the notion of the conditional PDF does not exist and there is no distinction between dependent and independent items. In what follows, we take a correlated multi-item continuous reiew inentory system as an illustration of the adantage of normalization. B. Model formulation and set of equations The illustratie model we consider in this paper is a correlated multi-item continuous reiew inentory system that includes two types of items: a anilla box and many optional components. Furthermore, as stated in Section II A, each optional component depends on the anilla box through a biariate Gaussian probability distribution. We use the following notations to formulate the model: A : fixed procuring cost, C : unit ariable procurement cost, D : expected annual demand, h : carrying cost, p : unit shortage cost, κ : serice cost rate, and β : aailable budget limit. Note that each of these terms can be used for both the anilla box and optional components wheneer possible. We distinguish the anilla box from the th optional component by the subscripts and o, respectiely. For instance, A and A o represent the fixed procuring costs of the anilla component and the th optional component, respectiely. The model is composed of the sum of the expected aerage annual cost (EAC) of the two types of items under budgetary constraint. The budgetary constraint, in turn, includes the serice costs of the anilla box and optional components. A detailed account of the model deelopment can be found in [5, 6]. The obectie of the model is to minimize the sum of EAC for the anilla box and m optional components under a stochastic budgetary constraint. That is, we would like to find (Q,r ) and ( Q o, r o ), where Q o = (Q o,q o2,,q om ) and r o = (r o,r o2,,r om ), that minimize EAC( Q, r) = EAC (Q,r )+EAC 2 ( Q o, r o ) (7) subect to Prob C (Q +r X )+ C o (Q o +r o X )+κ F X (r )+ κ o F X X (r o r ) β η (8) for = Q 0, r 0 and Q o 0, r o 0 for =,2,,m. (9) =

4 In addition, F X (r ) and F X X (r o r ) are the cumulatie density functions (CDFs) of f X (x ) and f X X (x x ), respectiely. EAC and EAC 2 are gien as EAC (Q,r ) = A ( ) D Q +C D +h Q 2 +r µ + p D (x r )f X (x)dx, (0) Q r EAC 2 ( Q o, r o ) = = [ A o D o +C o D o +h o Q o ( Qo 2 +r o µ o ) + p od o Q o r o (x r o )f X X (x r )dx where µ o is defined as in Eq. (5). It is readily shown that the constraint of Eq. (8) can be rewritten as C (Q +r )+ C o (Q o +r o )+κ F X (r )+ κ o F X X (r o r ) β +µ Y +z η σ Y, (2) = = ] 4,() where m µ Y C µ + C o µ o, σy 2 C2 σ2 + Co 2 σ2 o, (3) = = and z η = F ( η) with F ( η) being the inerse of the standard Gaussian CDF of the probability η. With the model, the Lagrangian function J using the Lagrangian relaxation can be written as J = EAC (Q,r )+EAC 2 ( Q o, r o )+λ C (Q +r )+ C o (Q o +r o ) +κ F X (r )+ = κ o F X X (r o r ) (β +µ Y +z η σ Y ). (4) = The first order necessary conditions can be achieed by differentiating J with respect to Q, r, Q o, r o, and λ: J = A od o Q o Q 2 + h o o 2 p od o Q 2 (x r o )f X X (x r )dx+λc o = 0 (5) o r o J = h o + p od o (x r o )f X X r o Q o r (x r )dx+λc o +λκ o F X X o r o r (r o r ) = 0 (6) o J = A D Q Q 2 + h 2 p D Q 2 (x r )f X (x)dx+λc = 0 (7) r J r = h + p D Q r (h o +λc o ) = J λ = C (Q +r )+ (x r )f X (x)dx+ r ( σ ρ = = p o D o Q o ) +λc +λκ d dr F X (r )+ r o (x r o ) r f X X (x r )dx = λκ o r F X X (r o r ) = 0 (8) [C o (Q o +r o )]+κ F X (r )+ κ o F X X (r o r ) (β +µ Y +z η σ Y ) = 0.(9) Note that Eqs. (5) (8) are simultaneous equations for Q o,r o,q, and r. = C. Simplification of equations using normalization Gien that Q o and r o are coupled with Q and r from Eqs. (5) (8), the equations are intractable to sole directly. For example, two equations [Eqs. (5) and (6)] contain three ariables Q o, r o, and r. It turns out, howeer, that the normalizations discussed in Section II A can not only simplify the arious expressions in Eqs. (5) (8), but also, more importantly, decouple the equations for the optional components [Eqs. (5) and (6)] from the equations for the anilla box [Eqs. (7) and (8)].

5 5 FIG. : The list of expressions that are used to simplify the first order necessary condition. Similarly to the normalization of Eq. (4), we further define the normalized reorder points as z r µ and z o r o µ o σ o = (r σ o µ ) ρ (r µ ) σ ρ 2, (20) where we hae used the definition of µ o of Eq. (5). Note that the normalized reorder point z o is a function of r o and r. With the normalization defined in Eqs. (4) and (20), the arious expressions in Eqs. (5) (8) can be simplified as shown in Fig.. We derie the simplified expressions in detail in VI. With the normalization, Eqs. (5) (9) can be re-expressed as Here, we define J = A od o Q o Q 2 + h o o 2 p od o Q 2 σ o L(z o )+λc o = 0 o (2) J = h o p od o G(z o )+λc o +λ κ o f(z o ) = 0 r o Q o σ o (22) J = A D Q Q 2 + h 2 p D Q 2 L(z )+λc = 0 (23) J = h p D G(z )+λc +λ κ f(z ) = 0 r Q (24) J λ = C (Q + z )+ C o (Q o +σ o z o )+κ F(z )+ κ o F(z o ) f(z) e z2 /2, G(z) =0 (β +µ Y +z η σ Y ) = 0. (25) z f(t) dt, F(z) G(z), and L(z) f(z) zg(z). (26) Note that, owing to the normalization, Eqs. (5) (8) are decoupled into two sets of equations: Eqs. (2) and (22) and Eqs. (23) and (24). Furthermore, each set of equations is identical and differs from the other only by the subscript. This implies that Eqs. (2) and (22) can be soled independently from Eqs. (23) and (24). This is expected because the PDF of the biariate Gaussian distribution [Eq. (3)] can be expressed as the product of the PDF of the normalized ariables Z and Z [Eq. (6)]. Thus, Eqs. (2) and (22) can be soled for Q o and z o directly. Similarly, Eqs. (23) and (24) can be soled for Q and z. Once z o and z are found, we can use Eqs. (4) and (20) to get r o and r. In the next section, we describe how to sole the set of Eqs. (2) (24) under the constraint of Eq. (25). =0 III. PROPOSED METHOD TO SOLVE THE MODEL SYSTEM Similar to the approach by [9], the proposed method to sole the model system consists of two parts. First, we regard Eqs. (2) (24) as a subproblem for a gien λ. Second, we repeatedly sole the subproblem until we find the solution λ to Eq. (25).

6 6 A. Procedures for soling the subproblem Foragienλ, weregardeqs.(2) (24) asasubproblem and solethem forq o, z o, Q, and z. The crucialpoint is that soling Eqs. (2) and (22) for z o and Q o is independent of soling Eqs. (23) and (24) for z and Q, respectiely. Gien that G(z o ) can be numerically ealuated for any z o, we note that Eqs. (2) and (22) are functions of Q o and z o only. By eliminating Q o from Eqs. (2) and (22), we obtain an equation for z o. In particular, Eqs. (2) and (22) can be rewritten, respectiely, in terms of Q o as Q 2 o = A od o +p o D o σ o L(z o ) h o /2+λC o and Q o = p o D o G(z o ) h o +λc o +λ(κ o /σ o )f(z o ). (27) We can eliminate Q o from Eq. (27), resulting in an equation for g o (z o ) in terms of z o only: g o (z o ) p o D o G(z o ) h o +λc o +λ(κ o /σ o )f(z o ) /2 Ao D o +p o D o σ o L(z o ) = 0. (28) h o /2+λC o There exists a unique solution z o = z o of Eq. (28) if the following three conditions are satisfied: (a) g o (z o ) is a continuous function in z o. (b) g o (z o ) is strictly monotonic in z o. (c) There exist two distinct alues z and z 2 such that g o (z ) g o (z 2 ) < 0. Condition (a) is immediately satisfied because f(z o ), G(z o ), and L(z o ) are continuous in z o. The following theorem satisfies condition (b). Theorem. g o (z o ) is strictly decreasing function in z o. The proof is gien in VII. Finally, condition (c) imposes a restriction on the alues of the parameters. Gien that g o (z o ) is a continuous and strictly decreasing function in z o from conditions (a) and (b), by assuming that z o > 0, the following two inequalities should be satisfied to meet condition (c): g o (0) > 0 and g o (+ ) < 0. It is easy to see that the condition g o (+ ) < 0 is satisfied by noting that G(+ ) = 0 and 0 < L(+ ) <. Therefore, the alues of the parameters hae to satisfy g o (0) > 0. That is, g o (0) = p o D o 2h o +λc o +λκ o /( σ o ) { A o D o +p o D o σ o / } /2 > 0. (29) h o /2+λC o The alues of the parameters should satisfy this inequality for Eq. (28) to hae a a unique solution. Gien that Eq. (28) is a function of z o only, one can use a simple search technique, such as a bisection method (see, for instance, [[]]), to sole it at least numerically if not analytically. Once the solution z o of g o(z o ) = 0 is obtained, we can substitute it back into Eq. (27) to get Q o, the solution of Q o. We repeat the same procedure for =,2,,m to obtain the normalized reorder point and order quantity for the optional components. The normalized reorder point z and order quantity Q for the anilla box can be obtained by applying a method similar to that used for the optional components. That is, Eqs. (23) and (24) can be rewritten respectiely in terms of Q o as The existence of a unique solution z of Q 2 = A D +p D L(z ) h /2+λC and Q = g (z ) p D G(z ) h +λc +λ(κ / )f(z ) p D G(z ) h +λc +λ(κ / )f(z ). (30) /2 A D +p D L(z ) = 0 (3) h /2+λC can be proen in a similar fashion to the case for zo of Eq. (28). Thus, we can sole numerically for z = z; subsequently, we can sole for Q = Q by substituting z into either Eq. (23) or (24).

7 7 TABLE I: Parameters for the anilla box and the two optional components with η = so that z η =.3. A C D h p κ µ β Vanilla box , ,000 Optional component A o C o D o h o p o κ o µ σ ρ B. Algorithm for soling the model system The solution scheme for the subproblem discussed in Section III A reduces the set of Eqs. (2) (25) to one equation [Eq. (32)] with one unknown λ, which is implicitly dependent on the ariables: g(λ) = C (Q + z )+ C o (Q o +σ o z o )+κ F(z )+ κ o F(z o ) (β +µ Y +z η σ Y ). (32) =0 If g(λ) > 0, then the constraint is iolated; otherwise (that is, g(λ) 0), the constraint is satisfied. By using the Lagrangian relaxation [9], the proposed algorithm to obtain the optimal reorder points and order quantities for the optional components and anilla box [r o, Q o, r, Q ] is as follows: Step : Find λ and λ 2 such that g(λ ) > 0 and g(λ 2 ) < 0. Step 2: For each λ and λ 2, sole the subproblem as follows: Step 2(a): For = to m, numerically sole Eq. (28) for z o to obtain zo, and substitute z o into Eq. (27) to get Q o. Step 2(b): Numerically sole Eq. (3) for z to obtain z, and substitute z into Eq. (30) to get Q. Step 3: Let λ new = (λ + λ 2 )/2 and find Q o, z o, Q, and z from Steps 2(a) and 2(b). If g(λ new) > 0, then let λ = λ new ; otherwise let λ 2 = λ new. Step 4: Repeat Steps 2 and 3 until g(λ ) < ǫ or g(λ 2 ) < ǫ, where ǫ is a predetermined error. Step 5: Use Eq. (20) to get r and r o from z and z o, respectiely. =0 IV. EXPERIMENTAL RESULTS AND DISCUSSION We illustrate the performance of the proposed method by an experiment that consists of one anilla box and two optional components (i.e., m = 2). The parameters for the anilla box and two optional components are listed in Table I. They are the same as the parameters used in [6] and ǫ = 0 7. Table II lists the solution to the model for a gien input from Table I together with the results of [6]. We also perform a sensitiity analysis of the order quantities and reorder points for the optional components with respect to the correlationcoefficient. It should be noted that Q and r are independent of ρ from Eqs. (30) and (3). For the sensitiity analysis, we first need to find the behaior of z o as ρ aries. Figure 2 shows that z o is almost constant with respect to ρ although z o decreases slightly for large alue of ρ. This implies that the normalized reorder points are insensitie to the correlation coefficient. Figure 3 shows the behaior of order quantities Q o of the optional components as the correlation coefficient ρ aries. From the first equation in Eq. (27), we see that the dependence of Q o on ρ stems from σ o L(z o ). Because z o is more or less insensitie to ρ [Fig. 2], so is L(z o ). Thus, considering ρ dependence only, we hae Q o σ o σ ρ 2. (33) This implies that Q o decreases as the absolute alue of ρ increases and Q o reaches its maximum when ρ = 0 as shown in Fig. 3. The maximum of Q o at ρ = 0 can also be proed as follows. From the first equation in Eq. (27), it can be readily shown that dq o p o D o L(zo )+z o G(z o ) σ 2 = ρ. (34) dρ {h o +2λC o } σ o Q o

8 8 TABLE II: Solutions for Q, r, Q o, and r o, together with λ and the total cost. Q r Q o r o Q o2 r o2 λ EAV( Q, r) Proposed method ,536,070 [6] ,536, optional optional 2 z o ρ FIG. 2: Plots of the normalized reorder points z o for the optional components ersus the correlation coefficient ρ. Thus, Q o is extreme when ρ = 0. Furthermore, d 2 Q o p o D o L(zo )+z o G(z o ) dρ 2 = σ < 0 (35) {h o +2λC o } Q o ρ=0 implies that Q o has its maximum when ρ = 0. Figure 4 shows the behaior of reorder points r o of the optional components as the correlation coefficient ρ aries. Unlike the order quantity, r o reaches its maximum at a positie alue of ρ o. For the behaior of r o, we can rewrite Eq. (20) as follows: σ r o = µ +ρ (r µ )+σ ρ 2 σ z o. (36) Because z o is almost independent of ρ, as ρ increases, the second term on the right-hand side of Eq. (36) also increases while the third term decreases. Thus, there is a trade-off between the second and the third terms, resulting in an optimum alue of r o. Furthermore, the maximum r o occurs when ρ max r µ σ 2 z2 o +σ2 (r µ ) 2. (37) This implies that ρ max depends on the safety stock r µ of the anilla box. Gien that the safety stock is a positie quantity, we hae 0 < ρ max < unlike the maximum of Q o. V. SUMMARY AND CONCLUSION In this paper, we presented an efficient method for finding an optimal solution (Q, r) to a correlated multi-item continuous reiew inentory model in which a biariate Gaussian probability distribution is used as a correlation

9 9 optional optional Q o ρ FIG. 3: Plots of the order quantities Q o for the optional components ersus the correlation coefficient ρ r o 8 6 r o ρ ρ FIG. 4: Plots of the reorder points r o for optional component ersus the correlation coefficient ρ while ρ 2 = 0.5. Inset: Same plot for optional component 2. between the anilla box and an optional component. By normalizations of the random ariables for the demands, we showed that the biariate Gaussian PDF can be expressed as a product of two independent Gaussian PDFs, which implies that the normalized random ariables are uncorrelated. To demonstrate the usefulness of the normalization, we soled a multi-item continuous reiew inentory(q, r) model in which the anilla box and optional components are correlated under stochastic payment and budget constraints. With normalization, we showed that the set of equations coupled with the anilla box and optional components was decoupled into sets of equations for the normalized quantities. Furthermore, each set of decoupled equations was reduced to a closed form and could be soled numerically without any approximation. We also performed the sensitiity analysis in terms of the correlation. We found that the order quantity and the reorder point of optional components depended on the strength of the correlation as we expected. In particular, we showed that the order quantity of an optional component reached its maximum when there was no correlation between the anilla box and the optional component. In addition, the reorder point of an optional component reached

10 a maximum that depended on the safety stock of the anilla box. The proposed method can be used as a dependable method for a generalized multi-item continuous reiew inentory model with complicated interactions among items. It would be interesting to inestigate how far the proposed method can be applied to other types of correlation, such as a multiariate Gaussian or other multiariate distributions. 0 Acknowledgments This research was supported by the Basic Science Research Program through the National Research Foundation of South Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF ). VI. APPENDIX I: SIMPLIFICATION OF VARIOUS FUNCTIONS In this appendix, we show how to simplify the expressions in Eqs. (5) (8). To this end, we define arious normalized ariables as follows: z x µ and z r µ, (A-) where z x µ o σ o In addition, from Eq. (3), we hae f X (x) = exp { (x µ ) 2 } σ and z o r o µ o σ o, (A-2) σ o σ ρ 2 and µ o µ +ρ σ (r µ ). 2σ 2 and f X X (x r ) = exp (x µ o) 2 σo 2σo 2 (A-3). (A-4) Finally, the G-function is defined as G(z) z e z2 /2 dz. (A-5) A. Ealuation of r (x r )f X (x)dx and d dr r (x r )f X (x)dx From Eqs. (A-) and (A-4), we hae (x r )f X (x)dx = σ (z z )e z 2 /2 dz = e z2 /2 z G(z ). (A-6) r z By using the aboe result, its deriatie becomes d dr r (x r )f X (x)dx = dz dr d dz e z2 /2 z G(z ) = G(z ). (A-7) B. Ealuation of r o (x r o)f X X (x r )dx and r o r o (x r o)f X X (x r )dx From Eqs. (A-2) (A-4), we hae (x r o )f X X (x r )dx = (x r o ) exp (x µ o) 2 r o r o σo 2σo 2 dx = σ ρ 2 (z z o ) e z2 /2 dz = σ o e z2 o /2 z o G(z o ). (A-8) z o

11 By using the aboe result, its deriatie becomes r o r o (x r o )f X X (x r )dx = σ o dz o dr o z o e z2 o /2 z o G(z o ) = G(z o ). (A-9) C. Ealuation of r o (x r o) r f X X (x r )dx By using Eq. (A-4), the partial deriatie becomes f X X r (x r ) = dµ o exp (x µ o) 2 dr µ o σo 2σo 2 = ( ) σ (x µo ) ρ exp (x µ o) 2 σ 3 o 2σo 2. Thus, by using Eqs. (A-2) and (A-4), we hae ( ) σ ρ (x r o ) f X X r o r (x r )dx = ( ) σ = ρ (z z o ) z e z2 /2 dz = z o { (x r o ) (x µ o) exp r o σ 3 o ) G(z o ). ( ρ σ (x µ o) 2 2σ 2 o } dx (A-0) D. Ealuation of d dr F X (r ), r F X X (r o r ), and r o F X X (r o r ) By using Eqs. (A-) and (A-4), we hae Thus, F X (r ) = r σ exp d F X (r ) = dz d z dr dr dz Similarly, by using Eqs. (A-2) (A-4), we hae F X X (r o r ) = ro σo exp { (x µ ) 2 { 2σ 2 e z2 /2 dz = (x µ o) 2 2σ 2 o } z dx = e z2 /2 dz. } σ e z2 /2. dx = zo e z2 /2 dz. (A-) Thus, r F X X (r o r ) = dz o dr z o zo ( ) e z2 /2 σ dz = ρ e z2 o /2. σo (A-2) Finally, by using Eqs. (A-2) and (A-4), we hae r o F X X (r o r ) = dz o dr o z o zo e z2 /2 dz = σo e z2 o /2. (A-3) VII. APPENDIX II: PROOF OF THEOREM For breity, we drop the subscript o for the optional component and assume λ > 0. To proe that g(z) of Eq. (28) is monotonically decreasing in z, it suffices to show that for any two alues of z such that z < z 2, g(z) satisfies g(z ) > g(z 2 ). We start with g(z) pdg(z) h+λc +λκ/σf(z) /2 AD +pdσl(z), (B-) h/2+λc

12 2 where the loss function L(z), the standard Gaussian PDF f(z), and G(z) are defined, respectiely, as L(z) σ e z2 /2 zg(z), f(z) e z2 /2, and G(z) e t2 /2 dt. From Eq. (B-), we hae { g(z ) g(z 2 ) = pdg(z ) h+λc +λκ/σf(z ) pdg(z 2 ) h+λc +λκ/σf(z 2 ) } + pdσ h/2+λc z L(z2 ) L(z ). (B-2) Since L(z ) < L(z 2 ) for z < z 2, the second term on the right hand side of Eq. (B-2) is positie. The first term can be rewritten as { } pdg(z ) h+λc +λκ/σf(z ) pdg(z 2 ) h+λc +λκ/σf(z 2 ) (h+λc){g(z ) G(z 2 )}+λκ/σ{g(z )f(z 2 ) G(z 2 )f(z )} = pd. (B-3) {h+λc +λκ/σf(z )}{h+λc +λκ/σf(z 2 )} From the definition of G(z), we hae G(z ) > G(z 2 ) when z < z 2. Thus, to show that g(z ) > g(z 2 ) when z < z 2, we are left with showing that G(z )f(z 2 ) G(z 2 )f(z ) is positie. To this end, consider the following expression: h(z ) h(z 2 ) G(z )f(z 2 ) G(z 2 )f(z ) f(z )f(z 2 ) = G(z ) f(z ) G(z 2) f(z 2 ). (B-4) Gien that f(z) > 0, to show g(z ) > g(z 2 ) when z < z 2 is equialent to proing that h(z ) > h(z 2 ) when z < z 2. That is, h(z) is a monotonically decreasing in z. Now, we express G(z)/f(z) in terms of an infinite series in z by using [2] where h(z) G(z) f(z) = ( ) n(2n )!! z 2n+ > 0, (B-5) n=0 By letting m = n, h(z) can be rewritten as h(z) = (2n )!! 3 5 (2n 3) (2n ). ( ) n(2n )!! = z + ( ) n(2n )!! n=0 z 2n+ n= z 2n+ = z + ( ) m+(2m+)!! z 2m+3 = z ( ) m(2m+)!! z 2m+3 > 0. (B-6) m=0 Since h(z) > 0 for all z, we hae, from Eq. (B-6), m=0 ( ) m(2m+)!! z 2m+3 < z. (B-7) m=0 To show that h(z) is monotonically decreasing in z, we will show that dh(z)/dz < 0. From Eq. (B-5), the deriatie becomes dh(z) = ( ) n (2n )!!( ) (2n+) dz z 2n+2 = ( z) ( ) n(2n+)!! z 2n+3. (B-8) n=0 By using Eq. (B-7), dh(z)/dz can be expressed as dh(z) dz = ( z) n=0 ( ) n(2n+)!! z 2n+3 < ( z) <. (B-9) z n=0

13 3 Since dh(z)/dz < 0 for all z, we hae h(z 2 ) < h(z ) when z 2 > z. That is, This completes the proof. h(z ) h(z 2 ) G(z ) f(z ) G(z 2) f(z 2 ) = G(z )f(z 2 ) G(z 2 )f(z ) f(z )f z2 > 0. (B-0) [] Abdel-Malek L, Areeratchakul N (2007) A quadratic programming approach to the multi-product newsendor problem with side constraints Euro. J. of Oper. Res. 76: [2] Abramowitz M, Stegun I, eds. (972) Handbook of Mathematical Functions, (Doer, New York, 972) p [3] Byrkett D (98) An Empirical Ealuation of Further Approximations to An Approximate Continuous Reiew Inentory Model Naal Research Logistics Quarterly 28: [4] Chen Y, Beaulieu N (2009) A simple polynomial approximation to the Gaussian Q-function and its application IEEE Communication Letters 3: [5] Tze-Ming Chen, Tai-Yue Wang, Jui-Ming Hu, Optimal ordering policy in simple inentory models with modularization and postponement, Computers and Industrial Engineering (CIE), th International Conference on, July 200, pp.-4 [6] Das C (975) Some Aids for Lot-Size Inentory Control under Normal Lead Time Demand AIEE Transactions 7: [7] Das C (985) On the Solution of Some Approximate Continuous Reiew Inentory Models Naal Research Logistics Quarterly 32: [8] Ernst R, Kamrad, B (2000) Ealuation of supply chain structures through modularization and postponement Euro. J. of Oper. Res. 24: [9] Ghalebsaz-Jeddi B, Shultes B, Hai R (2004) A multi-product continuous reiew inentory system with stochastic demand, backorders, and a budget constraint Euro. J. of Oper. Res. 58: [0] Herron D (967) Inentory Management for Minimum Cost Management Science 4:B [] Kincaid D, Cheney W (2002) Numerical Analysis: Mathematics of scientific computing, (Brooks/Cole, Pacific Groe, 2002). [2] Kundu A, Charkrabarti T (202) A multi-product continuous reiew inentory system in stochastic enironment with budget constraint Optim. Lett 6: [3] Swaminathan JM, Tayur S (998) Managing broader product lines through delayed differentiation using anilla boxes. Management Science 44:S6 S72. [4] Swaminathan JM, Tayur S (999) Managing design of assembly sequences for product lines that delay product differentiation. IIE Transactions 3: [5] Wang T, Hu J (2008) An inentory control system for products with optional components under serice leel and budget constraints. Euro. J. of Oper. Res. 89:4-58. [6] Wang T, Hu J (200) Heuristic method on soling an inentory model for products with optional components under stochastic payment and budget constraints, Exp. Sys. with App. 37: