In many response surface optimization problems for industrial processes, there are more than two responses of interest, and we want to find the optimal levels of the factors that influence the responses. This paper is to propose how to set up the desirability functions to find the optimum for a given set of data, and to propose how to analyse the data and the desirability functions to determine an optimal operating condition for the factors. To implement the proposed method in practice, a FORTRAN computer program was written and explained. Finally, an industrial example is illustrated to explain the proposed technique and the source list of the computer program is attached for the users.