The flux density has a linear relationship with the magnetic field intensity when the flux density of a magnetic material is in the millitesla (mT) region. However in the high-frequency region, the magnetic flux is concentrated near the surface of the magnetic material owing to the skin effect. Therefore, it is necessary to consider the nonlinear magnetic properties of the material. In this paper, we propose a modified frequency response analysis method that takes into account the nonlinear distribution of the complex magnetic permeability with respect to the magnetic flux density, in order to accurately compute the magnetic flux distribution in the high-frequency range. The validity of the proposed analysis method is verified by comparing the computed data with measurements obtained by using a prototype (toroidal coil).