The application of the differential operator method to the Blume-Emery-Griffiths model