In this paper, a fast H(M) magnetostatic hysteresis model is presented, tested, and used in Simulink and partial element equivalent circuit (PEEC) environments. The model is primarily designed to achieve full compatibility with the current 3-D PEEC formalism for nonlinear isotropic magnetic materials, compatibility not met by any direct model. The capabilities of the proposed model are shown in stand-alone configuration and as a core of circuit components within the Simulink environment. Its effectiveness in the 3-D PEEC formulation is then demonstrated.