This paper presents the shape optimization to minimize the BEMF(Back Electro-Motive Force) harmonics of PM type synchronous generators. RSM(Response Surface Methodology) is well adapted to make analytical model for a complex problem considering a lot of interaction of design variables. In this paper, RSM is used to find the optimal solution. The 2D-Finite Element Method is used to obtain the observer data of the BEMF and SQP(Sequential Quadratic Problem method) is used to solve the constrained nonlinear optimization problem.