Bayesian Integral-based Methods for Differential Equation Models