We investigate the dynamical behavior of a fractional order toxic-phytoplankton–phytoplankton–zooplankton system.
Stability analysis of biologically feasible equilibrium points is presented in terms of reproduction numbers.
We study the local stability of the toxic-phytoplankton-free equilibrium.
Some numerical experiments have been carried out in support of our analytic findings.