OPTIMAL POLICIES FOR INVENTORY MODEL WITH SHORTAGES, TIME-VARYING HOLDING AND ORDERING COSTS IN TRAPEZOIDAL FUZZY ENVIRONMENT

 

Pavan Kumar

Mathematics Division, School of Advanced Science & Language

VIT Bhopal University, India

E-mail: pavankmaths@gmail.com

 

Submission: 12/19/2019

Revision: 5/14/2020

Accept: 6/3/2020

 

ABSTRACT

This paper proposes the optimal policies for a fuzzy inventory model considering the holding cost and ordering cost as continuous functions of time. Shortages are allowed and partially backlogged. The demand rate is assumed in such to be linearly dependent on time during on-hand inventory, while during the shortage period, it remains constant. The inventory problem is formulated in crisp environment. Considering the demand rate, holding cost and ordering cost as trapezoidal fuzzy numbers, the proposed problem is transformed into fuzzy model. For this fuzzy model, the signed distance method of defuzzification is applied to determine the average total cost (ATC) in fuzzy environment. The objective is to optimize the ATC and the order quantity. One solved example is provided in order to show the applicability of the proposed model. The convexity of the cost function is verified with the help of 3D-graph.

Keywords: Inventory Model; Partial Backlogging; Time Dependent Demand Rate; Signed Distance Method.

Glossary

TrFN: Trapezoidal Fuzzy Number

EOQ: Economic Order Quantity

EPQ: Economic Production Quantity

ATC: Average Total Cost

pdf: Probability Density Function

 

1.       INTRODUCTION

            In inventory problems, the ordering cost is associated with ordering and receiving an order, and normally possesses the components- clerical/labor costs of processing orders, inspection and return of poor quality products, transport costs, handling costs. In practical life situations, all the above cost components depend upon time. Therefore, to consider the time dependant ordering cost will undoubtedly be a very close approximation to real case inventory problems.

            In literature, many inventory models with deteriorating items were studied by a large number of researchers (Dave & Patel, 1981; Abad, 1996; Chung & Ting, 1993; Goyal & Giri, 2001; Shah & Shukla, 2009). The optimal pricing and lot-sizing was proposed assuming the conditions of perishability and partial backordering (Abad, 1996). The concept of a reseller under partial backlogging was presented to determine the optimal price and order-size for the inventory model (Abad, 2001). The deterioration may be represented by a random variable with a known probability distribution function (pdf).

            An EOQ model for items with weibull distribution deterioration was presented with shortages and trended demand (Jalan; Giri; Chaudhary, 1996). Inventory models with ramp type demand rate, partial backlogging and Weibull deterioration rate was presented (Skouri & Konstantaras; Ganas, 2009). The ramp-type demand rate was studied by large number of researchers. An EOQ inventory model for Weibull distributed deteriorating items was proposed under ramp type demand and shortages (Mandal, 2010).

        In the real life business problems, the demand rate of any item is always variable. Several inventory models have been established by considering time-dependent demand. The time dependent demand was considered in some inventory models (Dave & Patel, 1981; Chung & Ting, 1993; Wu & Cheng, 2005; Mandal, 2010; Dutta & Kumar, 2015a). A deteriorating inventory model for (T, si)-policy was studied for time proportional demand function (Dave; Patel, 1981). By considering a linear trend in demand, a heuristic was studied for replenishment of deteriorating items for an inventory problem (Chung & Ting, 1993). By introducing the inflation when supplier credit linked to order quantity, an inventory model was proposed for time-dependent demand rate (Tripathi, 2011). 

        In inventory models, one important factor is shortages. In literature, two types of shortages are considered: partial backlogging and full backing. An EOQ model was proposed for deteriorating items with time varying demand and partial backlogging (CHANG; DYE, 1999). By considering the partial backlogging shortages, an inventory model for deteriorating items with exponential declining demand was formulated (Wu & Cheng, 2005). A comparison between two pricing and lot-sizing models with partial backlogging and deteriorated items was proposed by (Teng, Ouyang & Chen, 2007). A deteriorating inventory model for time dependent demand and holding cost with partial backlogging was introduced by (Mishra & Singh, 2011).

            The determination of optimal selling price as well as the lot size was presented for a varying rate of deterioration and exponential partial backlogging (Dye, 2007a). Afterwards, the joint pricing and ordering was proposed for a deteriorating inventory with partial backlogging (Dye, 2007b). In several cases, the demand rate depends on the selling price of the items. Assuming the price dependent demand rate, an inventory model was presented for deteriorating items and time varying holding cost (Roy, 2008).

            In inventory models, the concept of two-level trade credit was also introduced. An EOQ model with non instantaneous receipt and exponential deteriorating item was introduced under two-level trade credit (Liao, 2008). A deteriorating inventory model for waiting time partial backlogging was introduced by (Shah & Shukla, 2009). An inventory model was proposed by introducing the generalized type demand, deterioration and backorder rates (Hung, 2011). By introducing the multiple-market demand, an optimal production-inventory model was presented for deteriorating items (He, Wang & Lai, 2010). 

            In general, the parameters for holding cost, demand, ordering cost, deterioration rate, production cost and idle time cost are not specifically known. Due to uncertainty, we may assume the representation of some or all the parameters by fuzzy sets. The concept of fuzzy sets was introduced first time by (Zadeh, 1965). Some applications of fuzzy set theory to mathematical programming as well supply chain management were studied by  (Zimmermann, 1985), (Wang & Shu, 2005).

            The (Q, r) inventory model with backorder discount in fuzzy demand and fuzzy ordering cost was proposed by (Ho, et al., 2007). The fuzzy setup cost as well as fuzzy ordering cost was introduced for determining the single-vender and single-buyer inventory strategy (Yang, Tu & Wang, 2007). An interval-type uncertainty was introduced in inventory parameters. An interval number approach was applied to a partial backlogging inventory model for deteriorating items with time-varying demand and holding cost (Dutta & Kumar, 2015a).

            The global criteria method was suggested to solve the inventory model with variable holding cost and partial backlogging under interval uncertainty (Kumar & Keerthika, 2018). Uncertainty of neutrosophic type was also proposed. The optimal value of neutrosophic complex programming problem was determined by using lexicographic order (Khalifa, Kumar & Smarandache, 2020).

            The concepts of parabolic holding cost with salvage value were also introduced. An inventory planning problem was presented for time-varying linear demand and parabolic holding cost with salvage value (Kumar, 2019). Some applications of optimizations techniques were also suggested to various fields. In the field of a micro-company of the graphic sector, an application of a mathematical model for cost minimization was studied (Rodrigues, Marins & Souza, 2017).

            In the field of inventory, an application of fuzzy goal programming approach was implemented to multi- objective linear fractional inventory model (Dutta & Kumar, 2015b). In joint inventory-production problems, the production depends on demand. Moreover, the demand depends on the number of population of customers, selling price as well as advertisement. A fuzzy EPQ model for non-instantaneous deteriorating items was presented, where the production rate depends on demand which is proportional to population, selling price as well as advertisement (Biswas & Islam, 2019).

