Scientific Research

An Academic Publisher

Multi-Objective Mathematical Model for the Optimal Time to Harvest Sugarcane

**Author(s)**Leave a comment

^{1,2}, Wisanlaya Pornprakun

^{1,2}, Santipong Prasattong

^{1}, Chanasak Baitiang

^{1}

1. Introduction

Sugarcane is an important crop in Thailand because sugarcane is a raw material of the sugar industry, which is one of the top five agricultural products that makes revenue for the country each year. In the last ten years, the link between agriculture and industry and the structure of sugar markets have both changed. These changes have affected the amount and direction of the global sugar trade [1] . Thailand is an exporter with export volume ranked second in the world and and it also has a large role in the Asian market. Because Thailand has the advan- tage of location and an oversupply of sugar it can respond quickly to changes in demand in the Asian market. The forecast volume of exports in the next ten years is 11 million tons, which is an increase of 69 percent over the average volume of exports in the past ten years. Thus, sugarcane is an important eco- nomic crop in Thailand.

A survey by the Promotion Center of the Cane and Sugar Industry has re- ported that, in the crop year 2012/13 the average cost of production of sugarcane including the cost of transportation was 1196.31 baht per ton (transport cost was 143.61 baht per ton). This average cost of production was an increase of 246.31 baht per ton over the cost of 950 baht per ton in crop year 2011/12. This increase was because of increases in land rents and prices of labour, fertilizer and ma- terials. In crop year 2012/13, sugarcane production for the whole country was at the level of 94.64 million tons from sugarcane plantations of 9.33 million rai for an average yield of 10.14 tons per rai. This production represented a decrease from crop year 2011/12 when the annual production was 97.98 million tons. The decrease occurred because of drought during the crop year 2012/13 when the rains came later than usual.

Many researchers have discussed the renewable energy which can be pro- duced from the biomass of agricultural products. For example, the biomass can be used to produce ethanol which can then increase the profit and reduce the cost of the agricultural product. In 2008, Helenice de Oliveira Florentino [2] studied multi-objective optimization of the economics of sugarcane harvest bio- mass in Brazil. The aim of this work was to develop a model to optimize plant variety selection, to minimize cost of the residual biomass transfer process, to evaluate the economics of using this material, and to address sucrose production subject to planting area constraints, and considering distance from planting area to processing center. For this, 0 - 1 multiple objective linear programming tech- niques were used. The results showed the viability of the model when selecting varieties and the increased profit that could be obtained from residual biomass use. In 2010, Maximiliano Salles Scarpari and Edgar Gomes Ferreira de Beauclair [3] used linear programming to develop an optimized plan for sugarcane far- ming in Brazil. The program language used was General Algebraic Modelling System (GAMS) as this system was seen to be an excellent tool to allow profit maximization and harvesting time schedule optimization in the sugar mill stu- died. The results supported this optimized planning model as being a very useful tool for sugarcane management. In 2012, Francisco Regis Abreu Gomes [4] studied a bi-objective mathematical model for choosing sugarcane varieties with harvest residual biomass in energy co-generation. This study developed a bi- objective mathematical model for choosing sugarcane varieties that result in maximum revenue from electricity sales and minimum gathering cost of sugar- cane by harvesting residual biomass. The approach used to solve the proposed model was based on the $\epsilon $ -constraints method. Experiments were performed using real data from sugarcane varieties and costs and showed effectiveness of the model and method proposed. These experiments showed the possibility of increasing net revenue from electricity sale, i.e., discounting the cost increase with residual biomass gathering, by up to 98.44%.

At present, the Thai government determines the price of sugarcane each year depending on its quality (sweetness) and weight. In general, the sweetness in- creases over time while the weight decreases. In this paper, a model is developed to find the optimal harvest time of sugarcane in order to maximize the agri- culturists revenue and to minimize the cost. The mathematical model is for- mulated in a multi-objective optimization framework under the Office of Cane and Sugar Board (OCSB) conditions. The $\epsilon $ -constraints method is used to solve the multi-objective optimization problem. The results for production of fresh sugarcane (not burnt before cutting) and fired sugarcane (was burnt before cutting) in crop years 2012/13, 2013/14 and 2014/15 are compared and analyzed.

2. Multi-Objective Optimization

