Journal of Model Based Research

Journal of Model Based Research

Journal of Model Based Research

Current Issue Volume No: 1 Issue No: 3

Research Article Open Access Available online freely Peer Reviewed Citation

Determining the Efficiency of Fuzzy Logic EOQ Inventory Model with Varying Demand in Comparison with Lagrangian and Kuhn-Tucker Method Through Sensitivity Analysis

1Department of Mathematics, Cauvery College for Women, Bharathidasan University, Trichy, Tamilnadu -620018, India.

2Department of Mathematics, KhadirMohideen College, Bharathidasan University, Adhiramapattinam, Tamilnadu- 614701, India.

3Department of Mathematics, Saveetha Engineering College, Chennai, Tamilnadu -602105, India

4Department of Physics, Loyola College, Chennai, Tamilnadu- 600034, India.


This paper considers an EOQ inventory model with varying demand and holding costs. It suggests minimizing the total cost in a fuzzy related environment. The optimal policy for the nonlinear problem is determined by both Lagrangian and Kuhn-tucker methods and compared with varying price-dependent coefficient. All the input parameters related to inventory are fuzzified by using trapezoidal numbers. In the end, a numerical example discussed with sensitivity analysis is done to justify the solution procedure. This paper primarily focuses on the aspect of Economic Order Quantity (EOQ) for variable demand using Lagrangian, Kuhn-Tucker and fuzzy logic analysis. Comparative analysis of there methods are evaluated in this paper and the results showed the efficiency of fuzzy logic over the conventional methods. Here in this research trapezoidal fuzzy numbers are incorporated to study the price dependent coefficients with variable demand and unit purchase cost over variable demand. The results are very close to the crisp output. Sensitivity analysis also done to validate the model.

Author Contributions
Received 28 Jun 2020; Accepted 10 Jul 2020; Published 01 Aug 2020;

Academic Editor: Babak Mohammadi, College of Hydrology and Water Resources, China.

Checked for plagiarism: Yes

Review by: Single-blind

Copyright ©  2020 Kalaiarasi K,et al.

Creative Commons License     This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Competing interests

The authors have declared that no competing interests exist.


K. Kalaiarasi, M. Sumathi, H. Mary Henrietta, A. Stanley Raj (2020) Determining the Efficiency of Fuzzy Logic EOQ Inventory Model with Varying Demand in Comparison with Lagrangian and Kuhn-Tucker Method Through Sensitivity Analysis. Journal of Model Based Research - 1(3):1-12.

Download as RIS, BibTeX, Text (Include abstract )

DOI 10.14302/issn.2643-2811.jmbr-20-3465


In today’s competitive scenario, organizations face immense challenges for meeting the transitional consumer's demand, and maintaining the inventory plays a major role. Indulging the betterment of various promotional activities that yield a reputation for the concern. The classical economic order quantity by Harris 1 and Wilson 2 was developed for problems where demand remains constant. In recent study of inexplicit demand rates by Chen 3 and few more resulting in the time-dependent demand rates and varying holding cost. In 1965, Zadeh 4 introduced the concept of fuzzy that carries vagueness or unclear in sense. Hence fuzzy sets came into prominence in describing the vagueness and uncertainty that impressed many researchers. In this paper, the total cost has been taken and the proposed model economic order quantity (EOQ) and the input parameters such as ordering cost, purchase cost, order size, holding cost, unit purchase cost, and unit selling price are fuzzified using trapezoidal numbers. The optimization is carried out by both Lagrangian and Kuhn-Tucker methods and graded mean integration for defuzzification of the total cost.

In the early 1990’s the Economic order quantity (EOQ) and Economic Production Quantity (EPQ) played a vital role in the operations management area, contrastingly it failed in meeting with the real-world challenges. Since these models assumed that the items received or produced are of a perfect quality which is highly challenging. In inventory management, EOQ plays a vital role in minimizing holding and ordering costs. Ford W. Harris (1915) developed the EOQ model but it was further extended and extensively applied by R.H. Wilson. Wang5 examined an EOQ model where a proportion of defective items were represented as fuzzy variables. In the year 2000, Salameh et al 6 came up with a classical EOQ that assumed a random proportion of defective items, and the recognized imperfect items are sold in a single batch at an economical rate. Chen 3 in the year 2003, came up with an EOQ with a random demand that minimized the total cost. Firms to cope up with the current scenarios, have to adapt the diversity among inventory models due to the advancement of management strategies and varying production levels.