1.1.          Motivation:

      Several researchers studied the inventory problem with constant deterioration, time-sensitive demand and holding cost (Tripathi, 2011; Dutta & Kumar, 2015a; Kumar & Keerthika, 2018), etc. For instance, Dutta and Kumar (2015a) formulated such a model, and applied the interval number method to solve it.  Kumar and Keerthika (2018) also studied a similar model using uncertainty in model variables and parameters. They did not consider the trapezoidal fuzzy number to deal the uncertainties. Also, the ordering cost may be considered as time dependent along with the time-sensitive holding cost.

       In present research paper, we propose a fuzzy inventory model with shortages. We consider the holding and ordering costs as continuous functions of time. Our objective is to determine the optimal control policies for the proposed model.

1.2.          Main Contribution:

The main contributions of the present paper are summarized below:

·      Introducing the time-sensitive holding cost as well as ordering cost.

·      Applying the solution method based on defuzzification of trapezoidal fuzzy numbers.   

·      Convexity is justified with the help of a 3D-graph..

·      Based on the results, some managerial implications are proposed.

       The remaining outlay of the paper is organized as follows. In section 2, we describe some definitions. In section 3, the problem description and formulation is provided. In section 4, a fuzzy mathematical model of the proposed problem is developed. In section 5, the solution procedure is described. In section 6, a solved example is given to explain the procedure. In section 7, some special cases are discussed. In section 8, managerial implications and limitations are explored. In the last, the conclusions and future directions are reported in section 9.

2.       SOME DEFINITIONS

·       Definition 1: (Dutta & Kumar, 2015b)

         A Trapezoidal Fuzzy Number (TrFN), denoted by  is defined by the membership function as:

                                               =

