In this paper, we formulate the optimization problem of multi-output Reed-Muller (RM) circuits and suggest a non-exhaustive search method based on the genetic algorithm (GA) to solve it. Simple as well as randomly generated multi-output logic functions are used to test and compare the performance of GA and exhaustive-search methods. It is shown that the GA is well suited for this application, especially when the number of variables is more than ten.