Industrial optimization using three-factor Cobb-Douglas production function of non-linear programming with application

Shakoor Muhammad; Fazal Hanan; Sayyar Ali Shah; Aihua Yuan; Wahab Khan; Hua Sun; Shakoor Muhammad; Fazal Hanan; Sayyar Ali Shah; Aihua Yuan; Wahab Khan; Hua Sun

doi:10.3934/math.20231532

AIMS Mathematics

2023, Volume 8, Issue 12: 29956-29974. doi: 10.3934/math.20231532

Previous Article Next Article

Research article

Industrial optimization using three-factor Cobb-Douglas production function of non-linear programming with application

1.
School of Environmental and Chemical Engineering Jiangsu University of Science and Technology Zhenjiang 212003, China
2.
Department of Mathematics, Abdul Wali Khan University Mardan, Khyber Pakhtunkhwa, Pakistan
3.
Department of Computer Science, University of Science and Technology, Bannu, KP, Pakistan

Received: 22 July 2023 Revised: 09 October 2023 Accepted: 16 October 2023 Published: 06 November 2023
MSC : 90C30, 90C90

This paper is about the effectiveness of the Cobb-Douglas (C-D) production function in industrial optimization, estimating the number of factors used in the production process of the water industry, for instance, capital and human labor. Moreover, we have modeled a nonlinear optimization problem for a local water industry using two and three factors of production. For this purpose, we have taken into account the Cobb-Douglas production function with different production factors using the Lagrange multiplier method with the ordinary least squares method. In the course of the solution, a linear function is used to calculate the cost function, and the C-D production function is used to calculate the production function. The Lagrange multiplier method with the ordinary least squares method is then used to solve the constrained optimization problem for the product of production. Furthermore, we compared the outcomes from both examples of two- and three-factor C-D production functions in order to validate the Lagrange multiplier method for the C-D production function. Moreover, the three-factor C-D production function is solved by the Lagrange multiplier method with the ordinary least squares method, which provides optimal results as compared to previous studies in literature. The validity of the proposed methodology is explained by using the products of a local production industry in Pakistan.

Keywords:

Citation: Shakoor Muhammad, Fazal Hanan, Sayyar Ali Shah, Aihua Yuan, Wahab Khan, Hua Sun. Industrial optimization using three-factor Cobb-Douglas production function of non-linear programming with application[J]. AIMS Mathematics, 2023, 8(12): 29956-29974. doi: 10.3934/math.20231532

Related Papers:

[1]	Maolin Cheng, Bin Liu . A novel method for calculating the contribution rates of economic growth factors. AIMS Mathematics, 2023, 8(8): 18339-18353. doi: 10.3934/math.2023932
[2]	Chengjin Tang, Jiahao Guo, Yinghui Dong . Optimal investment based on performance measure with a stochastic benchmark. AIMS Mathematics, 2025, 10(2): 2750-2770. doi: 10.3934/math.2025129
[3]	Regina S. Burachik, Bethany I. Caldwell, C. Yalçın Kaya . Douglas–Rachford algorithm for control- and state-constrained optimal control problems. AIMS Mathematics, 2024, 9(6): 13874-13893. doi: 10.3934/math.2024675
[4]	Habibe Sadeghi, Fatemeh Moslemi . A multiple objective programming approach to linear bilevel multi-follower programming. AIMS Mathematics, 2019, 4(3): 763-778. doi: 10.3934/math.2019.3.763
[5]	Hsien-Chung Wu . Numerical method for solving the continuous-time linear programming problems with time-dependent matrices and piecewise continuous functions. AIMS Mathematics, 2020, 5(6): 5572-5627. doi: 10.3934/math.2020358
[6]	Abdullah Ali H. Ahmadini, Firoz Ahmad . Solving intuitionistic fuzzy multiobjective linear programming problem under neutrosophic environment. AIMS Mathematics, 2021, 6(5): 4556-4580. doi: 10.3934/math.2021269
[7]	M. E. Elbrolosy . Semi- $(E, F)$ -convexity in complex programming problems. AIMS Mathematics, 2022, 7(6): 11119-11131. doi: 10.3934/math.2022621
[8]	Nodari Vakhania . On preemptive scheduling on unrelated machines using linear programming. AIMS Mathematics, 2023, 8(3): 7061-7082. doi: 10.3934/math.2023356
[9]	Jiaxin Chen, Yuxuan Wu, Shuai Huang, Pei Wang . Multi-objective optimization for AGV energy efficient scheduling problem with customer satisfaction. AIMS Mathematics, 2023, 8(9): 20097-20124. doi: 10.3934/math.20231024
[10]	M Junaid Basha, S Nandhini . A fuzzy based solution to multiple objective LPP. AIMS Mathematics, 2023, 8(4): 7714-7730. doi: 10.3934/math.2023387

Abstract

1. Introduction

Optimization is a useful technique for determining the optimal solution to a problem. In other words, optimization is the problem of choosing suitable inputs under given circumstances in order to get the best possible output. For instance, optimization can be used in production models to adjust different inputs and make them more effective in order to get the best output for the production of a particular industry ^[1]. In this scenario, once we have modeled a problem, it can be solved using the available optimization techniques to find the optimal solution ^[2].

An optimization problem usually consists of three ingredients: an objective function, a set of constraints and a number of decision variables. There are two types of optimization problems: constrained optimization problems and unconstrained optimization problems. Constrained optimization problems have restriction(s) on the objective function, while unconstrained optimization problems have no restriction on the objective function ^[3,4].

