An Interactive Multiobjective Optimization Approach to Supplier Selection and Order Allocation Problems Using the Concept of Desirability

In supply chain management, selecting the right supplier is one of the most important decision-making processes for improving corporate competitiveness. In particular, when a buyer considers selecting multiple suppliers, one should consider the issue of order allocation with supplier selection. In this article, an interactive multiobjective optimization approach is proposed for the supplier selection and order allocation problem. Also, the concept of desirability is incorporated into the optimization model to take into account the principles of diminishing marginal utility. The results are presented by comparing them with the solutions from the weighting methods. This study shows the advantage of the proposed method in that the decision-maker directly checks the degree of desirability and learns his/her preference structure through improved solutions.


Introduction
With the development of the Internet and network technologies, the digital divide between businesses is being addressed.In this situation, companies are using strategies to enhance outsourcing to focus on core competences along with information.In particular, such rapid developments in information technology are pushing many companies into a race that transcends time and space.Furthermore, the global business environment surrounding corporations is changing from competition between individual companies to competition between supply chains.In order to ensure the competitive edge of the supply chain, companies have tried to find the right suppliers offering higher quality, reduced costs, and shorter lead times.Therefore, in supply chain management, selecting the right suppliers is one of the most important decision-making processes for improving corporate competitiveness [1][2][3].
In many cases, a single supplier may not be able to meet the buyer's requirements.In such cases, selecting multiple suppliers-multiple sourcing-would be a reasonable alternative [4].Whereas single sourcing significantly increases the disruption risk in the supply chain, multiple sourcing increases the fixed cost in terms of administrative and negotiating costs [5].However, multiple sourcing is preferred over a single sourcing, ensuring order flexibility [6].Therefore, multiple sourcing inevitably includes the problem of order allocation.Furthermore, the relationship between the buyer and supplier is influenced by order allocation decisions based on strategic purchase decisions [7].Therefore, the overall supplier selection problem should not only cover the selection of the right supplier but also the determination of the orders assigned to the selected supplier based on the given objectives and constraints [8].To date, however, only a few mathematical programming models to analyze such decisions have been published [9][10][11][12].In the field of supply chain management, it is necessary to develop a precise decision support model which simultaneously considers supplier selection and order allocation.In particular, researchers need to advance research in this direction to guide the decision-maker (DM)'s rational choice.
Supplier selection is inherently a multi-criteria decision-making (MCDM) problem since some conflicting performance criteria have an influence on the selection of suppliers [13,14].It is also considered one of the most familiar problems in MCDM [15].This problem has been studied from a variety of perspectives, such as green supplier selection [6,16] and global supplier selection [17].Generally, the performance criteria involves the factors such as cost, quality, and lead time.Each supplier has its own strengths and weaknesses, so it is very difficult to select a superior supplier at all dimensions of criteria.Thus, supplier selection problem has made MCDM a very challenging task.Ho et al. [18] provided comprehensive reviews for MCDM approaches to vendor selection problems between 2000 and 2008.In addition, a more recent review from Chai et al. [19] provides guidelines for a MCDM-based supplier selection model.
For solving a multiobjective optimization (MOO) problem in MCDM, the DM seeks a compromise solution which provides the greatest satisfaction in the presence of conflicting objectives.Thus, the DM's preference information plays a critical role in finding the solution.In the literature, the DM's preference information can enter the solving process of MOO problems in three different ways: (1) a priori; (2) a posteriori; and (3) progressive (interactive) articulation [20,21].For an a priori setting, multiple objective functions combined with preference convert into one single objective.For an a posteriori setting, the DM's preference information is articulated after optimization process by selecting the most preferred one from a set of non-dominated solutions, usually called a Pareto optimal set.For progressive optimization, the preference of the DM is incorporated into the solution search process.Iterative dialogues between the DM and optimization model contribute to find the most satisfactory solution in this optimization approach.For this reason, the progressive articulation approach is also referred to as the interactive approach.
Most studies in supplier selection and order allocation are categorized into the prior approach by the MOO categorization scheme.However, the prior approach may decrease the reliability of solutions because of the unrealistic assumption that the DM can specify the preference information in advance.Although the interactive approach has been broadly applied to various fields as an alternative to overcome this limitation on the DM's preference, it is rarely used to solve supplier selection and order allocation problems.To the best of the author's knowledge, Demirtas and Üstün [4]'s research is the only material related to using the interactive MOO method to solve supplier selection and order allocation problem.Demirtas and Üstün [4] proposed an interactive MOO method, and a reservation level driven Tchebycheff procedure.In this model, the solution process makes the DM express adjusting some or all of reservation values of objective functions, by generating candidate solutions with sampling weights.The interactive method that we propose in this paper has something in common with Demirtas and Üstün [4]'s research in terms of a kind of objective space reduction.Also, two methods share the Tchebycheff framework to find a solution by utilizing the concept of distance to the ideal vector.However, the proposed approach provides an integrated use of the desirability function approach in response to surface methodology and the step method (STEM).Our model has advantages in that it reflects the satisfaction of the DM more realistically by using the concept of desirability function, and it reduces the burden on the DM to express preference information.
Based on the above-mentioned background, the purpose of the research is summarized as follows.We aim to solve the multiple sourcing problem which deals with order allocation at the same time as supplier selection by using an interactive MOO method.The devised method, which progressively articulates the DM's preference information, is applied to a problem with three important criteria: cost, quality, and delivery.Also, to intuitively utilize the level of satisfaction, we borrow the concept of desirability from the research field of product and process design.We show that our method can be utilized effectively in the supplier selection and order allocation problem.
The organization of this paper is as follows.Section 2 presents previous works to be addressed in this study.The proposed interactive desirability function approach to supplier selection is presented in Section 3. Section 4 analyzes results and provides discussions.Conclusions are provided in Section 5.