A multi-objective optimization problem has a number of objective functions which are to be minimized or maximized. As in a single-objective optimization problem, the problem usually has a number of constraints which any feasible solution (including the optimal solution) must satisfy. In its general form, the mul- ti-objective optimization problem (MOOP) can be stated as follows (see, e.g., [5] ):

$\begin{array}{l}\text{Minimize}{f}_{m}\left(x\right),m=1,2,\cdots ,M\\ \text{subjectto}{g}_{j}\left(x\right)\ge 0,j=1,2,\cdots ,J\\ \text{}{h}_{k}\left(x\right)=0,k=1,2,\cdots ,K\\ \text{}{x}_{i}^{\left(L\right)}\le {x}_{i}\le {x}_{i}^{\left(U\right)},i=1,2,\cdots ,n\end{array}$ (1)

A solution $x$ is a vector of $n$ decision variables : $x={\left({x}_{1},{x}_{2},\cdots ,{x}_{n}\right)}^{\text{T}}$ . The last set of constraints are called variable bounds as they restrict each decision variable ${x}_{i}$ to take a value within a lower bound ${x}_{i}^{\left(L\right)}$ and an upper bound ${x}_{i}^{\left(U\right)}$ . The feasible region $D$ for the MOOP is the set of vectors $x$ that satisfy all constraints. If each objective function ${f}_{m}\left(x\right)$ is denoted by ${z}_{m}$ and the vector of all objective functions is denoted by $z$ , then the objective function space can be defined as:

$Z=\left\{z=\left({z}_{1},{z}_{2},\cdots ,{z}_{M}\right)|{z}_{m}={f}_{m}\left(x\right),\forall x\in D,m=1,2,\cdots ,M\right\}$ (2)

In this work, the $\epsilon $ -constraints method [5] was used to solve the MOOP. This method consists of reformulating a multi-objective problem by choosing the most important objective while maintaining other objectives constrained by upper bounds defined by a decision maker. For example, if ${f}_{\mu}\left(x\right)$ is selected as the most important objective, then the problem can be reformulated as follows:

$\begin{array}{l}\text{Minimize}{f}_{\mu}\left(x\right)\\ \text{subjectto}{f}_{m}\left(x\right)\le {\epsilon}_{m},m=1,2,\cdots ,M;m\ne \mu \\ \text{}{g}_{j}\left(x\right)\ge 0,j=1,2,\cdots ,J\\ \text{}{h}_{k}\left(x\right)=0,k=1,2,\cdots ,K\\ \text{}{x}_{i}^{\left(L\right)}\le {x}_{i}\le {x}_{i}^{\left(U\right)},i=1,2,\cdots ,n\end{array}$ (3)

where ${\epsilon}_{m}$ is an upper bound of objective $m$ , $m=1,2,\cdots ,M;m\ne \mu $ and $D$ is the set of feasible solutions to the $\epsilon $ -constrained problem.

3. Mathematical Model

In this section, the mathematical model is formulated in a multi-objective op- timization framework under the Office of the Cane and Sugar Board (OCSB) conditions. The main purpose is to find the optimal harvest time of sugarcane in order to maximize the agriculturists’ revenue and minimize the cost. This pro- blem is a bi-objective optimization problem to maximize revenue and minimize cost. For the $\epsilon $ -constraint formulation, the maximum of revenue from the sugarcane sale will be chosen as the objective function and the gathering cost of production of the sugarcane will be bounded above by an $\epsilon $ -constraint.

Revenue of selling: The government determines the sugarcane prices which are based on two main factors; 1) weight and 2) commercial cane sugar (C.C.S.).

1) Revenue from the weight of sugarcane: The OCSB classifies sugarcane going into factories into 2 types; a) fresh sugarcane and b) fired sugarcane. As determined by the government, agriculturists who sell fired sugarcane will have 20 baht/ton deducted from the price of sugarcane based on weight. The factory will share this amount of money between agriculturists who sell fresh sugarcane and increase their price based on weight by a maximum of 70 baht/ton. The price of fired sugarcane based on weight ${P}_{w}\left(B\right)$ is therefore defined by

${P}_{w}\left(B\right)={P}_{w}-20,$ (4)

where ${P}_{w}$ is the basic price of sugarcane based on weight set by the govern- ment (baht/ton). The price of fresh sugarcane based on weight ${P}_{w}\left(A\right)$ is de- fined by

