In this paper, a simple new approach for the grassroots and retrofit design of water utilization systems with multiple contaminants is presented. This approach uses a combination of mathematical programming and necessary conditions of optimality to automatically generate the optimal solution featuring minimum capital and operating costs. This paper presents the only existing method to solve this problem that can guarantee global optimality.