Literature Review
This section briefly describes the concept of desirability function that is the basis for the proposed methodology.Then we deal with the consideration of the DM's preference information in MOO problems.Also, the STEM is introduced.

Desirability Function Approach
The desirability function approach is in fact designed to solve the problem of product and process design (also called multiple response optimization).The desirability function approach, proposed by Harrington [22] and Derringer and Suich [23], is one of the most widely used methods in product and process design.This approach transforms an estimated response function into a scale-free function, called the desirability function, which ranges from zero to one.Thus, the value of the desirability function presents the degree of desirability or satisfaction for the corresponding response.If the larger-the-better (LTB) response is the objective that has to be maximized, the individual desirability function is defined as where ŷmin i is the minimum and ŷmax i is the maximum acceptable value of response i and r is the parameter (r > 0) that determines the shape of desirability functions.The function is convex if r > 1, and is concave if 0 < r < 1.If a response is of the smaller-the-better (STB) type, it is defined in a similar way to (1).Also, for a certain response, if the specified target value has to be attained, the response is referred to as the nominal-the-best (NTB) type.In this case, the individual desirability function is derived in a slightly different form from STB and LTB type.For more information, readers may refer to Derringer and Suich [23].
There are several ways to draw the optimal solution considering multiple individual desirability functions [22][23][24][25][26] but the most widely used method is to optimize by converting multiple desirability functions into single desirability.In most cases the concept of geometric means is often used for this aggregation.

Interactive Approach and Decision Maker's Preference Information
The supplier selection is inherently a MCDM problem, since some conflicting criteria have influence on the selection of suppliers [13,14].MOO is a mathematical optimization technique of MCDM, which involves more than one objective function to be optimized simultaneously.The ultimate goal of MOO is to find the best compromise solution from a set of the non-dominated solutions that satisfies multiple objectives simultaneously.Thus, it necessarily requires the involvement of the DM who provides the preference information among conflicting objectives.According to the timing of articulating the DM's preference information, the MOO is classified into three categories: the prior approach, the posterior approach, and the interactive approach [20,21].More detailed descriptions of the theories and applications of MOO studies can be found in the classical MOO textbooks: [20,27,28].
The interactive approach (also referred as the progressive approach) allows the DM to articulate his/her preference information progressively while solving the problem.The optimization process is repeated until the DM has found the most preferred solution to extracting preference information in an interactive manner.Specifically, at each iteration in an algorithm, the DM is asked to express some preference information given one or more solutions generated from the previous iteration.This information reflects the judgment based on the DM's implicit value function, and the process updates the DM's preference through adding the new information to an optimization model.Typically, a general interactive method has the following three steps [27].

•
Step 1: Find an initial solution.

•
Step 2: Interact with the DM.

