Estimation of Ordinary Differential Equations Solutions with Gaussian Processes and Polynomial Chaos Expansion