We establish analogues of the mean value theorem and Taylor’s theorem for fractional differential operators defined using a Mittag-Leffler kernel. We formulate a new model for the fractional Boussinesq equation by using this new Taylor series expansion.