Bayesian approximation of solutions to linear ordinary differential equations