Optimal Control Approaches to the Aggregate Production Planning Problem

In the area of production planning and control, the aggregate production planning (APP) problem represents a great challenge for decision makers in production-inventory systems. Tradeoff between inventory-capacity is known as the APP problem. To address it, static and dynamic models have been proposed, which in general have several shortcomings. It is the premise of this paper that the main drawback of these proposals is, that they do not take into account the dynamic nature of the APP. For this reason, we propose the use of an Optimal Control (OC) formulation via the approach of energy-based and Hamiltonian-present value. The main contribution of this paper is the mathematical model which integrates a second order dynamical system coupled with a first order system, incorporating production rate, inventory level, and capacity as well with the associated cost by work force in the same formulation. Also, a novel result in relation with the Hamiltonian-present value in the OC formulation is that it reduces the inventory level compared with the pure energy based approach for APP. A set of simulations are provided which verifies the theoretical contribution of this work.


Introduction
Nowadays one of the most important challenges faced by business is the adjustment of firm resources in order to satisfy market requirements subjected to fluctuations over time, namely costs, prices, existences, demands, etc. [1].In the case of facing fluctuating, hard-to-predict demand, several companies around the world have the conflicting goals of (1) limiting the buildup in finished goods inventory; and (2) minimizing changes in the capacity level.Traditionally, the most economical solution to absorb fluctuations in demand is a mix of two alternatives [2,3]: Adjust the capacity level (hiring/firing labor force, working overtime/undertime, subcontracting, etc.); this is known as the chase alternative [4,5]: track the expected monthly sales and compute the corresponding capacity requirements.

‚
Use of inventories (excess of SKUs, backlog of orders, or lost sales); this is known as the level plan: maintain a steady production rate over the entire year, using finished goods (smoothing/anticipation) stocks to absorb ongoing differences between output and sales.
This problem-known as the aggregate planning problem-can be stated as follows [6][7][8][9]: ‚ determine simultaneously the production rate P, inventory level β, and capacity levels C, ‚ for meeting a fluctuating demand (a set of forecasts), ‚ at each period of a finite planning time horizon in days (t = 1,2,..., N), ‚ for a given set of production resources, ‚ involving one product or a family of similar items (with small differences so that considering the problem from an aggregate viewpoint is justified), ‚ while minimizing total relevant costs (i.e., payroll, hiring/layoffs, overtime/undertime, inventory/shortage, etc.), ‚ and subject to non-constant, time varying constraints.
Aggregate production planning (APP, or workforce planning, production and employment smoothing, capacity and production planning), deals with matching capacity (via adjustment of production load, inventory, and employment levels) to changing demand, over a finite planning horizon, in order to achieve long-run profitability.By converting monthly sales forecasts, inventory levels, labor inputs, and production rates-of a single entity with characteristics representative of an entire product group-to a convenient aggregate load/capacity format (such as standard hours), a production plan is generated [10]: this involves a tradeoff between penalties for carrying inventory and varying the capacity level incurs the minimum total marginal cost over a calendar year.Because of this, the aggregate planning process has an economic importance due to the decisions involved (regarding the capacity and inventory levels necessary to meet anticipated demand over the planning period), as they impact the company's performance, i.e., profit maximization.This in turn requires complete and accurate information about: machine capacity, labor utilization, levels (inventory, safety stock, manpower adjustment, subcontract, storage), time (regular/overtime), costs (production, inventory, overtime/idle time, subcontracts, shortage, lost sales, break down, backorder, hiring/firing/training).With this idea on mind, the next section presents a discussion of APP in production-inventory systems, while Section 3 introduces the research statement of APP as a control engineering problem in the context of optimal control theory.In Section 4, the stability analysis and mathematical modeling of APP is present.In order to show the validity and usefulness of the proposed formulation, results are presented in Section 5. Finally, Section 6 shows the conclusions derived from the case study plus future research venues.