${P}_{w}\left(A\right)={P}_{w}+\frac{20{a}_{j,k}\left(B\right)}{{a}_{j,k}\left(A\right)},$ (5)

where ${a}_{j,k}\left(A\right)$ is the amount of fresh sugarcane (tons) from planted area $j$ in the harvest at time $k$ and ${a}_{j,k}\left(B\right)$ is the amount of fired sugarcane (tons) from planted area $j$ in the harvest at time $k$ . Reasonable values for the total planted areas of sugarcane were estimated from OCSB data. The actual values for ${a}_{j,k}$ were computed by the optimization program.

Therefore the revenue from weight of sugarcane is defined by

${P}_{1}\left(i\right)={P}_{w}\left(i\right){a}_{j,k}\left(i\right),\text{}i=A,B$ (6)

2) Revenue from C.C.S. of sugarcane. The price per ton based on C.C.S is defined by

${P}_{2}\left(i\right)=\left[{P}_{c}\left(1+0.06{y}_{j,k}\right)\right]{a}_{j,k}\left(i\right),\text{}i=A,B$ (7)

where ${P}_{c}$ is the price per ton determined by the government for sugarcane with 10 C.C.S.

${y}_{j,k}=\text{C}\text{.C}\text{.S}.-10$ , where C.C.S. is the average C.C.S. from sugarcane in planted area $j$ in the harvest at time $k$ . The factor 0.06 is the rate of change of price per 1 C.C.S. change from the base level of 10.

Therefore, revenue $\left[R{V}_{j,k}\left(i\right)\right]$ from the sale of sugarcane from planted area $j$ in the harvest at time $k$ is determined by adding Equation (6) and Equation (7) as shown in Equation (8)

$R{V}_{j,k}\left(i\right)={P}_{1}\left(i\right)+{P}_{2}\left(i\right)=\left\{{P}_{w}\left(i\right)+\left[{P}_{c}\left(1+0.06{y}_{j,k}\right)\right]\right\}{a}_{j,k}\left(i\right),\text{}i=A,B$ (8)

Gathering cost of production: The gathering cost of production can be sepa- rated into two parts: 1) average cost of production on the farm and 2) cost of transport. The total gathering cost $G{C}_{j,k}\left(i\right)$ baht/ton of sugarcane production from planted area $j$ in the harvest at time $k$ is given by

$G{C}_{j,k}\left(i\right)=\left({C}_{j}+{C}_{{T}_{j}}\right){a}_{j,k}\left(i\right),\text{}i=A,B$ (9)

where ${C}_{j}$ is the average total cost of production of sugarcane on farms in area $j$ and ${C}_{{T}_{j}}$ is cost of transport to factories for sugarcane produced in area $j$ . The average total cost of production ${C}_{j}$ consists of fixed costs such as farm rent and depreciation of equipment, and variable costs such as labor, materials and interest rates.

3.1. The Mathematical Model of Fresh Sugarcane

The multi-objective mathematical model of fresh sugarcane is described by Equa- tions (10)-(13). The objective function Equation (10) maximizes revenue from the fresh sugarcane selling. The constraint Equation (11) represents the second objective of the problem which is to minimize gathering cost of production, with upper bound given by ${\epsilon}_{e}$ . The decision variables ${X}_{j,k}$ are defined by ${X}_{j,k}=1$ means that planted area $j$ is harvested at time $k$ and ${X}_{j,k}=0$ means that $j$ is not harvested at $k$ . In the constraint set, Equation (12) and Equation (13) ensure that in each area the sugarcane is harvested only at one time $k$ .

$\text{Maximize}{\displaystyle \underset{j=1}{\overset{m}{\sum}}}{\displaystyle \underset{k=1}{\overset{n}{\sum}}}\mathrm{}R{V}_{j,k}\left(A\right)\cdot {X}_{j,k}$ (10)

$\text{subjectto}{\displaystyle \underset{j=1}{\overset{m}{\sum}}}{\displaystyle \underset{k=1}{\overset{n}{\sum}}}\mathrm{}G{C}_{j,k}\left(A\right)\cdot {X}_{j,k}\le {\epsilon}_{e}$ (11)

$\underset{k=1}{\overset{n}{\sum}}}{X}_{j,k}=1;\text{}\forall j=1,2,\cdots ,m$ (12)

