Yibo Xia

I notice that in the latest version of ```dynamicalsystems.jl``` ds becomes the integrator itself. It's a good improvement but I don't know how to get the derivative via get_du() now....