versão On-line ISSN 1012-277X
S. Afr. J. Ind. Eng. vol.25 no.3 Pretoria Nov. 2014
K. Adendorff; P.S. Kruger*
Department of Industrial and Systems Engineering University of Pretoria, South Africa. firstname.lastname@example.org
The study revisits the subject matter of inventory control, a continual part of the activities of wide-ranging organisations internationally. The mathematical model is presented of a particular situation that deals with the regular acquisition of a material required for a production process in a volatile environment of varying demand and fluctuating price. The usual process dynamics are demonstrated against a background of diverse choices of probability density function. The model makes use of Normal and Weibull distributions.
Die studie herbesoek die bakermat van voorraadbeheer om hernude klem te lê op die ontwikkeling van gepaste wiskundige modelle vir uiteenlopende voorraadsituasies wat tans internasionaal met welslae gebruik kan word. 'n Rudimentêre model word voorgehou wat desnieteenstaande eenvoud dog onmisbare gereedskap soos Dinamiese Programmering en keuse van kansdigtheidsfunksie voorhou vir stogastiese behandeling van die onbestendige gedrag met verloop van tyd van voorraadprys- en aanvraag. Die model gebruik Normaal- en Weibullverdelings.
"Almost any realistic decision problem has two fundamental characteristics: it is sequential and it is uncertain. An excellent example of a sequential decision problem under uncertainty is the holding of inventories." Arrow 
The present-day field of inventory control has reached such a stage that practising industrial and systems engineers may easily be overwhelmed by the extensive population of vagaries that may be encountered in attempting to select, design, and successfully operate a reliable inventory control system. The menu of opportunities is fraught, inter alia, with advanced concepts such as Dynamic Programming  and Stochastic Systems  that may be employed to good effect in practice.
The numeric case study presented here deals with an organisation in the agricultural sector that purchases a raw material to be used as an ingredient of a weed-killing preparation. The raw material price fluctuates with the passage of time, as does demand for the raw material.
To illustrate the creation of an inventory model suited to the problem, a two-period policy is initially considered without loss of applicability for n-period policies .
2 THE MODEL
The stochastic and dynamic nature of a two-period model, as proposed by , relates to achieving assurance of an outcome by using the applicable observations and constraints. To this end, the formulation of the model for the expected cost of the n-period optimal policy as provided by  is shown in equation (1)
n is the period,
s is the quantity of inventory on hand at the beginning of the period,
P is the market price of the commodity,
S is the stock level at the beginning of the first period,
I(S) is the holding and shortage cost in period one,
(r) is the distribution of demand (units/period),
(p) is the distribution of the price (Rand/unit), and
r and p are integration veriables.
A brief description of the particular modelling environment includes the following features:
• The inventory level is unconstrained.
• Inventory is replenished instantaneously if required.
• Surplus inventory is not sold.
• Inventory holding and shortage costs are considered.
• The distributions of demand and price are independent.
• The final model is created in two sequential stages by proceeding via a single period policy to an eventual two-period policy.
2.1 The Normal demand distribution
The Normal distribution is initially used to create the mathematical functions that are required for the single period model stage shown in Table 1. The values used are as follows:
Mean: 36 units/period,
Standard deviation: 10.5 units per period,
Inventory holding cost = c1 =10.00 Rand per unit period, and
Inventory shortage cost = c2 = 130.00 Rand per unit.
The obtained values of the terms in equation (1) are shown in Table 1. These values are provided primarily for purposes of illustration. The integral of the Normal distribution is obtained by fitting a logistic function to the cumulative function of the Normal distribution using non-linear regression. The same result may be obtained by using the built-in relevant function provided by, for example, Microsoft Excel. The values for x have been chosen to reflect the mean plus or minus three standard deviations of the demand distribution.
To determine the market price (Rand per unit) of the raw material at the beginning of the period, equation (2) is applied :
2.2 The Weibull price distribution
The Weibull distribution  is selected to serve as a tool in establishing the purchasing policy described below. The reason for this choice is that it is versatile, and is often used in analysing industrial and systems engineering problems.
To establish the correct two-period purchasing policy p2(s), the distribution given in Table 3 below is used to create functions to model the second period model stage, as shown in Table 4. The symbols and values are as follows:
Mean: 250 Rand/unit
Shape parameter: c = 3 (chosen for the sake of convenience)
Scale parameter: b = (Mean)/ [(c+1)/c)]=250/0.8930=280
pn = future price per unit (Rand)
φ(p) = probability density function of future price per unit.
The determination of p2(s) typically employs equation (3) over the range of s-values from 4.50 to 67.5 :
The value of p2(s) as a function of s is shown in Figure 3. The curve depicts the essence of the model that is used to operate the inventory system optimally, as follows:
• The current unit price p2(s) of the raw material must be determined at the beginning of the period.
• The optimum inventory level may be read from the curve in Figure 3.
• Should the actual inventory level, x, be less than s, an order size of (s - x) shall be purchased. Alternatively, if the actual inventory level is greater than s, replenishment should not take place.
The elementary case study of this paper emphasises features of a quantitative nature for purposes of demonstration, and does not cater for subjectivity. Its tutorial-like flavour is meant to encourage practitioners and the uninitiated to explore the trail of information, which is presented here in truncated form.
To have confidence in the optimisation dictated by equation (2), the operation of the system has been simulated and compared with non-optimal operation. The simulation exercise indicates a 9% saving in the cost of holding inventory.
 Arrow, K.J. 1957. Statistics and economic policy. Econometrica 25, 523-531. [ Links ]
 Bellmann, R. 1957. Dynamic programming. Princeton University Press. [ Links ]
 Fabian, T., Fisher, J.L., Sasieni, M.W. & Yardeni, A. 1958. Purchasing raw material on a fluctuating market. Cleveland, Ohio: Case Institute of Technology. [ Links ]
 Lai, C.D., Murthy, D. N. P. & Xie, M. 2011. Weibull distributions. Wiley Interdisciplinary Reviews: Computational Statistics, 3(3), 282-287. [ Links ]
 Naddor, E. 1966, Inventory systems. John Wiley & Sons, Inc. [ Links ]
* Corresponding author