BTW, I gather that you are using scipy function to generate ODE? You might find this interesting as well: Dynamical model in pymc3
BTW, I gather that you are using scipy function to generate ODE? You might find this interesting as well: Dynamical model in pymc3