Decision Science Approaches to the Aggregate Planning Problem
Models that have been used to solve the APP problem include: linear decision rule (for long-term strategic APP decisions), transportation model, dynamic programming, lot sizing model, linear programming (the most widely accepted method like in ( [11]), , heuristics [12], simulation [13], goal programming [14], micro spreadsheet analysis [15], multi-objective optimization [16], fuzzy formulation models [17,18], genetic algorithms (GAs) [19], and multiple criteria mixed integer programming [20].Based on Holt et al., in [2], in the area of models for planning/scheduling production-inventory systems, two types of models can be identified:

‚
Static models which are based on a finite planning horizon, considering deterministic demand.In this scenario, the receding planning horizon is bigger than the period by period plan.

‚
Dynamic models which consider a indefinite planning horizon, with a proper forecasting in demand; this based on the context that period by period decisions are based on a receding forecast over a planned time horizon.
where most of the existing aggregate planning models found in the literature [21]: ‚ Try to minimize an objective function representing "total relevant costs" (like production, inventory, shortage costs) over the fixed planning horizon.

‚
The usual constraints employed are inventory and capacity constraints, and are formulated to a single-objective function in linear programming [22].
A review of the literature over the past 20 years, in the area production planning modeling, shows the following: a taxonomy of the mathematical models used in aggregate planning; a literature review in aggregate, disaggregate, scheduling, and sequencing methods by feature, model, objectives, decision variables, and solution method; the achievements in optimal control methods-based on the maximum Pontryagin´s principle theory-that allowed analytical investigation of aggregate production planning systems to be carried out in order to gain an insight into their optimal behavior; and the literature in aggregate production planning models that have been introduced in the last four decades; the different planning models used in the production arena-from the supply chain level down to the manufacturing resource and capacity-as well the different approaches followed (namely conceptual, analytical, etc.) are summarized in [23].In recent years, it has become evident that many researchers and practitioners are increasingly aware of addressing:

‚
Real-life situations of management and decisions with the presence of multiple objectives; by developing goal programming optimization models.

‚
Variations over time; by developing multi-period optimization models.

‚
Novel optimization methods as convex optimization in [24].
Even though all of the APP models cited in the literature have shown good performance in the academic field that have not in practical life, according to [25] the following reasons that apply are:

‚
Current models do not capture the real world approaches in APP scenarios.

‚
The hypothesis in which items are homogeneous and are aggregated.

‚
The hypothesis that workforce has the same competence.

‚
In the mathematical modeling, the following areas are not considered, such as: human resources, marketing, and finance.

‚
The information provided by industry is not proper in the context of sales forecasting and cost information.

‚
Complex mathematical modeling and analysis of the APP process.
Furthermore, in [26] the following reasons apply:

‚
The hypothesis that uniform rates are assumed for different items, which do not capture the scenario for producing various items.

‚
There is an absence of interest from managers to adapt mathematical methods and proper techniques.

‚
The cost related to the collection and quantification of data in order to apply the above techniques.

‚
The real cost functions in organizations are not well established by actual techniques.
Each of the last mentioned issued can be grouped as follows: ‚ Regarding the modeling approach; based on [27], decision makers such as operations researchers and other mathematical decision builders and its application in real scenarios are far from business practice.

‚
Regarding the aggregation approach; even though it intends to (1) simplify/facilitate the myriad of calculations, and (2) to reduce the solving time of a large problem, it presents some issues [28]: aggregation does not result in a better optimal solution; it probably compromises optimality and may even result in an infeasible solution; it usually results in a better attained solution within a fixed solution time.

‚
Regarding the use of mathematical models; comments in journals indicate that because a worrying gap exists between theory and practice, as mathematical (optimizing) models have not had a significant impact on industry operations management practices managers find mathematical methods too daunting; the cost assumptions used in models are over-simplified and unrealistic; simplified/inflexible assumptions which limit their industrial applicability; cannot cope with calendar variations (i.e., holidays), and the revision of distant, and therefore, speculative forecasts cause instability in the schedules; broader concepts in the area of employment policy and inventory practices need to be introduced.
Summarizing: due to the highly complex constraints of the APP problem, exact optimal solutions provided by traditional optimization methods may very possibly be meaningless.It is true that a number of artificial intelligence approaches-combined with mathematical programming models-have been used to solve the APP problem, but little attention has been given to the consideration that marks the difference between a pure academic treatment of the APP problem and a result with real-life, practical implications: the simultaneous combination of many constraints affecting the quality of the APP.Based on this, in its simplest form the production planning paradigm describes firm operating in a market facing external demand which it tries to meet by utilizing a limited set of production resources that has limited ability to generate output in a given time period [29].In this context, the application of optimal control approaches presents a suitable mathematical tool with the intention to model and analyze the APP problem in a short time horizon, which is our case of analysis.Also, it is important to establish that optimal control is a dynamic optimization problem which can handle several decision variables and a control law in order to set up the conditions for a decision maker in the production planning context.

The Dynamic Nature of APP
The previous sections can be summarized as follows: an economical solution (for companies) to deal with a fluctuating demand is a mix between the chase alternative (hire/fire capacity) and level plan (inventory excess/shortages) strategies.This involves an inventory-capacity tradeoff, and is known as the aggregate planning problem: determine what must be the production rate, inventory and capacity levels necessary to minimize inventory holding, shortages, and production switching costs subject to time varying constraints.To address this problem, static and dynamic models-which try to minimize a linear programming objective function subject to inventory and capacity constraints-have been proposed.More recently, and in order to reflect more realistic situations, models with multiple objectives and variations over time, have been developed.As mentioned before, the proposed approaches have several shortcomings: they are over-simplified and unrealistic, attained solutions may be even impractical and infeasible, the use of speculative forecasts causes instability in the obtained schedules, etc.However, it is the premise of this paper that the main drawback of these proposals is, that they do not take into account the dynamic nature of the aggregate planning problem: the excessive inventory and production costs due to the use of speculative forecasts, is worsened by a poor understanding of the time lags between the ordering of goods and receipt into stock [30,31].In fact, a production system designed to deal with the APP problem should have as objectives:

‚
to buffer the production system from the customer; with a minimum reasonable inventory that absorbs the high frequency content in demand and allows having a level schedule.In this way, the variability of customer demand is reduced/avoided as switching production levels up and down frequently may be very expensive in practice [32].
‚ to buffer the customer from supply time lags; by selling goods straight off the shelf.In this way high customer service levels can be achieved, which can only be accomplished when the dynamic behavior of its constituent parts (i.e., materials/information flow, operations performed, resources/decision, rules/performance measures used) has been taken into account [33].

‚
Demand patterns; after a ramp increase there is a continuing freefall in inventory levels, after a step increase there is a permanent inventory deficit.So demand-without some form of averaging-results in excessive fluctuations in production rates (which are supposed to be absorbed by the inventory buffers).

‚
Lead times; the amount of inventory holding that is needed to satisfy a customer service level is dependent on the uncertainties in both demand and lead times.

‚
Inventory levels; trying to correct all the inventory discrepancy in a single time period, when in fact it may take many more time periods, provokes excessive (overshoots and undershoots around the target level).

Control Engineering & Production/Inventory Control
In [34] the author presents a series of reasons for using control engineering techniques in production/inventory control, i.e., the use of standard forms, the block diagram format, standard techniques that enable important performance metrics to be calculated without recourse to simulation, there are a number of techniques for transferring problems from one domain into another, etc. Authors like [35][36][37][38][39] present in-detail reviews of control engineering applications to production/inventory control.One of these applications that deserves to be mentioned apart-due to the number of research studies based on it-is the APIOBPCS concept.APIOBPCS stands for Automatic Pipeline, Inventory and Order Based Production Control System, and it is a well-established (both industrially and theoretically) production scheduling/control system model which operates on a knowledge of customer demand, inventory level, and unfilled orders, and that it provides an acceptable trade-off between production smoothing and a high level of stock turnover.The ordering policy/production algorithm of APIOBPCS is representative of work in [40].Most of the APIOBPCS research has been undertaken using both control theory mathematics and system dynamics (SD), the most representative works in SD are present in [41] where it is applied to APP via transfer function approach; in [42] applying SD for production-inventory systems; in [43] a production-inventory system for APP described by differential equations via control oriented approaches; and in [44] where the SD approach is extended in production-inventory in remanufacturing.Based on this, APIOBPCS models are usually expressed in a continuous control form, but there is a discrete version available as well.In any case, the APIOBPCS simulation models may be used confidently as a benchmark to demonstrate performance enhancement for a wide range of practical scenarios [45].

Order-Up-To (OUT) and Smoothing Policies
As inventories should have a stabilizing effect on material flow patterns, a minimum reasonable inventory is necessary to absorb the variations in demand and allow a level schedule.The level scheduling problem (LSP), or production smoothing problem (PSP), refers to the problem of finding level schedules where: ‚ production of a given product is constant over time or, ‚ the cumulative production amount of a product is proportional to time or, ‚ the items should be dispersed over the schedule as uniformly as possible, with a minimal total deviation from the final level of operations.
Setting a takt-paced production results in a leveled production/production smoothing (key to establish the strategic/market pull), which is a simple matter of buffering-with either a time backlog or inventory-the production line from demand variability.The tradeoff between stocks and schedule stability is reflected in the master production schedule.On the other hand, a replenishment strategy that strives to bring the inventory position up to a predetermined target level is called an OUT Sustainability 2015, 7, 16324-16339 policy.These kind of policies are very popular both in research as in practice since they are known to minimize inventory holding and shortage costs.Dejonckheere, et al. [39] presents two shortcomings of OUT policies are: 1.
Generation of a bullwhip effect.2.
If production is not flexible and costs are excessive, a not optimal option is present.
Regarding point 1; the bullwhip effect has largely been analyzed by the OR, system dynamics, and control theory communities, where two popular approaches are: the statistical inventory control approach and the control engineering approach.Most of the research developed analyzing bullwhip effect considered supply chain systems.However, from the point of view of production-inventory systems the bullwhip effect has been studied in [46] where it applies autoregressive models to achieve multiple steps demand; in [47] an adaptive base-stock policy to determine order quantities; and in [48] explores the inventory stability in the context of a seasonal supply chain.In the case of the first approach: steady-state models based on steady-state conditions-like deterministic, stochastic, economic game-theoretic, and simulation models-are insufficient and therefore, unable to describe, analyze, and find remedies for problems like the bullwhip effect.In the case of the second approach: the bullwhip effect can be avoided by smoothing the ordering pattern-a problem known in the literature as the "production smoothing problem"-which means that it is possible to dampen order fluctuations even in environments where decision makers have to rely on forecasts.
Regarding point 2; the total costs of a perfectly controlled system-defined as a system that faithfully tracks some reference or target signal [49], is composed by costs associated with perfectly tracking the target (i.e., traditional fixed and variable costs), and costs associated with not being in perfect control of a system (i.e., under produce/over produce costs, and excessive/insufficient inventory costs).As there is trade-off to be made between OUT policies' minimum inventory holding plus shortage costs, and smoothing policies' minimum production switching costs, when the cost structure is altered, there is a need to identify a set of values (related to the production rate, and inventory/capacity levels) that reduce the sum of total costs [50].

The Aggregate Planning Problem as an Optimal Control Problem
Fact #1: the aggregated planning problem has to do with finding an optimal set of parameters' values that minimizes a set of costs, facing a varying demand and within a finite planning time horizon; fact #2: the control engineering approach allows to analyze, design and simulate dynamic models; fact #3: in order to dampen order fluctuations and avoid the bullwhip effect (generated by OUT policies), a smoothing policy is required.When viewed together, these facts suggest (1) the understanding of the aggregated planning problem as an optimization problem; and (2) the use of control engineering-based tool capable of dealing with dynamic systems that allows the characterization of damping strategies.With this idea on mind, we state the research proposal of this paper as the formulation of APP as an optimization problem, applying OC techniques via an energy-based formulation for the dynamical system which describes the behavior and nature of the problem.

Mathematical Modeling of APP
Based on a previous work [51], and [52] which applies second-order differential equation analogizing production system with mechanical vibration systems, here the approach is to propose an energy based analogy to obtain the dynamic equation for APP in production-inventory systems.In order to develop a consistent energy based analogy for APP, the following energy based function for APP, where E APP is the total energy of the APP system which is proposed: , and the potential energy is calculated as: V " 1 2 WP 2 , where in our analysis W pPq " WP.The Lagrangian is: Calculating the Lagrangian equation is: From where, for stability purposes, the damping is of the form: φ " ´pβ ptq ´αq ‚ P. Finally, for purposes of APP, the equation of the system is the following model: In order to achieve an understanding about the α parameter (EOQ) in the formulation, please refer to the Appendix.The Damping factor in Equation ( 3) has the characteristic to be a function over time, Lipschitz and it requires convergence in a scalar value.Regarding to a Production-Inventory system, the inventory level βptq presents an ODE form, as presented in [53], in this paper the approach is of the form:

Optimal Control Basic Concepts and Notation
Optimal control theory has as its objective the maximization of the return from, or the minimization of the cost of, the operation of physical, social, and economic processes [54].Based on this, in this work the interest is to apply Optimal Control approaches to problems in APP with dynamics in production-inventory level taking into account a contribution of the EOQ.

Optimal Control Formulation for APP: Continuous Inventory Policy
An optimal control is defined as an admissible control which minimizes an objective function.Given a dynamic system with initial condition x 0 , and which evolves in time according to ‚ x " f px, u, tq, the objective is to find control vector which is admissible and achieves a minimum for the cost functional.The optimal control problem is: x " f px, u, tq , x pt 0 q " x 0 , u ptq P m (5b) In economics, several problems in OC present a discount factor e ´δt .In this work, this discount factor is applied to production-inventory systems, as the ratio between the EOQ and Capacity.Based on this, the optimal control problem is of the form: From where the Hamiltonian is: Defining the Hamiltonian "Present value" by: Considering the relation: m ptq " λ ptq e δt , after some manipulation, a detailed discussion is on [55], the Hamiltonian "Present value" is: The interest is to apply both Optimal Control formulations to the dynamic system from Equations ( 3) and ( 4), which results are provided in the following section.

Theorem 1. (LaSalle's Invariance Principle)
Let f pxq be a locally Lipschitz function defined over a domain D Ă n and Ω Ă D be a compact set that is positively invariant with respect to ‚ x " f pxq.Let V pxq be a continuously differentiable function defined over D such that ‚ V pxq ď 0 in Ω.Let E be the set of all points in Ω where ‚ V pxq " 0, and M be the largest invariant set in E. Then every solution starting in Ω approaches M as t Ñ 8 .
In order to apply Theorem 1, for stability purposes, in the dynamical system conformed by Equations ( 3) and ( 4), the input (demand function) is of the form: d ptq " 0. From which, the dynamical system has the form: Considering that: x 1 " P, x 2 " ‚ P and x 3 " β, with an equilibrium point in px 1 , x 2 , x 3 q " p0, 0, 0q from Equations ( 10)- (12).To probe stability, a Lyapunov candidate function is proposed such as: Corollary 2. ( [56]) Let x = 0 be an equilibrium point for ‚ x " f pxq.Let V : D Ñ R be a continuously differentiable positive definite function on a domain D containing the origin x = 0, such that ‚ V pxq ď 0 in D. Let S "

"
x P D| ‚ V pxq " 0 * and suppose that no solution can stay identically in S, other than the trivial solution.Then, the origin is asymptotically stable.Applying the condition ‚ V ď 0 to Equation ( 13): Sustainability 2015, 7, 16324-16339 After some algebraic manipulation Equation ( 14) is of the form: Substituting Equations ( 10)- (12) in Equation ( 15): Grouping terms and factorizing in Equation( 16): The following condition must be satisfied in Equation ( 17) for stability purposes: α " K 2 .From which we have: Considering that x 1 , x 2 , x 3 ą 0 and From Equation (19), in order to satisfy ‚ V ď 0, a necessary and sufficient condition is to achieve:

Energy Based Optimal Control Formulation: Continuous Inventory Policy Approach
Consider the following energy based-Optimal Control problem: Calculating the Hamiltonian from: Applying first condition from Pontryagin Maximum Principle: We achieve the following set of ODE, which are the co-states: Sustainability 2015, 7, 16324-16339 Applying the second condition from Pontryagin Maximum Principle: and B 2 H Bu 2 " 1 ą 0 thus we have a minimum.Finally, substituting u " λ 3 ´λ2 C in the set of states and co-states the following simulations are achieved.
In order to present a discussion over the previous results, Figure 1 establishes the production rate behavior over time which has the characteristic of a maximum and minimum level over the horizon presented.Figure 2 presents a maximum level in the same time horizon this is based on the formulation of the first order differential equation which characterizes it.Figure 3 relates that the demand level increase conforms with the time horizon increases.

10
( ) Applying the second condition from Pontryagin Maximum Principle: In order to present a discussion over the previous results, Figure 1 establishes the production rate behavior over time which has the characteristic of a maximum and minimum level over the horizon presented.Figure 2 presents a maximum level in the same time horizon this is based on the formulation of the first order differential equation which characterizes it.Figure 3 relates that the demand level increase conforms with the time horizon increases.

Energy Based-Optimal Control with a Discount Factor: A Hamiltonian "Present Value"
Based on a previous work, in [57] which applies a discount factor in the performance index, this approach considers the following problem with a discount factor ( ) The first condition for the Hamiltonian-Present value, which are the co-states:

Energy Based-Optimal Control with a Discount Factor: A Hamiltonian "Present Value"
Based on a previous work, in [57] which applies a discount factor in the performance index, this approach considers the following problem with a discount factor The first condition for the Hamiltonian-Present value, which are the co-states: where the Hamiltonian as the form: The following set of ODE, are for the co-states: Applying the second condition from Pontryagin Maximum Principle: and B 2 H Bu 2 " 1 ą 0 thus we have a minimum.Substituting u " m 3 ´m2 C in the set of states and co-states, after some manipulation, the following simulations are achieved.
The introduction of the Hamiltonian present value produces a lower a higher production rate level in the dynamics, based on Figure 4 and a lower inventory level (almost half of the energy based inventory level from Section 5.1) which is present in Figure 5. Finally, the demand level increases conform with time horizon increases as is shown in Figure 6.
Applying the second condition from Pontryagin Maximum Principle: The introduction of the Hamiltonian present value produces a lower a higher production rate level in the dynamics, based on Figure 4 and a lower inventory level (almost half of the energy based inventory level from Section 5.1) which is present in Figure 5. Finally, the demand level increases conform with time horizon increases as is shown in Figure 6.

Conclusions
This research work presents an Optimal Control formulation for APP problems via the approach of Energy-Based and Hamiltonian-Present value.Stability analysis via LaSalle invariance principle establishes a condition in the inventory level which considers the EOQ parameter.The main contribution of this paper is the mathematical model which integrates a second order dynamical system coupled with a first order dynamical system which incorporates production rate, inventory level, and capacity as well with the associated cost of the work force in the same formulation.
Simulations show conforming with the increased production rate, inventory level achieves a maximum level when that demand level grows.A novel result in relation with the Hamiltonian-present value in the Optimal Control formulation is that it reduces the inventory level compared with the

Conclusions
This research work presents an Optimal Control formulation for APP problems via the approach of Energy-Based and Hamiltonian-Present value.Stability analysis via LaSalle invariance principle establishes a condition in the inventory level which considers the EOQ parameter.The main contribution of this paper is the mathematical model which integrates a second order dynamical system coupled with a first order dynamical system which incorporates production rate, inventory level, and capacity as well with the associated cost of the work force in the same formulation.
Simulations show conforming with the increased production rate, inventory level achieves a maximum level when that demand level grows.A novel result in relation with the Hamiltonian-present value in the Optimal Control formulation is that it reduces the inventory level compared with the pure energy based approach for APP.Simulations show conforming with the increased production rate, inventory level achieves a maximum level when that demand level grows.
Further work presents the idea of integrating the associated cost with dynamics, and to extend the case studies in the cost functional.Also the interest is to extend the mathematical model, in discrete formulation, to apply Model Predictive Control (MPC).Future research will consider the use of robust optimal control, another way of dealing with uncertainty, where a deterministic uncertain-but-bounded quantity is used (i.e., future demand can be bounded between lower and upper limits, without needing to define the probability of occurrence of each possible event within these limits), and the constraints regarding the operation of the system.Our intention is to extend this work with the application of MPC strategies via a suitable dynamic system which integrates dynamics in the inventory level.A practical real life problem which addresses this approach is our interest as well.
in the set of states and co-states the following simulations are achieved.

Figure 1 .
Figure 1.Production Rate for APP (time scale in days).

Figure 2 .
Figure 2. Inventory level for APP (time scale in days).

Figure 2 .
Figure 2. Inventory level for APP (time scale in days).

Figure 3 .
Figure 3. Demand level for APP (time scale in days).

Figure 3 .
Figure 3. Demand level for APP (time scale in days).
of states and co-states, after some manipulation, the following simulations are achieved.

Figure 4 .
Figure 4. Production rate for APP (Hamiltonian-present value, time scale in days).

Figure 5 .
Figure 5. Inventory level for APP (Hamiltonian-present value, time scale in days).

Figure 5 .
Figure 5. Inventory level for APP (Hamiltonian-present value, time scale in days).

Figure 5 .
Figure 5. Inventory level for APP (Hamiltonian-present value, time scale in days).

Figure 6 .
Figure 6.Demand level for APP (Hamiltonian-present value, time scale in days).

Figure 6 .
Figure 6.Demand level for APP (Hamiltonian-present value, time scale in days).