Processes for optimizing the coil shape of deflection yoke are proposed A very accurate and practical winding modeler is developed and volume integral equation method (VIEM) is used for field calculation. Two steps of optimizations are done by using (1+1) evolution strategy. Those are dimensional optimization and pin-position optimization Various techniques are applied for reducing computational time for the optimization.