•
Step 3: Generate one or some new solution(s).If the DM is satisfied with the current solution, stop.Otherwise, go to step 2.
Over the years, a number of various interactive methods have been developed.The most well-known interactive methods are the step method (STEM) [29], the Geoffrion-Dyer-Feinberg (GDF) method [30], the Zionts-Wallenius (ZW) method [31], the reference direction approach (RDA) [32], and the Nondifferentiable Interactive Multiobjective BUndle-based optimization System (NIMBUS) [27].Each of interactive methods has a different scheme from others in terms of which type of preference information is asked from the DM.No matter what kind of preference information is asked, an iterative process is valuable for the DM in that he/she can learn about the structure of the problem and expect the result of the interaction.Furthermore, each employs different computational algorithms.Therefore, a particular interactive method may be suitable for a specified MOO problem given mathematical properties or assumptions.In other words, there is no the sole interactive method that is decidedly superior to the others, in a universal sense.Even though a lot of extensions and variants of the interactive methods can be applied to the supplier selection problem, in this paper we propose an interactive approach based on the STEM with some modifications.

STEM
The STEM is proposed by Benayoun et al. [29] and it is known as the first interactive MOO method.Originally STEM is designed for the linear MOO problems, and it also can be applied to integer and nonlinear MOO problems [28].The interactive process in the STEM begins with generating an initial feasible solution and identifying the ideal objective function values.The ideal objective function values are obtained by optimizing each individual objective function over the initial feasible region.At each iteration, the DM examines the current solution with the ideal objective function value, and presents his/her preference depending on the concept of level of satisfaction.In particular, the DM provides an amount of relaxation by which the objective function value can be sacrificed, in order to improve some other unsatisfied objectives.From this numerical information presented from the DM, it generates a new solution.If the DM satisfied with this new solution, this solution becomes a final solution as a compromise and the interactive process terminates.Otherwise this process is repeated until no further relaxation is accepted.The mathematical procedure for the STEM has been given in the following steps.
Step 0: Construct the pay-off table.In this step, the ideal objective function vector is obtained by maximizing each objective function individually.Let f * be the ideal solution vector of the following K problems, Step 1: Let h be iteration counter and set be zero (h = 0).Calculate π i values for use in weighting the objectives: where c ij are the coefficients of the ith objective.The first term is to place the most weight on the objectives with the greatest relative ranges.The second term normalizes the gradients of the objective functions.
Step 2: Let h = h + 1, and S h+1 = S. S 1 = S means the solution process begin with the original feasible region.At this point, the relaxation index set be null (J = ∅).Compute relative weights: Step 3: Solve the weighted minimax program: If all objective function values are satisfactory, then stop; the current solution may become the final solution.
Step 4: Determine satisfactory objective and then specify the amount of relaxation, ∆ i (i ∈ J).
Step 5: Form reduced feasible region: where h + 1 is increased iteration counter.Thus, S h+1 is the feasible region in the (h + 1)th iteration.f i (x) h is the objective function value in the hth iteration and ∆ i is the amount of relaxation of particular objective function.Go to Step 2 with the reduced feasible region.
The STEM is widely used to solve not only linear but also nonlinear MOO problems for the reason that it is simple and straightforward to generate the solution by reducing feasible regions in the objective space.Despite these advantages, the STEM shows two major drawbacks.First, the STEM does not take into account different degrees of satisfaction within the acceptable interval of a relaxed objective.For example, suppose that one objective f e is determined by the DM as a satisfactory function, and that he/she decides to relax as much as the ∆ e .In this case, the weight of objective function e becomes zero since e / ∈ J, so this function does not serve as an objective function in the next iteration.Therefore, the solution of the next iteration depends on the assumption that the satisfaction of the DM is indifferent within the interval of [f e − ∆ e , f e * ].This disadvantage has been pointed out by Jeong and Kim [33].The second drawback is in the dialogue scheme between the DM and the optimization model.The dialogue is made mainly by asking the DM to acceptable amount of relaxation of a satisfactory objective.However, it increases the burden on the DM because he/she must provide specific numerical information at each iteration.In this paper, we intend to utilize the STEM to supplier selection problem by modifying it in such a way that we enjoy its merits and complement the shortcomings mentioned above.More specifically, we try to overcome the first limitation by applying the concept of desirability function and the second one by proposing a modification to the DM's preference information representation.

