“…Such system can be very difficult to model and analyze with classical differential equations, but non – locality gives fractional derivative built-in ability to incorporate memory effects [11] . Fractional differential equations appear naturally in numerous fields of study including physics, polymer rheology, regular variation in thermodynamics, biophysics, blood flow phenomena, aerodynamics, electrodynamics of complex medium, viscoelasticity, capacitor theory, electrical circuits, electron-analytical chemistry, biology, control theory, and fitting of experimental data [ [12] , [13] , [14] , [23] , [24] , [25] , [26] ]. Recently there are many studies on epidemiological disease modeling using fractional order differential equations [27] , [28] , [29] , [30] , [31] , [32] , [33] .…”