The play model is a simple and efficient hysteresis model that is mathematically equivalent to the classical Preisach model. The play model requires only a one-dimensional integral, whereas the Preisach model is basically described by a two-dimensional integral. The relation between the Preisach distribution function and the shape function of the play model is explained. Based on the Everett integral for the Preisach model, a method for identification of the play model is presented, where the shape function is determined from given symmetric hysteresis loops. For the application to the magnetic field analysis using the magnetic vector potential, the B-input play model is also explained. Isotropic and anisotropic vector play models are presented, which are required for 2D and 3D magnetic field analyses.