If at least one of the objective functions or the constraint function is nonlinear, then the problem is known as a nonlinear optimization problem. There are many techniques that have been used for the maximization or minimization of nonlinear optimization problems. Sometimes, a problem cannot be modeled correctly using linear programming; therefore, one can use nonlinear programming approaches ^[5,6] to model the problem. For the constrained optimization problem under consideration in this paper, the objective function, known as the cost function, is linear while the constraint function, known as the Cobb-Douglas (C-D) production function, is nonlinear. In the course of solution, the constrained optimization problem is converted into an unconstrained optimization problem and then solved by the Lagrange multiplier method using the ordinary least squares approach.

There have been many problems in literature in which the C-D production function has been used. A mixed-integer linear programming (MILP) model was established for the optimization of production scheduling ^[7]. A two-factor C-D production function is carried out with the aim of picking the suitable C-D production model for calculating the production process of the selected manufacturing industries in Bangladesh ^[8]. By combining the allometric scaling concept, which is used to estimate the parameters of the C-D function, with the application in transportation problems in China, a novel algorithm for creating geographical C-D models is developed ^[9].

For the improvement of different factors in the Polish metallurgical industry, the power regression C-D function was used with the aim of developing a number of production factors, for instance, net production, production sold and volume of steel production ^[10]. The proper management of a country's resources is an important issue for its economic development. Along similar lines, the optimization of water management for three industries that rely on water demand prediction, subject to a number of ambiguities, is handled by the use of the C-D production function ^[11]. A two-factor production function was used to model sustainable economic development with the goal of labor production in relation to the commodity production system's capital-labor ratio ^[12]. An economic model has been presented by the application of proper Inada conditions to the C-D production function, which converges to or diverges from per capita product and the steady state of capital ^[13]. Moreover, a two-factor C-D production function was presented in which the effects of labor force and capital on agricultural heritage systems are carried out to maximize profit as well as the sustainability of agricultural heritage through the application of a two-factor C-D production function in which they examine the impacts of major factors on agricultural productivity ^[14]. An algorithmic or analytical procedure was conducted to handle the issue of optimal utilization of resources towards a feasible and profitable model via the C-D production function ^[15]. An application of the C-D production function model was used to find the role of land in urban economic growth ^[16]. A general oligopolistic market equilibrium with nonlinear programming was considered, in which each firm's factor contributes to the system, and then solved by tensor variational inequality ^[17,18].

It becomes necessary in a real-world optimization problem to adopt a set of nonlinear terms in a mathematical model in order to get particular operational features of the decision-making problem. On the other hand, when nonlinear terms exist in the course of a solution, they add to the computational complexity of the problem. For this purpose, the researchers have developed proper transformation as well as linearization methods for the optimization problems that consist of nonlinear terms.

In this paper, a special type of production function, which was founded by Cobb and Douglas in 1928 and is known as the C-D production function, is under consideration. This function is based on empirical studies that have been applied to the economy for optimal production ^[19]. This function supplies a number of different inputs to the problem and, as a result, produces a unique output for the problem. These inputs may be two or more than two in number, depending on the factors of production used in the industry. Also, the function having more than two factors is known as the generalized C-D production function ^[20,21]. There have been many results in literature in which two-factor C-D function was presented for production. This work presents the optimal solution technique of the C-D function for three factors of production using the Lagrange multiplier method and the ordinary least squares method, with applications. Moreover, Nervole's C-D production function with three inputs using the ordinary least squares method has been presented in ^[22]. Furthermore, this work also developed an environment to transform and linearize an optimization problem with nonlinear objectives or nonlinear constraints by using existing techniques ^[23].

The proposed work is summarized as follows:

(1) To develop a model in order to transform as well as linearize the nonlinear terms.

(2) To choose a suitable model ^[1] as an optimization problem and to solve it using a novel methodology known as C-D production function using the Lagrange multiplier method with the ordinary least squares method.

(3) To find whether the best possible solution to the problem is obtained or not ^[2].

(4) To apply the C-D production function using the Lagrange multiplier method with the ordinary least squares method for the solution of the two-factor C-D production function as well as for the three-factor C-D production function.

(5) To apply the case study to industrial optimization in order to minimize costs more efficiently.

First we solved the two-factor C-D production function ^[24] using the Lagrange multiplier method with the ordinary least squares method. Second we used three factors of the C-D production function with the same technique. The obtained results show that cost minimization in cases of three factors of production is more efficient as compared to cost minimization in cases of two factors of production.

2. Cost minimization

The problem under consideration depends on three basic factors:

● Cost of the company for manufacturing production.

● Quantity of production.

● Income from the sales of production according to market prices.

For instance, the role of the water industry is to filter water in different purifying tanks, then add chemicals, pack them properly, and then sell them in the market. In this task, the main objective of the industry is to minimize costs or maximize production. There have been many approaches in literature in which the total cost of an industry is presented as a linear function ^[24,25].

$\begin{equation} K = \sum\limits_{i = 1}^{n}q_i x_i+R_s = q^T x+R_s \end{equation}$

(2.1)

where n represents production factors.

$x = \left[{\begin{array}{cc} x_1 \\ x_2 \\ . \\ . \\ . \\ x_n \end{array}}\right]$

is a vector of production factors.

$q = \left[{\begin{array}{cc} q_1 \\ q_2 \\ . \\ . \\ . \\ q_n \end{array}}\right]$

is the vector of prices of production factors.

The following relationships present the production volume ^[26]:

$\begin{equation} y = p(x) = p(x_1,x_2...x_n). \end{equation}$

(2.2)

The income of the company from production sales is given by:

$\begin{equation} S(y) = b^T y \end{equation}$

(2.3)

where,

