Add plot with variance around points for GP's and models that have a `return_std()` argument

[mastoffel]

This might need a change in the plotting interface to make custom plotting for specific models possible. So potentially a bigger issue.

[mastoffel]

now addressed in #238