The most significant assumptions in the subdomain technique (i.e., based on the formal resolution of Maxwell's equations applied in subdomain) is defined by: "The iron parts (i.e., the teeth and the back-iron) are considered to be infinitely permeable, i.e., µ iron → +∞, so that the saturation effect is neglected". In this paper, the author presents a new scientific contribution on improving of this method in two-dimensional (2-D) and in Cartesian coordinates by focusing on the consideration of iron. The subdomains connection is carried out in the two directions (i.e., xand y-edges). The improvement was performed by solving magnetostatic Maxwell's equations for an air-or iron-core coil supplied by a direct current. To evaluate the efficacy of the proposed technique, the magnetic flux density distributions have been compared with those obtained by the 2-D finite-element analysis (FEA). The semi-analytical results are in quite satisfying agreement with those obtained by the 2-D FEA, considering both amplitude and waveform.3 of 37 has been applied to surface-mounted/-inset magnets machines [32][33][34][35], surface-inset magnet machines [33], and others electrical machines.-Concept wave impedance: They are based on a direct solution of Maxwell's field equations in homogeneous multi-layers of magnetic material properties, viz., the magnetic permeability and the electrical conductivity. This approach, developed by Mishkin (1953) [36], was first applied to squirrel-cage induction machine in Cartesian coordinates with three-layers (i.e., stator slotting, air-gap, and rotor slotting). It was used and enhanced by many authors, viz., * simplification of the electromagnetic theory [37]; * extended with an infinite number of layers [38]; * converted into equivalent circuits and terminal impedance [39]; * included the curvature effect with the magnetizing current [40]; * incorporated spatial harmonics in the multi-layers theory by considering isotropic and anisotropic (e.g., laminated, composite, and toothed) regions [41,42]; * introduced the nonlinear B(H) curve in homogenous layers by an iterative procedure [43,44]; * taking account of the effect of slot openings [45], i.e., the multi-layers model is combined with the subdomain technique for slotted structures by assuming infinitely permeable tooth tips; * included the current penetration effect in conductive layers [43,46]. The analytical solution for the electromagnetic field in conductive layers is then defined by Bessel functions.