$y = \left[{\begin{array}{cc} y_1 \\ y_2 \\ . \\ . \\ . \\ y_n \end{array}}\right]$

is a vector of the quantity of produced goods.

$b = \left[{\begin{array}{cc} b_1 \\ b_2 \\ . \\ . \\ . \\ b_n \end{array}}\right]$

is the vector prices of produced goods.

The profit of the company is the difference between its income and its cost of production.

$\begin{equation} Z(x,y) = S(y)-K(x). \end{equation}$

(2.4)

The main focus of this paper is to minimize the cost at a specific level of production and, consequently, to maximize production with more factors of production.

The optimization problem originating from our data analysis consists of a linear function and a nonlinear function. We have used the linear function as a cost function, while the nonlinear C-D production function is used as production output. Let us consider a linear cost function as an objective function and the nonlinear C-D production function as a constrained function. Then, this becomes a constrained optimization problem with an equality constraint. For the sake of simplicity, we convert this constrained optimization problem to an unconstrained optimization problem by using the Lagrange multiplier method, and then we solve this constrained optimization problem for stationary points. Let us start with a function $Z$ of two independent variables that is subject to an equality constraint function $g$ . The objective and constraint function are given as follows:

$\begin{equation*} min\,\,z(x, y). \end{equation*}$

Subject to:

$\begin{equation*} g(x, y) = 0. \end{equation*}$

We have, using the Lagrange multiplier,

$\begin{equation} L(x, y, \lambda) = z(x, y) + \lambda g(x, y). \end{equation}$

(2.5)

In order to solve the problem, we must first determine the values of $x$ , $y$ , and $\lambda$ . Note that we consider the cost function to be an objective function and the C-D production function as a constraint function, and solve it for different factors of production. Furthermore, we extend this work as an application to optimize production factors in a specific industry.

The production function for three inputs is as follows:

$\begin{equation} P = A X^{a_1} Y^{a_2} Z^{a_3}. \end{equation}$

(2.6)

After linearizing, we get

$\begin{equation} \ln P = \ln A +a_1 \ln X+ a_2 \ln Y+ a_3 \ln Z. \end{equation}$

(2.7)

After this, we will use least squares linear regression to find out the structural parameters, which are $A, a_1, a_2, a_3$ . The corresponding cost functions is given by

$\begin{equation} C(X,Y,Z) = \alpha_1+ \alpha_2 X+\alpha_3 Y+\alpha_4 Z. \end{equation}$

(2.8)

Along similar lines as Eq (2.5), we can get the corresponding Lagrange function for three-factor C-D production function.

In other words, the overall methodology can also be summarized as follows:

(1) First, model an objective and constraint function using two inputs for the given industry.

(2) Transform a constrained optimization problem into an unconstrained optimization problem.

(3) Solve the unconstrained problem using the Lagrange multiplier method with ordinary least squares.

(4) Repeat the above procedure for three factors of production for the given industry.

(5) Compare the outcomes in both cases.

3. Production function and its structural parameters

Given input prices, a cost function shows how much it costs to produce various output levels. In the course of solving such problems, one or both productions and factors of production may be stated by using their values. It is good practice to present the products of an industry in proper units that have a number of production components. In a similar way, human labor can be calculated when needed. When the availability of aggregated data is smaller, it can be measured with headcount or work time. When a higher level of aggregation is available, then the value of human labor seems more suitable. The most challenging task is the quantitative description of the capital used in the industry. In the majority of analyses, it becomes challenging due to the use of a number of factors of production. The greater a company's asset base, the less it is associated with high productivity ^[27]. In the following, we analyzed a local water industry with a two-factor C-D production function and a three-factor C-D production function using the Lagrange multiplier method with the ordinary least squares method and compared their results.

3.1. Cost minimization in the water industry using the two-factor C-D production function

For this problem, the data is taken from a local water industry Abysin Water Industry, which is one of the local registered branches of Chemtronics Water Services in Lahore, Pakistan. First, the problem is solved for two factors, and then it is generalized. In this case, the cost function consists of two factors (labor and capital) of production, which is given by:

$\begin{equation} C(X,Y) = \alpha_1+\alpha_2 X_ +\alpha_3 Y \end{equation}$

(3.1)

where,

C(X, Y) represents the cost of the industry in rupees,

$\alpha_1$ represents the fixed cost of the industry,

$\alpha_2$ represents the unit price of labor per hour,

$\alpha_3$ represents the unit price of capital per kg,

$X$ represents the number of labor hours,

$Y$ represents the amount of capital in kg.

All the prices in the industry are in Pakistani currency (the rupee). The raw materials used as capital are taken in kilograms, i.e., the unit for capital is kg.

Here,

$\begin{equation*} \alpha_1 = 75000,\; \; \alpha_2 = 1400,\; \; \alpha_3 = 100. \end{equation*}$

Putting these values in Eq (3.1), we get,

$\begin{equation} C(X,Y ) = 75000+ 1400 X +100 Y. \end{equation}$

(3.2)

For production, we use the C-D production function, which is given by:

$\begin{equation} P = A X^a Y^b \end{equation}$

(3.3)

where

$P$ represents the amount of production in liters.

After Eq (3.3) we have,

$\begin{equation} \ln P = \ln A +a \ln X+ b \ln Y. \end{equation}$

(3.4)

The analyzed data for the water industry with two inputs, i.e., human labor and capital, and production as an output for the year 2022 is given in Table 1.

Table 1. Inputs and output values of the water industry with two inputs.

Months	Production( $P$ )	Labors( $X$ )	Capital( $Y$ )	$Ln(P)$	$Ln(X)$	$Ln(Y)$
	in litres	in hours	in kg
