I’m working on model using Differential Equation that has two state variables, current and temperature. So when using `DifferentialEquations.jl`

I define something like

`f(du, u, p, t) = ...`

`u`

is a vector with two elements, each of which I’d like to have different units and therefore different types. And du has yet different units from `u`

since we add `1/s`

when we take the time derivative.

Any advice or packages on how to deal with this in a type stable way? For now I strip all the units off and pass floats to `f`

, but I’d rather propagate the units all the way through if possible.