Proposed Model
Across fifty years of evaluating suppliers, many researchers have proposed different sets of criteria.
The first set of criteria was proposed by Dickson [14], who identified 23 different criteria evaluated in supplier selection.Evans [34] and Shipley [35] agreed that price, quality and delivery are the most important criteria for evaluating suppliers.Ellram [36] proposed that the quality dimension should be divided into product quality and service quality, and it is suggested to use them with price and delivery time to select suppliers.Weber et al. [37] surveyed based on Dickson's 23 criteria and concluded that price, delivery, quality, production capacity, and localization are the most important criteria.Pi and Low [38] proposed quality, delivery, price and service for supplier evaluation.Amid et al. [39] uses price, quality, and service, and Jadidi et al. [11] utilized price, quality, and lead time to the supplier selection and order allocation model.As shown by the literature, the most important criteria for supplier selection problems are cost, quality, and delivery.We also use those criteria by using the measures: total purchasing costs, the number of rejects, and the number of late delivery for cost, quality, and delivery, respectively.We assume the buyer considers a single item that should be purchased under known total demand.Also, information about criteria and production capacity is already known for a set of potential suppliers.The notations of the proposed model are presented in Table 1.There are three original objective functions: f 1 (minimizing total purchasing cost), f 2 (minimizing total number of defects), f 3 (minimizing the number of late delivery).Each objective function can be shown as follows: Using the functions defined above, MOO problem considering the supplier selection with order allocation can be formulated as follows [12]: The first and second constraint present demand satisfaction and capacity restriction, respectively, and the last constraint ensures non-negativity of the decision variables.Now we transform the original objective functions into the desirability functions.Because all original objective functions are to be minimized (STB type), the individual desirability is defined as where f i max are the maximum values of the objective functions as obtained from the payoff table and f i * (x) are the ideal points from maximizing the objective functions individually.We represent the desirability-based supplier selection problem in the form of a MOO problem: The individual desirability functions d i (i = 1, . . ., k) are objective functions to be maximized simultaneously.
The optimization process consists of five major steps.The overall procedure is presented in Figure 1.In addition, a pseudo code is described in Algorithm 1 to aid readers understanding.

Case r = 1
We apply a modified STEM to supplier selection problem.In order to present how the method works, we deal with the linear type desirability function transformation.
Step 1: Calculate π i values for use in weighting the objectives: The coefficients c ij , which reflect the gradients of the objective functions in the weight calculation is the coefficient of desirability function, not the original objective function.If the desirability functions are nonlinear, i.e., r = 1, the weights calculation is somewhat different from that of the linear case.As mentioned in Section 2.2, originally the STEM is developed for multiobjective linear problems, but nonlinear extensions have been proposed, for example, in Vanderpooten and Vincke [40], Eschenauer et al. [41], and Sawaragi et al. [42].A weight calculation scheme for the nonlinear case (r = 2) is presented in following Section 3.2.At this point, calculated π i values are shown as Step 2: Let h be iteration counter (set h as zero for the first iteration).For h = 0, the relaxation index set be null (J = ∅).
Step 4: Determine satisfactory objective by asking the DM to express satisfactory objectives.If h = 0, go to Step 2 by setting h = h + 1.We assume the DM satisfied with d 2 , so relax d 2 .
Step 4: Present new desirability function vectors to the DM as shown in Table 4, who is asked to select an acceptable rate of relaxation, δ.Unlike original STEM, the proposed approach does not compel the DM to ask specified amount of relaxation from satisfactory level.Then, the optimal solution corresponding selected rate of relaxation becomes the new vector of desirability function values.
Step 4: As shown in Table 5, present new desirability function vectors as δ increases, then ask to the DM an acceptable rate of relaxation, δ. ).This interactive process can be repeated until the final solution satisfies the DM.At this stage, we assume that the DM is satisfied with the current solution and terminate the iterative procedure.

