Machine learning of linear differential equations using Gaussian processes