A novel method to solve variable-order fractional delay differential equations based in lagrange interpolations.

In this work, we present a novel numerical method based on the fundamental theorem of fractional calculus and the Lagrange polynomial interpolation to solve numerically fractional delay differential equations. We focus on the fractional derivative with power-law, exponential decay and Mittag-Leffler kernel of Liouville-Caputo type with constant and variable-order. The numerical methods were applied to simulate the Duffing attractor, El-Niño/Southern-Oscillation, and Ikeda systems. [ABSTRACT FROM AUTHOR]