Case r = 2
All three original objective functions considered in the supplier selection model are linear, and the defined desirability functions in Section 3.1 are also linear.This sub-section discusses the need for nonlinear desirability functions and presents a problem-solving process in r = 2 case.
The weighted sum method assumes that the DM's overall desirability (utility) is determined by the weighted sum of multiple objective functions.Thus, the marginal contribution to utility of each linear objective function is constant (see Figure 2a).Therefore, the optimization result is highly dependent upon the predetermined weights, although the objective functions are normalized by using the ideal points and the nadir points.On the other hand, the STEM modifies the concept of satisfaction by accepting the DM's preference information in the problem-solving process.In particular, the STEM assumes that the desirability of the DM is the same up to the acceptable level, i.e., f i (x) + ∆ i for a STB type objective function (see Figure 2b).Namely, the STEM considers that the DM's satisfaction is indifferent within the interval of [f i * , f i (x) + ∆ i ].With this concept the STEM improves the other objective functions at a level that does not impair the satisfaction of the satisfactory objective.However, we need to consider the principles of diminishing marginal utility, overlooked in the STEM.The principles of diminishing marginal utility explain that the closer the objective function value is to the optimum value, the less contribution it will make.In this study, the desirability function enables the proposed model to consider the principles of diminishing marginal utility.If the parameter r of desirability functions is not equal to 1, the functions become nonlinear.The parameter r exhibits a diminishing marginal contribution to the maximum cumulative desirability.The nonlinearity of desirability function provides a clue to successfully deal with the principles of diminishing marginal utility in the proposed model.This characteristic is also highlighted in Jeong and Kim [33]'s STEM based on the desirability function.Figure 2c shows an example of a nonlinear desirability function.
Step 1: Calculate π i values for use in weighting the objectives.The several weight schemes used to extend the STEM to nonlinear MOO problems were suggested by various researchers (Eschenauer et al. [41]; Sawaragi et al. [42]; Vanderpooten and Vincke [40]).In this study, we utilize Vanderpooten and Vincke [40]'s calculation (assume that the denominators are not equal to zero) as follows: By (14), the determined weights are π 1 = π 2 = 1 and π 3 = 0.986.
Step 4: We assume the DM satisfied with d 2 , so relax d 2 .Go to Step 2.
Step 3 and Step 4: Solve the weighted minimax program and present new desirability function vectors as δ increases, then ask to the DM an acceptable rate of relaxation, δ.Use the computation results as shown in Table 7.

Summary of Results and Discussion
This study successfully applied the concept of desirability function in response surface methodology, which is used in product and process design, to supplier selection problems.The desirability function diverts the DM's recognition system for each objective function from linear coupling and realistically reflects the degree of satisfaction for each objective function.Since the proposed method assumes that the DM's preference information is not completely known, it is difficult to discuss the superiority of solutions by directly comparing the results with other methods.However, we explain the advantages of the proposed method by comparing them with the solutions from two weighting methods: the weighted sum method and weighted geometric method.Several weighting vectors are assumed for the purpose of comparing solutions.We adopted three sets of weighting parameter w, proposed be Jadidi et al. [11] for the same supplier selection problem.The results are shown in Table 8.
The results of the weighted sum method show that the more weights are assigned to the first objective function (minimizing total purchasing costs), the larger the contribution is made to the value of the first desirability function.However, when w 1 is larger than 0.6, the individual desirability value for f 2 equals to 0. Namely, the second objective function, the number of late delivery, has the worst value.Furthermore, there is no differences between w = (0.6, 0.2, 0.2) and w = (0.8, 0.1, 0.1), although the overall objective function values ∑w i d i differ.These extreme results show that the DM may not satisfied with the results.Next, we also tested the weighted geometric model using the same set of w.In fact, the weighted geometric mean is a popular method for unifying individual desirability functions to a single function in desirability function approach.The weighted geometric method results in a somewhat balanced solution avoiding extreme values in one objective function, even if the weight is to one side.This method, however, also may result in a controversial solution.In Table 8, we found the third desirability value for late delivery increases as w 3 decreases.Thus, the results do not necessarily guarantee that the DM can find the most satisfactory solution even if he/she can decide preference information in advance at one time.The proposed method can complement the shortcomings of these methods because the DM directly checks the degree of desirability and learns a preference structure through improved solutions.In other words, the advantage of the proposed method is that the solution changes in the object function space can be detected through the DM's preference information.
The last row section in Table 8 summarizes the solutions from the proposed method.The results show that the proposed method prevents the emergence of extreme values in two ways.First, the initial solution d 0 describes this characteristic.In each case, for r, the initial solution is a neutral compromise solution because it is calculated by using the Tchebycheff metric without preference information.This setting eliminates the unnecessary iteration that causes a desirability function value to deviate from an extreme value.The results in the 12th and 15th row show reasonable levels of desirability functions, for both cases, r = 1 and r = 2. Second, the proposed method also prevents the extreme value in the problem-solving process.The original STEM allows the DM to relax one objective function, which may lead to significantly improved values than expected.In such a case, the DM might want to relax an excessively improved value for the purpose of adjustment.However, the proposed method does not ask the DM directly for the amount of relaxation, so it induces the satisfactory solution to find the change of desirability functions according to the rate of relaxation.Accordingly, the DM can choose the most satisfactory solution among the various possible solutions.Therefore, it helps to avoid occurrence of an extreme value by showing the alternatives that the DM can choose.
The proposed method presents the changes of the desirability function values to the DM according to the rate of relaxation of a particular desirability function in the form of a table.If a MOO method requires preference information that is difficult for the DM to express-for example, specific numerical information-it may be difficult to find the satisfactory solution.In this regard, the information presentation of the proposed method provides an additional advantage in that it can ease the burden of the DM by showing a set of expected candidate solutions and selecting the most satisfactory solution rather than requiring specific values.If the number of desirability functions increases, it is recommended to use the graph form.The graph form may also be more useful because the scales of the desirability are all the same.An alternative interactive MOO approach to supplier selection and order allocation is Demirtas and Üstün [4].In this interactive method, the DM can control objective values directly in a similar way to our research.The authors uses the reservation level that represents an objective function value which must be equaled or exceeded to be considered acceptable, in the maximization context.The method repeatedly reduces the objective space by adjusting reservations level from the DM's preference information.Although there are differences in, for example, how local weights are used and how the DM expresses preference information, both studies show that interactive MOO approaches can help supplier selection and order allocation problems.We expect that a variety of interactive MOO methods can be used for supplier selection and order allocation problems.