${X}_{j,k}\in \left\{0,1\right\};\text{}\forall j=1,2,\cdots ,m;\text{}\forall k=1,2,\cdots ,n$ (13)

3.2. The Mathematical Model of Fired Sugarcane

The multi-objective mathematical model of fired sugarcane is described by Equations (14)-(17). The objective function Equation (14) maximizes revenue from the fired sugarcane selling. The constraint Equation (15) represents the second objective of the problem which is to minimize gathering cost of pro- duction, with upper bound given by ${\epsilon}_{e}$ . As for the fresh sugarcane model, the decision variables ${X}_{j,k}$ are defined by ${X}_{j,k}=1$ means that planted area $j$ is harvested at time $k$ and ${X}_{j,k}=0$ means that $j$ is not harvested at $k$ . In the constraint set, Equation (16) and Equation (17) ensure that in each area the sugarcane is harvested only at one time $k$ .

$\text{Maximize}{\displaystyle \underset{j=1}{\overset{m}{\sum}}}{\displaystyle \underset{k=1}{\overset{n}{\sum}}}\mathrm{}R{V}_{j,k}\left(B\right)\cdot {X}_{j,k}$ (14)

$\text{subjectto}{\displaystyle \underset{j=1}{\overset{m}{\sum}}}{\displaystyle \underset{k=1}{\overset{n}{\sum}}}\mathrm{}G{C}_{j,k}\left(B\right)\cdot {X}_{j,k}\le {\epsilon}_{e}$ (15)

$\underset{k=1}{\overset{n}{\sum}}}{X}_{j,k}=1;\text{}\forall j=1,2,\cdots ,m$ (16)

${X}_{j,k}\in \left\{0,1\right\};\text{}\forall j=1,2,\cdots ,m;\text{}\forall k=1,2,\cdots ,n$ (17)

Method of determining value of ${\epsilon}_{e}$ . The model will be solved for $p$ values of ${\epsilon}_{e}$ defined as follows:

${\epsilon}_{e}={\epsilon}_{e-1}+\Delta \epsilon ;\text{}e=1,2,\cdots ,p,$ (18)

where $\Delta \epsilon \mathrm{}=\mathrm{}\frac{UB-LB}{p-1}$ and $LB$ and $UB$ are lower and upper bounds defined

as follows.

#Math_93# (19)

3.3. Technologies and Equipment Used

The multi-objective mathematical models for fresh and fired sugarcane were solved using the linear programming and mixed-integer programming package GLPK (GNU Linear Programming Kit) version 4.52.

3.4. Data Used in Experiments

In this work, the optimal times to harvest sugarcane were determined for crop years 2012/13, 2013/14 and 2014/15. Examples of the real data used in the ex- periments for crop year 2012/13 are presented in Table 1 & Table 2, respectively

Table 1. The quantity of sugarcane ( $\times {10}^{5}$ tons) delivered into factories in four regions of Thailand for crop year 2012/13.

Table 2. The price of sugarcane and average total cost of production (baht/ton) in crop year 2012/13.

[6] . Table 1 shows quantity of fresh and fired sugarcane delivered to factories in the four regions of Thailand, Table 2 the price of the sugarcane and the average total cost of production, Table 3 shows the average C.C.S. of the sugarcane. The data in Table 3 were used to determine the values of the ${y}_{j,k}$ in the C.C.S. price ${P}_{2}\left(i\right)$ for the sugarcane in Equation (7).

A comparison of the quantities of fresh and fired sugarcane delivered to the sugar mills in the three crop years 2012/2013, 2013/2014 and 2014/2015 are shown in Figure 1. It can be seen that the percentages of fresh and fired sugar- cane delivered are approximately 30% and 70%, respectively.

4. Results and Discussion

The results of crop years 2012/13, 2013/14 and 2014/15 are presented in this section. The results are given for both fresh sugarcane and fired sugarcane for the four regions of Thailand.

4.1. Results for Crop Year 2012/13

The mathematical model for fresh sugarcane is described by objective function Equation (10) subject to constraints Equations (11)-(13) and the mathematical model for fired sugarcane is described by objective function Equation (14) subject to constraints Equations (15)-(17). The data from Tables 1-3 were used in the optimization. As stated earlier, the data in Table 3 were used to determine the values of the ${y}_{j,k}$ and the price ${P}_{2}\left(i\right)$ in Equation (7). The results are as follows.

