The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations. Therefore, we shall introduce a modified fractional differential operator Dα proposed by M. Caputo in his work on the theory of viscoelasticity [11].