Jan.	36,000	300	900	10.49127	5.703782	6.802395
Feb.	34,000	272	833	10.43412	5.605802	6.725034
March	36,700	315	924	10.51053	5.752573	6.828712
April	37,000	330	942	10.51867	5.799093	6.848005
May	35,700	278.5	874	10.48291	5.629418	6.77308
June	36,400	321	922	10.50232	5.771441	6.826545
July	33,600	269	860	10.42228	5.594711	6.756932
Aug.	35,500	308	950	10.47729	5.7301	6.856462
Sep.	36,875	338	957	10.51529	5.823046	6.863803
Oct.	37,000	321	914	10.51867	5.771441	6.817831
Nov.	33,600	266	875	10.42228	5.583496	6.774224
Dec.	36,500	336	890	10.50507	5.817111	6.791221

| Show Table

DownLoad: CSV

We have used ordinary least squares regression in order to find the structural parameters of the given C-D production function. The data analysis has been done using Microsoft Excel, which is given in Regressions 1 and 2 for two and three inputs, respectively.

Figure Regression 1. Linear regression for two inputs.

DownLoad: Full-Size Img PowerPoint

Figure Regression 2. Linear regression for three inputs.

DownLoad: Full-Size Img PowerPoint

Structural parameters for two inputs are as follow:

$\begin{equation*} \begin{split} & A = 3212.468, \\ & a = 0.3568,\\ & b = 0.0542,\\ \end{split} \end{equation*}$

$\begin{equation} L(X,Y,\lambda) = C(X,Y)+ \lambda H(X,Y), \end{equation}$

(3.5)

where,

$\begin{equation*} H(X,Y) = 36000-3212.468 X^{0.3568} Y^{0.0542}. \end{equation*}$

Using Eq (3.2) and H(x, y) in Eq (3.5), we have,

$\begin{equation} L = 75000+1400 X+100 Y+\lambda(36000-3212.468 X^{0.3568} Y^{0.0542}). \end{equation}$

(3.6)

Now taking partial derivatives of L with respect to $X, Y$ and $\lambda$ respectively and equating to zero, i.e.,

$\begin{equation*} \frac{\partial L}{\partial X} = 0, \end{equation*}$

we get,

$\begin{equation*} 1400 -1146.2085 \lambda X^{-0.6432} Y^{0.0542} = 0, \end{equation*}$

which implies that

$\begin{equation} \lambda = \frac{1.2214}{X^{-0.6432} Y^{0.0542}}. \end{equation}$

(3.7)

Now, when

$\begin{equation*} \frac{\partial L}{\partial Y} = 0, \end{equation*}$

we get

$\begin{equation*} 100- 174.1157 \lambda X^{0.3568} Y^{-0.9458} = 0, \end{equation*}$

which implies

$\begin{equation} \lambda = \frac{100}{X^{0.3568} Y^{-0.9458}}. \end{equation}$

(3.8)

Similarly, when

$\begin{equation*} \frac{\partial L}{\partial \lambda} = 0, \end{equation*}$

we get,

$\begin{equation} 36000-3212.468 X^{0.3568} Y^{0.0542} = 0. \end{equation}$

(3.9)

Comparing (3.7) and (3.8),

$\begin{equation*} \frac{1.2214}{X^{-0.6432} Y^{0.0542}} = \frac{100}{X^{0.3568} Y^{-0.9458}}, \end{equation*}$

we get

$\begin{equation} Y = 2.1266 X. \end{equation}$

(3.10)

Putting Eq (3.10) in Eq (3.9), we have

$\begin{equation*} 36000-3212.468 X^{0.3568} (2.1266 X )^{0.0542} = 0. \end{equation*}$

Solving for $X$ , we have

$\begin{equation} X = 323.7554. \end{equation}$

(3.11)

From (3.10), we have

$\begin{equation} Y = 688.4982, \end{equation}$

(3.12)

and using the values of X and Y in Eq (3.2), we have,

$\begin{equation} C(X,Y) = 597107. \end{equation}$

(3.13)

Table 2 compares the actual cost with the estimated cost.

Table 2. Actual and estimated cost values for water industry using two inputs.

Months	Estimated cost	Actual cost	Error (R)	Square of the error $(R^{2})$
	in rupees	in rupees
Jan.	585,000	608,000	23000	529000000
Feb.	539,100	570,000	30900	954810000
March	608,400	588,000	-20400	416160000
April	631,200	608,000	-23200	538240000
May	552,300	598,000	45700	2088490000
June	616,600	621,000	4400	19360000
July	537,600	574,000	36400	1324960000
Aug.	601,200	564,000	-37200	1383840000
Sep.	643,900	622,000	-21900	479610000
Oct.	615,800	598,000	-17800	316840000
Nov.	534,900	582,000	47100	2218410000
Dec.	634,400	632,000	-2400	5760000
			$\sum R^{2}$	10275480000
			Mean square error	856290000

| Show Table

DownLoad: CSV

Figure 1 represents the actual and theoretical costs of the industry. For example, in January 2022, we can see that the actual cost was greater than the estimated cost, but from the data analysis, our calculated cost is less than both the actual and estimated cost. Similarly, if we compare the cost from the table with our calculations, we can see the difference. It means that there is a sufficient difference in both costs, meaning that costs are minimized to a great extent.

Figure 1. Estimated cost versus actual cost of the industry using two inputs.

DownLoad: Full-Size Img PowerPoint

3.2. Cost minimization in the water industry using the three-factor C-D function

In this case, the cost function consists of three factors: capital, chemicals, and labor. The cost function is given by:

$\begin{equation} C(X,Y,Z) = \alpha_1+ \alpha_2 X+\alpha_3 Y+\alpha_4 Z \end{equation}$

