I am not sure that I understand the structure of `sol.k`

.

For an ODE with 3 state variable, `sol.k[1:end][1]`

which should contain the first derivatives at all time is only of Array of Array of 3 elements, I guess it contains the derivatives at time t_0.

And `sol.k[:][:]`

is an Array of Array of Array. I tried to play around with `VectorOfArray()`

and `convert(Array,..)`

but I can only reduce it to a Array of Array

Is there a simple way to get an Array of type Array{Float64,2} of all state derivative at different time steps.

```
function lorenz(du,u,p,t)
du[1] = 10.0*(u[2]-u[1])
du[2] = u[1]*(28.0-u[3]) - u[2]
du[3] = u[1]*u[2] - (8/3)*u[3]
end
u0 = [1.0;0.0;0.0]
tspan = (0.0,100.0)
prob = ODEProblem(lorenz,u0,tspan)
sol = solve(prob)
sol.k[:][1]
1-element Array{Array{Float64,1},1}:
[4.94066e-324, 1.97626e-323, 2.47033e-323]
sol.k[:][:]
1250-element Array{Array{Array{Float64,1},1},1}:
[[4.94066e-324, 1.97626e-323, 2.47033e-323]]
[[-10.0, 28.0, 0.0], [-9.99782, 27.9982, 0.00016083], [-9.99557, ...blablabla
```

Thank you for your help,