def toy_rbf_1d(optimize=True, plot=True):
try:import pods<EOL>except ImportError:<EOL><INDENT>print('<STR_LIT>')<EOL>return<EOL><DEDENT>data = pods.datasets.toy_rbf_1d()<EOL>m = GPy.models.GPRegression(data['<STR_LIT:X>'], data['<STR_LIT:Y>'])<EOL>if optimize:<EOL><INDENT>m.optimize('<STR_LIT>')<EOL><DEDENT>if plot:<EOL><INDENT>m.plot()<EOL><DEDENT>return m<EOL>
Run a simple demonstration of a standard Gaussian process fitting it to data sampled from an RBF covariance.