(3.14)

where

$\alpha_1$ represents the fixed cost of the industry,

$\alpha_2$ represents the unit price of labor per hour,

$\alpha_3$ represents the unit price of chemicals per kg,

$\alpha_4$ represents the unit price of capital per kg,

$X$ represents the number of labor hours,

$Y$ represents the amount of chemicals in kg,

$Z$ represents the amount of capitals in kg,

where

$\begin{equation*} \alpha_1 = 75000,\; \; \alpha_2 = 1250,\; \; \alpha_3 = 2000 ,\; \; \alpha_4 = 100. \end{equation*}$

Putting the above values in Eq (3.14), we get

$\begin{equation} C(X,Y,Z) = 75000+ 1250 X +2000 Y +100 Z. \end{equation}$

(3.15)

Now, the production function, as in ^[26], is given by

$\begin{equation} P = A X^{a_1} Y^{a_2} Z^{a_3} \end{equation}$

(3.16)

where

$P$ represents amount of production in liters.

After linearizing, we get

$\begin{equation} \ln P = \ln A +a_1 \ln X+ a_3 \ln Y+ a_3 \ln Z. \end{equation}$

(3.17)

A data analysis of the water industry consisting of three inputs for the year 2022 is given in Table 3. The three inputs are human labor, capital, and chemicals, respectively.

Table 3. Inputs and output values of the water industry with three inputs.

Months	Labors( $X$ )	Chemicals( $Y$ )	Capitals( $Z$ )	Production( $P$ )	$Ln(X)$	$Ln(Y)$	$Ln(Z)$	$Ln(P)$
	in hours	in kg	in kg	in litres
Jan.	240	60	900	36,000	5.4806	4.0943	6.8023	10.4912
Feb.	220	52	833	34,000	5.3936	3.9512	6.7250	10.4341
March	252	63	924	36,700	5.5294	4.1431	6.8287	10.5105
April	264	66	942	37,000	5.5759	4.1896	6.8480	10.5186
May	222	56.5	874	35,700	5.4026	4.0342	6.773	10.4829
June	260	61	922	36,400	5.5606	4.1108	6.8265	10.5023
July	215	54	860	33,600	5.3706	3.9889	6.7569	10.4222
Aug.	250	58	950	35,500	5.5214	4.0604	6.8564	10.4772
Sep.	268	70	957	36,875	5.5909	4.2484	6.8638	10.5152
Oct.	256	65	914	37,000	5.5451	4.1743	6.8178	10.5186
Nov.	212	54	875	33,600	5.3565	3.9889	6.7742	10.4222
Dec.	262	74	890	36,500	5.5683	4.3040	6.7912	10.5050

| Show Table

DownLoad: CSV

This data is collected from the water industry for the year 2022, and all results are calculated on a monthly basis.

From these analyses, we have,

$\begin{equation*} A = 3938.188,\; \; a_1 = 0.2944,\; \; a_2 = 0.0673,\; \; a_3 = 0.0457. \end{equation*}$

Putting these values in (3.16), we get

$\begin{equation} P = 3938.188 X^{0.2944}Y^{0.0673} Z^{0.0457}. \end{equation}$

(3.18)

Using the Lagrange multipliers method, we have,

$\begin{equation} L(X,Y,Z,\lambda) = C(X,Y,Z )+ \lambda H(X,Y,Z ), \end{equation}$

(3.19)

where

$\begin{equation} H(X,Y,Z) = 36000-3938.188 X^{0.2944} Y^{0.0673} Z^{0.0457}. \end{equation}$

(3.20)

Putting Eqs (3.15) and (3.20) in Eq (3.19), we obtain

$\begin{equation} L = 75000+ 1250 X +2000 Y +100 Z + \\ \lambda(36000-3938.188 X^{0.2944} Y^{0.0673} Z^{0.0457}). \end{equation}$

(3.21)

Now, taking partial derivatives of L with respect to $X$ , $Y$ , $Z$ and $\lambda$ and setting all partial derivatives equal to zero, we have,

$\begin{equation*} \frac{\partial L}{\partial X} = 0. \end{equation*}$

we then get

$\begin{equation*} 1250-1159.4925 \lambda X^{-0.7056} Y^{0.0673} Z^{0.0457} = 0. \end{equation*}$

Solving the above equation for $\lambda$ , we get

$\begin{equation} \lambda = \frac{1250}{1159.4925 X^{-0.7056} Y^{0.0673} Z^{0.0457}}. \end{equation}$

(3.22)

also setting

$\begin{equation*} \frac{\partial L}{\partial Y} = 0, \end{equation*}$

we get

$\begin{equation*} 2000-265.04 \lambda X^{0.2944} Y^{-0.9327} Z^{0.0457} = 0. \end{equation*}$

Solving above equation for $\lambda$ , we get

$\begin{equation} \lambda = \frac{7.5460}{ X^{0.2944} Y^{-0.9327} Z^{0.0457}}. \end{equation}$

(3.23)

Finally setting

$\begin{equation*} \frac{\partial L}{\partial Z} = 0, \end{equation*}$

we get

$\begin{equation*} 100-179.9751 \lambda X^{0.2944} Y^{0.0673} Z^{-0.9543} = 0. \end{equation*}$

Solving above equation for $\lambda$ , we get

$\begin{equation} \lambda = \frac{100}{179.9751 X^{0.2944} Y^{0.0673} Z^{-0.9543}}. \end{equation}$

(3.24)

Similarly, when

$\begin{equation*} \frac{\partial L}{\partial \lambda} = 0, \end{equation*}$

we get

