This paper describes the general and efficient design method of an axial and a radial shim coils to correct field impurities of various harmonic orders in the imaging volume of Magnetic Resonance Imaging magnet. Shim coils are optimized by BCLSF subroutine of IMSL, which is the well-known commercial package for optimization, aiming at maximizing the magnitude of the desired field component as well as minimizing other field components. In order to evaluate their effect, the developed method was applied to the MRI magnet constructed in KERI.