Aim of this paper is to compare the features of integral and differential formulations for the solution of transient eddy current problems in nonlinear media. The integral method used in this paper reduces the nonlinear eddy current problem to the analysis of an equivalent network. Two dual edge clement differential formulations are also used. Their field estimates automatically verify Faraday's and Ampere's laws and furnish the distribution of the constitutive error in solution domain, which provides a useful tool for the refinement of the discretization. The methods are applied to analyze a massive magnetic circuit characterized by the presence of soft and hard magnetic materials with narrow gaps.