Numerical simulations study of a spin-1 Blume–Emery–Griffiths model on a square lattice