ov_exp_quad(x, 1.0, 1.0);
  vector[N] mu = rep_vector(0, N);
  for (n in 1:N)
    K[n, n] = K[n, n] + 0.1;
}
parameters {
  vector[N] y;
}
model {
  y ~ multi_normal(mu, K);
}</code></pre>
<p>The input data are just the vector of inputs <code>x</code> and its size
<code>N</code>. Such a model can be used with values of <code>x</code> evenly
spaced over some interval in order to plot sample draws of functions
from a Gaussian process.</p>
<div id="multivariate-inputs" class="section level3">
<h3 class="unnumb