Briefly, the input parameters of inventory models are often taken as crisp values due to variability in nature.H.J. Zimmerman 7 studied the fuzziness in operational research. In a classical EOQ, Park 8 fuzzified the decision variables ordering cost and holding cost into trapezoidal fuzzy numbers that gave rise to study the fuzzy EOQ models, solving non-linear programming to obtain the optimal solution for economic order quantity. In 2002, Hsieh 9 initiated two production inventory models that applied the graded-mean representation method for the defuzzification process. In traditional inventory models when minimizing the annual costs, the demand rate was always assumed to be independent. Despite the promotional setups and the deterioration items, a wavering demand arises. In 2003, Sujit 10 came up with an inventory model that involves a fuzzy demand rate fuzzy deterioration rate. In 2013, Dutta and Pavan Kumar 11 proposed an inventory model without shortfall with fuzziness in demand, holding cost, and ordering cost. The study of inexplicit demand rates was doneby Chan 12 and few more resulting in the time-dependent demand rates and varying holding cost. A detailed study of Lagrangian method to solve the non-linear programming of the total cost was done by Kalaiarasi et al13. The Kuhn-Tucker optimization technique was for the NPP(non-linear programming) problem was carried out by Kalaiarasi et al14.

Considering these inputs, the non-linear programming problem is solved for total cost function using Lagrangian and Kuhn-Tucker method and concluded with a sensitivity analysis, which exhibits the variations between the fuzzy and crisp values.

Considering these inputs, non-linear programming is solved for total cost function using Lagrangian and Kuhn-Tucker method and done a sensitivity analysis between the slight variations between the fuzzy and crisp values.


Definition 1:

A fuzzy set à defined on R (∞,-∞), if the membership function of à is defined by

Definition 2:

A trapezoidal fuzzy number Ã=(a, b, c, d) with a membership function μà is defined by

The Fuzzy Arithmetical Operations

Function principle 12 is proposed to be as the fuzzy arithmetical operations by trapezoidal fuzzy numbers. Defining some fundamental fuzzy arithmetical operations under function principle as follows

Suppose Ã=(x1, x2, x3, x4) and = (y1, y2, y3, y4) be two trapezoidal fuzzy numbers. Then

The addition of is

à ⊕ B̅ =(x1+y1, x2+y2, x3+y3, X4+y4)

Where x1, x2, x3, x4,y1, y2, y3, y4 are any real numbers.

(2) The multiplication of is

à ⊕ = (c1, c2, c3, c4)

Where Z1={x1y1, x1y4, x4y1, x4y4}, Z2 = ={x2y2, x2y3, x3y2, x3y4}, C1 = min Z1, C2 = min Z1, C1 = max Z2, C1 = max Z2.

If x1, x2, x3, x4,y1, y2, y3 and y4 are all zero positive real numbers then

à ⊗ = (x1y1, x2y2, x3y3, x4y4).

(3) The subtraction of is

Where also x1, x2, x3, x4,y1, y2, y3 and y4 are any real numbers.

(4) The division of is

where y1, y2, y3 and y4 are positive real numbers. Also x1, x2, x3, x4,y1, y2, y3 and y4are nonzero positive numbers.

(5) For any ∝ ∈ R

  Extension of the Lagrangian Method.

Solving a nonlinear programming problem by obtaining the optimum solution was discussed by Taha13 with equality constraints, also by solving those inequality constraints using Lagrangian method.

Suppose if the problem is given as

Minimize y = f(x)

Sub to gi(x) ≥ 0, i = 1, 2, · · ·, m.

The constraints are non-negative say x ≥ 0 if included in the m constraints. Then the procedure of the extension of the Lagrangian method will involve the following steps.

Step 1:

Solve the unconstrained problem

Min y = f(x)

If the resulting optimum satisfies all the constraints, then stop since all the constraints are inessential. Or else set K = 1 and move to step 2.

Step 2:

Activate any K constraints (i.e., convert them into equalities) and optimize f(x) subject to the K active constraints by the Lagrangian method. If the resulting solution is feasible with respect to the remaining constraints, the steps have to be repeated. If all sets of active constraints taken K at a time are considered without confront a feasible solution, go to step 3.

