I have a solution to a DAE via DifferentialEquations.jl, which uses Hermite interpolation as can be seen here:

```
julia> sol
retcode: Success
Interpolation: 3rd order Hermite
...
```

(the example in the tutorial page linked above can be used as a MWE if you need one)

I can use the interpolation to obtain solutions to my system via:

```
julia> sol(1987.4)
5-element Array{Float64,1}:
1877.24
913.053
3522.07
0.512887
731.781
```

which are interpolated values of `sol.u`

at time `t = 1987.4`

.

`sol.u`

is an `Array{Array{Float64,1},1}`

with a corresponding array `sol.t`

calculated at specific times.

I’d also like to access the interpolation for my derivatives `sol.du`

, which again is an `Array{Array{Float64,1},1}`

corresponding to `sol.t`

.

Is it possible to do this? Neither `sol(1987.4).du`

or `sol.du(1987.4)`

work.

I’d be happy to set up my own linear interpolation of the `sol.du`

array after the fact, but at face value it seems like accessing this data should be possible already since `struct HermiteInterpolation`

uses `t, u, du`

, and I’ve just missed that section of the docs.