Three‐dimensional MR reconstruction of high‐contrast magnetic susceptibility by the variational born iterative method based on the magnetic field volume integral equation