·       Definition 2: (Dutta & Kumar, 2015b)

         A fuzzy set is called in LR-Form, if there exist reference functions L (for left), R (for right), and scalars m > 0 and n > 0 with membership function

                                =

Where 𝜎 is a real number called the mean value of  m and n are called the left and right spreads, respectively. The functions L and R map  and are decreasing. A LR-Type fuzzy number can be represented as =

·       Definition 3: (Dutta & Kumar, 2015b)

        The fuzzy arithmetic operations under the functional principle for trapezoidal fuzzy numbers are described below. Let   () be two trapezoidal fuzzy numbers, and let k be any scalar. Then arithmetical operations are defined as:

(i)     Fuzzy Addition:   = ().

(ii)   Fuzzy Multiplication:   = ().

(iii) Fuzzy Subtraction:   = ()

(iv)  Fuzzy Division:   = ()

(v)    Fuzzy Scalar Multiplication:  k =  

·       Definition 4: (Dutta & Kumar, 2015b)

         Let  be a fuzzy set defined on R. Then the Signed Distance of  is defined as:

                   () =             

where [] = [a + ( b a)α, d (d c)α],  α [0, 1], is α cut of fuzzy set, which is a close interval.

 

3.       PROBLEM DESCRIPTION AND FORMULATION

            The main objective of this paper is to minimize the total cost of the proposed inventory model, that possesses the following features: (i) the demand rate is a time dependent function, (ii) the holding cost is an increasing linear function of the storage time, and (iii) the ordering cost depends on the total time elapsed up to the beginning of each cycle.

3.1.          Notations

The following notations are used in this paper.

Holding cost per unit per time unit (Time dependant).

Purchase cost per unit

Shortage cost per unit per time unit

Cost of lost sales per unit

           : Deterioration rate

           : Cycle time (decision variable) 

Time at which shortage starts, i.e., inventory exhausted time, (decision variable), 0 ≤

Length of waiting time

      : Maximum inventory level during a cycle of length T

         : Maximum amount of demand backlogged during a cycle of length T

           : (=  + ) order quantity during a cycle of length T

                          : Time dependant ordering cost per order.

          : Inventory holding cost per cycle

          : Deterioration cost per cycle

                            : Shortage cost per cycle

                             : Lost sales cost per cycle

                   C(, T)  : Average total cost per time unit per cycle

         : Optimal value of X, where X is any variable or parameter.

        : Inventory level at time, 0 ≤

       : Inventory level during the period [0, ]

       : Inventory level during the period []

3.2.          Assumptions

The proposed model is based on the following assumptions:

       I.          Inventory system involves only one item.

     II.          Planning horizon is infinite. Each replenishment cycle is of same length, T. Only one order is placed for the period [0, T]. No lead time between placing an order and arrival of the order.

   III.          Deteriorating rate, q (0 < q < 1), is constant and occurs as soon as items are received into inventory. There is no replacement or repair of deteriorated units.

  IV.          Demand during the on-hand inventory, the demand rate, R(t) increases linearly with respect to time, and during the shortages period, the demand rate becomes constant:                  

R(t) = , where  are arbitrary  constants.

    V.          Shortages are allowed and partially backlogged with backordered rate

B(t) =   ,                                                         

            The backlogging parameter δ>0, and  For the special case with d = 0, B(t) = 1, that is, the fully

backlogged case. In the proposed model, we assume δ<1 to approximate by Maclaurin’s series.

  VI.          Holding cost,  per unit per time unit is time dependant and assume its function form is

(t) = ht, where h > 0 is the holding cost scale parameter.

VII.          The ordering cost, depends on the total time elapsed up to the beginning of each cycle and is taken as

                

                      where   ,  and   is the ordering cost at time zero.