$\begin{equation} 36000-3938.188 X^{0.2944} Y^{0.0673} Z^{0.0457} = 0. \end{equation}$

(3.25)

Comparing (3.22) and (3.23), we get

$\begin{equation*} \frac{1250}{1159.4925 X^{-0.7056} Y^{0.0673} Z^{0.0457}} = \frac{7.5460}{X^{0.2944} Y^{-0.9327} Z^{0.0457}} \end{equation*}$

which implies

$\begin{equation} X = 6.9996 Y. \end{equation}$

(3.26)

Now, comparing (3.23) and (3.24), we get

$\begin{equation*} \frac{7.5460}{ X^{0.2944} Y^{-0.9327} Z^{0.0457}} = \frac{100}{179.9751 X^{0.2944} Y^{0.0673} Z^{-0.9543}}, \end{equation*}$

which implies

$\begin{equation} Z = 13.5809 Y. \end{equation}$

(3.27)

Using (3.25) and (3.26) in (3.27)

$\begin{equation*} (6.9996 Y)^{0.2944} Y^{0.0673} (13.5809 Y)^{0.0457} = 9.1412, \end{equation*}$

yields,

$\begin{equation*} Y^{0.4074} = 4.5755. \end{equation*}$

Thus,

$\begin{equation} Y = 41.7932. \end{equation}$

(3.28)

Putting (3.28) in (3.26) and (3.27) respectively, we get

$\begin{equation} \begin{split} & X = 292.5356,\\ & Z = 567.5892, \\ \end{split} \end{equation}$

(3.29)

and using the values of X, Y and Z in Eq (3.15), we have

$\begin{equation} C(X,Y,Z) = 581014. \end{equation}$

(3.30)

Table 4 presents actual and estimated costs on a monthly basis for the water industry for the year 2022.

Table 4. Actual and estimated cost values for the water industry using three inputs.

Months	Estimated cost	Actual cost	Error (R)	Square of the error $(R^{2}$ )
	in rupees	in rupees
Jan.	585,000	608,000	23000	529000000
Feb.	537,300	570,000	32700	1069290000
March	608,400	588,000	-20400	416160000
April	631,200	608,000	-23200	538240000
May	552,900	598,000	45100	2034010000
June	614,200	621,000	6800	46240000
July	537,750	574,000	36250	1314062500
Aug.	598,500	564,000	-34500	1190250000
Sep.	645,700	622,000	-23700	561690000
Oct.	616,400	598,000	-18400	338560000
Nov.	535,500	582,000	46500	2162250000
Dec.	639,500	632,000	-7500	56250000
			$\sum R^{2}$	10256002500
			Mean square error	854666875

| Show Table

DownLoad: CSV

Figure 2 represents the actual and estimated cost of the industry for the three-factors C-D production function.

Figure 2. Estimated cost versus actual cost for the industry for three inputs.

DownLoad: Full-Size Img PowerPoint

4. Results and comparative analysis

We have solved the C-D production function with two and three factors of production using the Lagrange multiplier method with the ordinary least squares method. Moreover, this is an optimal solution approach for the C-D production function with three factors of production using the Lagrange multiplier method with the ordinary least squares method. Despite the fact that the proposed approach is a different solution technique as compared to the existing solution techniques in the literature, we still compared the general features of the proposed methodology for the water industry with different factors to Nervole's approach. In the following, we compared the presented solution approach with Nervole's approach ^[22]. In Table 5, we have differentiated Nerlove's C-D function and the C-D production function with the Lagrange multiplier method.

Table 5. Comparative analysis between Nerlove's C-D function and the C-D production function with the Lagrange multiplier method.

No.	Nerlove's C-D cost function	C-D production function
1	Nerlove used the C-D cost model	We used the C-D production function
2	Nerlove used the function to estimate the cost	We used the given function as an output
3	Nerlove approach needs much algebra for evaluation	It has a simple implementation
4	Computational complexity is too much	Having less computational complexity

| Show Table

DownLoad: CSV

In addition, the paired t-test is used to determine if there is a significant difference between the means of two related data sets as well as to find the mean square error from the findings of both the two- and three-factor C-D production functions. In the case of the findings of the two-factor C-D function, the outcomes of a paired t-test on two data sets result in a P-value of 0.547 and degrees of freedom (df) of 11. The paired t-test compares the means of the two variables to determine if there is a significant difference between them. The null hypothesis is that there is no significant difference between the means of the two variables. This means that there is no significant difference between the means of the two variables at the 5% level of significance. This showed that the P-value is greater than the critical value of 0.05, indicating that we cannot reject the null hypothesis.

In the case of the findings from the three-factor C-D function, the findings of a paired t-test result in a P-value of 0.559 and degrees of freedom (df) of 11. Moreover, the findings from the three-factor C-D using a paired t-test were also conducted, which also showed that there is no significant difference between the means of the two variables with the same 5% level of significance.

Based on the findings of the paired t-test, we conclude that there is no significant difference between the means of the two related variables from the findings of two and three factors in the C-D production functions. This shows that the C-D production function plays a key role in the production problem when using the Lagrange multiplier method with the ordinary least squares method. Furthermore, we have worked on the cost comparison of two and three factors in the C-D production function. In the case of two factors of production, the cost value is $597107$ per unit, while in the case of three factors of production, the cost value is $581014$ per unit. In both cases, our calculated cost is less than the actual cost of the industry. Besides, the cost calculation for the three-factor C-D production function is less than that of the two-factor C-D production function. Clearly, we can see the differences in Tables 6–9.

Table 6. Data set for the two-factor C-D production function.

Paired samples statistics (Pair 1)
	Mean	N	Std. deviation	Std. error mean
