A shimming method for correcting inhomogeneity of a static magnetic field generated by a magnet of a nuclear magnetic resonance imaging machine includes: measuring the magnetic field at a plurality of points over a reference surface; generating a polynomial that solves Laplace's equation with boundary conditions given on the reference surface, the polynomial representing the magnetic field on the reference surface and having a plurality of harmonic terms, each associated with a coefficient; determining the coefficients from the field sampling values; defining a grid for positioning a plurality of correction elements and relating it to the field structure; and calculating the position and magnitude parameters of the correction elements, such that the correction elements affect the coefficients of the magnetic field to obtain the desired field characteristics, wherein the reference surface is a superquadric surface, such that the magnetic field is corrected in a volume delimited by the superquadric surface.