A perturbation approach is used to derive surface impedance boundary conditions (SIBCs) for nonlinear magnetic conductors described by means of Brauer's model [1]. A new technique is proposed to solve the nonlinear diffusion equation of the magnetic field in conductive media and then its solution is coupled with a BEM formulation for the case of parallel magnetic conductors in terms of E and B fields.