Estimated cost	591700.0000	12	40670.96351	11740.69586
Actual cost	597083.3333	12	21997.76159	6350.20679
Paired samples correlations (Pair 1)
		N	Correlation	Sig.
Estimated cost and Actual cost		12	0.690	0.013

| Show Table

DownLoad: CSV

Table 7. Paired sample test.

	Paird differences
	Mean	Std. derivation	Std. eror mean	95s% confidence interval of the difference		t	df	Sig. (2-tailed)
	Mean	Std. derivation	Std. eror mean	Lower	Upper	t	df	Sig. (2-tailed)
Pair1 estimated cost-actual cost	5383.33333	30041.96560	8672.36846	-24471.08762	13704.42096	-0.621	11	0.547

| Show Table

DownLoad: CSV

Table 8. Data set for the three-factor C-D production function.

Paired samples statistics (Pair 1)
	Mean	N	Std. deviation	Std. error mean
Estimated cost	591862.5000	12	41320.23016	11928.12300
Actual cost	597083.3333	12	21997.76159	6350.20679
Paired samples correlations (Pair 1)
		N	Correlation	Sig.
Estimated cost and Actual cost		12	0.709	0.010

| Show Table

DownLoad: CSV

Table 9. Paired sample test.

	Paird differences
	Mean	Std. derivation	Std. error mean	95s% Confidence interval of the difference		t	df	Sig. (2-tailed)
	Mean	Std. derivation	Std. error mean	Lower	Upper	t	df	Sig. (2-tailed)
Pair1 estimated cost-Actual cost	5220.83333	30043.78213	8672.89265	-24309.74179	13866.07513	-0.602	11	0.559

| Show Table

DownLoad: CSV

5. Conclusions

In this paper, the optimal solution developed by C-D is carried out with different production factors. From these analyses, we concluded that the C-D production function plays a key role in the production problem when using the Lagrange multiplier method with the ordinary least squares method. Moreover, we solved the constrained optimization problem with a two-factor and three-factor C-D production function using the Lagrange multiplier with the ordinary least squares method. In the case of two factors of production, the cost value is $597107$ per unit, while in the case of three factors of production, the cost value is $581014$ per unit. This showed that, with more production factors in the C-D production function, the cost value is minimized to a high extent. This validates that the C-D production function with more factors using the Lagrange multiplier is more effective than previous approaches in literature. Moreover, the presented solution methodology is compared to Nervole's C-D production function. This means that the individual expression of each factor as an input has a key role in obtaining the best optimized results. Moreover, we optimized the overall cost of the water industry using the three-factor C-D production function as an application of C-D production using the Lagrange multiplier method with the ordinary least squares method.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

We are thankful to the National Natural Science Foundation of China (under research grant no. 22150410332) and the start-up foundation for the introduction of talent at Jiangsu University of Science and Technology, China.

Conflict of interest

The authors declare that they have no competing interest.

References

