Numerical discretization for fractional differential equations with feedback control

In this paper, we introduce a numerical scheme for fractional differential equations with feedback control. Due to the possibility of dealing with a feedback control as a functional delay, we construct a numerical method based on Euler method accompanied with piecewise constant interpolation. The method is based on the idea of separating the current state and the prehistory function. The convergence of the method is stated and proved. Numerical experiments are given to clarify the good agreement between numerical and theoretical results.