Abstract
Fractional derivative operators of non-integer order can be utilized as a powerful tool to model nonlinear fractional differential equations. In this paper, we propose numerical solutions for simulating fractional-order derivative operators with the power-law and exponential-law kernels. We construct the numerical schemes with the help the fundamental theorem of fractional calculus and the Lagrange polynomial interpolation. These schemes are applied to simulate the dynamical fractional-order model of the immune response (FMIR) to the uncomplicated influenza A virus (IAV) infection, which focuses on the control of the infection by the innate and adaptive immunity. Numerical results are then presented to show the applicability and efficiency on the FMIR.