2D arrays with ModelingToolkit

This actually works:

# initial state
u0 = Dict(pos=>POS0, vel=>VEL0)

But I could not find that in the documentation of ODEProblem…