Step 3:

If K = m, stop; there’s no feasible solution.

Otherwise set K = K + 1 and go to step 2.

Graded Mean Integration Representation Method

Graded mean Integration Representation Method was introduced by Hsieh et al 14 based on the integral value of graded mean h-level of generalized fuzzy number for the defuzzification of generalized fuzzy number. First, a generalized fuzzy number is defined as follows: Ã-=(α1, α2, α3, α4,)LR. By graded mean integration are the inverses of L and R are L-1 and R-1 respectively. The graded mean h-level value of the generalized fuzzy number Ã-=(α1, α2, α3, α4,)LR is given by h/2[L-1(h)+R-1(h)]. Then the graded mean integration representation of P(Ã) with grade then

In this paper trapezoidal fuzzy numbers is used as fuzzy parameters for the production inventory model. Let Then the graded mean integration representation is given by the formula as




 Inventory Model For Eoq

The input parameters for the corresponding model are

K - Ordering cost

a - constant demand rate coefficient

b - price-dependent demand rate coefficient

P - selling price

Q - Order size

c - unit purchasing cost

g - constant holding cost coefficient

Let’s consider the total cost 15 and modify in a fuzzy environment

The total cost per cycle is given by


Partially differentiating w.r.t Q,



Equating the economic order quantity in crisp values obtained as

Inventory Model For Fuzzy Order Quantity

By using trapezoidal numbers, fuzzified input parameters are as follows

The optimal order quantity

Partially differentiating w.r.t ‘Q’ and equating to zero,

Hence the optimal economic order quantity for crisp values is derived,

Applying the graded mean representation

Now partially differentiating w.r.t Q1, Q2, Q3, Q4 and equating to zero,