[1]	J. A. Snyman, New gradient-based trajectory and approximation methods, In: Practical mathematical optimization, Boston: Springer, 97 (2005), 97–150. https://doi.org/10.1007/0-387-24349-6_4
[2]	J. Nocedal, S. J. Wright, Numerical optimization, New York: Springer, 2006. https://doi.org/10.1007/978-0-387-40065-5
[3]	J. E. Dennis, R. B. Schnabel, Numerical methods for unconstrained optimization and nonlinear equations, Society for Industrial and Applied Mathematics, 1996. https://doi.org/10.1137/1.9781611971200
[4]	P. C. Pendharkar, J. A. Rodger, Nonlinear programming and genetic search application for production scheduling in coal mines, An. Oper. Res., 95 (2000), 251–267. https://doi.org/10.1023/A:1018958209290 doi: 10.1023/A:1018958209290
[5]	D. P. Bertsekas, Nonlinear programming, J. Oper. Res. Soc., 48 (1997), 334. https://doi.org/10.1057/palgrave.jors.2600425 doi: 10.1057/palgrave.jors.2600425
[6]	D. G. Luenberger, Y. Ye, Linear and nonlinear programming, 2 Eds., Addison-Wesley Publishing Co., 1984.
[7]	Z. Liu, Z. Liu, M. Liu, J. Wang, Optimization of flow shop scheduling in precast concrete component production via mixed-integer linear programming, Adv. Civ. Eng., 2021 (2021), 6637248. https://doi.org/10.1155/2021/6637248 doi: 10.1155/2021/6637248
[8]	M. M. Hossain, A. K. Majumder, T. Basak, An application of non-linear Cobb-Douglas production function to selected manufacturing industries in Bangladesh, Open J. Stat., 2 (2012), 460–468. https://doi.org/10.4236/ojs.2012.24058 doi: 10.4236/ojs.2012.24058
[9]	H. Luo, Y. Chen, An allometric algorithm for fractal-based Cobb-Douglas function of geographical systems, Discrete Dyn. Nat. Soc., 2014 (2014), 910457. https://doi.org/10.1155/2014/910457 doi: 10.1155/2014/910457
[10]	B. Gajdzik, R. Gawlik, Choosing the production function model for an optimal measurement of the restructuring efficiency of the Polish metallurgical sector in years 2000–2015, Metals, 8 (2018), 23. https://doi.org/10.3390/met8010023 doi: 10.3390/met8010023
[11]	F. Zhang, Q. Tan, C. Zhang, S. Guo, P. Guo, A regional water optimal allocation model based on the Cobb-Douglas production function under multiple uncertainties, Water, 9 (2017), 923. https://doi.org/10.3390/w9120923 doi: 10.3390/w9120923
[12]	O. Yankovyi, V. Koval, L. Lazorenko, O. Poberezhets, M. Novikova, V. Gonchar, Modeling sustainable economic development using production functions, Studies Appl. Econ., 39 (2021), 1–14. https://doi.org/10.25115/eea.v39i5.5090 doi: 10.25115/eea.v39i5.5090
[13]	L. A. Quezada-Téllez, G. Fernández-Anaya, D. Brun-Battistini, B. Nuñez-Zavala, J. E. Macías-Díaz, An economic model for OECD economies with truncated $M$ -derivatives: exact solutions and simulations, Mathematics, 9 (2021), 1780. https://doi.org/10.3390/math9151780 doi: 10.3390/math9151780
[14]	S. Liu, W. Jiao, Q. Min, J. Yin, The influences of production factors with profit on agricultural heritage systems: a case study of the rice-fish system, Sustainability, 9 (2017), 1842. https://doi.org/10.3390/su9101842 doi: 10.3390/su9101842
[15]	S. Saha, J. Sarkar, A. Dwivedi, N. Dwivedi, A. M. Narasimhamurthy, R. Roy, A novel revenue optimization model to address the operation and maintenance cost of a data center, J. Cloud Comp., 5 (2016), 1. https://doi.org/10.1186/s13677-015-0050-8 doi: 10.1186/s13677-015-0050-8
[16]	G. Xu, X. Yin, G. Wu, N. Gao, Rethinking the contribution of land element to urban economic growth: evidence from 30 provinces in China, Land, 11 (2022), 801. https://doi.org/10.3390/land11060801 doi: 10.3390/land11060801
[17]	A. Barbagallo, S. G. Bianco, G. Tensor, Variational inequalities: theoretical results, numerical methods and application to an equilibrium Model, J. Nonlinear Var. Anal., 4 (2020), 87–105. https://doi.org/10.23952/jnva.4.2020.1.07 doi: 10.23952/jnva.4.2020.1.07
[18]	M. Zhang, J. Geng, S. Wu, A new infeasible interior-point algorithm with full-Newton steps for linear programming based on a kernel function, J. Nonlinear Funct. Anal., 2021 (2021), 31. https://doi.org/10.23952/jnfa.2021.31 doi: 10.23952/jnfa.2021.31
[19]	X. Wang, Y. Fu, Some characterizations of the Cobb-Douglas and CES production functions in microeconomics, Abstr. Appl. Anal., 2013 (2013), 761832. https://doi.org/10.1155/2013/761832 doi: 10.1155/2013/761832
[20]	M. M. Hossain, T. Basak, A. K. Majumder, Application of non-linear Cobb-Douglas production function with autocorrelation problem to selected manufacturing industries in Bangladesh, Open J. Stat., 3 (2013), 173–178. https://doi.org/10.4236/ojs.2013.33019 doi: 10.4236/ojs.2013.33019
[21]	K. V. Murthy, Arguing a case for Cobb-Douglas production function, Rev. Commer. Studies, 20-21 (2002), 1–17.
[22]	M. Green, Electric cost modeling calculations: regulations, technology, and the role of renewable energy, 2 Eds., 2021. https://doi.org/10.1016/C2019-0-03162-7
[23]	M. Asghari, A. M. Fathollahi-Far, S. M. Mirzapour Al-e-hashem, M. A. Dulebenets, Transformation and linearization techniques in optimization: a state-of-the-art survey, Mathematics, 10 (2022), 283. https://doi.org/10.3390/math10020283 doi: 10.3390/math10020283
[24]	T. Franik, E. Franik, Application of nonLinear programming for optimization of factors of production in mining industry, Fribourg: Alcon Pharmacetutical Ltd anovartis Company, 2009.
[25]	T. Franik, Productivity of hard coal mining during the reform period against the background of changes in the mining and quarrying section, Econ. Miner. Mater., 21 (2005), 47–61.
[26]	G. E. Vîlcu, A geometric perspective on the generalized Cobb-Douglas production functions, Appl. Math. Lett., 24 (2011), 777–783. https://doi.org/10.1016/j.aml.2010.12.038 doi: 10.1016/j.aml.2010.12.038
[27]	T. Franik, The analysis of productiveness of branch of mining hard coal in Poland with using the function of production, Gospod. Surowcami Miner., 23 (2007), 77–91.

This article has been cited by:

1.	Dmitry B. Rokhlin, Gennady A. Ougolnitsky, A Simple Model for Targeting Industrial Investments with Subsidies and Taxes, 2024, 12, 2227-7390, 822, 10.3390/math12060822
2.	Shakoor Muhammad, Taimur Ali, Yousaf Khan, Saleem Abdullah, Optimal photovoltaic location selection systems through TOPSIS, AHP, and GIS techniques: a case study in Pakistan, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05860-2

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(1695) PDF downloads(87) Cited by(2)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(4) / Tables(9)

AIMS Mathematics

Industrial optimization using three-factor Cobb-Douglas production function of non-linear programming with application

Related Papers:

Abstract

1. Introduction

2. Cost minimization

3. Production function and its structural parameters

3.1. Cost minimization in the water industry using the two-factor C-D production function

3.2. Cost minimization in the water industry using the three-factor C-D function

4. Results and comparative analysis

5. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Mathematics

Industrial optimization using three-factor Cobb-Douglas production function of non-linear programming with application

Related Papers:

Abstract

1. Introduction

2. Cost minimization

3. Production function and its structural parameters

3.1. Cost minimization in the water industry using the two-factor C-D production function

3.2. Cost minimization in the water industry using the three-factor C-D function

4. Results and comparative analysis

5. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog