prob = ODEProblem(myODE,eltype(prior).(x0),tspan,prior)
https://docs.juliadiffeq.org/latest/basics/faq/#Native-Julia-solvers-compatibility-with-autodifferentiation-1