The above derived results depict that Q1 > Q2 > Q3 > Q4 failing to satisfy the constraints 0 ≤ Q1 ≤ Q2 ≤ Q3 ≤ Q4 . So, converting the inequality constraint Q2 - Q1 ≥ 0 into equality constraint Q2 - Q1 = 0 Optimizing P(TC(Q,P) subject to Q2 - Q1 = 0 by Lagrangian method.



From equations (1) and (2) the results are,





Since Q3 > Q4 which does not satisfy the constraint 0 ≤ Q1 ≤ Q2 ≤ Q3 ≤ Q4. Now converting the inequality constraints Q2 - Q1 ≥ 0, Q3 - Q2 ≥ 0 into equality constraints Q2 - Q1 = 0 and Q3 - Q2 = 0. Optimizing, 



From equations (1’) , (2’) and (3’),

In the above-mentioned results Since Q1 > Q4 which does not satisfy the constraint 0 ≤ Q1 ≤ Q2 ≤ Q3 ≤ Q4. Converting the inequality constraints Q2 - Q1 ≥ 0, Q3 - Q2 ≥ 0 and Q4 - Q3 ≥ 0 into equality constraints Q2 - Q1 = 0, Q3 - Q2 = 0 and Q4 - Q3 = 0 . Optimizing,


After Partially differentiation (Appendix)

satisfies the required inequality constraints.

Kuhn-Tucker Optimization Method

By applying Kuhn Tucker conditions, the total cost is minimized by finding the solution of Q1, Q2, Q3, Q4 with Q1, Q2, Q3, Q4

The above conditions simplify to the following 𝞴1, 𝞴2, 𝞴3, 𝞴4.

It is known that, Q1 > 0 then in


𝞴1Q1 = 0 arrive at 𝞴4 = 0 In a similar fashion, if 𝞴1= 𝞴2 = 𝞴3 = 0, Q4 ≤ Q3 ≤ Q2 ≤ Q1 does not satisfy 0 ≤ Q1 ≤ Q2 ≤ Q3 ≤ Q4. Therefore the conclusion is Q2 = Q1, Q3 = Q2 and Q4 = Q3

i.e., Q1 = Q2 = Q3 = Q4 = Q*


Numerical Examples

Let us consider an integrated inventory system having the following statistics with crisp parameters having following values 15, the ordering cost K=520 units, the optimal selling price P=36.52 units, unit purchasing cost c=4.75, g=0.2, constant demand rate coefficient a=100 units, price-dependent demand rate coefficient b=1.5.

As mentioned earlier, the trapezoidal numbers

yielding the below results. Equation 1 was the optimal order quantity for crisp (equation 1) values and optimization to the EOQ is done by applying Lagrangian (equation 2) and Kuhn-Tucker (equation 3) methods under graded-mean defuzzification method. Figure 1 represents the crisp and fuzzified output varying on demand.

Figure 1.Variable demand for crisp and fuzzy comparison for EOQ
 Variable demand for crisp and fuzzy comparison for EOQ


The fuzzified output varies at a larger rate than crisp output while varying the demand. This shows that the optimization results can be obtained only by restricting the controlling parameters. Figure 2 shows the variations in unit purchase cost for crisp and fuzzy output. An analysis of the total cost in a fuzzy environment is studied herewith by applying the Lagrangian and Kuhn-Tucker methods for optimization and using trapezoidal numbers. The value does not show much variation in optimal solutions. While varying the price-dependent demand rate coefficient among both the methods, fuzzified values decrease in Lagrangian and increases in the Kuhn-Tucker method. In Figure 3, the curves of crisp and fuzzy EOQ remain closer. Table 2 shows the comparison of graded mean values between the two optimization methods in operations research 13 Lagrangian and Kuhn-Tucker methods. The values are varied to an increase and decrease of 25% and 40% which clearly shows that the crisp values of EOQ stay perfectly aligned between both the methods. Figure 4 collates the three parameters viz., price dependent coefficient, purchase cost, and varying demand which shows the variations in 3D model. Table 1 & Table 2.

Figure 2.Variable demand for unit purchase cost for EOQ
 Variable demand for unit purchase cost for EOQ

Figure 3.Variable demand for price dependent coefficient for EOQ
 Variable demand for price dependent coefficient for EOQ

Figure 4.Three-dimensional variation of price dependent coefficient, demand and purchase cost.
 Three-dimensional variation of price dependent coefficient, demand and purchase cost.

Table 1. Comparison of Lagrangian and Kuhn Tucker method for the variation done in price dependent coefficient
Price-dependentco-efficient Lagrangian method Kuhn-Tucker method
b=1 EOQ = 263.6169 EOQ = 263.6169
Fuzzy = 260.2202 Fuzzy = 321.2251
 b=1.5 EOQ = 222.4949 EOQ = 222.4949
Fuzzy = 207.7356 Fuzzy = 258.333
 b=2 EOQ = 171.7967 EOQ = 171.7967
Fuzzy = 147.1273 Fuzzy = 194.5951

Table 2. Sensitivity analysis for Lagrangian and Kuhn Tucker method
Parameters % change parameters Graded-mean values  CrispEOQ LagrangianMethod Kuhn-Tucker Method
 K(520) +40% 728(708,718,738,748) 263.2595 246.140 305.902
+25% 650(630,640,660,670) 248.7569 232.451 288.931
-25% 390(370,380,400,410) 192.6863 179.430 223.231
-40% 312(292,302,322,332) 172.3438 160.136 199.342
 P(36.52) +40% 51.128(31.128,41.128,61.128,71.128) 159.7376 106.502 128.4246
+25% 45.65(25.65,35.65,55.65,65.65) 185.7729 147.128 178.9012
-25% 27.39(7.39,17.39,37.39,47.39) 253.9615 236.628 289.3353
-40% 21.912(1.912,11.912,31.912,41.912) 271.093 257.475 315.0044
 a(100) +40% 140(120,130,150,160) 305.4398 301.807 374.9838
+25% 125(105,115,135,145) 277.2588 271.688 338.6781
-25% 75(55,65,85,95) 148.7803 127.422 168.3332
-40% 60(40,50,70,80) 75.5944 52.1485 49.2549
 b(1.5) +40% 2.1(1.6,1.8,2.2,3) 159.7376 134.7605 187.338
+25% 1.9(1.6,1.8,2,2.2) 185.7729 180.8325 227.4581
-25% 1.1(0.5,0.8,1.4,1.7) 253.9615 250.2417 310.0121
-40% 0.9(0.6,0.8,1,1.2) 271.0938 271.8465 334.0561
 c(4.75) +40% 6.7(6.4,6.5,6.9,7) 188.0425 169.2837 210.5158
+25% 5.9(5.5,5.7,6,6.5) 199.0055 179.1926 222.8382
-25% 3.6(3,3.2,3.8,4.6) 256.9150 223.8899 278.4223
-40% 2.8(2.2,2.7,3,3.2) 287.2397 258.3873 321.3222
 g(0.2) +40% 0.28(0.18,0.2,0.3,0.5) 188.0425 170.5264 212.0612
+25% 0.25(0.14,0.15,0.26,0.54) 199.0055 176.1179 219.0145
-25% 0.15(0.1,0.14,0.16,0.2) 256.9150 241.4086 300.208
-40% 0.12(0.1,0.11,0.13,0.14) 287.2397 271.2692 337.3417


  1. 1.Harris F. (1915) Operations and cost:AW Shaw Co.Chicago.
  1. 2.Wilson R. (1934) A scientific routine for stock control:. , Harvard Business Review 13, 116-128.
  1. 3.C Y. (2003) A probabilistic approach for traditional EOQ model: Journal of information optimization sciences v 24. 249-253.
  1. 4.L A Zadeh. (1965) Fuzzy sets: Information Control v 8. 338-353.
  1. 5.Wang X, Tang W, Zaho R. (2007) Random fuzzy EOQ model with imperfect quality items. Fuzzy Optimization and Decision-Making v 6. 139-153.
  1. 6.M K Salameh, M Y Jaber. (2000) Economic production quantity model for items with imperfect quality:. , International Journal of Production Economics 64, 59-64.
  1. 7.Hans J Zimmermann. (1996) Fuzzy set theory and its applications: 3rdEd.Dordrecht Kluwer, Academic Publishers v 19.
  1. 8.K S Park. (1987) Fuzzy Set Theoretical Interpretation of economic order quantity. , IEEE Trans. Systems Man. CybernetSMCv 17, 1082-10.
  1. 9.C H Hsieh. (2002) Optimization of fuzzy production inventory models: Information Sciences. v 146. 29-40.
  1. 10.Sujit De Kumar, P K Kundu, Goswami A. (2003) An Economic production quantity inventory model involving fuzzy demand rate and fuzzy deterioration rate:. , Journal of Applied Mathematics and Computing.12(1-2): 251-260.
  1. 11.Dutta D, Kumar Pavan. (2013) Optimal policy for an inventory model without shortages considering fuzziness in demand, holding cost and ordering cost. , International Journal of Advanced Innovation and Research 2(3), 320-325.
  1. 12.Chan.Wang.(1996) Backorder fuzzy inventory model under function principle: Information Science,95,pgs1-2,71-79.
  1. 13.M K Sumathi.Mary HenriettaH.(2019) Optimization of fuzzy inventory model for EOQ using Lagrangian method:. , Malaya Journal of Matematik.v 7-3.
  1. 14.M K Sumathi, Mary Henrietta H. (2019) Optimization of fuzzy inventory EOQ model using Kuhn-Tucker method:Journal of Advanced Research and Dynamical control systems. v 11,09 Special issue.
  1. 15.H A Taha. (1997) . , Operations Research:Prentice Hall,New Jersey,USA
  1. 16.S H Chen, C H Hsieh. (1999) Graded Mean Integration Representation of Generalized fuzzy numbers:. , Journal of the Chinese Fuzzy Systems Association 5, 1-7.
  1. 17.TotanGarai Dipankar.Chakraborty, Tapan Kumar Roy.(2019) Fully fuzzy inventory model with price-dependent demand and time varying holding cost under fuzzy decision variables:. , Journal of Intelligent and Fuzzy Systems v 36, 3725-3738.

Cited by (4)

  1. 1.Raj A. Stanley, Henrietta H. Mary, 2023, , , 647(), 21, 10.1007/978-3-031-27409-1_3
  1. 2.Stanley Raj A., Mary Henrietta H., Kalaiarasi K., Sumathi M., 2021, , , 1374(), 123, 10.1007/978-981-16-0708-0_11
  1. 3.Kalaiarasi K., Mary Henrietta H., Sumathi M., 2021, WITHDRAWN: Comparison of hexagonal fuzzy numbers with varied defuzzification methods in optimization of EOQ using GP technique, Materials Today: Proceedings, (), 10.1016/j.matpr.2020.12.209
  1. 4.Gribanova Ekaterina, 2022, Elaboration of an Algorithm for Solving Hierarchical Inverse Problems in Applied Economics, Mathematics, 10(15), 2779, 10.3390/math10152779