1) Fresh Sugarcane :

The value of ${\epsilon}_{e}$ was calculated from the formula in Equation (18). The values of the lower bound $LB=1.46\times {10}^{8}$ baht and upper bound $UB=5.84\times {10}^{9}$ baht were estimated for each region. In the computations, 100 values of ${\epsilon}_{e}$ were used equally spaced between $LB$ and $UB$ in order to determine how the maximum revenue and the optimal cutting time changed as the upper bound on the gathering cost changed.

The results of revenue from the fresh sugarcane sales, gathering cost and the optimal harvest time for the four regions in Thailand; Northern, Central, Eastern and North-Eastern in each ${\epsilon}_{e}$ are presented in Table 4. It was found that the revenue increased at a rate greater than the increase of cost as shown in Figure 2(a).

For example, for ${\epsilon}_{2}$ , the revenue would be $0.199\times {10}^{9}$ baht, the gathering cost would be $0.196\times {10}^{9}$ baht and the optimal harvest times in each area would be as follows: 1) Northern should be harvested in 16-30 April 2013, 2)

Table 3. Average C.C.S. of sugarcane in crop year 2012/13.

Table 4. The optimal harvest time for fresh sugarcane in 2012/13 with $p=100$ ex- periments.

Figure 1. Comparison between quantities of fresh and fired sugarcane delivered to sugar mills.

Figure 2. Comparison between revenue from sugarcane sales and ga- thering cost of production for 2012/2013. (a) Fresh Sugarcane. (b) Fired Sugarcane.

Central should be harvested in 15-30 November 2012, 3) Eastern should be harvested in 1-15 December 2012 and 4) North-Eastern should be harvested in 1-16 May 2013.

2) Fired Sugarcane :

The values of the $LB$ and $UB$ for gathering cost of the fired sugarcane were estimated as described for the fresh sugarcane. With these values, 100 values of ${\epsilon}_{e}$ were used in the range ${\epsilon}_{e};4.69\times {10}^{8}\le {\epsilon}_{e}\le 9.87\times {10}^{9}$ . The cal- culated revenues from the fired sugarcane sales, gathering costs and the optimal harvest time for the four regions for each ${\epsilon}_{e}$ are shown in Table 5. It was again found that the revenue increased at a rate greater than the increase of cost as shown in Figure 2(b).

For example, for ${\epsilon}_{5}$ , It was found that the optimal harvest times in each area were as follows: 1) Northern should be harvested in 16-30 April 2013, 2) Central should be harvested in 15-30 November 2012, 3) Eastern should be harvested in 1-16 May 2013 and 4) North-Eastern should be harvested in 16-30 April 2013. The maximum values of revenue from the fired sugarcane selling would be $9.45\times {10}^{8}$ baht and the minimum values of gathering cost would be $8.05\times {10}^{8}$ baht.

4.2. Results for Crop Year 2013/14

The mathematical model for fresh sugarcane is described by objective function (10) subject to constraints (11)-(13) and the mathematical model of fired sugarcane is described by objective function (14) subject to constraints (15)-(17). Real data for crop year 2013/14 were used for the C.C.S. and the prices ${P}_{1}\left(i\right)$ for the weight, ${P}_{2}\left(i\right)$ for the C.C.S. and the gathering cost.

The results of revenues from the fresh and fired sugarcane sales, gathering costs and the optimal harvest time for the four regions in Thailand; Northern, Central, Eastern and North-Eastern for each ${\epsilon}_{e}$ are presented in Table 6 and Table 7, respectively. It was found that in both cases the revenue increased at a rate greater than the increase of cost as shown in Figure 3(a) and Figure 3(b).

4.3. Results for Crop Year 2014/15

The mathematical model for fresh sugarcane is described by objective function (10) subject to constraints (11)-(13) and the mathematical model for fired sugarcane is described by objective function (14) subject to constraints (15)-(17). Real data for crop year 2014/15 were used for the C.C.S. and the prices ${P}_{1}\left(i\right)$ for the weight, ${P}_{2}\left(i\right)$ for the C.C.S. and the gathering cost.

Table 5. The optimal harvest time for fired sugarcane in 2012/13 with $p=100$ ex- periments.