3.3.          MATHEMATICAL MODEL

            The inventory is replenished at time t = 0, when the inventory level is at its maximum, . Afterwards, the inventory level begins to decrease during the period [0, ], and eventually becomes zero, when t = . Further, during the period [, T], there is partially backlogging. The governing differential equations are:

              + =,   0 ≤ t ≤                                                                                    (1)

             =  ,             ≤ t ≤ T                                                                                     (2)

            With boundary conditions:                                                                                    

            , the inventory level decreases due to the demand as well as deterioration, and the inventory level is governed by (1). Using the boundary conditions, the solution of (1) is given by

             ) =  [(t) ] ,  0 ≤ t ≤                                                       (3)

            Therefore, the maximum level of the proposed inventory model, for each cycle, can be written as follows:  

              =  =]                                                                           (4)                                

            , the inventory level is sensitive with the demand. When the demand is more, the inventory level starts to decreases, and vice-versa. However, we assume that a fraction of demand is backlogged. The level of the inventory is governed by (2). Applying the boundary conditions, the solution of (2) can be written as:

              ) =  [log{1+ d(T t)} log{1+ d(T)}],  ≤ t ≤ T                                        (5)

            As per the calculations illustrated in Appendix-A, we can write the order quantity as follows:

             Q = + b[T ]                                                                                       (6)

            Hence, the average total cost per unit time per cycle is 

            Þ C(, T) =   { + b()[]}                       

            Applying the Macraurin series, and ignoring higher power terms, we get

           ÞC(,T)                 (7)                    

