A three-phase method for systematically and simultaneously computing multiple optimal solutions