A shape optimization method has been developed regarding magnetic materials in 2-dimensional geometry. The method combines (i) complex expressions of magnetic fields. (ii) computation of complex magnetizations based on integral equations and (iii) the Quasi-Newton method. Details of the computational procedure and calculation results are reported.