4.       FUZZY MATHEMATICAL MODEL

            In real-life situations, the decision maker has to face the ambiguity in model parameters. In this research paper, to face with the fuzzy type uncertainty and ambiguity in parameters of the inventory model, we consider the demand rate, holding cost and ordering cost as represented by trapezoidal fuzzy numbers (TrFN) as follows:

                                                                 (8)                            

            Therefore, we can the average total cost function per unit time in fuzzy environment case, as follows:

 () =

                 =

                     

                        

                 = {(   

                        

       = (), (say),                                                                                                                         (9)

    where   

             =  {(+ + }

             =  {(+ + }

             =  {(+ + }

             =  {(+ +                                                                                   

            (α) and (α), of trapezoidal fuzzy number (), can be written as follows:

(α) =   

=

            And    (α)  

 =

Þ=   (10)

And     (11)

            By applying the Signed Distance Method, the defuzzified value of average total cost, which is represented by the fuzzy number (), can be written as follows:

            (, T) = ,

=

=       (12)

            Hence, the optimization problem can be written as:

                                                                                               (13)

5.       SOLUTION PROCEDURE

            In this section, the solution procedure for the fuzzy model is presented. Since the objective is to minimize the total cost (, T) with respect to the decision variables  and T; the necessary conditions for minimizing the total cost are:

 = 0, and    = 0                                                                    (14)

            The sufficient conditions for maximizing the cost function (, T) are:

 > 0,   > 0, and   > 0                                                                                            (15)

            The satisfaction of these conditions confirms the convexity of the cost function ().

6.       SOLVED EXAMPLE

            Consider an inventory situation with crisp parameters having the following values:

            R(t) = ,

            where ,   = 4, and = 80, (t) == ,  q = 0.07,  δ = 0.5,   = 2.5, = 10,  = 5.

6.1.          Crisp Model. The optimal solutions are as follows:

 ,     1.5871 ,      

6.2.          Fuzzy Model. We set some trapezoidal fuzzy numbers as follows:

                                                                                  (16)

            For each of these input parameters, the variations in the values are arranged arbitrary and their defuzzified values are computed by applying the signed distance method. Applying the proposed solution procedure on package- MATHEMATICA-8.0, the optimal solutions are as follows:

 ,    1.550 ,     ,

6.3.          Convexity.

            The convexity of the cost function  is shown graphically in Figure 1 and Figure 2.

Figure 1: Convexity in Crisp Model

Figure 2: Convexity in Fuzzy Model

            The obtained results are summarized in Table-1.  The comparison for both models between the Time parameters is illustrated in Figure 3, while the comparison between the cost function values is illustrated in Figure 4.                               

Table 1: Optimal Results

 

Crisp Model

1.5871

Fuzzy Model

1.550

 

Figure 3: Comparison of Time parameter

Figure 4: Comparison of Cost parameter

7.       SPECIAL CASES

            In this section, we consider the following two special cases:

(i)         For  = 0, the proposed model changes from partially backlogged to fully backlogged. In fully backlogged inventory model, the total cost function is given by

            C(, T)                          (17)

(ii)       For = 0, the proposed model changes from linearly increasing type ordering cost to constant type ordering cost model. Subsequently, in constant ordering cost inventory model, the total cost function is given by                                               

            C(, T)                          (18)

8.       MANAGERIAL IMPLICATIONS AND LIMITATIONS

8.1.          Managerial Implications

            The following managerial implications are explored from the proposed research work: which is supportive of any operational managers in the manufacturing unit:

·       The model developed in this research work is supportive for the inventory management, where the holding cost and ordering cost are is continuously increasing with time, and where the decision maker targets to optimize the average total inventory costs.

·       The model considered the demand rate, ordering cost, and holding cost of the items with some ambiguity which is a natural case, and therefore, provides an approach to the inventory managers to optimize the various inventory costs with ambiguity.

8.2.          Limitations

            There exist some limitations of the proposed model which are listed below:

·       In the derivation of the proposed model, we used the approximations of exponential term. Therefore, this model will provide the approximate results only.

·       The examples, which we used to validate the proposed model, are hypothetical.

·       The proposed problem is limited to vagueness only. However, in real-world applications, the DM has to encounter the randomness and multiple choice cases.

·       The role of deterioration is assumed constant. However, in real-world applications, the DM has to face the time dependent deterioration.

9.       CONCLUSIONS AND FUTURE DIRECTIONS

            In this paper, a deteriorating inventory model with time-dependant demand rate and varying holding cost is presented. The shortages are permitted and partially backlogging case is assumed. The classical optimization technique is implemented over the cost function to obtain the optimal order quantity and optimal total cost. For the practical application of this model, a solved example is considered by using the  package. The proposed model can assist the manufacturer and retailer in accurately determining the economic order quantity, cycle time and total cost, to be used in inventory control of seasonal items.

            There are several directions for future research. We can discuss the sensitivity test to identify the most sensitive parameters in the model. Another direction of research is to consider the randomness in inventory parameter. Also, the demand and holding cost may be considered as quadratic time dependent. 

REFERENCES

Abad, P.  (1996). Optimal pricing and lot-sizing under conditions of perishability and partial backordering. Management Science, 42, 1093–1104.

Abad, P.  (2001). Optimal price and order-size for a reseller under partial backlogging. Computers and Operation Research, 28, 53–65.

Alamri, A., & Balkhi, Z. (2007). The effects of learning and forgetting on the optimal production lot size for deteriorating items with time varying demand and deterioration rates. International Journal of Production Economics, 107, 125–138.

Biswas, A. K., & Islam, S. (2019). A fuzzy EPQ model for non-instantaneous deteriorating items where production depends on demand which is proportional to population, selling price as well as advertisement. Independent Journal of Management & Production, 10( 5), 1679-1703. 

Chang, H., & Dye, C. (1999). An EOQ model for deteriorating items with time varying demand and partial backlogging. Journal of the Operational Research Society, 50, 1176–1182.

Chung, K. J., & Ting, P.  S. (1993). A heuristic for replenishment of deteriorating items with a linear trend in demand.  Journal of the Operational Research Society, 44, 1235-1241.

Dave, U., &Patel, L.  (1981).  (T, si)-policy inventory model for deteriorating items with time proportional demand. Journal of the Operational Research Society, 32, 137–142. 

Dutta, D., & Kumar, p.  (2015a). A partial backlogging inventory model for deteriorating items with time-varying demand and holding cost: An interval number approach. Croatian Operational Research Review, 6(2), 321-334. DOI:10.17535/crorr.2015.0025

Dutta, D., & Kumar, P.  (2015b). Application of fuzzy goal programming approach to multi- objective linear fractional inventory model, International Journal of Systems Science, 46(12), 2269-2278. 

Dye, C. (2007a).  Determining optimal selling price and lot size with a varying rate of deterioration and exponential partial backlogging. European Journal of Operational Research, 181, 668–678.

Dye, C. (2007b). Joint pricing and ordering for a deteriorating inventory with partial backlogging. Omega, 35(2), 184–189. 

Goyal, S., & Giri, B. (2001). Recent trends in modeling of deteriorating inventory. European Journal of Operational Research, 134, 1–16. 

He, Y., Wang, S., & Lai, K. (2010). An optimal production-inventory model for deteriorating items with multiple-market demand. European Journal of Operational Research, 203(3), 593–600. 

Ho, L-.H., Cheng, C-.L., Yang, M-.F., & Lo, M-.C. (2007). (Q, r) inventory model with backorder discount in fuzzy demand and fuzzy ordering cost. Journal of Information and Optimization Sciences, 28(4),  561-572.

Hung, K. (2011). An inventory model with generalized type demand, deterioration and backorder rates. European Journal of Operational Research, 208(3), 239–242. 

Jalan, A., Giri, R., & Chaudhary, K. (1996). EOQ model for items with weibull distribution deterioration shortages and trended demand. International Journal of System Science, 27, 851–855. 

Khalifa, H. A. E., Kumar, P., & Smarandache, F. (2020). On optimizing neutrosophic complex programming using lexicographic order. Neutrosophic Sets and Systems, 32, 330-343.

Kumar, P.  (2019). An inventory planning problem for time-varying linear demand and parabolic holding cost with salvage value. Croatian Operational Research Review, 10(2), 187-199.  DOI:10.17535/crorr.2019.0017

Kumar, P., & Keerthika, P.  S. (2018). An inventory model with variable holding cost and partial backlogging under interval uncertainty: Global criteria method. International Journal of Mechanical Engineering & Technology, 9(11),  1567-1578.

Liao, J. (2008). An EOQ model with non instantaneous receipt and exponential deteriorating item under two-level trade credit. International Journal of Production Economics, 113, 852–861.

Mandal, B. (2010). An EOQ inventory model for Weibull distributed deteriorating items under ramp type demand and shortages. OPSEARCH, 47(2), 158–165.

Yang, M-.F., Tu, H-.J., & Wang, C-.M. (2007). Determining the single-vender and single-buyer inventory strategy with fuzzy setup cost and fuzzy ordering cost. Journal of Statistics and Management Systems, 10(4),  499-510.

Mishra, V. K., Singh, L. S. (2011). Deteriorating inventory model for time dependent demand and holding cost with partial backlogging. International Journal of Management Science and Engineering Management, 6(4), 267-271.

Rodrigues, P. C. C., Marins, F. A. S., & Souza, F. B. D. (2017). Application of a mathematical model for the minimization of costs in a micro-company of the graphic sector. Independent Journal of Management & Production, 8(5), 676-692.

Roy, A. (2008). An inventory model for deteriorating items with price dependent demand and time varying holding cost. Advanced Modeling and Optimization, 10, 25–37.

Shah, N., & Shukla, K. (2009). Deteriorating inventory model for waiting time partial backlogging. Applied Mathematical Sciences, 3, 421–428.

Skouri, K., Konstantaras, S., & Ganas, I. (2009). Inventory models with ramp type demand rate, partial backlogging and weibull deterioration rate. European Journal of Operational Research, 192, 79–92.

Teng, J., Ouyang, L., & Chen. L. (2007). A comparison between two pricing and lot-sizing models with partial backlogging and deteriorated items. International Journal of Production Economics, 105, 190–203.

Tripathi, R. P. (2011).  Inventory model with time dependent demand rate under inflation when supplier credit linked to order quantity. International Journal of Business and Information Technology, 1(3),  174-183.

Wang, J., & Shu, Y. F. (2005). Fuzzy decision modelling for supply chain management. Fuzzy Sets and System, (150), 107–127.

Wu, O., & Cheng (2005). An inventory model for deteriorating items with exponential declining demand and partial backlogging. Yugoslav Journal of Operations Research, 15(2), 277–288. 

Zadeh, L. A. (1965). Fuzzy sets, Information and Control, (8), 338-353.

Zimmermann, H. J. (1985). Application of fuzzy set theory to mathematical programming, Information Science, (36), 29-58.                                                                        

 

Appendix-A

The maximum amount of demand backlogged per cycle is as follows:

 =  ) =  [log{1+ d(T)}]                                                             (A1)

So, the order quantity per cycle is given by

Q =  +  =  ] +  [log{1+ d(T)}]                                          

=  ] +  [d(T) ], using Maclaurin series,         

=  + b[T ], taking terms upto power two, as q<1.          

Q = + b[T ]                                                                                (A2)

The inventory holding cost per cycle is

 = dt

= dt

= dt

=  [   + ()( + )]

Expanding  into ascending powers of q, taking the terms up to 2nd power in q, and  neglecting all terms, as q<1. After simplification we have

 =                                                                                                          (A3)

The deterioration cost per cycle is

 =

=

= 

=                

=               

=   +    + +   ]   

=   ()                                                                                                                 (A4)

The total ordering cost per order for the cycle [0, T] is given by

    for                                 

The shortage cost per cycle is

= [ ]

=   dt

=  b[ ]                                                                         (A5)

The lost sale cost per cycle is

 = b dt = b[(T)             (A6)