Table 6. The optimal harvest time for fresh sugarcane in 2013/14 with $p=100$ ex- periments.

Table 7. The optimal harvest time for fired sugarcane in 2013/14 with $p=100$ ex- periments.

1) Fresh Sugarcane.

The results of revenues from the fresh sugarcane sales, gathering costs, values of ${\epsilon}_{e}$ and the optimal harvest time for the four regions in Thailand: Northern, Central, Eastern and North-Eastern for each ${\epsilon}_{e}$ are presented in Table 8. As shown in Figure 4(a), it was found that for this crop year only 20% of the ${\epsilon}_{e}$ values gave revenue greater than the gathering cost.

2) Fired Sugarcane.

For this crop year it was found, as shown in Figure 4(b), that the revenue was less than the gathering cost for all ${\epsilon}_{e}$ . Therefore there was no optimal harvest time in this case.

4.4. Summary of Results

A summary of the optimal harvest times for fresh and fired sugarcane for the crop years 2012/13, 2013/14 and 2014/15 are shown in Table 9 and Table 10.

Figure 3. Comparison between revenue from sugarcane sales and gathering cost of production for 2013/2014. (a) Fresh Sugarcane; (b) Fired Sugarcane.

5. Conclusions

In this paper, a multi-objective mathematical model has been presented for the calculation of the optimal times to harvest fresh and fired sugarcane in Thailand for the crop years 2012/13, 2013/14 and 2014/15. The two objectives in the model were to maximize revenue and minimize cost. The $\epsilon $ -constraint method was used to change the multi-objective optimization model into a single-ob- jective model to maximize revenue subject to an $\epsilon $ -constraint of an upper limit on the cost. The GNU Linear Programming Kit version 4.52 was used to compute

Figure 4. Comparison between revenue from sugarcane sales and gathering cost of production for 2014/2015. (a) Fresh Sugarcane; (b) Fired Sugarcane.

the optimal harvest time, maximum revenue and minimum gathering cost for a range of 100 values of $\epsilon $ -constraint. The results of the computations are sum- marized in Table 9 and Table 10.

The results show that the methods used in this paper can be used to predict optimal harvesting times. In this model, the main factor in determining the optimal harvest time is the change in C.C.S. with time as shown in Table 3. However, there are a number of other factors that are not included in the model that could change the optimal harvest times. These factors include: 1) changes in fixed and variable production costs during a year, 2) changes in delivery costs

Table 8. The optimal harvest time for fresh sugarcane in 2014/15 with $p=100$ ex- periments.

Table 9. Optimal harvesting times for fresh sugarcane.

Table 10. Optimal harvesting times for fired sugarcane.

during a year, 3) decisions by the sugar mills on their requirements for sugar cane.

Acknowledgements

This research was partially funded by King Mongkut’s University of Technology North Bangkok (Contract No. KMUTNB-GOV-59-23) and Centre of Excellence in Mathematics, The Commission on Higher Education, Thailand.

Conflicts of Interest

The authors declare no conflicts of interest.

Cite this paper

*Applied Mathematics*,

**8**, 329-343. doi: 10.4236/am.2017.83028.

[1] | Naranong, V., et al. (2013) A Study of Reform the Structure of Thailand’s Sugar and Cane Industry. Thailand Development Research Institute. |

[2] |
Florentino, H.O., Moreno, E.V. and Sartori, M.P. (2008) Multiobjective Optimization of Economic Balances of Sugarcane Harvest Biomass. Science Agricultural, 65, 561-564. https://doi.org/10.1590/s0103-90162008000500018 |

[3] | Scarpari, M.S. and Ferreira, E.G. (2010) Optimized Agricultural Planning of Sugarcane Using Linear Programming. Revista Investigacion Operacional, 31, 126-132. |

[4] | Gomes, F.R.A. (2012) Bi-Objective Mathematical Model for Choosing Sugarcane Varieties with Harvest Residual Biomass in Energy Cogeneration. International Journal of Agricultural and Biological Engineering, 5, 50-58. |

[5] | Deb, K. (2001) Multi-Objective Optimization Using Evolutionary Algorithms. John Wiley and Sons, Chichester. |

[6] |
Office of the Cane and Sugar Board. http://www.ocsb.go.th |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.