DifferentialEquations.jl look at the ODEsystem

General Usage
7

du isn’t a matrix. But if you’re looking for the operator, in this case it’s equivalent to the Jacobian, so just calculate the Jacobian.