In the present work, a modified mixed-integer linear programming model was implemented in Microsoft Excel® and minimized using the Solver tool to obtain information to devise a heat exchanger network with multiple utilities from a set of hot and cold streams and selected utilities by hand. Regarding the mixed-integer linear programming problem, the summation of utility energy was added to the model, and this energy was equal to that from the Temperature Interval method and the Grand Composite Curve. Moreover, feasible temperature ranges for heat exchange were considered according to the second law of thermodynamics. Also, the last two temperature intervals from the rank-ordered ones were assigned for water energy balances. An energy balance was introduced into the algorithm for each interval between its temperature and the process pinch temperature in the case of the boiler-feed water. Seven stream sets collected from the literature were used for the mixed-integer linear programming formulation testing, and six of them are presented in this article. Because of boiler-feed water generation and the low cost of utilities, the annualized cost of a heat exchanger network with multiple utilities can be lower than that of a network without multiple utilities.