The paper presents an algorithm for solving the equations simulating the steady-state behavior of the multi-stage flash(MSF) desalination process. The model can be used either for designing new plants or analysis and optimization of existing units. The set of equations relating the large number of design and operating variables for any stage (i-th stage) in MSF plants is broken down into three subsets. The first subset contains equations describing processes in the tubes of the preheater inside the flashing chamber. The second subset deals with the processes inside the flashing chamber. These two subsets of equations are subsequently solved iteratively involving only a single variable for each subset. These variables are: (1) saturation temperature T-vi of the flashed-off vapor and (2) temperature of the unevaporated brine flowing from the chamber. The third subset of equations considers the existing interactions between the first two subsets. This subset defines the mass of vapor formed by flashing D-i in the i-th stage. All equations in the above-mentioned subsets are solved by a reliable and efficient one-dimensional fixed-point iteration. The main advantages of this method are less sensitivity to initial guesses, fewer iterations to obtain the required solution, and no need for calculating derivatives. The algorithm is implemented using the computer-aided design (CAD) interative package L-A-S (Linear Algebra and Systems). Detailed results are presented to show the dependence of the important factors controlling the fresh water cost, which are plant performance ratio, specific heat transfer area, specific brine flow rate, and specific cooling water now rate, on the most significant two design variables, namely the total number of stages and the top brine temperature. The predicted data hem the model are compared with published data of a typical MSF plant in operation at Kuwait. The agreement was found to be very good.