I was curious about Gaussian processes (GPs) so I implemented one and used it to approximate a Van Der Pol oscillator. Here’s a sample plot:

The GP's predictions for position given 10 training points. Ribbon plot shows predicted uncertainty.

Shout out to Jie Wang and his tutorial on GP regression, it had everything I needed to implement this.

Find the code here.