l
the <a href="../gp.m">gp</a> function.</p>

<center><img src="f1.gif" alt="f1.gif"></center><br>

<p>Let's ask the model to compute the (joint) negative log probability
(density) <tt>nlml</tt> (also called marginal likelihood or evidence)
and to generalize from the training data to other (test) inputs
<tt>z</tt>:</p>

<pre>
  nlml = gp(hyp, @infExact, meanfunc, covfunc, likfunc, x, y)

  z = linspace(-1.9, 1.9, 101)';
  [m s2] = gp(hyp, @infExact, meanfunc, covfunc, likfunc, x, y, z);

  f = [m+2*sqrt(s2)