The shimming work assistance unit performs singular value decomposition of a response matrix, which represents the relationship between an error magnetic field distribution and an adjusted magnetic moment placement distribution. From the multiple eigenmodes obtained, the eigenmodes are selected and added one by one in order from the eigenmode with the highest singular value, and the residual magnetic field error, which represents the fluctuation range of the difference between the magnetic field distribution, generated by the placement of the shimming magnetic moments corresponding to said eigenmode, and the error magnetic field distribution, is displayed on a display unit as a function graph of eigenmode order (line graph (1)).