To calculate the position and number of passive shimming ferromagnetic elements, a simulation algorithm is designed for magnetic resonance imaging(MRI)permanent magnet system.