How to implement fractional-order LIF model using Brian2?

I am trying to implement a fractional-order leaky integrate and fire model according to this paper: Fractional-order leaky integrate-and-fire model with long-term memory and power law dynamics

The fractional-order LIF model involves a fractional differential equation to describe the evolution of membrane potentials. I am wondering if Brian2 supports simulation of such fractional-order models (fractional differential equations).