In this work, we employ the Bayesian inference framework to solve problem of estimating solution and particularly, its derivatives, which satisfy a known differential equation, from given noisy scarce observations data only. To address key issue accuracy robustness derivative estimation, use Gaussian processes jointly model solution, equation. By regarding linear equation as constraint, process...