An effective optimum design method is presented. The method combines the three-dimensional boundary element method with the (1+1) evolution strategy. The reduced scalar potential formulation with the magnetic surface charge as the unknown variable is used under the assumption that the yokes are not saturated. It is found, through the numerical example, that the global optimum shape of the magnet can be easily found within a reasonable number of generations. More elaborate design can be achieved by increasing the number of design variables.