 times input-output covariance</li><li>Compute predictive covariance matrix, non-central moments</li><li>Centralize moments</li><li>Vectorize derivatives</li></ol></div><pre class="codeinput"><span class="keyword">function</span> [M, S, V, dMdm, dSdm, dVdm, dMds, dSds, dVds, dMdi, dSdi, dVdi, <span class="keyword">...</span>
  dMdt, dSdt, dVdt, dMdX, dSdX, dVdX] = gp2d(gpmodel, m, s)
</pre><h2>Code<a name="2"></a></h2><pre class="codeinput">input = gpmodel.inputs;  target = gpmodel.targets; X = gpmodel.hyp;