Today I was asked whether it was possible to solve in Julia differential equations involving numbers with uncertainties. Of course the answer is yes. What I find really amazing about Julia is that the two packages don’t know anything about each other, yet they can work together without any effort. Here is an short example based on this tutorial: https://nbviewer.jupyter.org/gist/giordano/e82a3959d8f64301129d64d004e10b4e

The only think I found weird was that the extrema of the time interval had to be `Measurement`

objects, otherwise `solve`

would complain (because it can’t convert `Measurement`

type to `Float64`

). Why it’s that? I can see why `u0`

should have the same type as that returned by `f`

, but not why also the time should match the same type.