A numerical method for solving fractional delay differential equations based on the operational matrix method. (June 2021)