Conclusions
In this study, a new interactive MOO approach is proposed for supplier selection and order allocation problems with the concept of desirability.The multiple objectives for this problem are defined as purchasing costs, quality, and lead time.Then, we applied an interactive MOO method based on STEM to this problem.We presented how to solve the problem step by step to deliver the problem-solving procedure.It is also shown that the principles of diminishing marginal utility can be used to determine the level of the satisfaction through the use of desirability function.Two cases involving linear and nonlinear desirability function are described to explain the detailed procedure of the proposed method.The obtained results are presented along with results from a prioiri methods and have been discussed accordingly.
This study has following salient features that contribute to the research stream on supplier selection and order allocation.

1.
The use of desirability concept was shown to be an excellent tool for reflecting the level of satisfaction and has advantages in that the sensitivity of satisfaction can be adjusted by using desirability parameter.

2.
The proposed method alleviates the appearance of extreme values that can be derived when the model uses pre-determined weights.

3.
The interactive MOO method, which progressively articulates the DM's preference information to deal with supplier selection problem, may be an alternative.4.
Modification of the original STEM allowed us to relieve the DM's burden in terms of presenting preference information and reduce unnecessary iteration in the optimization process.
This study tried to deal with the concept of desirability more intuitively by repeatedly articulating the preference information of the DM.However, our model will also require the parameter value of the desirability function, in advance, to determine its shape.Therefore, it is necessary to develop a method that progressively reflects the information of this parameter, since the shape of the desirability function plays a critical role in deriving the solution.It is possible to complement this shortcoming if a method is developed that can estimate this parameter for each objective function utilized for the supplier selection and order allocation problems.Also, as stated in Section 4, we could not directly compare the solution quality of the developed method with solutions of other methods.Nonetheless, in order to compare the solution quality, we suggest that setting the proxy on the premise that the DM's preference information is known and comparing the posteriorly derived results is an alternative that can also improve the method's validation.Finally, developing an interactive method under incomplete information and applying it to supplier selection problems would also be a challenging work for practical application.Recent research considering incomplete weights presented by Liao and Xu [43] is expected to be useful in these respects.

Figure 1 .
Figure 1.The overall procedure of the proposed method.
Table 2 presents the payoff table that includes f i max and f i * (x).For a STB type function, f i max are recognized as the nadir points which are maximum values from each column in the payoff table, and the ideal values are on the diagonal of the table.

Algorithm 1 .
Pseudo code of the proposed method.

Table 4 .
Changes of desirability function values according to the rate of relaxation (h = 1).If the DM is satisfied with the new solution with selected rate of relaxation, δ*, then stop.Otherwise, set h = h + 1 and set the new feasible reason S h+1 .Go to Step 2. Assume the DM choose 20% relaxation of d 2 .Therefore, a new solution is d 1* = (0.589, 0.460, 0.615) with f = (61.167,5.083, 3.582).Also, we assume the DM satisfied with d 1 , and relax d 1 .

Table 5 .
Changes of desirability function values according to the rate of relaxation (h = 2).

Table 7 .
Changes of desirability function values according to the rate of relaxation (h = 1).

Table 8 .
Comparison with some weighting methods.
1A special case of weighted geometric mean with equal weights.